The Class for CSI Geometry Parameters. More...

#include <CsiGeometryPar.h>

Public Member Functions
	CsiGeometryPar ()
	Constructor.

virtual	~CsiGeometryPar ()
	Destructor.

void	clear ()
	Clears.

void	Print (const int cid, int debuglevel=80)
	Print crystal information.

void	PrintAll (int debuglevel=80)
	Print all crystals information.

void	read ()
	Gets geometry parameters from gearbox.

int	CsiVolNameToCellID (const G4String VolumeName)
	Get Cell Id.

G4Material *	GetMaterial (int cid)
	Get pointer to the Geant4 Material.

ROOT::Math::XYZVector	GetPosition (int cid)
	Get the position of the crystal.

ROOT::Math::XYZVector	GetOrientation (int cid)
	Get the orientation of the crystal.

int	GetEnclosureID (int cid)
	Get Enclosure ID from cell ID.

int	GetSlotID (int cid)
	Get Slot ID in the Enclosure from cell ID.

double	GetMaterialProperty (int cid, const char *propertyname)
	Get material property.

double	GetTauFast (int cid)
	Get crystal fast time constant.

double	GetTauSlow (int cid)
	Get crystal slow time constant.

Static Public Member Functions
static CsiGeometryPar *	Instance ()
	Static method to get a reference to the CsiGeometryPar instance.

Private Attributes
int	m_cellID
	The Cell ID information.

std::vector< int >	m_thetaID
	The Theta ID information.

std::vector< int >	m_phiID
	The Phi ID information.

std::vector< int >	m_BoxID
	The index of the enclosure.

std::vector< int >	m_SlotID
	The slot index of the crystal in the enclosure.

std::vector< ROOT::Math::XYZVector >	m_Position
	Position of the nominal centre of the crystal.

std::vector< ROOT::Math::XYZVector >	m_Orientation
	Orientation of the crystal.

Static Private Attributes
static CsiGeometryPar *	m_B4CsiGeometryParDB = 0
	Pointer that saves the instance of this class.

Detailed Description

The Class for CSI Geometry Parameters.

This class provides CSI gemetry paramters for simulation, reconstruction and so on. These parameters are gotten from gearbox.

Definition at line 48 of file CsiGeometryPar.h.

Constructor & Destructor Documentation

◆ CsiGeometryPar()

CsiGeometryPar ( )

Constructor.

Definition at line 32 of file CsiGeometryPar.cc.

{
  clear();
  read();
 
  PrintAll();
}

◆ ~CsiGeometryPar()

~CsiGeometryPar ( )

virtual

Destructor.

Definition at line 40 of file CsiGeometryPar.cc.

{
  if (m_B4CsiGeometryParDB) {
    delete m_B4CsiGeometryParDB;
    B2INFO("m_B4CsiGeometryParDB deleted ");
  }
}

Member Function Documentation

◆ clear()

void clear ( )

Clears.

Definition at line 48 of file CsiGeometryPar.cc.

{
  m_cellID = 0;
 
  m_Position.clear();
  m_Orientation.clear();
  m_BoxID.clear();
  m_SlotID.clear();
  m_thetaID.clear();
  m_phiID.clear();
 
}

◆ CsiVolNameToCellID()

int CsiVolNameToCellID ( const G4String VolumeName )

Get Cell Id.

Definition at line 134 of file CsiGeometryPar.cc.

{
  int cellID = 0;
 
  vector< string > partName;
  split(partName, VolumeName, is_any_of("_"));
 
  int iEnclosure = -1;
  int iCrystal   = -1;
  for (std::vector<string>::iterator it = partName.begin() ; it != partName.end(); ++it) {
    if (equals(*it, "Enclosure")) iEnclosure = boost::lexical_cast<int>(*(it + 1)) - 1;
    else if (equals(*it, "Crystal")) iCrystal = boost::lexical_cast<int>(*(it + 1)) - 1;
  }
 
  cellID = 3 * iEnclosure + iCrystal;
 
  if (cellID < 0) B2WARNING("CsiGeometryPar: volume " << VolumeName << " is not a crystal");
 
  return cellID;
}

◆ GetEnclosureID()

int GetEnclosureID ( int cid )

inline

Get Enclosure ID from cell ID.

Definition at line 90 of file CsiGeometryPar.h.

90{ return m_BoxID.at(cid); };

◆ GetMaterial()

G4Material * GetMaterial ( int cid )

Get pointer to the Geant4 Material.

Definition at line 156 of file CsiGeometryPar.cc.

{
  int iEnclosure = GetEnclosureID(cid) + 1;
  int iSlot      = GetSlotID(cid) + 1;
 
  GearDir content = GearDir("/Detector/DetectorComponent[@name=\"CSI\"]/Content/");
 
  GearDir enclosureContent(content);
  string gearPath = "Enclosures/Enclosure";
  string enclosurePath = (format("/%1%[%2%]") % gearPath % iEnclosure).str();
  enclosureContent.append(enclosurePath);
 
  string slotName = (format("CrystalInSlot[%1%]") % iSlot).str();
  int iCry = enclosureContent.getInt(slotName);
 
 
  GearDir crystalContent(content);
  crystalContent.append((format("/EndCapCrystals/EndCapCrystal[%1%]/") % (iCry)).str());
  string strMatCrystal = crystalContent.getString("Material", "Air");
 
  return geometry::Materials::get(strMatCrystal);
 
}

◆ GetMaterialProperty()

double GetMaterialProperty	(	int	cid,
		const char *	propertyname
	)

Get material property.

Definition at line 181 of file CsiGeometryPar.cc.

{
  G4Material* mat = GetMaterial(cid);
  G4MaterialPropertiesTable* properties = mat->GetMaterialPropertiesTable();
  G4MaterialPropertyVector* property = properties->GetProperty(propertyname);
 
  return property->Value(0);
}

◆ GetOrientation()

ROOT::Math::XYZVector GetOrientation ( int cid )

inline

Get the orientation of the crystal.

Definition at line 87 of file CsiGeometryPar.h.

87{ return m_Orientation.at(cid); };

◆ GetPosition()

ROOT::Math::XYZVector GetPosition ( int cid )

inline

Get the position of the crystal.

Definition at line 84 of file CsiGeometryPar.h.

84{ return m_Position.at(cid); };

◆ GetSlotID()

int GetSlotID ( int cid )

inline

Get Slot ID in the Enclosure from cell ID.

Definition at line 93 of file CsiGeometryPar.h.

93{ return m_SlotID.at(cid); };

◆ GetTauFast()

double GetTauFast ( int cid )

inline

Get crystal fast time constant.

Definition at line 99 of file CsiGeometryPar.h.

99{ return GetMaterialProperty(cid, "FASTTIMECONSTANT"); };

Belle2::csi::CsiGeometryPar::GetMaterialProperty

double GetMaterialProperty(int cid, const char *propertyname)

Get material property.

Definition: CsiGeometryPar.cc:181

◆ GetTauSlow()

double GetTauSlow ( int cid )

inline

Get crystal slow time constant.

Definition at line 102 of file CsiGeometryPar.h.

102{ return GetMaterialProperty(cid, "SLOWTIMECONSTANT"); };

◆ Instance()

CsiGeometryPar * Instance ( )

static

Static method to get a reference to the CsiGeometryPar instance.

Returns: A reference to an instance of this class.

Definition at line 26 of file CsiGeometryPar.cc.

{
  if (!m_B4CsiGeometryParDB) m_B4CsiGeometryParDB = new CsiGeometryPar();
  return m_B4CsiGeometryParDB;
}

◆ Print()

void Print	(	const int	cid,
		int	debuglevel = `80`
	)

Print crystal information.

Definition at line 190 of file CsiGeometryPar.cc.

{
  B2DEBUG(debuglevel, "Cell ID : " << cid);
 
  B2DEBUG(debuglevel, "   Position  x : " << GetPosition(cid).x());
  B2DEBUG(debuglevel, "   Position  y : " << GetPosition(cid).y());
  B2DEBUG(debuglevel, "   Position  z : " << GetPosition(cid).z());
 
  B2DEBUG(debuglevel, "   Orientation  x : " << GetOrientation(cid).x());
  B2DEBUG(debuglevel, "   Orientation  y : " << GetOrientation(cid).y());
  B2DEBUG(debuglevel, "   Orientation  z : " << GetOrientation(cid).z());
 
  B2DEBUG(debuglevel, "   Material : " << GetMaterial(cid)->GetName());
 
  B2DEBUG(debuglevel, "   Slow time constatnt : " << GetMaterialProperty(cid, "SLOWTIMECONSTANT"));
  B2DEBUG(debuglevel, "   Fast time constatnt : " << GetMaterialProperty(cid, "FASTTIMECONSTANT"));
  B2DEBUG(debuglevel, "   Light yield : "         << GetMaterialProperty(cid, "SCINTILLATIONYIELD"));
 
  //GetMaterial(cid)->GetMaterialPropertiesTable()->DumpTable();
}

◆ PrintAll()

void PrintAll ( int debuglevel = 80 )

Print all crystals information.

Definition at line 211 of file CsiGeometryPar.cc.

{
  for (uint i = 0; i < m_thetaID.size(); i++)
    Print(i, debuglevel);
}

◆ read()

void read ( )

Gets geometry parameters from gearbox.

Position of the nominal centre of crystals in the box

Definition at line 61 of file CsiGeometryPar.cc.

{
 
  GearDir content = GearDir("/Detector/DetectorComponent[@name=\"CSI\"]/Content/");
  string gearPath = "Enclosures/Enclosure";
  int nEnc = content.getNumberNodes(gearPath);
 
  int iCell = 0;
 
  for (int iEnc = 1; iEnc <= nEnc; iEnc++) {
    // Build the box (same for all)
    //double length = content.getLength("Enclosures/Length") * CLHEP::cm;
    //double thk    = content.getLength("Enclosures/Thickness") * CLHEP::cm;
    //double halflength = 15.0 * CLHEP::cm;
    //double zshift = 0.5 * length - thk - halflength; /*< Shift of the box along z-axis (cry touches panel) **/
 
    string enclosurePath = (boost::format("/%1%[%2%]") % gearPath % iEnc).str();
 
    // Connect the appropriate Gearbox path
    GearDir enclosureContent(content);
    enclosureContent.append(enclosurePath);
 
    // Read position
    double PosZ  = enclosureContent.getLength("PosZ") * CLHEP::cm;
    double PosR  = enclosureContent.getLength("PosR") * CLHEP::cm;
    double PosT  = enclosureContent.getAngle("PosT") ;
 
    // Read Orientation
    double Phi1  = enclosureContent.getAngle("AngPhi1") ;
    double Theta = enclosureContent.getAngle("AngTheta") ;
    double Phi2  = enclosureContent.getAngle("AngPhi2") ;
 
    //Transform3D zsh = Translate3D(0, 0, zshift);
    Transform3D m1 = RotateZ3D(Phi1);
    Transform3D m2 = RotateY3D(Theta);
    Transform3D m3 = RotateZ3D(Phi2);
    Transform3D position = Translate3D(PosR * cos(PosT), PosR * sin(PosT), PosZ);
 
    Transform3D Tr    = position * m3 * m2 * m1;
 
    int nSlots = enclosureContent.getNumberNodes("CrystalInSlot");
    for (int iSlot = 1; iSlot <= nSlots; iSlot++) {
      iCell++;
 
      //Thread the strings
      string slotPath = (boost::format("/Enclosures/Slot[%1%]") % iSlot).str();
 
      GearDir slotContent(content);
      slotContent.append(slotPath);
 
      double SlotX = slotContent.getLength("PosX") * CLHEP::cm;
      double SlotY = slotContent.getLength("PosY") * CLHEP::cm;
      double SlotZ = slotContent.getLength("PosZ") * CLHEP::cm;
      Transform3D Pos = Translate3D(SlotX, SlotY, SlotZ);
 
      Transform3D CrystalPos = Tr * Pos;
      RotationMatrix CrystalRot = CrystalPos.getRotation();
 
      m_Position.push_back(ROOT::Math::XYZVector(CrystalPos.getTranslation()) * 1.0 / CLHEP::cm);
      m_Orientation.push_back(ROOT::Math::XYZVector(CrystalRot.colZ()));
 
      m_thetaID.push_back(CrystalPos.getTranslation().z() > 0 ? 0 : 1);
      m_phiID.push_back(iCell - 9 * m_thetaID.back());
 
      m_BoxID.push_back(iEnc - 1);
      m_SlotID.push_back(iSlot - 1);
    }
    //
  }
}

Member Data Documentation

◆ m_B4CsiGeometryParDB

CsiGeometryPar * m_B4CsiGeometryParDB = 0

staticprivate

Pointer that saves the instance of this class.

Definition at line 130 of file CsiGeometryPar.h.

◆ m_BoxID

std::vector<int> m_BoxID

private

The index of the enclosure.

Definition at line 118 of file CsiGeometryPar.h.

◆ m_cellID

int m_cellID

private

The Cell ID information.

Definition at line 109 of file CsiGeometryPar.h.

◆ m_Orientation

std::vector<ROOT::Math::XYZVector> m_Orientation

private

Orientation of the crystal.

Definition at line 127 of file CsiGeometryPar.h.

◆ m_phiID

std::vector<int> m_phiID

private

The Phi ID information.

Definition at line 115 of file CsiGeometryPar.h.

◆ m_Position

std::vector<ROOT::Math::XYZVector> m_Position

private

Position of the nominal centre of the crystal.

Definition at line 124 of file CsiGeometryPar.h.

◆ m_SlotID

std::vector<int> m_SlotID

private

The slot index of the crystal in the enclosure.

Definition at line 121 of file CsiGeometryPar.h.

◆ m_thetaID

std::vector<int> m_thetaID

private

The Theta ID information.

Definition at line 112 of file CsiGeometryPar.h.

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

beast/csi/geometry/include/CsiGeometryPar.h
beast/csi/geometry/src/CsiGeometryPar.cc

Public Member Functions

Static Public Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CsiGeometryPar()

◆ ~CsiGeometryPar()

Member Function Documentation

◆ clear()

◆ CsiVolNameToCellID()

◆ GetEnclosureID()

◆ GetMaterial()

◆ GetMaterialProperty()

◆ GetOrientation()

◆ GetPosition()

◆ GetSlotID()

◆ GetTauFast()

◆ GetTauSlow()

◆ Instance()

◆ Print()

◆ PrintAll()

◆ read()

Member Data Documentation

◆ m_B4CsiGeometryParDB

◆ m_BoxID

◆ m_cellID

◆ m_Orientation

◆ m_phiID

◆ m_Position

◆ m_SlotID

◆ m_thetaID