70 fLogicalVolumeVector(),fPhysicalVolumeVector()
77 fLogicalVolumeVector.clear();
78 fPhysicalVolumeVector.clear();
129 temperature = 2.73*
kelvin;
134 density = 2.03*
g/
cm3;
136 Concrete->
AddElement(elH , fractionmass= 0.01);
137 Concrete->
AddElement(elO , fractionmass= 0.529);
138 Concrete->
AddElement(elNa , fractionmass= 0.016);
139 Concrete->
AddElement(elHg , fractionmass= 0.002);
140 Concrete->
AddElement(elAl , fractionmass= 0.034);
141 Concrete->
AddElement(elSi , fractionmass= 0.337);
142 Concrete->
AddElement(elK , fractionmass= 0.013);
143 Concrete->
AddElement(elCa , fractionmass= 0.044);
144 Concrete->
AddElement(elFe , fractionmass= 0.014);
145 Concrete->
AddElement(elC , fractionmass= 0.001);
164 spanningAngleCylinder);
170 fLogicalVolumeVector.push_back(worldCylinder_log);
172 name =
"shieldWorld";
178 fPhysicalVolumeVector.push_back(fWorldVolume);
193 spanningAngleShield);
199 fLogicalVolumeVector.push_back(aShield_log);
209 for (i=1; i<=18; i++)
214 pos_z = startz + (i-1) * (2*heightShield);
223 fPhysicalVolumeVector.push_back(pvol);
230 innerRadiusShield = 0*
cm;
231 outerRadiusShield = 100*
cm;
233 startAngleShield = 0*
deg;
234 spanningAngleShield = 360*
deg;
241 spanningAngleShield);
245 fLogicalVolumeVector.push_back(aRest_log);
260 fPhysicalVolumeVector.push_back(pvol_rest);
270 G4cout <<
" B01DetectorConstruction:: Creating Importance Store " <<
G4endl;
271 if (!fPhysicalVolumeVector.size())
273 G4Exception(
"B01DetectorConstruction::CreateImportanceStore"
275 ,
"no physical volumes created yet!");
278 fWorldVolume = fPhysicalVolumeVector[0];
287 for (std::vector<G4VPhysicalVolume *>::iterator
288 it = fPhysicalVolumeVector.begin();
289 it != fPhysicalVolumeVector.end() - 1; it++)
291 if (*it != fWorldVolume)
293 imp = std::pow(2., n++);
294 G4cout <<
"Going to assign importance: " << imp <<
", to volume: "
295 << (*it)->GetName() <<
G4endl;
304 *(fPhysicalVolumeVector[fPhysicalVolumeVector.size()-1]),++n);
313 if (!fPhysicalVolumeVector.size())
315 G4Exception(
"B01DetectorConstruction::CreateWeightWindowStore"
317 ,
"no physical volumes created yet!");
320 fWorldVolume = fPhysicalVolumeVector[0];
328 std::set<G4double, std::less<G4double> > enBounds;
329 enBounds.insert(1 *
GeV);
334 std::vector<G4double> lowerWeights;
336 lowerWeights.push_back(1);
340 for (std::vector<G4VPhysicalVolume *>::iterator
341 it = fPhysicalVolumeVector.begin();
342 it != fPhysicalVolumeVector.end() - 1; it++)
344 if (*it != fWorldVolume)
346 lowerWeight = 1./std::pow(2., n++);
347 G4cout <<
"Going to assign lower weight: " << lowerWeight
349 << (*it)->GetName() <<
G4endl;
351 lowerWeights.clear();
352 lowerWeights.push_back(lowerWeight);
361 gRestCell(*(fPhysicalVolumeVector[fPhysicalVolumeVector.size()-1]), ++n);
371 std::ostringstream os;
415 G4String concreteSDname =
"ConcreteSD";
434 for (std::vector<G4LogicalVolume *>::iterator it =
435 fLogicalVolumeVector.begin();
436 it != fLogicalVolumeVector.end(); it++){
G4bool RegisterPrimitive(G4VPrimitiveScorer *)
G4VPhysicalVolume * GetWorldVolume()
CLHEP::Hep3Vector G4ThreeVector
int universe_mean_density
void AddLowerWeights(const G4GeometryCell &gCell, const std::vector< G4double > &lowerWeights)
virtual G4VPhysicalVolume * Construct()
void SetForceSolid(G4bool)
Definition of the B01DetectorConstruction class.
B01DetectorConstruction()
void Weighted(G4bool flg=true)
void Weighted(G4bool flg=true)
G4GLOB_DLL std::ostream G4cout
G4VWeightWindowStore * CreateWeightWindowStore()
void AddImportanceGeometryCell(G4double importance, const G4GeometryCell &gCell)
~B01DetectorConstruction()
void SetFilter(G4VSDFilter *value)
static G4WeightWindowStore * GetInstance()
void SetSensitiveDetector(const G4String &logVolName, G4VSensitiveDetector *aSD, G4bool multi=false)
G4String GetCellName(G4int i)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
void AddNewDetector(G4VSensitiveDetector *aSD)
G4VIStore * CreateImportanceStore()
static G4IStore * GetInstance()
static G4SDManager * GetSDMpointer()
void MultiplyKineticEnergy(G4bool flg=true)
void SetGeneralUpperEnergyBounds(const std::set< G4double, std::less< G4double > > &enBounds)
virtual void ConstructSDandField()
void AddElement(G4Element *element, G4int nAtoms)
void DivideByVelocity(G4bool flg=true)
void SetVisAttributes(const G4VisAttributes *pVA)