| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 /// file: DamageModel.hh
28 /// brief: Parameters related to the DNA Damage model
29 #ifndef MOLECULAR_DAMAGE_MODEL
30 #define MOLECULAR_DAMAGE_MODEL
31
32 #include "DamageModelMessenger.hh"
33
34 #include "G4Electron_aq.hh"
35 #include "G4Hydrogen.hh"
36 #include "G4MolecularConfiguration.hh"
37 #include "G4OH.hh"
38 #include "G4SystemOfUnits.hh"
39 #include "Randomize.hh"
40 #include "globals.hh"
41
42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43
44 class DamageModel
45 {
46 public:
47 DamageModel();
48
49 virtual ~DamageModel() { delete fpDamageMessenger; };
50
51 // methods
52 G4bool IsDirectStrandBreak(const G4double&) const;
53
54 G4bool IsIndirectStrandDamage(const G4MolecularConfiguration*) const;
55
56 G4bool IsIndirectBaseDamage(const G4MolecularConfiguration*) const;
57
58 G4bool IsInducedStrandBreak(const G4MolecularConfiguration*) const;
59
60 // setters
61 void SetDirectDamageLower(G4double d)
62 {
63 fDirectDmgLower = d;
64 if (fDirectDmgLower < fDirectDmgUpper) fDirectDmgUpper = fDirectDmgLower;
65 };
66
67 void SetDirectDamageUpper(G4double d)
68 {
69 fDirectDmgUpper = d;
70 if (fDirectDmgLower > fDirectDmgUpper) fDirectDmgLower = fDirectDmgUpper;
71 };
72
73 // Indirect OH
74 void SetIndirectOHStrandChance(G4double d)
75 {
76 CheckValidProbability("DamageStrandOH", d);
77 fStrandOH = d;
78 }
79
80 void SetIndirectOHBaseChance(G4double d)
81 {
82 CheckValidProbability("DamageBaseOH", d);
83 fBaseOH = d;
84 }
85
86 void SetInductionOHChance(G4double d)
87 {
88 CheckValidProbability("InductionOH", d);
89 fInductionOH = d;
90 };
91
92 // Indirect Eaq
93 void SetIndirectEaqStrandChance(G4double d)
94 {
95 CheckValidProbability("DamageStrandEaq", d);
96 fStrandEaq = d;
97 }
98
99 void SetIndirectEaqBaseChance(G4double d)
100 {
101 CheckValidProbability("DamageBaseEaq", d);
102 fBaseEaq = d;
103 }
104
105 void SetInductionEaqChance(G4double d)
106 {
107 CheckValidProbability("InductionEaq", d);
108 fInductionEaq = d;
109 };
110
111 // Indirect H
112 void SetIndirectHStrandChance(G4double d)
113 {
114 CheckValidProbability("DamageStrandH", d);
115 fStrandH = d;
116 }
117
118 void SetIndirectHBaseChance(G4double d)
119 {
120 CheckValidProbability("DamageBaseH", d);
121 fBaseH = d;
122 }
123
124 void SetInductionHChance(G4double d)
125 {
126 CheckValidProbability("InductionH", d);
127 fInductionH = d;
128 };
129
130 protected:
131 static void CheckValidProbability(const G4String& str, G4double p);
132
133 private:
134 DamageModelMessenger* fpDamageMessenger;
135 G4double fDirectDmgLower = 17.5 * eV, fDirectDmgUpper = 17.5 * eV;
136 G4double fInductionOH = 0, fInductionEaq = 0, fInductionH = 0;
137 G4double fStrandOH = 1.0, fStrandEaq = 1.0, fStrandH = 1.0;
138 G4double fBaseOH = 1.0, fBaseEaq = 1.0, fBaseH = 1.0;
139 const G4MoleculeDefinition* fOH = G4OH::Definition();
140 const G4MoleculeDefinition* fe_aq = G4Electron_aq::Definition();
141 const G4MoleculeDefinition* fH = G4Hydrogen::Definition();
142 };
143
144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
145
146 #endif // MOLECULAR_DAMAGE_MODEL
147