| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // Based on the work of M. Terrissol and M. C.
27 //
28 // Users are requested to cite the following p
29 // - M. Terrissol, A. Baudre, Radiat. Prot. Do
30 // - M.C. Bordage, J. Bordes, S. Edel, M. Terr
31 // M. Bardies, N. Lampe, S. Incerti, Phys. M
32 //
33 // Authors of this class:
34 // M.C. Bordage, M. Terrissol, S. Edel, J. Bor
35 //
36 // 15.01.2014: creation
37 //
38 //
39 // Modified for Adenine material by S. Zein on
40 // Based on the study by S. Zein et. al. Nucl.
41
42 #include "G4DNACPA100ExcitationStructure.hh"
43
44 #include "G4Material.hh"
45 #include "G4SystemOfUnits.hh"
46 G4DNACPA100ExcitationStructure::G4DNACPA100Exc
47 {
48 fpGuanine = G4Material::GetMaterial("G4_GUAN
49 fpG4_WATER = G4Material::GetMaterial("G4_WAT
50 fpDeoxyribose = G4Material::GetMaterial("G4_
51 fpCytosine = G4Material::GetMaterial("G4_CYT
52 fpThymine = G4Material::GetMaterial("G4_THYM
53 fpAdenine = G4Material::GetMaterial("G4_ADEN
54 fpPhosphate = G4Material::GetMaterial("G4_PH
55
56 if (fpGuanine != nullptr) {
57 InitialiseGuanine();
58 }
59
60 if (fpG4_WATER != nullptr) {
61 InitialiseWater();
62 }
63
64 if (fpDeoxyribose != nullptr) {
65 InitialiseDeoxyribose();
66 }
67
68 if (fpCytosine != nullptr) {
69 InitialiseCytosine();
70 }
71
72 if (fpThymine != nullptr) {
73 InitialiseThymine();
74 }
75
76 if (fpAdenine != nullptr) {
77 InitialiseAdenine();
78 }
79
80 if (fpPhosphate != nullptr) {
81 InitialisePhosphate();
82 }
83 }
84
85 void G4DNACPA100ExcitationStructure::Initialis
86 {
87 auto index = fpGuanine->GetIndex();
88 fEnergyConstant[index].push_back(8.2315 * eV
89 fEnergyConstant[index].push_back(11.0928 * e
90 fEnergyConstant[index].push_back(11.5984 * e
91 fEnergyConstant[index].push_back(11.7906 * e
92 fEnergyConstant[index].push_back(11.9382 * e
93 fEnergyConstant[index].push_back(12.4424 * e
94 fEnergyConstant[index].push_back(13.3581 * e
95 fEnergyConstant[index].push_back(15.1381 * e
96 fEnergyConstant[index].push_back(16.4059 * e
97 fEnergyConstant[index].push_back(16.496 * eV
98 fEnergyConstant[index].push_back(16.8457 * e
99 fEnergyConstant[index].push_back(17.297 * eV
100 fEnergyConstant[index].push_back(18.0608 * e
101 fEnergyConstant[index].push_back(18.441 * eV
102 fEnergyConstant[index].push_back(19.2414 * e
103
104 fnLevels[index] = fEnergyConstant[index].siz
105 }
106 void G4DNACPA100ExcitationStructure::Initialis
107 {
108 auto index = fpG4_WATER->GetIndex();
109 // The following values are extracted from t
110 // Paul Sabatier University, Toulouse, Franc
111 // Page 36
112 fEnergyConstant[index].push_back(8.17 * eV);
113 fEnergyConstant[index].push_back(10.13 * eV)
114 fEnergyConstant[index].push_back(11.31 * eV)
115 fEnergyConstant[index].push_back(12.91 * eV)
116 fEnergyConstant[index].push_back(14.50 * eV)
117
118 fUConstant[index].push_back(61.91 * eV);
119 fUConstant[index].push_back(59.52 * eV);
120 fUConstant[index].push_back(48.36 * eV);
121 fUConstant[index].push_back(70.71 * eV);
122 fUConstant[index].push_back(796.2 * eV);
123
124 fnLevels[index] = fEnergyConstant[index].siz
125 }
126 void G4DNACPA100ExcitationStructure::Initialis
127 {
128 auto index = fpDeoxyribose->GetIndex();
129 fEnergyConstant[index].push_back(11.2410 * e
130 fEnergyConstant[index].push_back(11.7927 * e
131 fEnergyConstant[index].push_back(12.6579 * e
132 fEnergyConstant[index].push_back(12.8163 * e
133 fEnergyConstant[index].push_back(13.3238 * e
134 fEnergyConstant[index].push_back(13.9487 * e
135 fEnergyConstant[index].push_back(14.4374 * e
136 fEnergyConstant[index].push_back(14.7433 * e
137 fEnergyConstant[index].push_back(15.0818 * e
138 fEnergyConstant[index].push_back(15.6112 * e
139 fEnergyConstant[index].push_back(16.0547 * e
140 fEnergyConstant[index].push_back(16.8319 * e
141 fEnergyConstant[index].push_back(17.4294 * e
142 fEnergyConstant[index].push_back(18.0000 * e
143 fEnergyConstant[index].push_back(18.2696 * e
144 fEnergyConstant[index].push_back(18.6049 * e
145 fEnergyConstant[index].push_back(19.8378 * e
146
147 fnLevels[index] = fEnergyConstant[index].siz
148 }
149 void G4DNACPA100ExcitationStructure::Initialis
150 {
151 auto index = fpCytosine->GetIndex();
152 // The following values are extracted from t
153 fEnergyConstant[index].push_back(9.3222 * eV
154 fEnergyConstant[index].push_back(10.4601 * e
155 fEnergyConstant[index].push_back(11.3044 * e
156 fEnergyConstant[index].push_back(11.9986 * e
157 fEnergyConstant[index].push_back(13.4528 * e
158 fEnergyConstant[index].push_back(14.7371 * e
159 fEnergyConstant[index].push_back(16.2286 * e
160 fEnergyConstant[index].push_back(16.5877 * e
161 fEnergyConstant[index].push_back(17.0741 * e
162 fEnergyConstant[index].push_back(17.1875 * e
163 fEnergyConstant[index].push_back(18.638 * eV
164 fEnergyConstant[index].push_back(19.6884 * e
165 fnLevels[index] = fEnergyConstant[index].siz
166 }
167 void G4DNACPA100ExcitationStructure::Initialis
168 {
169 // The following values are extracted from t
170 auto index = fpThymine->GetIndex();
171
172 fEnergyConstant[index].push_back(9.639 * eV)
173 fEnergyConstant[index].push_back(11.8278 * e
174 fEnergyConstant[index].push_back(12.0876 * e
175 fEnergyConstant[index].push_back(12.9656 * e
176 fEnergyConstant[index].push_back(13.9555 * e
177 fEnergyConstant[index].push_back(15.0774 * e
178 fEnergyConstant[index].push_back(15.4078 * e
179 fEnergyConstant[index].push_back(15.4689 * e
180 fEnergyConstant[index].push_back(16.1964 * e
181 fEnergyConstant[index].push_back(16.8955 * e
182 fEnergyConstant[index].push_back(17.5018 * e
183 fEnergyConstant[index].push_back(18.2979 * e
184 fEnergyConstant[index].push_back(18.4495 * e
185 fEnergyConstant[index].push_back(19.3186 * e
186 fnLevels[index] = fEnergyConstant[index].siz
187 }
188 void G4DNACPA100ExcitationStructure::Initialis
189 {
190 auto index = fpAdenine->GetIndex();
191 // The following values are extracted from t
192
193 fEnergyConstant[index].push_back(8.5114 * eV
194 fEnergyConstant[index].push_back(10.1294 * e
195 fEnergyConstant[index].push_back(11.0606 * e
196 fEnergyConstant[index].push_back(11.5849 * e
197 fEnergyConstant[index].push_back(12.1533 * e
198 fEnergyConstant[index].push_back(13.356 * eV
199 fEnergyConstant[index].push_back(13.6554 * e
200 fEnergyConstant[index].push_back(15.3296 * e
201 fEnergyConstant[index].push_back(16.179 * eV
202 fEnergyConstant[index].push_back(16.7676 * e
203 fEnergyConstant[index].push_back(17.3489 * e
204 fEnergyConstant[index].push_back(17.5568 * e
205 fEnergyConstant[index].push_back(18.554 * eV
206 fEnergyConstant[index].push_back(19.0866 * e
207 fnLevels[index] = fEnergyConstant[index].siz
208 }
209 void G4DNACPA100ExcitationStructure::Initialis
210 {
211 // The following values are extracted from t
212 auto index = fpPhosphate->GetIndex();
213 fEnergyConstant[index].push_back(12.9963 * e
214 fEnergyConstant[index].push_back(12.9972 * e
215 fEnergyConstant[index].push_back(14.3109 * e
216 fEnergyConstant[index].push_back(15.2221 * e
217 fEnergyConstant[index].push_back(16.0591 * e
218 fEnergyConstant[index].push_back(16.0622 * e
219 fEnergyConstant[index].push_back(17.6365 * e
220 fEnergyConstant[index].push_back(17.6401 * e
221 fEnergyConstant[index].push_back(18.8803 * e
222 fnLevels[index] = fEnergyConstant[index].siz
223 }
224
225 //....oooOO0OOooo........oooOO0OOooo........oo
226
227 G4double G4DNACPA100ExcitationStructure::Excit
228
229 {
230 G4double ionisation = 0.;
231
232 if (level < fnLevels[MatID]) {
233 ionisation = fEnergyConstant[MatID][level]
234 }
235 else {
236 std::ostringstream oss;
237 oss << " material was not found. ";
238 G4Exception("G4DNACPA100ExcitationStructur
239 oss.str().c_str());
240 }
241
242 return ionisation;
243 }
244
245 //....oooOO0OOooo........oooOO0OOooo........oo
246
247 G4double G4DNACPA100ExcitationStructure::UEner
248 {
249 G4double UEnergy = 0.;
250 if (level < fnLevels[MatID]) {
251 UEnergy = fUConstant[MatID][level];
252 }
253 else {
254 std::ostringstream oss;
255 oss << " material was not found. ";
256 G4Exception("G4DNACPA100ExcitationStructur
257 oss.str().c_str());
258 }
259
260 return UEnergy;
261 }
262