34 #include "SensitiveDetector.hh"
35 #include "Materials.hh"
55 : fRadiatorDescription(0)
82 G4double foilGasRatio = radThickness/(radThickness+gasGap);
109 G4double totDensity = foilDensity*foilGasRatio
110 + gasDensity*(1.0-foilGasRatio);
112 G4double fractionFoil = foilDensity*foilGasRatio/totDensity;
113 G4double fractionGas = gasDensity*(1.0-foilGasRatio)/totDensity;
115 radiatorMat->AddMaterial(ch2, fractionFoil);
116 radiatorMat->AddMaterial(air, fractionGas);
127 foilGasRatio = 0.99999;
128 foilDensity = 1.2928*
mg/
cm3;
129 gasDensity = 0.178*
mg/
cm3;
130 totDensity = foilDensity*foilGasRatio + gasDensity*(1.0-foilGasRatio);
132 fractionFoil = foilDensity*foilGasRatio/totDensity;
133 fractionGas = gasDensity*(1.0-foilGasRatio)/totDensity;
136 pipeMat->AddMaterial(air, fractionFoil);
137 pipeMat->AddMaterial(he, fractionGas);
147 =
new G4Box(
"World", worldSizeR, worldSizeR, worldSizeZ/2.);
157 G4double radThick = foilNumber*(radThickness + gasGap) - gasGap + detGap;
158 G4double radZ = startZ + 0.5*radThick;
161 =
new G4Box(
"Radiator", 1.1*absorberRadius, 1.1*absorberRadius, 0.5*radThick);
167 "Radiator", logicRadiator, physicsWorld,
false, 0 );
179 G4double zElectrode1 = radZ - radThick/2. - electrodeThick/2.;
180 G4double zElectrode2 = radZ + radThick/2. + electrodeThick/2.;
184 G4cout <<
"electrodeThick = " << electrodeThick/
mm <<
" mm" << G4endl <<
G4endl;
190 G4double zPipe = zElectrode2 + electrodeThick/2. + pipeLength/2. + pipeDist/2.;
193 G4cout <<
"pipeLength = " << pipeLength/
mm <<
" mm" << G4endl <<
G4endl;
198 G4double zMylar1 = zPipe - pipeLength/2. - mylarThick/2 - 0.01*
mm;
199 G4double zMylar2 = zPipe + pipeLength/2. + mylarThick/2 + 0.01*
mm;
203 G4cout <<
"fMylarThick = " << mylarThick/
mm <<
" mm" << G4endl <<
G4endl;
208 G4double zMylar = zElectrode2 + electrodeThick/2. + mylarThick/2. + 1.0*
mm;
209 zMylar += ( pipeLength + pipeDist );
212 G4cout <<
"mylarThick = " << mylarThick/
mm <<
" mm" << G4endl <<
G4endl;
216 G4double absorberZ = zMylar + mylarThick/2. + absorberThickness/2.;
219 =
new G4Box(
"Absorber", 10.*
mm, 10.*
mm, absorberThickness/2.);
225 "Absorber", logicAbsorber, physicsWorld,
false, 0);
240 G4cout <<
"\n The WORLD is made of "
241 << worldSizeZ/
mm <<
"mm of " << worldMaterial->
GetName();
242 G4cout <<
", the transverse size (R) of the world is "
243 << worldSizeR/
mm <<
" mm. " <<
G4endl;
244 G4cout <<
" The ABSORBER is made of "
245 << absorberThickness/
mm <<
"mm of " << absorberMaterial->
GetName();
246 G4cout <<
", the transverse size (R) is "
247 << absorberRadius/
mm <<
" mm. " <<
G4endl;
248 G4cout <<
" Z position of the (middle of the) absorber "
249 << absorberZ/
mm <<
" mm." <<
G4endl;
256 G4cout <<
"fRadiatorMat = " << radiatorMat->GetName() <<
G4endl;