source: trunk/examples/advanced/xray_fluorescence/src/XrayFluoNistMaterials.cc@ 809

Last change on this file since 809 was 807, checked in by garnier, 17 years ago

update
File size: 16.8 KB
RevLine 
[807]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// $Id: XrayFluoDetectorConstruction.hh
28// GEANT4 tag $Name: xray_fluo-V03-02-00
29//
30// Author: Alfonso Mantero (Alfonso.Mantero@ge.infn.it)
31//
32// History:
33// -----------
34// 20 Aug 2001 Alfonso Mantero Created
35//
36// -------------------------------------------------------------------
37
38#include "XrayFluoNistMaterials.hh"
39
40XrayFluoNistMaterials::XrayFluoNistMaterials()
41{ CreateMaterials();}
42
43XrayFluoNistMaterials::~XrayFluoNistMaterials()
44{
45 delete dolorite;
46 delete HPGe;
47 delete mars1;
48 delete galactic;
49 delete madaBasalt;
50 delete icelandicBasalt;
51 delete GaAs;
52}
53XrayFluoNistMaterials* XrayFluoNistMaterials::instance = 0;
54
55XrayFluoNistMaterials* XrayFluoNistMaterials::GetInstance()
56{
57 if (instance == 0)
58 {
59 instance = new XrayFluoNistMaterials;
60
61 }
62 return instance;
63}
64
65G4Material* XrayFluoNistMaterials::GetMaterial(G4String material)
66{
67
68 //instancing G4NistManager
69 nistMan = G4NistManager::Instance();
70 nistMan->SetVerbose(0);
71
72 G4Material* mat = nistMan->FindOrBuildMaterial(material);
73 if (!mat) {
74 mat = G4Material::GetMaterial(material);
75 }
76 if (!mat) {G4cout << material << "Not Found, Please Retry"<< G4endl;}
77 return mat;
78}
79
80
81void XrayFluoNistMaterials::CreateMaterials()
82{
83
84 G4double density;
85 std::vector<G4int> natoms;
86 std::vector<G4double> fractionMass;
87 std::vector<G4String> elements;
88
89
90 //instancing G4NistManager
91 nistMan = G4NistManager::Instance();
92 nistMan->SetVerbose(1);
93
94 // Materials Definitions
95
96
97 ///////////////////////
98 // Madagascar Basalt //
99 ///////////////////////
100
101
102 // Define Madagascar Basalt main components 0054.PP.0044 sample
103 density = 3*g/cm3;
104 elements.push_back("Si"); fractionMass.push_back(0.1992); // 0.007093 mol/g(mat)
105 elements.push_back("Ti"); fractionMass.push_back(0.02027); // 4.235e-4
106 elements.push_back("Al"); fractionMass.push_back(0.04758); // 0.001763
107 elements.push_back("Fe"); fractionMass.push_back(0.1303); // 0.002333
108 elements.push_back("Mn"); fractionMass.push_back(0.001549);// 2.820e-5
109 elements.push_back("Mg"); fractionMass.push_back(0.08141); // 0.003350
110 elements.push_back("Ca"); fractionMass.push_back(0.06468); // 0.001614
111 elements.push_back("Na"); fractionMass.push_back(0.01692); // 7.360e-4
112 elements.push_back("K"); fractionMass.push_back(0.008576);// 2.193e-4
113 elements.push_back("P"); fractionMass.push_back(0.001977);// 6.383e-5
114 elements.push_back("O"); fractionMass.push_back(0.427538);// 0.02672
115
116 // sum is 0.04434383 total number of moles of atoms in one gram of material
117 // 2.248766e8 g per 10.000.000 moles.
118
119 G4Material* madaBasaltMain= nistMan->ConstructNewMaterial("MadaBasaltMain",elements, fractionMass, density);
120 elements.clear();
121 fractionMass.clear();
122
123// Define Madagascar Basalt traces components 0054.PP.0044 sample
124 density = 3*g/cm3;
125
126
127 elements.push_back("Ti"); natoms.push_back(33);
128 elements.push_back("Ba"); natoms.push_back(4131);
129 elements.push_back("Ce"); natoms.push_back(694);
130 elements.push_back("Co"); natoms.push_back(965);
131 elements.push_back("Cr"); natoms.push_back(5584);
132 elements.push_back("La"); natoms.push_back(269);
133 elements.push_back("Nb"); natoms.push_back(259);
134 elements.push_back("Nd"); natoms.push_back(410);
135 elements.push_back("Ni"); natoms.push_back(389);
136 elements.push_back("Rb"); natoms.push_back(227);
137 elements.push_back("Sc"); natoms.push_back(212);
138 elements.push_back("Sr"); natoms.push_back(8686);
139 elements.push_back("V"); natoms.push_back(4203);
140 elements.push_back("Y"); natoms.push_back(272);
141 elements.push_back("Zn"); natoms.push_back(1440);
142 elements.push_back("Th"); natoms.push_back(19);
143 elements.push_back("Sm"); natoms.push_back(93);
144 elements.push_back("Eu"); natoms.push_back(32);
145 elements.push_back("Gd"); natoms.push_back(89);
146 elements.push_back("Tb"); natoms.push_back(13);
147 elements.push_back("Yb"); natoms.push_back(15);
148 elements.push_back("Lu"); natoms.push_back(2);
149 elements.push_back("Ta"); natoms.push_back(15);
150 elements.push_back("Hf"); natoms.push_back(62);
151
152 //tot 28114/10e7 weight: 2335253.28 g per 10e6 moles
153
154 G4Material* madaBasaltTraces= nistMan->ConstructNewMaterial("MadaBasaltTraces", elements, natoms, density);
155 elements.clear();
156 natoms.clear();
157
158 // Define Madacagascar Basalt complete material 0054.PP.0044 sample
159 density = 3*g/cm3;
160
161 madaBasalt= new G4Material("MadaBasalt", density, 2);
162 madaBasalt->AddMaterial(madaBasaltMain, 0.9897);
163 madaBasalt->AddMaterial(madaBasaltTraces, 0.0103);
164
165
166
167 ///////////////////////
168 // Iceland Basalt //
169 ///////////////////////
170
171 elements.push_back("Si"); fractionMass.push_back(0.2313);
172 elements.push_back("Ti"); fractionMass.push_back(0.0127);
173 elements.push_back("Al"); fractionMass.push_back(0.0702);
174 elements.push_back("Fe"); fractionMass.push_back(0.1134);
175 elements.push_back("Mn"); fractionMass.push_back(0.0019);
176 elements.push_back("Mg"); fractionMass.push_back(0.0349);
177 elements.push_back("Ca"); fractionMass.push_back(0.0756);
178 elements.push_back("Na"); fractionMass.push_back(0.0892);
179 elements.push_back("K"); fractionMass.push_back(0.0032);
180 elements.push_back("P"); fractionMass.push_back(0.00096);
181 elements.push_back("S"); fractionMass.push_back(0.0004);
182 elements.push_back("O"); fractionMass.push_back(0.36624);
183
184 // Define Icelandic Basalt main components 0029.PP.0035 sample
185 density = 3*g/cm3;
186 G4Material* icelandicBasaltMain= nistMan->ConstructNewMaterial("IceBasaltMain",elements, fractionMass, density);
187
188 // Define Icelandic Basalt traces components 0029.PP.0035 sample
189 density = 3*g/cm3;
190
191 elements.push_back("Ba"); natoms.push_back(756);
192 elements.push_back("Ce"); natoms.push_back(328);
193 elements.push_back("Co"); natoms.push_back(643);
194 elements.push_back("Cr"); natoms.push_back(1000);
195 elements.push_back("Cu"); natoms.push_back(1396);
196 elements.push_back("Ga"); natoms.push_back(190);
197 elements.push_back("La"); natoms.push_back(103);
198 elements.push_back("Mo"); natoms.push_back(9);
199 elements.push_back("Nb"); natoms.push_back(114);
200 elements.push_back("Nd"); natoms.push_back(104);
201 elements.push_back("Ni"); natoms.push_back(544);
202 elements.push_back("Rb"); natoms.push_back(78);
203 elements.push_back("S"); natoms.push_back(5550);
204 elements.push_back("Sc"); natoms.push_back(531);
205 elements.push_back("Sr"); natoms.push_back(1353);
206 elements.push_back("U"); natoms.push_back(22);
207 elements.push_back("V"); natoms.push_back(4533);
208 elements.push_back("Y"); natoms.push_back(408);
209 elements.push_back("Zn"); natoms.push_back(1259);
210 elements.push_back("Zr"); natoms.push_back(1274);
211
212 G4Material* icelandicBasaltTraces= nistMan->ConstructNewMaterial("IceBasaltTraces", elements, natoms, density);
213
214 elements.clear();
215 natoms.clear();
216
217 // Define Icelandic Basalt complete material 0029.PP.0035 sample
218 density = 3*g/cm3;
219 icelandicBasalt= new G4Material("IceBasalt", density, 2);
220 icelandicBasalt->AddMaterial(icelandicBasaltMain, 0.9978);
221 icelandicBasalt->AddMaterial(icelandicBasaltTraces, 0.0022);
222
223
224 ///////////////////////
225 // Dolorite //
226 ///////////////////////
227
228 // Define dolorite main components 0055.PP.0038 sample
229
230 density = 3*g/cm3;
231
232 elements.push_back("Fe"); fractionMass.push_back(0.1750);
233 elements.push_back("Ti"); fractionMass.push_back(0.0082);
234 elements.push_back("Ca"); fractionMass.push_back(0.0753);
235 elements.push_back("Si"); fractionMass.push_back(0.2188);
236 elements.push_back("Al"); fractionMass.push_back(0.0676);
237 elements.push_back("Mg"); fractionMass.push_back(0.0008);
238 elements.push_back("O"); fractionMass.push_back(0.4377);
239 elements.push_back("Mn"); fractionMass.push_back(0.0015);
240 elements.push_back("Na"); fractionMass.push_back(0.0134);
241 elements.push_back("K"); fractionMass.push_back(0.0011);
242 elements.push_back("P"); fractionMass.push_back(0.0006);
243
244
245 G4Material* dolorite = nistMan->ConstructNewMaterial("Dolorite", elements, fractionMass, density);
246
247 elements.clear();
248 natoms.clear();
249
250 // define traces in dolorite 0055.PP.0038 sample
251
252 density = 3*g/cm3;
253
254 elements.push_back("Nb"); natoms.push_back(5);
255 elements.push_back("Zr"); natoms.push_back(91);
256 elements.push_back("Y"); natoms.push_back(29);
257 elements.push_back("Sr"); natoms.push_back(140);
258 elements.push_back("Rb"); natoms.push_back(3);
259 elements.push_back("Ga"); natoms.push_back(20);
260 elements.push_back("Zn"); natoms.push_back(99);
261 elements.push_back("Ni"); natoms.push_back(77);
262 elements.push_back("Sc"); natoms.push_back(32);
263 elements.push_back("V"); natoms.push_back(314);
264 elements.push_back("Cr"); natoms.push_back(130);
265 elements.push_back("Co"); natoms.push_back(56);
266 elements.push_back("Cu"); natoms.push_back(119);
267 elements.push_back("Ba"); natoms.push_back(38);
268 elements.push_back("Ce"); natoms.push_back(15);
269 elements.push_back("Nd"); natoms.push_back(9);
270
271 G4Material* tracesOfDolorite= nistMan->ConstructNewMaterial("TracesOfDolorite", elements, natoms, density);
272
273 elements.clear();
274 natoms.clear();
275
276 // define dolorite (full) -- 0055.PP.0038 sample
277
278 density = 3*g/cm3;
279 dolorite = new G4Material("Dolorite", density, 2);
280 dolorite->AddMaterial(tracesOfDolorite, 0.0027842352);
281 dolorite->AddMaterial(dolorite, 0.9972157648);
282
283 ///////////////////////
284 // Mars1 //
285 ///////////////////////
286
287
288 // define mars1 -- 01.PP.0030 sample
289
290 density = 3*g/cm3;
291
292 elements.push_back("Fe"); fractionMass.push_back(0.100916);
293 elements.push_back("Ti"); fractionMass.push_back(0.0186804);
294 elements.push_back("Ca"); fractionMass.push_back(0.0404091);
295 elements.push_back("Si"); fractionMass.push_back(0.196378);
296 elements.push_back("Al"); fractionMass.push_back(0.103282);
297 elements.push_back("Mg"); fractionMass.push_back(0.0241622);
298 elements.push_back("Mn"); fractionMass.push_back(0.00184331);
299 elements.push_back("Na"); fractionMass.push_back(0.0177908);
300 elements.push_back("K"); fractionMass.push_back(0.00574498);
301 elements.push_back("P"); fractionMass.push_back(0.00280169);
302 elements.push_back("O"); fractionMass.push_back(0.48799152);
303
304
305 G4Material* mars1Main = nistMan->ConstructNewMaterial("Mars1 Main components", elements, fractionMass, density);
306
307 elements.clear();
308 fractionMass.clear();
309
310 elements.push_back("Nb"); natoms.push_back(55);
311 elements.push_back("Zr"); natoms.push_back(433);
312 elements.push_back("Y"); natoms.push_back(58);
313 elements.push_back("Sr"); natoms.push_back(968);
314 elements.push_back("Rb"); natoms.push_back(16);
315 elements.push_back("Ga"); natoms.push_back(24);
316 elements.push_back("Zn"); natoms.push_back(109);
317 elements.push_back("Ni"); natoms.push_back(70);
318 elements.push_back("Sc"); natoms.push_back(21);
319 elements.push_back("V"); natoms.push_back(134);
320 elements.push_back("Cr"); natoms.push_back(141);
321 elements.push_back("Co"); natoms.push_back(30);
322 elements.push_back("Cu"); natoms.push_back(19);
323 elements.push_back("Ba"); natoms.push_back(580);
324 elements.push_back("Pb"); natoms.push_back(4);
325 elements.push_back("S"); natoms.push_back(444);
326 elements.push_back("U"); natoms.push_back(2);
327
328 density = 3*g/cm3;
329 G4Material* tracesOfMars1 = nistMan->ConstructNewMaterial("TracesOfMars1", elements, natoms, density);
330
331 elements.clear();
332 natoms.clear();
333
334 density = 3*g/cm3;
335 mars1 = new G4Material("Mars1", density, 2);
336 mars1->AddMaterial(tracesOfMars1, 0.0044963163);
337 mars1->AddMaterial(mars1Main, 0.9955036837);
338
339 ///////////////////////
340 // Anorthosite //
341 ///////////////////////
342
343
344 density = 2.8*g/cm3;
345
346 elements.push_back("Fe"); fractionMass.push_back(0.095283);
347 elements.push_back("Mn"); fractionMass.push_back(0.00137086);
348 elements.push_back("Ni"); fractionMass.push_back(5e-5);
349 elements.push_back("Cu"); fractionMass.push_back(5.2e-4);
350 elements.push_back("Na"); fractionMass.push_back(0.017635);
351 elements.push_back("Mg"); fractionMass.push_back(0.0245361);
352 elements.push_back("Al"); fractionMass.push_back(0.0800355);
353 elements.push_back("Si"); fractionMass.push_back(0.232204);
354 elements.push_back("Ca"); fractionMass.push_back(0.0635368);
355 elements.push_back("K"); fractionMass.push_back(0.00464912);
356 elements.push_back("C"); fractionMass.push_back(0.000837803);
357 elements.push_back("P"); fractionMass.push_back(0.00176742);
358 elements.push_back("Ti"); fractionMass.push_back(0.0240879);
359 elements.push_back("Cl"); fractionMass.push_back(0.00014);
360 elements.push_back("Pd"); fractionMass.push_back(0.00001);
361 elements.push_back("Cd"); fractionMass.push_back(0.00018);
362 elements.push_back("Ag"); fractionMass.push_back(0.00048);
363 elements.push_back("S"); fractionMass.push_back(0.00144);
364 elements.push_back("V"); fractionMass.push_back(0.00228);
365 elements.push_back("Ba"); fractionMass.push_back(0.00151);
366 elements.push_back("O"); fractionMass.push_back(0.447026);
367
368 anorthosite = nistMan->ConstructNewMaterial("Anorthosite", elements, fractionMass, density);
369
370 elements.clear();
371 fractionMass.clear();
372
373 //define gallium arsenide
374
375 elements.push_back("Ga"); natoms.push_back(1);
376 elements.push_back("As"); natoms.push_back(1);
377
378 density = 5.32 * g/cm3;
379 GaAs = nistMan->ConstructNewMaterial("gallium arsenide", elements, natoms, density);
380
381 elements.clear();
382 natoms.clear();
383
384
385 // define germanium
386
387 density = 5.32 * g/cm3;
388
389 elements.push_back("Ge"); natoms.push_back(1);
390 HPGe = nistMan->ConstructNewMaterial("HPGe",elements, natoms, density);
391
392 elements.clear();
393 natoms.clear();
394
395 //define scintillator
396
397 elements.push_back("C"); natoms.push_back(9);
398 elements.push_back("H"); natoms.push_back(10);
399
400 density = 1.032*g/cm3;
401 Sci = nistMan->ConstructNewMaterial("Scintillator", elements, natoms, density);
402
403 elements.clear();
404 natoms.clear();
405
406 //define vacuum
407
408 density = universe_mean_density; //from PhysicalConstants.h
409 G4double pressure = 3.e-18*pascal;
410 G4double temperature = 2.73*kelvin;
411 Vacuum = new G4Material("Galactic", 1., 1.01*g/mole, density,
412 kStateGas,temperature,pressure);
413
414 //define basalt
415 density = 3.*g/cm3;
416
417 elements.push_back("Fe"); fractionMass.push_back(0.1200);
418 elements.push_back("Ti"); fractionMass.push_back(0.0160);
419 elements.push_back("Ca"); fractionMass.push_back(0.0750);
420 elements.push_back("Si"); fractionMass.push_back(0.2160);
421 elements.push_back("Al"); fractionMass.push_back(0.0710);
422 elements.push_back("Mg"); fractionMass.push_back(0.0590);
423 elements.push_back("O"); fractionMass.push_back(0.4430);
424
425
426
427 basalt = nistMan->ConstructNewMaterial("Basalt", elements, fractionMass, density);
428
429 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
430}
431
Note: See TracBrowser for help on using the repository browser.