65 G4Element* elH =
new G4Element(
"Hydrogen",
"H", Z=1., A=1.00794*g/mole);
66 G4Element* elC =
new G4Element(
"Carbon",
"C", Z=6., A= 12.011 *g/mole);
67 G4Element* elN =
new G4Element(
"Nitrogen",
"N", Z=7., A= 14.00674*g/mole);
68 G4Element* elO =
new G4Element(
"Oxygen",
"O", Z=8., A= 15.9994*g/mole);
69 G4Element* elNa =
new G4Element(
"Sodium",
"Na", Z=11., A= 22.989768*g/mole);
70 G4Element* elSi =
new G4Element(
"Silicon",
"Si", Z=14., A= 28.0855*g/mole);
71 G4Element* elAr =
new G4Element(
"Argon",
"Ar", Z=18., A= 39.948*g/mole);
72 G4Element* elI =
new G4Element(
"Iodine",
"I", Z=53., A= 126.90447*g/mole);
73 G4Element* elCs =
new G4Element(
"Cesium",
"Cs", Z=55., A= 132.90543*g/mole);
78 G4double density, massfraction;
82 const G4double expTemp= STP_Temperature+20.*kelvin;
85 density= universe_mean_density;
86 G4Material* Vacuum=
new G4Material(
"Vacuum", density, nel=2);
87 Vacuum-> AddElement(elN, .7);
88 Vacuum-> AddElement(elO, .3);
91 density= 1.2929e-03 *g/cm3;
92 G4Material* Air=
new G4Material(
"Air", density, nel=3,
94 G4double ttt= 75.47+23.20+1.28;
95 Air-> AddElement(elN, massfraction= 75.47/ttt);
96 Air-> AddElement(elO, massfraction= 23.20/ttt);
97 Air-> AddElement(elAr, massfraction= 1.28/ttt);
101 const G4double denAr= 1.782e-03 *g/cm3 * STP_Temperature/expTemp;
102 G4Material* Ar=
new G4Material(
"ArgonGas", Z=18., A, denAr,
106 const G4double denEthane= 1.356e-3 *g/cm3 * STP_Temperature/expTemp;
107 G4Material* Ethane=
new G4Material(
"Ethane", denEthane, nel=2,
109 Ethane-> AddElement(elC, natoms=2);
110 Ethane-> AddElement(elH, natoms=6);
113 density= (denAr+denEthane)/2.;
114 G4Material* ArEthane=
new G4Material(
"ArEthane", density, nel=2,
116 ArEthane-> AddMaterial(Ar, massfraction= denAr/density/2.);
117 ArEthane-> AddMaterial(Ethane, massfraction= denEthane/density/2.);
121 density= 2.33 *g/cm3;
122 G4Material* Si=
new G4Material(
"SiliconWafer", Z=14., A, density);
126 density= 2.70 *g/cm3;
127 G4Material* Al=
new G4Material(
"Al", Z=13., A, density);
131 density= 7.87 *g/cm3;
132 G4Material* Fe=
new G4Material(
"Iron", Z=26., A, density);
136 density= 11.35 *g/cm3;
137 G4Material* Pb=
new G4Material(
"Lead", Z=82., A, density);
140 density= 1.032 *g/cm3;
141 G4Material* Scinti=
new G4Material(
"Scinti", density, nel=2);
142 Scinti-> AddElement(elC, natoms=8);
143 Scinti-> AddElement(elH, natoms=8);
146 density= 2.64 *g/cm3;
147 G4Material* Quartz=
new G4Material(
"Quartz", density, nel= 2);
148 Quartz-> AddElement(elSi, natoms=1);
149 Quartz-> AddElement(elO, natoms=2);
152 density= 3.67 *g/cm3;
153 G4Material* NaI=
new G4Material(
"NaI", density, nel= 2);
154 NaI-> AddElement(elNa, natoms=1);
155 NaI-> AddElement(elI, natoms=1);
158 density= 4.51 *g/cm3;
159 G4Material* CsI=
new G4Material(
"CsI", density, nel= 2);
160 CsI-> AddElement(elCs, natoms=1);
161 CsI-> AddElement(elI, natoms=1);