Belle II Software  release-06-01-15
Public Member Functions | Protected Attributes | List of all members
BeamabortCreator Class Reference

The creator for the BEAMABORT geometry. More...

#include <BeamabortCreator.h>

Inheritance diagram for BeamabortCreator:
Inheritance graph
[legend]
Collaboration diagram for BeamabortCreator:
Collaboration graph
[legend]

Public Member Functions

 BeamabortCreator ()
 Constructor.
 
virtual ~BeamabortCreator ()
 Destructor.
 
virtual void create (const GearDir &content, G4LogicalVolume &topVolume, geometry::GeometryTypes type)
 Creation of the detector geometry from Gearbox (XML). More...
 
 BELLE2_DEFINE_EXCEPTION (DBNotImplemented, "Cannot create geometry from Database.")
 Exception that will be thrown in createFromDB if member is not yet implemented by creator.
 
virtual void createFromDB (const std::string &name, G4LogicalVolume &topVolume, GeometryTypes type)
 Function to create the geometry from the Database. More...
 
virtual void createPayloads (const GearDir &content, const IntervalOfValidity &iov)
 Function to create the geometry database. More...
 

Protected Attributes

SensitiveDetector * m_sensitive
 SensitiveDetector BEAMABORT.
 

Detailed Description

The creator for the BEAMABORT geometry.

Definition at line 27 of file BeamabortCreator.h.

Member Function Documentation

◆ create()

void create ( const GearDir &  content,
G4LogicalVolume &  topVolume,
geometry::GeometryTypes  type 
)
virtual

Creation of the detector geometry from Gearbox (XML).

Parameters
[in]contentXML data directory.
[in]topVolumeGeant world volume.
[in]typeGeometry type.

Implements CreatorBase.

Definition at line 60 of file BeamabortCreator.cc.

62  {
63 
64  m_sensitive = new SensitiveDetector();
65 
66  B2INFO("BeamabortCreator phase2");
67  //Visualization Attributes
68  G4VisAttributes* orange = new G4VisAttributes(G4Colour(1, 2, 0));
69  orange->SetForceAuxEdgeVisible(true);
70  G4VisAttributes* magenta = new G4VisAttributes(G4Colour(1, 0, 1));
71  magenta->SetForceAuxEdgeVisible(true);
72  //lets get the stepsize parameter with a default value of 5 µm
73  double stepSize = content.getLength("stepSize", 5 * CLHEP::um);
74 
75  //Lets loop over all the Active nodes
76  BOOST_FOREACH(
77  const GearDir & activeParams, content.getNodes("Active")) {
78 
79  int phase = activeParams.getInt("phase");
80  //Positioned PIN diodes
81  std::vector<double> x_pos;
82  std::vector<double> y_pos;
83  std::vector<double> z_pos;
84  std::vector<double> phi;
85  std::vector<double> thetaX;
86  std::vector<double> thetaY;
87  std::vector<double> thetaZ;
88  std::vector<double> svdAngle;
89  std::vector<double> r;
90  std::vector<double> deltaX;
91  unsigned int dimz = 0;
92 
93  for (double z : activeParams.getArray("z", {0})) {
94  z *= CLHEP::cm;
95  z_pos.push_back(z);
96  dimz++;
97  }
98 
99  for (double ThetaZ : activeParams.getArray("ThetaZ", {0})) {
100  thetaZ.push_back(ThetaZ);
101  }
102 
103  if (thetaZ.size() != dimz) { B2ERROR("Diamond data not consistent (i.e. not same number of all position parmeters)"); return;}
104 
105  if (phase == 2 || phase == 3) {
106 
107  for (double Phi : activeParams.getArray("Phi", {0})) {
108  phi.push_back(Phi);
109  }
110 
111  for (double r_dia : activeParams.getArray("r_dia", {0})) {
112  r_dia *= CLHEP::cm;
113  r.push_back(r_dia);
114  }
115  for (double dX : activeParams.getArray("deltaX", {0})) {
116  dX *= CLHEP::cm;
117  deltaX.push_back(dX);
118  }
119  for (double addAngle : activeParams.getArray("addAngle", {0})) {
120  svdAngle.push_back(addAngle);
121  }
122  if (phi.size() != dimz || r.size() != dimz || deltaX.size() != dimz || svdAngle.size() != dimz) { B2ERROR("Diamond data not consistent (i.e. not same number of all position parmeters)"); return;}
123  }
124  if (phase == 1) {
125  for (double x : activeParams.getArray("x", {0})) {
126  x *= CLHEP::cm;
127  x_pos.push_back(x);
128  }
129  for (double y : activeParams.getArray("y", {0})) {
130  y *= CLHEP::cm;
131  y_pos.push_back(y);
132  }
133 
134  for (double ThetaX : activeParams.getArray("ThetaX", {0})) {
135  thetaX.push_back(ThetaX);
136  }
137  for (double ThetaY : activeParams.getArray("ThetaY", {0})) {
138  thetaY.push_back(ThetaY);
139  }
140  if (x_pos.size() != dimz || y_pos.size() != dimz || thetaX.size() != dimz || thetaY.size() != dimz) { B2ERROR("Diamond data not consistent (i.e. not same number of all position parmeters)"); return;}
141  }
142 
143 
144  //create beamabort package
145  G4double dx_opa = 12. / 2. * CLHEP::mm;
146  G4double dy_opa = 18. / 2. * CLHEP::mm;
147  G4double dz_opa = 3.1 / 2. * CLHEP::mm;
148  G4VSolid* s_pa = new G4Box("s_opa", dx_opa, dy_opa, dz_opa);
149  G4double dx_ipa = 6. / 2. * CLHEP::mm;
150  G4double dy_ipa = 6. / 2. * CLHEP::mm;
151  G4double dz_ipa = 5.6 / 2. * CLHEP::mm;
152  G4VSolid* s_ipa = new G4Box("s_ipa", dx_ipa, dy_ipa, dz_ipa);
153  s_pa = new G4SubtractionSolid("s_pa", s_pa, s_ipa, 0, G4ThreeVector(0., 4.0, 0.));
154  G4LogicalVolume* l_pa = new G4LogicalVolume(s_pa, G4Material::GetMaterial("Al6061"), "l_pa");
155  l_pa->SetVisAttributes(magenta);
156  G4Transform3D transform;
157  for (unsigned int i = 0; i < dimz; i++) {
158  B2INFO("DIA-" << i << "RotateZ3D phi: " << phi[i]);
159 
160  if (phase == 1) {
161  transform = G4Translate3D(x_pos[i], y_pos[i], z_pos[i]) * G4RotateX3D(thetaX[i]) *
162  G4RotateY3D(thetaY[i]) * G4RotateZ3D(thetaZ[i]);
163  } else if (phase == 2 || phase == 3) {
164  transform = G4Translate3D(deltaX[i], 0, 0) * G4RotateZ3D(phi[i]) *
165  G4Translate3D(r[i], 0, z_pos[i]) * G4RotateX3D(-M_PI / 2 - thetaZ[i]) *
166  G4RotateZ3D(svdAngle[i]) * G4RotateY3D(M_PI / 2);
167  }
168  new G4PVPlacement(transform, l_pa, TString::Format("p_dia_pa_%d", i).Data(), &topVolume, false, 0);
169  B2INFO("DIA-" << i << " placed at: " << transform.getTranslation() << " mm ");
170  }
171 
172  //create beamabort volumes
173  G4double dx_ba = 4.5 / 2. * CLHEP::mm;
174  G4double dy_ba = 4.5 / 2. * CLHEP::mm;
175  G4double dz_ba = 0.5 / 2. * CLHEP::mm;
176  G4Box* s_BEAMABORT = new G4Box("s_BEAMABORT", dx_ba, dy_ba, dz_ba);
177  G4LogicalVolume* l_BEAMABORT = new G4LogicalVolume(s_BEAMABORT, geometry::Materials::get("Diamond"),
178  "l_BEAMABORT", 0, m_sensitive);
179  l_BEAMABORT->SetVisAttributes(orange);
180 
181  //Lets limit the Geant4 stepsize inside the volume
182  l_BEAMABORT->SetUserLimits(new G4UserLimits(stepSize));
183  l_BEAMABORT->SetVisAttributes(orange);
184 
185  for (unsigned int i = 0; i < dimz; i++) {
186 
187  if (phase == 1)
188  transform = G4Translate3D(x_pos[i], y_pos[i],
189  z_pos[i]) * G4RotateX3D(thetaX[i]) * G4RotateY3D(thetaY[i]) * G4RotateZ3D(thetaZ[i]);
190  if (phase == 2 || phase == 3)
191  transform = G4Translate3D(deltaX[i], 0, 0) * G4RotateZ3D(phi[i]) *
192  G4Translate3D(r[i], 0, z_pos[i]) * G4RotateX3D(-M_PI / 2 - thetaZ[i]) *
193  G4RotateZ3D(svdAngle[i]) * G4RotateY3D(M_PI / 2) * G4Translate3D(0, 4.2, 0);
194 
195  new G4PVPlacement(transform, l_BEAMABORT, TString::Format("p_dia_%d", i).Data(), &topVolume, false, i);
196 
197  B2INFO("DIA-sensitive volume-" << i << " placed at: " << transform.getTranslation() << " mm "
198  << " at phi angle = " << phi[i] << " at theta angle = "
199  << thetaZ[i]);
200  B2INFO("DIA-sensitive volume-" << i << " G4RotateZ3D of: phi= " << phi[i] << " G4RotateX3D = "
201  << (-M_PI / 2 - thetaZ[i]));
202  }
203  }
204  }
Belle2::beamabort::BeamabortCreator::m_sensitive
SensitiveDetector * m_sensitive
SensitiveDetector BEAMABORT.
Definition: BeamabortCreator.h:47
Belle2::geometry::Materials::get
static G4Material * get(const std::string &name)
Find given material.
Definition: Materials.h:63
Belle2::PXD::SensitiveDetector
VXD::SensitiveDetector< PXDSimHit, PXDTrueHit > SensitiveDetector
The PXD Sensitive Detector class.
Definition: SensitiveDetector.h:23

◆ createFromDB()

void createFromDB ( const std::string &  name,
G4LogicalVolume &  topVolume,
GeometryTypes  type 
)
virtualinherited

Function to create the geometry from the Database.

Parameters
namename of the component in the database, could be used to disambiguate multiple components created with the same creator
topVolumeTop volume in which the geometry has to be placed
typeType of geometry to be build

Reimplemented in GeoVXDServiceCreator, GeoHeavyMetalShieldCreator, GeoTOPCreator, GeoSVDCreator, GeoSTRCreator, GeoServiceMaterialCreator, GeoCOILCreator, GeoPXDCreator, GeoKLMCreator, GeoEKLMCreator, GeoBKLMCreator, GeoFarBeamLineCreator, GeoCryostatCreator, GeoBeamPipeCreator, MyDBCreator, GeoECLCreator, GeoCDCCreator, GeoARICHCreator, and GeoMagneticField.

Definition at line 17 of file CreatorBase.cc.

◆ createPayloads()

void createPayloads ( const GearDir &  content,
const IntervalOfValidity &  iov 
)
virtualinherited

Function to create the geometry database.

This function should be implemented to convert Gearbox parameters to one ore more database payloads

Parameters
contentGearDir pointing to the parameters which should be used for construction
iovinterval of validity to use when generating payloads

Reimplemented in GeoVXDServiceCreator, GeoHeavyMetalShieldCreator, GeoTOPCreator, GeoSVDCreator, GeoSTRCreator, GeoServiceMaterialCreator, GeoCOILCreator, GeoPXDCreator, GeoKLMCreator, GeoEKLMCreator, GeoBKLMCreator, GeoFarBeamLineCreator, GeoCryostatCreator, GeoBeamPipeCreator, MyDBCreator, GeoMagneticField, GeoECLCreator, GeoCDCCreator, and GeoARICHCreator.

Definition at line 24 of file CreatorBase.cc.

The documentation for this class was generated from the following files: