9#include <vxd/geometry/GeoVXDCreator.h>
10#include <vxd/geometry/GeoCache.h>
11#include <vxd/simulation/SensitiveDetector.h>
12#include <vxd/simulation/SensitiveDetectorBase.h>
13#include <simulation/background/BkgSensitiveDetector.h>
15#include <geometry/CreatorFactory.h>
16#include <geometry/Materials.h>
17#include <geometry/utilities.h>
18#include <framework/gearbox/GearDir.h>
20#include <boost/algorithm/string.hpp>
21#include <boost/format.hpp>
23#include <G4ReflectionFactory.hh>
24#include <G4LogicalVolume.hh>
29#include <G4Polycone.hh>
30#include <G4PVPlacement.hh>
31#include <G4AssemblyVolume.hh>
32#include <G4UserLimits.hh>
33#include <G4Point3D.hh>
35#include <G4TessellatedSolid.hh>
36#include <G4QuadrangularFacet.hh>
37#include <G4TriangularFacet.hh>
39#include <CLHEP/Units/PhysicalConstants.h>
40#include <CLHEP/Units/SystemOfUnits.h>
70 for (G4UserLimits* userLimit :
m_UserLimits)
delete userLimit;
76 vector<VXDGeoPlacement> placements,
bool originCenter,
bool allowOutside)
79 B2DEBUG(100,
"Creating component " << name);
80 vector<VXDGeoComponent> subComponents;
81 subComponents.reserve(placements.size());
84 bool widthResize = component.getWidth() <= 0;
85 bool lengthResize = component.getLength() <= 0;
86 bool heightResize = component.getHeight() <= 0;
91 B2FATAL(
"A component is requested that was not created before!");
95 B2DEBUG(100,
"SubComponent " << p.getName());
96 B2DEBUG(100, boost::format(
"Placement: u:%1% cm, v:%2% cm, w:%3% + %4% cm") % p.getU() % p.getV() % p.getW() % p.getWOffset());
97 B2DEBUG(100, boost::format(
"Dimensions: %1%x%2%x%3% cm") % sub.getWidth() % sub.getLength() % sub.getHeight());
101 if (!allowOutside) B2FATAL(
"Cannot place component " << p.getName() <<
" outside of component " << name);
102 }
else if (sub.getHeight() + p.getWOffset() > component.getHeight()) {
105 B2FATAL(
"Subcomponent " << p.getName() <<
" does not fit into volume: "
106 <<
"height " << sub.getHeight() <<
" > " << component.getHeight());
108 component.getHeight() = sub.getHeight() + p.getWOffset();
112 double minWidth = max(abs(p.getU() + sub.getWidth() / 2.0), abs(p.getU() - sub.getWidth() / 2.0));
113 double minLength = max(abs(p.getV() + sub.getLength() / 2.0), abs(p.getV() - sub.getLength() / 2.0));
114 if (minWidth > component.getWidth() + component.getWidth() * numeric_limits<double>::epsilon()) {
116 B2FATAL(
"Subcomponent " << p.getName() <<
" does not fit into volume: "
117 <<
"minWidth " << minWidth <<
" > " << component.getWidth());
119 component.setWidth(minWidth * 2.0);
121 if (minLength > component.getLength() + component.getLength() * numeric_limits<double>::epsilon()) {
123 B2FATAL(
"Subcomponent " << p.getName() <<
" does not fit into volume: "
124 <<
"minLength " << minLength <<
" > " << component.getLength());
126 component.setLength(minLength * 2.0);
128 subComponents.push_back(sub);
132 if (component.getWidth() <= 0 || component.getLength() <= 0 || component.getHeight() <= 0) {
133 B2FATAL(
"At least one dimension of component " << name <<
" is zero which does not make sense");
137 if (!component.getVolume()) {
138 G4VSolid* componentShape =
createTrapezoidal(name, component.getWidth(), component.getWidth2(), component.getLength(),
139 component.getHeight());
140 component.setVolume(
new G4LogicalVolume(componentShape,
Materials::get(component.getMaterial()), name));
143 B2DEBUG(100, boost::format(
"Component %1% dimensions: %2%x%3%x%4% cm") % name % component.getWidth() % component.getLength() %
144 component.getHeight());
147 for (
size_t i = 0; i < placements.size(); ++i) {
151 G4Transform3D transform =
getPosition(component, s, p, originCenter);
154 assembly.
add(s.getVolume(), transform);
156 new G4PVPlacement(transform, s.getVolume(), name +
"." + p.getName(), component.getVolume(),
false, i);
161 if (component.getColor().empty()) {
162 B2DEBUG(200,
"Component " << name <<
" is an Air volume, setting invisible");
165 B2DEBUG(200,
"Component " << name <<
" color: " << component.getColor());
166 setColor(*component.getVolume(), component.getColor());
168 B2DEBUG(100,
"--> Created component " << name);
175 G4RotationMatrix rotation(params.getAlpha(), params.getBeta(), params.getGamma());
177 return G4Transform3D(rotation, translation);
183 double u(placement.
getU()), v(placement.
getV()), w(0);
184 switch (placement.
getW()) {
186 w = - mother.
getHeight() / 2.0 - daughter.getHeight() / 2.0;
189 w = - mother.
getHeight() / 2.0 + daughter.getHeight() / 2.0;
195 w = mother.
getHeight() / 2.0 - daughter.getHeight() / 2.0;
198 w = mother.
getHeight() / 2.0 + daughter.getHeight() / 2.0;
205 return G4Translate3D(u, v, w + placement.
getWOffset());
211 G4LogicalVolume& envelopeVolume)
214 G4LogicalVolume* top = &topVolume;
215 if (params.getInsideEnvelope()) {
216 top = &envelopeVolume;
220 const double width = params.getWidth();
221 const double length = params.getLength();
222 const double height = params.getHeight();
227 const std::vector<VXDGeoRadiationSensorsPositionPar>& Positions = params.getPositions();
230 const double r = position.getRadius();
231 const double z = position.getZ();
232 const double theta = position.getTheta();
234 const std::map<int, double>& Sensors = position.getSensors();
235 for (
const std::pair<const int, double>& sensor : Sensors) {
238 const double phi = sensor.second;
239 const int id = sensor.first;
241 const std::string name = params.getSubDetector() +
".DiamondSensor." + std::to_string(
id);
243 if (not shape) shape =
new G4Box(
"radiationSensorDiamond", width / 2 * CLHEP::cm, length / 2 * CLHEP::cm, height / 2 * CLHEP::cm);
244 G4LogicalVolume* volume =
new G4LogicalVolume(shape, material, name);
247 volume->SetSensitiveDetector(sensitive);
249 G4Transform3D transform = G4RotateZ3D(phi - M_PI / 2) * G4Translate3D(0, r * CLHEP::cm,
250 z * CLHEP::cm) * G4RotateX3D(-M_PI / 2 - theta);
251 new G4PVPlacement(transform, volume, name, top,
false, 1);
261 const double tana =
tan(angle);
262 height = min(tana * length, min(tana * width, height));
263 offset = height / tana;
265 const double hwidth = width / 2.0;
266 const double hwidth2 = width2 / 2.0;
267 const double hlength = length / 2.0;
268 const double hheight = height / 2.0;
270 if (width2 <= 0 || width == width2) {
272 return new G4Box(name, hwidth, hlength, hheight);
274 return new G4Trd(name, hwidth, hwidth - offset, hlength, hlength - offset, hheight);
279 return new G4Trap(name, hheight, 0, 0, hlength, hwidth, hwidth2, 0, hlength - offset, hwidth - offset, hwidth2 - offset, 0);
283 const G4Transform3D& placement,
289 G4Transform3D ladderPlacement = placement * G4RotateZ3D(phi) * ladderPos *
getAlignment(parameters.getAlignment(ladder));
294 vector<G4Point3D> lastSensorEdge;
296 VxdID sensorID(ladder);
300 std::map<string, VXDGeoSensor>::iterator it =
m_sensorMap.find(p.getSensorTypeID());
302 B2FATAL(
"Invalid SensorTypeID " << p.getSensorTypeID() <<
", please check the definition of " << sensorID);
305 string name =
m_prefix +
"." + (string)sensorID;
310 G4Transform3D reflection;
311 if (p.getFlipU()) reflection = reflection * G4ReflectX3D();
312 if (p.getFlipV()) reflection = reflection * G4ReflectY3D();
313 if (p.getFlipW()) reflection = reflection * G4ReflectZ3D();
315 G4VSolid* sensorShape =
createTrapezoidal(name, s.getWidth(), s.getWidth2(), s.getLength(),
321 s.setVolume(
new G4LogicalVolume(sensorShape, sensorMaterial, name));
328 G4VSolid* activeShape =
createTrapezoidal(name +
".Active", s.getActiveArea().getWidth(), s.getActiveArea().getWidth2(),
329 s.getActiveArea().getLength(), s.getActiveArea().getHeight());
337 G4LogicalVolume* active =
new G4LogicalVolume(activeShape, sensorMaterial, name +
".Active",
342 setColor(*active, s.getActiveArea().getColor());
346 G4Transform3D activePosition = G4Translate3D(s.getActiveArea().getWidth() / 2.0, s.getActiveArea().getLength() / 2.0, 0) *
347 getPosition(s, s.getActiveArea(), s.getActivePlacement(),
false);
349 G4ReflectionFactory::Instance()->Place(activePosition * reflection, name +
".Active", active, s.getVolume(),
350 false, (
int)sensorID,
false);
356 G4RotationMatrix rotation(0, -M_PI / 2.0, -M_PI / 2.0);
357 G4Transform3D sensorAlign =
getAlignment(parameters.getAlignment(sensorID));
358 G4Transform3D sensorPlacement = G4Rotate3D(rotation) * sensorAlign * reflection;
360 if (s.getSlanted()) {
361 sensorPlacement = G4TranslateX3D(
m_ladder.getSlantedRadius() -
m_ladder.getRadius()) * G4RotateY3D(
362 -
m_ladder.getSlantedAngle()) * sensorPlacement;
364 sensorPlacement = G4Translate3D(0.0, 0.0, p.getZ()) * sensorPlacement;
367 sensorPlacement = ladderPlacement * sensorPlacement;
369 assembly.
add(s.getVolume());
370 assembly.
place(volume, sensorPlacement);
374 double u = s.getWidth() / 2.0 +
m_ladder.getGlueSize();
375 double v = s.getLength() / 2.0;
376 double w = s.getHeight() / 2.0 +
m_ladder.getGlueSize();
377 std::vector<G4Point3D> curSensorEdge(4);
379 curSensorEdge[0] = sensorPlacement * reflection * G4Point3D(u, v, + w);
380 curSensorEdge[1] = sensorPlacement * reflection * G4Point3D(u, v, - w);
381 curSensorEdge[2] = sensorPlacement * reflection * G4Point3D(-u, v, - w);
382 curSensorEdge[3] = sensorPlacement * reflection * G4Point3D(-u, v, + w);
384 if (lastSensorEdge.size()) {
387 for (
int i = 0; i < 4; ++i) glueOK &= curSensorEdge[i].z() <= lastSensorEdge[i].z();
389 B2WARNING(
"Cannot place Glue at sensor " + (
string)sensorID +
390 " since it overlaps with the last module in z");
393 G4TessellatedSolid* solidTarget =
new G4TessellatedSolid(
m_prefix +
".Glue." + (
string)sensorID);
396 solidTarget->AddFacet(
new G4QuadrangularFacet(
397 curSensorEdge[3], curSensorEdge[2], curSensorEdge[1], curSensorEdge[0], ABSOLUTE));
399 solidTarget->AddFacet(
new G4QuadrangularFacet(
400 lastSensorEdge[0], lastSensorEdge[1], lastSensorEdge[2], lastSensorEdge[3], ABSOLUTE));
403 solidTarget->AddFacet(
new G4TriangularFacet(
404 curSensorEdge[3], curSensorEdge[0], lastSensorEdge[0], ABSOLUTE));
405 solidTarget->AddFacet(
new G4TriangularFacet(
406 lastSensorEdge[0], lastSensorEdge[3], curSensorEdge[3], ABSOLUTE));
408 solidTarget->AddFacet(
new G4TriangularFacet(
409 curSensorEdge[1], curSensorEdge[2], lastSensorEdge[2], ABSOLUTE));
410 solidTarget->AddFacet(
new G4TriangularFacet(
411 lastSensorEdge[2], lastSensorEdge[1], curSensorEdge[1], ABSOLUTE));
413 solidTarget->AddFacet(
new G4TriangularFacet(
414 curSensorEdge[0], curSensorEdge[1], lastSensorEdge[1], ABSOLUTE));
415 solidTarget->AddFacet(
new G4TriangularFacet(
416 lastSensorEdge[1], lastSensorEdge[0], curSensorEdge[0], ABSOLUTE));
418 solidTarget->AddFacet(
new G4TriangularFacet(
419 curSensorEdge[2], curSensorEdge[3], lastSensorEdge[3], ABSOLUTE));
420 solidTarget->AddFacet(
new G4TriangularFacet(
421 lastSensorEdge[3], lastSensorEdge[2], curSensorEdge[2], ABSOLUTE));
423 solidTarget->SetSolidClosed(
true);
426 m_prefix +
".Glue." + (
string)sensorID);
428 new G4PVPlacement(G4Transform3D(), glue,
m_prefix +
".Glue." + (
string)sensorID, volume,
false, 1);
432 lastSensorEdge.resize(4);
433 lastSensorEdge[0] = sensorPlacement * reflection * G4Point3D(u, -v, + w);
434 lastSensorEdge[1] = sensorPlacement * reflection * G4Point3D(u, -v, - w);
435 lastSensorEdge[2] = sensorPlacement * reflection * G4Point3D(-u, -v, - w);
436 lastSensorEdge[3] = sensorPlacement * reflection * G4Point3D(-u, -v, + w);
440 return ladderPlacement;
460 sensorInfo.getSensorID(),
461 sensorInfo.getSensorTypeID(),
463 sensorInfo.getFlipU(),
464 sensorInfo.getFlipV(),
465 sensorInfo.getFlipW()
474 VxdID ladder(layerID, ladderID, 0);
477 string path = (boost::format(
"Align[@component='%1%']/") % ladder).str();
480 B2WARNING(
"Could not find alignment parameters for ladder " << ladder);
484 params.getLength(
"dv"),
485 params.getLength(
"dw"),
486 params.getAngle(
"alpha"),
487 params.getAngle(
"beta"),
488 params.getAngle(
"gamma")
494 VxdID sensorID(ladder);
497 std::map<string, VXDGeoSensorPar>::iterator it = vxdGeometryPar.
getSensorMap().find(p.getSensorTypeID());
499 B2FATAL(
"Invalid SensorTypeID " << p.getSensorTypeID() <<
", please check the definition of " << sensorID);
508 string pathSensor = (boost::format(
"Align[@component='%1%']/") % sensorID).str();
511 B2WARNING(
"Could not find alignment parameters for sensorID " << sensorID);
544 string path = (boost::format(
"descendant::Component[@name='%1%']/") % name).str();
545 GearDir params(componentsDir, path);
547 B2FATAL(
"Could not find definition for component " << name);
553 params.getString(
"Color",
""),
554 params.getLength(
"width", 0),
555 params.getLength(
"width2", 0),
556 params.getLength(
"length", 0),
557 params.getLength(
"height", 0),
558 params.getAngle(
"angle", 0)
561 if (c.getWidth() <= 0 || c.getLength() <= 0 || c.getHeight() <= 0) {
562 B2DEBUG(100,
"One dimension empty, using auto resize for component");
569 int chipID = params.getInt(
"activeChipID");
578 string path = (boost::format(
"Ladder[@layer=%1%]/") % layer).str();
579 GearDir paramsLadder(components, path);
581 B2FATAL(
"Could not find Ladder definition for layer " << layer);
588 paramsLadder.
getAngle(
"slantedAngle", 0),
589 paramsLadder.
getLength(
"slantedRadius", 0),
590 paramsLadder.
getLength(
"Glue/oversize", 0),
591 paramsLadder.
getString(
"Glue/Material",
"")
597 sensorInfo.getInt(
"@id"),
598 sensorInfo.getString(
"@type"),
599 sensorInfo.getLength(
"."),
600 sensorInfo.getBool(
"@flipU",
false),
601 sensorInfo.getBool(
"@flipV",
false),
602 sensorInfo.getBool(
"@flipW",
false)
609 vector<VXDGeoPlacementPar> result;
610 for (
const GearDir& component : path.getNodes(
"Component")) {
612 if (!component.exists(
"@type")) {
613 type = component.getString(
"@name");
615 type = component.getString(
"@type");
617 int nPos = max(component.getNumberNodes(
"u"), component.getNumberNodes(
"v"));
618 nPos = max(nPos, component.getNumberNodes(
"w"));
619 nPos = max(nPos, component.getNumberNodes(
"woffset"));
620 for (
int iPos = 1; iPos <= nPos; ++iPos) {
621 string index = (boost::format(
"[%1%]") % iPos).str();
624 component.getLength(
"u" + index, 0),
625 component.getLength(
"v" + index, 0),
626 component.getString(
"w" + index,
"bottom"),
627 component.getLength(
"woffset" + index, 0)
The Class for BeamBackground Sensitive Detector.
GearDir is the basic class used for accessing the parameter store.
virtual std::string getString(const std::string &path="") const noexcept(false) override
Get the parameter path as a string.
Base class for all Sensitive Detectors to create hits during simulation.
static const double mm
[millimeters]
The Class for VXD Alignment payload.
The Class for VXD geometry component.
Class holding all parameters for an VXD geometry component.
double getWidth() const
get the width of the component
double & getHeight()
get the height of the component
double getLength() const
get the length of the component
The Class for VXD Ladder payload.
const std::string & getGlueMaterial() const
get the glue material
double getSlantedAngle() const
get the slant angle for slanted sensors
const std::vector< VXDGeoSensorPlacementPar > & getSensors() const
get list of sensors
double getRadius() const
get the radius of all sensors except slanted ones
double getShift() const
get the shift along the u coordinate for all sensors in the ladder
double getSlantedRadius() const
get the radius for slanted sensors
double getGlueSize() const
get the additional glue size, e.g.
Struct containing all parameters of one ladder.
The Class for VXD placement payload.
Class holding all parameters to place a VXD geometry subcomponent.
double getU() const
get local u coordinate where to place the component
double getV() const
get local v coordinate where to place the component
EPosW getW() const
get local w position where to place the component
@ c_bottom
Place the component at the bottom of the mother.
@ c_above
Place the component above the mother.
@ c_center
Place the component at the center of the mother.
@ c_top
Place the component at the top of the mother.
@ c_below
Place the component below the mother.
double getWOffset() const
get offset to local w position where to place the component
The Class for VXD Radiation Sensor parameters.
The Class for VXD Radiation Sensor Position parameters.
The Class for VXD Sensor payload.
The Class for VXD Sensor Placement payload.
Struct holding the information where a sensor should be placed inside the ladder.
Struct holding all parameters for a completeVXD Sensor.
The Class for VXD geometry.
const std::map< std::string, VXDGeoComponentPar > & getComponentMap() const
get component maps
const std::map< std::string, VXDGeoSensorPar > & getSensorMap() const
get sensor map
const VXDGlobalPar & getGlobalParams() const
get global parameters
std::map< std::string, VXDAlignmentPar > & getAlignmentMap()
get alignment map
const std::vector< std::string > & getComponentInsertOder() const
get component insert order
std::map< int, VXDGeoLadderPar > & getLadderMap()
get ladder map
std::map< std::string, int > & getSensitiveChipIdMap()
get sensitive chip id map
bool getOnlyActiveMaterial() const
Get whether only active materials will be placed for tracking studies.
bool getActiveChips() const
Get whether chips are sensitive.
std::string getDefaultMaterial() const
Get default material.
void addSensorPlacement(VxdID ladder, VxdID sensor, const G4Transform3D &placement)
Remember how sensor is placed into ladder.
static GeoCache & getInstance()
Return a reference to the singleton instance.
void addLadderPlacement(VxdID halfShell, VxdID ladder, const G4Transform3D &placement)
Remember how ladder is placed into half-shell.
Class to group some Geant4 volumes and place them all at once with a given transformation matrix.
void place(G4LogicalVolume *mother, const G4Transform3D &transform)
Place all the volumes already added to the assembly in the given mother.
void add(G4LogicalVolume *volume, const G4Transform3D &transform=G4Transform3D())
Add a volume to the assembly.
std::map< std::string, VXDGeoSensor > m_sensorMap
Map containing Information about all defined sensor types.
float m_minimumElectrons
minimum number of electrons to be deposited by a particle to be saved
G4VSolid * createTrapezoidal(const std::string &name, double width, double width2, double length, double &height, double angle=0)
Create a trapezoidal solid.
virtual void readLadderComponents(int layerID, int ladderID, GearDir content, VXDGeometryPar &vxdGeometryPar)
Read parameters for ladder components and their alignment corresponding to the given ladder id.
void createDiamonds(const VXDGeoRadiationSensorsPar ¶ms, G4LogicalVolume &topVolume, G4LogicalVolume &envelopeVolume)
Create diamond radiation sensors.
GeoVXDRadiationSensors m_radiationsensors
Diamond radiation sensor "sub creator".
std::vector< Simulation::SensitiveDetectorBase * > m_sensitive
List to all created sensitive detector instances.
std::string m_prefix
Prefix to prepend to all volume names.
bool m_onlyActiveMaterial
If this is true, only active Materials will be placed for tracking studies.
double m_activeStepSize
Stepsize to be used inside active volumes.
float m_distanceTolerance
tolerance for Geant4 steps to be merged to a single step
virtual void readLadder(int layer, GearDir components, VXDGeometryPar &geoparameters)
Read parameters for a ladder in layer with given ID from gearbox and layer store them in payload.
VXDGeoLadder m_ladder
Parameters of the currently active ladder.
std::map< std::string, Belle2::VxdID > m_halfShellVxdIDs
Used for translation of half-shell name into a VxdID to consistently handle it in hierarchy.
bool m_onlyPrimaryTrueHits
If true only create TrueHits from primary particles and ignore secondaries.
G4Transform3D getPosition(const VXDGeoComponent &mother, const VXDGeoComponent &daughter, const VXDGeoPlacement &placement, bool originCenter)
Return the position where a daughter component is to be placed.
float m_electronTolerance
tolerance for the energy deposition in electrons to be merged in a single step
std::vector< G4UserLimits * > m_UserLimits
Vector of G4UserLimit pointers.
G4Transform3D getAlignment(const VXDAlignmentPar ¶ms)
Get Alignment from paylead object.
void readComponent(const std::string &name, GearDir components, VXDGeometryPar &vxdGeometryPar)
Read parameters for component name from Gearbox into geometry payload.
virtual ~GeoVXDCreator()
The destructor of the GeoVXDCreator class.
GeoVXDAssembly createSubComponents(const std::string &name, VXDGeoComponent &component, std::vector< VXDGeoPlacement > placements, bool originCenter=true, bool allowOutside=false)
Place a list of subcomponents into an component.
void readSubComponents(const std::vector< VXDGeoPlacementPar > &placements, const GearDir &componentsDir, VXDGeometryPar &vxdGeometryPar)
Read parameters for all components in placement container from Gearbox into payload.
std::map< std::string, VXDGeoComponent > m_componentCache
Cache of all previously created components.
G4Transform3D placeLadder(int ladderID, double phi, G4LogicalVolume *volume, const G4Transform3D &placement, const VXDGeometryPar ¶meters)
Place ladder corresponding to the given ladder id into volume setLayer has to be called first to set ...
virtual void setCurrentLayer(int layer, const VXDGeometryPar ¶meters)
Read parameters for given layer and store in m_ladder.
std::vector< VXDGeoPlacementPar > getSubComponents(const GearDir &path)
Return vector of VXDGeoPlacements with all the components defined inside a given path.
virtual SensitiveDetectorBase * createSensitiveDetector(VxdID sensorID, const VXDGeoSensor &sensor, const VXDGeoSensorPlacement &placement)=0
Return a SensitiveDetector implementation for a given sensor.
bool m_seeNeutrons
Make sensitive detectors also see neutrons.
std::string m_defaultMaterial
Name of the Material to be used for Air.
std::string m_currentHalfShell
Current half-shell being processed (need to know ladder parent for hierarchy)
GeoVXDCreator(const std::string &prefix)
Constructor of the GeoVXDCreator class.
Base class for Sensitive Detector implementation of PXD and SVD.
void setOptions(bool seeNeutrons, bool onlyPrimaryTrueHits, float distanceTolerance, float electronTolerance, float minimumElectrons)
Set all common options.
Class to uniquely identify a any structure of the PXD and SVD.
void setSensorNumber(baseType sensor)
Set the sensor id.
double getAngle(const std::string &path="") const noexcept(false)
Get the parameter path as a double converted to the standard angle unit.
double getLength(const std::string &path="") const noexcept(false)
Get the parameter path as a double converted to the standard length unit.
std::vector< GearDir > getNodes(const std::string &path="") const
Get vector of GearDirs which point to all the nodes the given path evaluates to.
static G4Material * get(const std::string &name)
Find given material.
double tan(double a)
tan for double
Namespace to provide code needed by both Vertex Detectors, PXD and SVD, and also testbeam telescopes.
Common code concerning the geometry representation of the detector.
void setVisibility(G4LogicalVolume &volume, bool visible)
Helper function to quickly set the visibility of a given volume.
void setColor(G4LogicalVolume &volume, const std::string &color)
Set the color of a logical volume.
Abstract base class for different kinds of events.