89 G4double
a, z, density;
92 G4double fractionmass, pressure, temperature;
99 auto elH =
new G4Element(name=
"Hydrogen",symbol=
"H" , z= 1.,
a);
102 auto elC =
new G4Element(name=
"Carbon", symbol=
"C", z=6.,
a);
105 auto elN =
new G4Element(name=
"Nitrogen",symbol=
"N" , z= 7.,
a);
108 auto elO =
new G4Element(name=
"Oxygen" ,symbol=
"O" , z= 8.,
a);
111 auto elAr =
new G4Element(name=
"Argon", symbol=
"Ar", z=18.,
a);
119 density = 1.39*g/cm3;
120 auto mylar =
new G4Material(name=
"Mylar", density, nel=3);
121 mylar->AddElement(elO,2);
122 mylar->AddElement(elC,5);
123 mylar->AddElement(elH,4);
127 auto CH2 =
new G4Material (
"Polypropelene" , 0.91*g/cm3, 2);
128 CH2->AddElement(elH,2);
129 CH2->AddElement(elC,1);
133 density = 3.700*mg/cm3;
135 auto Kr =
new G4Material(name=
"Kr",z=36.,
a, density );
139 density = 1.7836*mg/cm3;
140 auto argon =
new G4Material(name=
"Argon" , density, ncomponents=1);
141 argon->AddElement(elAr, 1);
143 density = 1.25053*mg/cm3;
144 auto nitrogen =
new G4Material(name=
"N2" , density, ncomponents=1);
145 nitrogen->AddElement(elN, 2);
147 density = 1.4289*mg/cm3;
148 auto oxygen =
new G4Material(name=
"O2" , density, ncomponents=1);
149 oxygen->AddElement(elO, 2);
151 density = 1.2928*mg/cm3;
153 temperature = STP_Temperature;
154 pressure = 1.0e-8*STP_Pressure;
156 auto air =
new G4Material(name=
"Air" , density, ncomponents=3,
157 kStateGas,temperature,pressure);
158 air->AddMaterial( nitrogen, fractionmass = 0.7557 );
159 air->AddMaterial( oxygen, fractionmass = 0.2315 );
160 air->AddMaterial( argon, fractionmass = 0.0128 );
164 density = 5.858*mg/cm3;
166 auto Xe =
new G4Material(name=
"Xenon",z=54.,
a, density );
170 density = 1.842*mg/cm3;
171 auto CarbonDioxide =
new G4Material(name=
"CO2", density, nel=2);
172 CarbonDioxide->AddElement(elC,1);
173 CarbonDioxide->AddElement(elO,2);
177 density = 5.0818*mg/cm3;
178 auto Xe20CO2 =
new G4Material(name=
"Xe20CO2", density, ncomponents=2);
179 Xe20CO2->AddMaterial( Xe, fractionmass = 0.922 );
180 Xe20CO2->AddMaterial( CarbonDioxide, fractionmass = 0.078 );
184 density = 3.601*mg/cm3;
185 auto Kr20CO2 =
new G4Material(name=
"Kr20CO2", density, ncomponents=2);
186 Kr20CO2->AddMaterial( Kr, fractionmass = 0.89 );
187 Kr20CO2->AddMaterial( CarbonDioxide, fractionmass = 0.11 );
189 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
208 G4GeometryManager::GetInstance()->OpenGeometry();
209 G4PhysicalVolumeStore::GetInstance()->Clean();
210 G4LogicalVolumeStore::GetInstance()->Clean();
211 G4SolidStore::GetInstance()->Clean();
219 G4bool checkOverlaps =
true;
241 G4cout <<
"zRad = " << zRad/mm <<
" mm" << G4endl;
242 G4cout <<
"radThick = " << radThick/mm <<
" mm" << G4endl;
243 G4cout <<
"fFoilNumber = " <<
fFoilNumber << G4endl;
244 G4cout <<
"fRadiatorMat = " <<
fRadiatorMat->GetName() << G4endl;
245 G4cout <<
"WorldMaterial = " <<
fWorldMaterial->GetName() << G4endl;
265 G4double zStart = 0.5*(-radThick + radSliceThick) +
fDetGap;
270 G4double zSlice = zStart + j*radSliceThick;
271 G4cout << zSlice/mm <<
" mm" <<
"\t";