Belle II Software development
Public Member Functions | Protected Attributes | Private Member Functions | List of all members
BeamabortCreator Class Reference

The creator for the BEAMABORT geometry. More...

#include <BeamabortCreator.h>

Inheritance diagram for BeamabortCreator:
CreatorBase

Public Member Functions

 BeamabortCreator ()
 Constructor.
 
virtual ~BeamabortCreator ()
 Destructor.
 
virtual void create (const GearDir &content, G4LogicalVolume &topVolume, geometry::GeometryTypes type) override
 Creates the ROOT Objects for the BeamPipe geometry.
 
virtual void createPayloads (const GearDir &content, const IntervalOfValidity &iov) override
 creates DB payload for BeamabortGeo class
 
virtual void createFromDB (const std::string &name, G4LogicalVolume &topVolume, geometry::GeometryTypes type) override
 Create the geometry from the Database.
 
 BELLE2_DEFINE_EXCEPTION (DBNotImplemented, "Cannot create geometry from Database.")
 Exception that will be thrown in createFromDB if member is not yet implemented by creator.
 

Protected Attributes

BeamabortGeo m_config
 geometry parameters object
 

Private Member Functions

BeamabortGeo createConfiguration (const GearDir &param)
 Reads beamabort geometry parameters from the xml files and createst DB class BeamabortGeo.
 
void createGeometry (G4LogicalVolume &topVolume, geometry::GeometryTypes type)
 Create detector geometry.
 

Detailed Description

The creator for the BEAMABORT geometry.

Definition at line 32 of file BeamabortCreator.h.

Constructor & Destructor Documentation

◆ BeamabortCreator()

BeamabortCreator ( )

Constructor.

Definition at line 52 of file BeamabortCreator.cc.

53 {}

◆ ~BeamabortCreator()

~BeamabortCreator ( )
virtual

Destructor.

Definition at line 55 of file BeamabortCreator.cc.

56 {}

Member Function Documentation

◆ create()

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

Creates the ROOT Objects for the BeamPipe geometry.

Parameters
contentA reference to the content part of the parameter description, which should to be used to create the ROOT objects.
topVolumeGeant4 logical top volume.
typeGeometry type.

Implements CreatorBase.

Definition at line 67 of file BeamabortCreator.h.

68 {
69 m_config = createConfiguration(content);
70
71 // override geometry configuration from the DB
72 DBStore::Instance().addConstantOverride("BeamabortGeo", new BeamabortGeo(m_config));
73
74 createGeometry(topVolume, type);
75 }
Belle2::beamabort::BeamabortCreator::createGeometry
void createGeometry(G4LogicalVolume &topVolume, geometry::GeometryTypes type)
Create detector geometry.
Definition: BeamabortCreator.cc:58
Belle2::beamabort::BeamabortCreator::m_config
BeamabortGeo m_config
geometry parameters object
Definition: BeamabortCreator.h:98
Belle2::beamabort::BeamabortCreator::createConfiguration
BeamabortGeo createConfiguration(const GearDir &param)
Reads beamabort geometry parameters from the xml files and createst DB class BeamabortGeo.
Definition: BeamabortCreator.h:35
Belle2::DBStore::Instance
static DBStore & Instance()
Instance of a singleton DBStore.
Definition: DBStore.cc:28
Belle2::DBStore::addConstantOverride
void addConstantOverride(const std::string &name, TObject *obj, bool oneRun=false)
Add constant override payload.
Definition: DBStore.cc:204

◆ createConfiguration()

BeamabortGeo createConfiguration ( const GearDir &  param)
inlineprivate

Reads beamabort geometry parameters from the xml files and createst DB class BeamabortGeo.

Definition at line 35 of file BeamabortCreator.h.

36 {
37 BeamabortGeo BeamabortGeoConfig;
38 BeamabortGeoConfig.initialize(param);
39 return BeamabortGeoConfig;
40 }

◆ createFromDB()

virtual void createFromDB ( const std::string &  name,
G4LogicalVolume &  topVolume,
geometry::GeometryTypes  type 
)
inlineoverridevirtual

Create the geometry from the Database.

Reimplemented from CreatorBase.

Definition at line 86 of file BeamabortCreator.h.

87 {
88 DBObjPtr<BeamabortGeo> dbObj;
89 if (!dbObj) {
90 // Check that we found the object and if not report the problem
91 B2FATAL("No configuration for " << name << " found.");
92 }
93 m_config = *dbObj;
94 createGeometry(topVolume, type);
95 }

◆ createGeometry()

void createGeometry ( G4LogicalVolume &  topVolume,
geometry::GeometryTypes  type 
)
private

Create detector geometry.

Parameters
topVolumeGeant4 logical top volume.
typeGeometry type.

Definition at line 58 of file BeamabortCreator.cc.

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

◆ createPayloads()

virtual void createPayloads ( const GearDir &  content,
const IntervalOfValidity &  iov 
)
inlineoverridevirtual

creates DB payload for BeamabortGeo class

Reimplemented from CreatorBase.

Definition at line 78 of file BeamabortCreator.h.

79 {
80 DBImportObjPtr<BeamabortGeo> importObj;
81 importObj.construct(createConfiguration(content));
82 importObj.import(iov);
83 }

Member Data Documentation

◆ m_config

BeamabortGeo m_config
protected

geometry parameters object

Definition at line 98 of file BeamabortCreator.h.

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