Belle II Software development
Public Member Functions | Private Member Functions | List of all members
GeoHeavyMetalShieldCreator Class Reference

The creator for the HeavyMetalShield geometry of the Belle II detector. More...

#include <GeoHeavyMetalShieldCreator.h>

Inheritance diagram for GeoHeavyMetalShieldCreator:
CreatorBase

Public Member Functions

virtual void create (const GearDir &content, G4LogicalVolume &topVolume, geometry::GeometryTypes type) override
 The old create member: create the configuration object(s) on the fly and call the geometry creation routine.
 
virtual void createPayloads (const GearDir &content, const IntervalOfValidity &iov) override
 Create the configuration objects and save them in the Database.
 
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.
 

Private Member Functions

HeavyMetalShieldGeometryPar createConfiguration (const GearDir &param)
 Create a parameter object from the Gearbox XML parameters.
 
void createGeometry (const HeavyMetalShieldGeometryPar &parameters, G4LogicalVolume &topVolume, geometry::GeometryTypes type)
 Create the geometry from a parameter object.
 

Detailed Description

The creator for the HeavyMetalShield geometry of the Belle II detector.


Definition at line 31 of file GeoHeavyMetalShieldCreator.h.

Member Function Documentation

◆ create()

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

The old create member: create the configuration object(s) on the fly and call the geometry creation routine.

Implements CreatorBase.

Definition at line 42 of file GeoHeavyMetalShieldCreator.h.

43 {
44 HeavyMetalShieldGeometryPar config = createConfiguration(content);
45 createGeometry(config, topVolume, type);
46 }
Belle2::VXD::GeoHeavyMetalShieldCreator::createConfiguration
HeavyMetalShieldGeometryPar createConfiguration(const GearDir &param)
Create a parameter object from the Gearbox XML parameters.
Definition: GeoHeavyMetalShieldCreator.cc:50
Belle2::VXD::GeoHeavyMetalShieldCreator::createGeometry
void createGeometry(const HeavyMetalShieldGeometryPar &parameters, G4LogicalVolume &topVolume, geometry::GeometryTypes type)
Create the geometry from a parameter object.
Definition: GeoHeavyMetalShieldCreator.cc:82

◆ createConfiguration()

HeavyMetalShieldGeometryPar createConfiguration ( const GearDir &  param)
private

Create a parameter object from the Gearbox XML parameters.

Definition at line 50 of file GeoHeavyMetalShieldCreator.cc.

51 {
52 HeavyMetalShieldGeometryPar heavyMetalShieldGeometryPar;
53 //Read the definition of all shields
54 for (const GearDir& shield : param.getNodes("Shield")) {
55 VXDPolyConePar shieldPar(
56 shield.getString("@name"),
57 shield.getString("Material", "Air"),
58 shield.getAngle("minPhi", 0),
59 shield.getAngle("maxPhi", 2 * M_PI),
60 (shield.getNodes("Cutout").size() > 0),
61 //shield.getLength("Cutout/width", 0.),
62 shield.getLength("Cutout/width1", 0.),
63 shield.getLength("Cutout/width2", 0.),
64 shield.getLength("Cutout/height", 0.),
65 shield.getLength("Cutout/depth", 0.)
66 );
67
68 for (const GearDir& plane : shield.getNodes("Plane")) {
69 VXDPolyConePlanePar planePar(
70 plane.getLength("posZ"),
71 plane.getLength("innerRadius"),
72 plane.getLength("outerRadius")
73 );
74 shieldPar.getPlanes().push_back(planePar);
75 }
76
77 heavyMetalShieldGeometryPar.getShields().push_back(shieldPar);
78 }
79 return heavyMetalShieldGeometryPar;
80 }

◆ 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 58 of file GeoHeavyMetalShieldCreator.h.

59 {
60 DBObjPtr<HeavyMetalShieldGeometryPar> dbObj;
61 if (!dbObj) {
62 // Check that we found the object and if not report the problem
63 B2FATAL("No configuration for " << name << " found.");
64 }
65 createGeometry(*dbObj, topVolume, type);
66 }

◆ createGeometry()

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

Create the geometry from a parameter object.

Definition at line 82 of file GeoHeavyMetalShieldCreator.cc.

84 {
85
86 // Create the shields
87 for (const VXDPolyConePar& shield : parameters.getShields()) {
88
89 string name = shield.getName();
90 double minZ(0), maxZ(0);
91
92 // Create a polycone
93 double minPhi = shield.getMinPhi();
94 double dPhi = shield.getMaxPhi() - minPhi;
95 int nPlanes = shield.getPlanes().size();
96 if (nPlanes < 2) {
97 B2ERROR("Polycone needs at least two planes");
98 return ;
99 }
100 std::vector<double> z(nPlanes, 0);
101 std::vector<double> rMin(nPlanes, 0);
102 std::vector<double> rMax(nPlanes, 0);
103 int index(0);
104 minZ = numeric_limits<double>::infinity();
105 maxZ = -numeric_limits<double>::infinity();
106
107
108 for (const VXDPolyConePlanePar& plane : shield.getPlanes()) {
109 z[index] = plane.getPosZ() / Unit::mm;
110 minZ = min(minZ, z[index]);
111 maxZ = max(maxZ, z[index]);
112 rMin[index] = plane.getInnerRadius() / Unit::mm;
113 rMax[index] = plane.getOuterRadius() / Unit::mm;
114 ++index;
115 }
116
117 G4VSolid* geoShield = new G4Polycone(name + " IR Shield", minPhi, dPhi, nPlanes, z.data(), rMin.data(), rMax.data());
118
119 // Cutouts (if present)
120 if (shield.getDoCutOut()) {
121 //double sizeX = shield.getCutOutWidth() / Unit::mm / 2.;
122 double sizeX1 = shield.getCutOutWidth1() / Unit::mm / 2.;
123 double sizeX2 = shield.getCutOutWidth2() / Unit::mm / 2.;
124 double sizeY = shield.getCutOutHeight() / Unit::mm / 2.;
125 double depth2 = shield.getCutOutDepth() / Unit::mm / 2.;
126 double sizeZ = (maxZ - minZ) / 2.;
127 double sign = math::sign<double>(minZ);
128 double minAbsZ = min(fabs(minZ), fabs(maxZ));
129
130 G4ThreeVector origin1(0, 0, sign * (minAbsZ + sizeZ));
131 G4ThreeVector origin2(0, 0, sign * (minAbsZ + depth2));
132
133 //G4Box* box1 = new G4Box("Cutout", sizeX, sizeY, sizeZ);
134 G4Trd* box1 = new G4Trd("Cutout1", sizeX1, sizeX2, sizeY, sizeY, sizeZ);
135 G4Box* box2 = new G4Box("Cutout2", 100 / Unit::mm, sizeY, depth2);
136
137 geoShield = new G4SubtractionSolid(name + " IR Shield", geoShield, box1, G4Translate3D(origin1));
138 geoShield = new G4SubtractionSolid(name + " IR Shield", geoShield, box2, G4Translate3D(origin2));
139 }
140
141 string materialName = shield.getMaterial();
142 G4Material* material = Materials::get(materialName);
143 if (!material) B2FATAL("Material '" << materialName << "', required by " << name << " IR Shield could not be found");
144
145 G4LogicalVolume* volume = new G4LogicalVolume(geoShield, material, name + " IR Shield");
146 setColor(*volume, "#cc0000");
147 //setVisibility(*volume, false);
148 new G4PVPlacement(0, G4ThreeVector(0, 0, 0), volume, name + " IR Shield", &topVolume, false, 0);
149 }
150 }
Belle2::Unit::mm
static const double mm
[millimeters]
Definition: Unit.h:70
Belle2::geometry::Materials::get
static G4Material * get(const std::string &name)
Find given material.
Definition: Materials.h:63
Belle2::geometry::setColor
void setColor(G4LogicalVolume &volume, const std::string &color)
Set the color of a logical volume.
Definition: utilities.cc:100

◆ createPayloads()

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

Create the configuration objects and save them in the Database.

If more than one object is needed adjust accordingly

Reimplemented from CreatorBase.

Definition at line 50 of file GeoHeavyMetalShieldCreator.h.

51 {
52 DBImportObjPtr<HeavyMetalShieldGeometryPar> importObj;
53 importObj.construct(createConfiguration(content));
54 importObj.import(iov);
55 }
The documentation for this class was generated from the following files: