61 fDetectorMessenger(0),
62 fSolidWorld(0), fLogicWorld(0), fPhysiWorld(0),
63 fSolidAbsorber(0), fLogicAbsorber(0), fPhysiAbsorber(0),
64 fAbsorberMaterial(0), fAbsorberThickness(0.), fAbsorberRadius(0.),
65 fZAbsorber(0.), fZStartAbs(0.), fZEndAbs(0.),
66 fWorldMaterial(0), fWorldSizeR(0.), fWorldSizeZ(0.)
112 G4double fractionmass, pressure, temperature;
139 density = 1.39*
g/
cm3;
153 density = 3.700*
mg/
cm3;
159 density = 1.7836*
mg/
cm3;
163 density = 1.25053*
mg/
cm3;
167 density = 1.4289*
mg/
cm3;
171 density = 1.2928*
mg/
cm3;
179 air->
AddMaterial( nitrogen, fractionmass = 0.7557 );
186 density = 5.858*
mg/
cm3;
192 density = 1.842*
mg/
cm3;
199 density = 5.0818*
mg/
cm3;
202 Xe20CO2->
AddMaterial( CarbonDioxide, fractionmass = 0.078 );
206 density = 3.601*
mg/
cm3;
209 Kr20CO2->
AddMaterial( CarbonDioxide, fractionmass = 0.11 );
281 G4cout <<
"\n The WORLD is made of "
283 G4cout <<
", the transverse size (R) of the world is "
285 G4cout <<
" The ABSORBER is made of "
289 G4cout <<
" Z position of the (middle of the) absorber "
303 for (
size_t j=0 ; j<theMaterialTable->size() ; j++)
304 { material = (*theMaterialTable)[j];
305 if (material->
GetName() == materialChoice)
323 for (
size_t j=0 ; j<theMaterialTable->size() ; j++)
324 { material = (*theMaterialTable)[j];
325 if(material->
GetName() == materialChoice)