61 fDetectorMessenger(0),
62 fSolidWorld(0), fLogicWorld(0), fPhysiWorld(0),
63 fSolidAbsorber(0), fLogicAbsorber(0), fPhysiAbsorber(0),
64 fSolidRadSlice(0), fLogicRadSlice(0), fPhysiRadSlice(0),
65 fSolidRadiator(0), fLogicRadiator(0), fPhysiRadiator(0),
66 fWorldMaterial(0), fAbsorberMaterial(0), fRadiatorMat(0),
68 fWorldSizeR( 22000.*
mm),
69 fWorldSizeZ( 44000.*
mm),
70 fAbsorberThickness( 1.*
mm),
71 fAbsorberRadius( 20000.*
mm),
72 fZAbsorber( 21990.*
mm),
75 fRadThickness( 100.*
mm),
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;
178 air->
AddMaterial( nitrogen, fractionmass = 0.7557 );
184 density = 5.858*
mg/
cm3;
190 density = 1.842*
mg/
cm3;
197 density = 5.0818*
mg/
cm3;
200 Xe20CO2->
AddMaterial( CarbonDioxide, fractionmass = 0.078 );
204 density = 3.601*
mg/
cm3;
207 Kr20CO2->
AddMaterial( CarbonDioxide, fractionmass = 0.11 );
239 G4bool checkOverlaps =
true;
290 G4double zSlice = zStart + j*radSliceThick;
291 G4cout << zSlice/
mm <<
" mm" <<
"\t";
326 G4cout <<
"\n The WORLD is made of "
328 G4cout <<
", the transverse size (R) of the world is "
330 G4cout <<
" The ABSORBER is made of "
334 G4cout <<
" Z position of the (middle of the) absorber "
348 for (
size_t j=0 ; j<theMaterialTable->size() ; j++)
349 { material = (*theMaterialTable)[j];
350 if (material->
GetName() == materialChoice)
368 for (
size_t j=0 ; j<theMaterialTable->size() ; j++)
369 { material = (*theMaterialTable)[j];
370 if(material->
GetName() == materialChoice)