17 #include "GeoModelKernel/GeoMaterial.h"
22 #include "G4LogicalVolume.hh"
23 #include "G4Material.hh"
24 #include "G4PVPlacement.hh"
25 #include "G4ThreeVector.hh"
28 #include "CLHEP/Units/SystemOfUnits.h"
43 <<
" side muon wall: "
49 const double zScintillator =
m_zLength / 2.;
50 const double yScintillator =
m_yLength / 2.;
51 const double xScintillator =
m_xLength / 2.;
52 const double dzmuonwall = 750. *
CLHEP::mm;
53 const double dymuonwall = 425. *
CLHEP::mm;
54 const double dxmuonwall = xScintillator;
56 G4Box *wall =
new G4Box(
"MuonWall", dxmuonwall, dymuonwall, dzmuonwall);
60 if (StatusCode::SUCCESS !=
detStore()->
retrieve(materialManager, std::string(
"MATERIALS"))) {
65 const GeoMaterial *geoAir = materialManager->
getMaterial(
"tile::Air");
66 const GeoMaterial *geoScintillator = materialManager->
getMaterial(
"tile::Scintillator");
68 G4Material *airMaterial = theMaterialFactory.
Build(geoAir);
69 G4Material *scintillatorMaterial = theMaterialFactory.
Build(geoScintillator);
71 G4LogicalVolume *wallV =
new G4LogicalVolume(wall, airMaterial,
"MuonWall");
74 double zLayer, yLayer, xLayer = 0.;
78 const int nrOfLayers = 12;
80 for (
int j = 0; j < nrOfLayers; j++) {
83 G4Box *scintillatorLayer =
new G4Box(
"ScintillatorLayer", xScintillator, yScintillator, zScintillator);
85 G4LogicalVolume *scintillatorLayerV =
new G4LogicalVolume(scintillatorLayer, scintillatorMaterial,
86 "MuScintillatorLayer");
92 zLayer = -450. *
CLHEP::mm + (2 * j + 1) * zScintillator;
94 yLayer = -yScintillator + 25. *
CLHEP::mm;
95 zLayer = -450. *
CLHEP::mm + (2 * (j - 6) + 1) * zScintillator;
100 G4ThreeVector(xLayer,yLayer,zLayer),
102 "MuScintillatorLayer",
112 const double yScintillator1 = yScintillator + 50. *
CLHEP::mm;
113 const double zScintillator1 = zScintillator - 50. *
CLHEP::mm;
115 const int nScintLayers = 3;
117 for (
int j = 0; j < nScintLayers; j++) {
119 G4Box* scintillatorLayer =
new G4Box(
"SideScintiLayer", xScintillator, yScintillator1, zScintillator1);
121 G4LogicalVolume* scintillatorLayerV =
new G4LogicalVolume(scintillatorLayer, scintillatorMaterial,
122 "SideMuScintiLayer");
125 zLayer = -750. *
CLHEP::mm + (2 * j + 1) * zScintillator1;
129 G4ThreeVector(xLayer,yLayer,zLayer),