GeoSTRCreator Class Reference

The GeoSTRCreator class to create the Belle2 structure geometry. More...

#include <GeoSTRCreator.h>

Inheritance diagram for GeoSTRCreator:

Public Member Functions
	GeoSTRCreator ()
	The Constructor of the GeoSTRCreator class.
virtual	~GeoSTRCreator ()
	The destructor of the GeoSTRCreator class.
virtual void	create (const GearDir &content, G4LogicalVolume &topVolume, geometry::GeometryTypes type) override
	Creates the Geant4 objects for the structure geometry from Gearbox.
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
STRGeometryPar	createConfiguration (const GearDir &param)
	Create a parameter object from the Gearbox XML parameters.
void	createGeometry (const STRGeometryPar &parameters, G4LogicalVolume &topVolume, geometry::GeometryTypes type)
	Create the geometry from a parameter object.
void	readShield (const GearDir &content, STRGeometryPar &parameters, std::string side)
	Read xml files from GearDir for one of the ECL shields.
void	readPole (const GearDir &content, STRGeometryPar &parameters, std::string side)
	Read xml files from GearDir for one of the Pole Pieces.

Detailed Description

The GeoSTRCreator class to create the Belle2 structure geometry.

The creator for the Belle II structural elements such as radiation shields and pole piece

Definition at line 38 of file GeoSTRCreator.h.

Constructor & Destructor Documentation

GeoSTRCreator()

GeoSTRCreator

(

)

The Constructor of the GeoSTRCreator class.

Definition at line 49 of file GeoSTRCreator.cc.

    {
 
    }

~GeoSTRCreator()

~GeoSTRCreator ( )

virtual

The destructor of the GeoSTRCreator class.

Definition at line 55 of file GeoSTRCreator.cc.

    {
 
    }

Member Function Documentation

create()

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

inlineoverridevirtual

Creates the Geant4 objects for the structure geometry from Gearbox.

This is the 'old' way of building the geometry from Gearbox (xml files). Expected to be deprecated in the future.

Parameters

content	A reference to the content part of the parameter description, which should to be used to create the ROOT objects.
topVolume	Geant4 logical top volume.
type	Geometry type.

Implements CreatorBase.

Definition at line 55 of file GeoSTRCreator.h.

      {
        STRGeometryPar config = createConfiguration(content);
        createGeometry(config, topVolume, type);
      }

createConfiguration()

STRGeometryPar createConfiguration ( const GearDir &  param )

private

Create a parameter object from the Gearbox XML parameters.

If more than one object is created these could be assigned to members or you could return a tuple.

Definition at line 148 of file GeoSTRCreator.cc.

    {
      STRGeometryPar strGeometryPar;
 
      // Get STR geometry parameters from Gearbox (no calculations here)
      readShield(param, strGeometryPar, "FWD_Shield");
      readShield(param, strGeometryPar, "BWD_Shield");
      readPole(param, strGeometryPar, "PolePieceL");
      readPole(param, strGeometryPar, "PolePieceR");
 
 
      return strGeometryPar;
    };

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 71 of file GeoSTRCreator.h.

      {
        DBObjPtr<STRGeometryPar> dbObj;
        if (!dbObj) {
          // Check that we found the object and if not report the problem
          B2FATAL("No configuration for " << name << " found.");
        }
        createGeometry(*dbObj, topVolume, type);
      }

createGeometry()

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

private

Create the geometry from a parameter object.

If more than one object is needed these could instead be taken from member variables or the number of parameters could be increased.

Definition at line 62 of file GeoSTRCreator.cc.

    {
 
 
 
 
      for (int iShield = 0; iShield < parameters.NECLSHIELDS; iShield++) {
        float  _mass = 0;
        string side;
        if (iShield == parameters.FWD_ECLSHIELD) {
          side = "FWD_Shield";
        } else if (iShield == parameters.BWD_ECLSHIELD) {
          side = "BWD_Shield";
        } else {
          B2FATAL("Only 2 ECL shields should be defined. Can't retrieve info for shield #" << iShield);
        }
 
        int nLayers = parameters.getNLayers(iShield);
 
        for (int iLayer = 0; iLayer < nLayers; iLayer++) {
          G4Material* LayerMat = Materials::get(parameters.getLayerMaterial(iShield, iLayer));
          int nPlanes = parameters.getLayerNPlanes(iShield, iLayer);
 
          //Thread the strings
          string shapeName   = (boost::format("%1%Layer_%2%") % side % (iLayer + 1)).str();
          string logiVolName = (boost::format("logi%1%Layer_%2%") % side % (iLayer + 1)).str();
          string physVolName = (boost::format("phys%1%Layer_%2%") % side % (iLayer + 1)).str();
 
          G4Polycone* LayerShape = new G4Polycone(shapeName.c_str(), 0, 2 * M_PI, nPlanes,
                                                  parameters.getLayerPlaneZ(iShield, iLayer),
                                                  parameters.getLayerPlaneInnerRadius(iShield, iLayer),
                                                  parameters.getLayerPlaneOuterRadius(iShield, iLayer));
 
          // Create logical volume
          G4LogicalVolume* logiShieldLayer = new G4LogicalVolume(LayerShape, LayerMat, logiVolName, 0, 0, 0);
 
          //Place physical volume
          new G4PVPlacement(0, G4ThreeVector(0, 0, 0), logiShieldLayer, physVolName, &topVolume, false, 0);
 
          B2DEBUG(1, "Mass of " << side << " layer " << iLayer
                  << " = " << logiShieldLayer->GetMass() / CLHEP::kg << ".kg.");
 
          _mass += (logiShieldLayer->GetMass() / CLHEP::kg);
        }
 
        B2DEBUG(1, "Total mass of side " << side << " = " << _mass << " kg");
      }
 
      //Place pole pieces
      for (int iPole = 0; iPole < parameters.NPOLEPIECES; iPole++) {
 
        string side;
        if (iPole == parameters.FWD_POLEPIECE) {
          side = "PolePieceR";
        } else if (iPole == parameters.BWD_POLEPIECE) {
          side = "PolePieceL";
        } else {
          B2FATAL("Only 2 pole pieces should be defined. Can't retrieve info for pole #" << iPole);
        }
 
 
        G4Material* PoleMat = Materials::get(parameters.getPoleMaterial(iPole));
        int nPlanes = parameters.getPoleNPlanes(iPole);
 
        //Thread the strings
        string shapeName   = (boost::format("%1%") % side).str();
        string logiVolName = (boost::format("logi%1%") % side).str();
        string physVolName = (boost::format("phys%1%") % side).str();
 
        G4Polycone* PoleShape = new G4Polycone(shapeName.c_str(), 0, 2 * M_PI, nPlanes,
                                               parameters.getPolePlaneZ(iPole),
                                               parameters.getPolePlaneInnerRadius(iPole),
                                               parameters.getPolePlaneOuterRadius(iPole));
 
        // Create logical volume
        G4LogicalVolume* logiPole = new G4LogicalVolume(PoleShape, PoleMat, logiVolName, 0, 0, 0);
 
        //Place physical volume
        new G4PVPlacement(0, G4ThreeVector(0, 0, 0), logiPole, physVolName, &topVolume, false, 0);
 
        B2DEBUG(1, "Total mass of " << side << " = " << logiPole->GetMass() / CLHEP::kg  << " kg");
      }
 
 
    }

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 63 of file GeoSTRCreator.h.

      {
        DBImportObjPtr<STRGeometryPar> importObj;
        importObj.construct(createConfiguration(content));
        importObj.import(iov);
      }

readPole()

void readPole ( const GearDir &  content, STRGeometryPar &  parameters, std::string  side  )

private

Read xml files from GearDir for one of the Pole Pieces.

Definition at line 162 of file GeoSTRCreator.cc.

    {
      // Check if method was called using the correct name for the shields
      std::size_t foundF = side.find("PolePieceR");
      std::size_t foundB = side.find("PolePieceL");
 
      int iPole = 0;
      if (foundF != std::string::npos) { iPole = parameters.FWD_POLEPIECE; }
      else if (foundB != std::string::npos) { iPole = parameters.BWD_POLEPIECE; }
      else { B2FATAL("No data for the Pole Piece requested " << side << "(not found)");}
 
 
      //Thread the strings
      std::string polePath   = (boost::format("/%1%/") % side).str();
 
      // Connect the appropriate Gearbox path
      GearDir poleContent(content);
      poleContent.append(polePath);
 
      // Retrieve material material
      parameters.setPoleMaterial(iPole, poleContent.getString("Material", "Air"));
 
      // Read the shape parameters
      const std::vector<GearDir> planes = poleContent.getNodes("Plane");
      parameters.setPoleNPlanes(iPole, planes.size());
      B2DEBUG(1, "Number of planes on side  " << side << " : " << planes.size());
 
      for (unsigned int iPlane = 0; iPlane < planes.size(); iPlane++) {
        parameters.setPolePlaneZ(iPole, iPlane, planes.at(iPlane).getLength("posZ") / Unit::mm);
        parameters.setPolePlaneInnerRadius(iPole, iPlane, planes.at(iPlane).getLength("innerRadius") / Unit::mm);
        parameters.setPolePlaneOuterRadius(iPole, iPlane, planes.at(iPlane).getLength("outerRadius") / Unit::mm);
 
      }
    }

readShield()

void readShield ( const GearDir &  content, STRGeometryPar &  parameters, std::string  side  )

private

Read xml files from GearDir for one of the ECL shields.

Definition at line 197 of file GeoSTRCreator.cc.

    {
 
      // Check if method was called using the correct name for the shields
      std::size_t foundF = side.find("FWD_Shield");
      std::size_t foundB = side.find("BWD_Shield");
 
      int iShield = 0;
      if (foundF != std::string::npos) { iShield = parameters.FWD_ECLSHIELD; }
      else if (foundB != std::string::npos) { iShield = parameters.BWD_ECLSHIELD; }
      else { B2FATAL("No data for the ECL shield called " << side << "(not found)");}
 
 
      std::string gearPath = (boost::format("%1%/Layers/Layer") % side).str();
 
      // Retrieve the number of layers in the shield
      int nLayers = content.getNumberNodes(gearPath);
      parameters.setNLayers(iShield, nLayers);
 
 
      for (int iLayer = 0 ; iLayer < nLayers ; ++iLayer) {
        //Thread the strings
        std::string layerPath   = (boost::format("/%1%[%2%]/") % gearPath % (iLayer + 1)).str();
 
        // Connect the appropriate Gearbox path
        GearDir layerContent(content);
        layerContent.append(layerPath);
 
        // Retrieve material material
        parameters.setLayerMaterial(iShield, iLayer, layerContent.getString("Material", "Air"));
 
        // Read the shape parameters
        const std::vector<GearDir> planes = layerContent.getNodes("Plane");
        parameters.setLayerNPlanes(iShield, iLayer, planes.size());
        B2DEBUG(1, "Number of planes on side  " << side << " layer " << iLayer
                << " : " << planes.size());
 
        for (unsigned int iPlane = 0; iPlane < planes.size(); iPlane++) {
          parameters.setLayerPlaneZ(iShield, iLayer, iPlane, planes.at(iPlane).getLength("posZ") / Unit::mm);
          parameters.setLayerPlaneInnerRadius(iShield, iLayer, iPlane, planes.at(iPlane).getLength("innerRadius") / Unit::mm);
          parameters.setLayerPlaneOuterRadius(iShield, iLayer, iPlane, planes.at(iPlane).getLength("outerRadius") / Unit::mm);
        }
      }
    }

---

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

• structure/geometry/include/GeoSTRCreator.h
• structure/geometry/src/GeoSTRCreator.cc