GeometryData Class Reference

EKLM geometry data. More...

#include <GeometryData.h>

Inheritance diagram for GeometryData:

Public Types
enum	DataSource {   c_Gearbox ,   c_Database }
	Geometry data source. More...

Public Member Functions
void	saveToDatabase (const IntervalOfValidity &iov) const
	Save geometry data to database.
double	getStripLength (int strip) const
	Get strip length.
double	getMaximalStripLength () const
	Get maximal strip length.
int	getNStripsDifferentLength () const
	Get number of strips with different lengths.
int	getStripLengthIndex (int positionIndex) const
	Get index in length-based array.
int	getStripPositionIndex (int lengthIndex) const
	Get index in position-based array.
bool	hitInEKLM (double z) const
	Check if z coordinate may be in EKLM.
void	getSectionTransform (HepGeom::Transform3D *t, int n) const
	Get section transformation.
void	getLayerTransform (HepGeom::Transform3D *t, int n) const
	Get layer transformation.
void	getSectorTransform (HepGeom::Transform3D *t, int n) const
	Get sector transformation.
void	getPlaneTransform (HepGeom::Transform3D *t, int n) const
	Get plane transformation.
void	getStripTransform (HepGeom::Transform3D *t, int n) const
	Get strip transformation.
void	getSheetTransform (HepGeom::Transform3D *t, int n) const
	Get plastic sheet element transformation.
bool	beamBackgroundStudy () const
	Whether to perform beam-background study.
int	getNSections () const
	Get number of sections.
int	getNLayers () const
	Get number of layers.
int	getNDetectorLayers (int section) const
	Get number of detector layers.
int	getNSectors () const
	Get number of sectors.
int	getNPlanes () const
	Get number of planes.
int	getNSegments () const
	Get number of segments.
int	getNSegmentSupportElementsSector () const
	Get number of segment support elements *in one sector).
int	getNStrips () const
	Get number of strips.
void	checkDetectorLayerNumber (int section, int layer) const
	Check if number of detector layers is correct (fatal error if not).
void	checkDetectorLayer (int section, int layer) const
	Check if detector layer number is correct (fatal error if not).
void	checkSegmentSupport (int support) const
	Check if segment support number is correct (fatal error if not).
void	checkStripSegment (int strip) const
	Check if number of strip in a segment is correct (fatal error if not).
double	getSolenoidZ () const
	Get solenoid center Z coordinate.
const EndcapStructureGeometry *	getEndcapStructureGeometry () const
	Get endcap structure geometry data.
const ElementPosition *	getSectionPosition () const
	Get position data for sections.
const ElementPosition *	getLayerPosition () const
	Get position data for layers.
double	getLayerShiftZ () const
	Get Z distance between two layers.
const ElementPosition *	getSectorPosition () const
	Get position data for sectors.
const ElementPosition *	getSectorSupportPosition () const
	Get position data for sector support structure.
const SectorSupportGeometry *	getSectorSupportGeometry () const
	Get sector support geometry data.
const ElementPosition *	getPlanePosition () const
	Get position data for planes.
const PlasticSheetGeometry *	getPlasticSheetGeometry () const
	Get plastic sheet geometry data.
const SegmentSupportGeometry *	getSegmentSupportGeometry () const
	Get segment support geometry data.
const SegmentSupportPosition *	getSegmentSupportPosition (int plane, int support) const
	Get position data for segment support structure.
const StripGeometry *	getStripGeometry () const
	Get strip geometry data.
const ElementPosition *	getStripPosition (int strip) const
	Get position data for strips.
const ShieldGeometry *	getShieldGeometry () const
	Get shield layer details geometry data.

Static Public Member Functions
static const GeometryData &	Instance (enum DataSource dataSource=c_Database, const GearDir *gearDir=nullptr)
	Instantiation.

Protected Member Functions
	ClassDef (Belle2::EKLMGeometry, 6)
	Class version.

Protected Attributes
const EKLMElementNumbers *	m_ElementNumbers
	Element numbers.
bool	m_BeamBackgroundStudy = false
	ROOT streamer.
int	m_NSections
	Number of sections.
int	m_NLayers
	Number of layers in one section.
int *	m_NDetectorLayers
	Number of detector layers.
int	m_NSectors
	Number of sectors in one layer.
int	m_NPlanes
	Number of planes in one sector.
int	m_NSegments
	Number of segments in one plane.
int	m_NSegmentSupportElementsSector
	Number of segment support elements in one sector.
int	m_NStrips
	Number of strips in one plane.
double	m_SolenoidZ
	Solenoid center Z coordinate.
EndcapStructureGeometry	m_EndcapStructureGeometry
	Section structure geometry data.
ElementPosition	m_SectionPosition
	Position data for sections.
ElementPosition	m_LayerPosition
	Position data for layers.
double	m_LayerShiftZ
	Z distance between two layers.
ElementPosition	m_SectorPosition
	Position data for sectors.
ElementPosition	m_SectorSupportPosition
	Position data for sector support structure.
SectorSupportGeometry	m_SectorSupportGeometry
	Sector support geometry data.
ElementPosition	m_PlanePosition
	Position data for planes.
PlasticSheetGeometry	m_PlasticSheetGeometry
	Plastic sheet geometry data.
SegmentSupportGeometry	m_SegmentSupportGeometry
	Segment support geometry data.
SegmentSupportPosition *	m_SegmentSupportPosition
	Position data for segment support structure.
StripGeometry	m_StripGeometry
	Strip geometry data.
ElementPosition *	m_StripPosition
	Position data for strips.
ShieldGeometry	m_ShieldGeometry
	Shield layer details geometry data.

Private Member Functions
	GeometryData (enum DataSource dataSource, const GearDir *gearDir)
	Constructor.
	~GeometryData ()
	Destructor.
void	readEndcapStructureGeometry (const GearDir &gd)
	Read section structure geometry data.
void	initializeFromGearbox (const GearDir *gearDir)
	Initialize from Gearbox (XML).
void	initializeFromDatabase ()
	Initialize from database.
void	calculateSectorSupportGeometry ()
	Calculate sector support geometry data.
void	fillStripIndexArrays ()
	Fill strip index arrays.
void	readXMLDataStrips (const GearDir &gd)
	Read strip parameters from XML database.
void	calculateShieldGeometry ()
	Calculate shield geometry data.

Private Attributes
EKLMGeometry *	m_Geometry
	Copy of data in this class used to write it to database.
int	m_nStripDifferent
	Number of strips with different lengths in one plane.
int *	m_StripAllToLen
	Number of strip in length-based array.
int *	m_StripLenToAll
	Number of strip in position-based array.
double	m_MinZForward
	Minimal z coordinate of the forward section.
double	m_MaxZBackward
	Maximal z coordinate of the backward section.

Detailed Description

EKLM geometry data.

All data and function results are in CLHEP units unless noted otherwise.

Definition at line 38 of file GeometryData.h.

Member Enumeration Documentation

DataSource

enum DataSource

Geometry data source.

Enumerator
c_Gearbox	Gearbox (XML).
c_Database	Database.

Definition at line 43 of file GeometryData.h.

                      {
 
        c_Gearbox,
 
        c_Database,
 
      };

Constructor & Destructor Documentation

GeometryData()

GeometryData ( enum DataSource  dataSource, const GearDir *  gearDir  )

explicitprivate

Constructor.

Parameters

[in]	dataSource	Data source.
[in]	gearDir	XML data directory.

Definition at line 605 of file GeometryData.cc.

{
  m_Geometry = nullptr;
  switch (dataSource) {
    case c_Gearbox:
      initializeFromGearbox(gearDir);
      break;
    case c_Database:
      initializeFromDatabase();
      break;
  }
  m_MinZForward = m_SolenoidZ + m_SectionPosition.getZ() -
                  0.5 * m_SectionPosition.getLength();
  m_MaxZBackward = m_SolenoidZ - m_SectionPosition.getZ() +
                   0.5 * m_SectionPosition.getLength();
  calculateSectorSupportGeometry();
  fillStripIndexArrays();
  calculateShieldGeometry();
}

~GeometryData()

~GeometryData ( )

private

Destructor.

Definition at line 626 of file GeometryData.cc.

{
  if (m_Geometry != nullptr)
    delete m_Geometry;
  free(m_StripLenToAll);
  free(m_StripAllToLen);
}

Member Function Documentation

beamBackgroundStudy()

bool beamBackgroundStudy ( ) const

inlineinherited

Whether to perform beam-background study.

Definition at line 1677 of file EKLMGeometry.h.

    {
      return m_BeamBackgroundStudy;
    }

calculateSectorSupportGeometry()

void calculateSectorSupportGeometry ( )

private

Calculate sector support geometry data.

Definition at line 117 of file GeometryData.cc.

{
  HepGeom::Point3D<double> p;
  Line2D line23Outer(0, m_SectorSupportPosition.getY(), 1, 0);
  Line2D line23Inner(0, m_SectorSupportPosition.getY() +
                     m_SectorSupportGeometry.getThickness(), 1, 0);
  Line2D line23Prism(0, m_SectorSupportPosition.getY() +
                     m_SectorSupportGeometry.getThickness() +
                     m_SectorSupportGeometry.getCorner3LY(), 1, 0);
  Line2D line41Outer(m_SectorSupportPosition.getX(), 0, 0, 1);
  Line2D line41Inner(m_SectorSupportPosition.getX() +
                     m_SectorSupportGeometry.getThickness(), 0, 0, 1);
  Line2D line41Prism(m_SectorSupportPosition.getX() +
                     m_SectorSupportGeometry.getThickness() +
                     m_SectorSupportGeometry.getCorner4LX(), 0, 0, 1);
  Line2D line41Corner1B(m_SectorSupportPosition.getX() +
                        m_SectorSupportGeometry.getCornerX(), 0, 0, 1);
  Circle2D circleInnerOuter(0, 0, m_SectorSupportPosition.getInnerR());
  Circle2D circleInnerInner(0, 0, m_SectorSupportPosition.getInnerR() +
                            m_SectorSupportGeometry.getThickness());
  Circle2D circleOuterInner(0, 0, m_SectorSupportPosition.getOuterR() -
                            m_SectorSupportGeometry.getThickness());
  Circle2D circleOuterOuter(0, 0, m_SectorSupportPosition.getOuterR());
  HepGeom::Point3D<double> intersections[2];
  /* Corner 1. */
  p.setX(m_SectorSupportPosition.getX());
  p.setY(m_SectorSupportPosition.getOuterR() -
         m_SectorSupportGeometry.getDeltaLY());
  p.setZ(0);
  m_SectorSupportGeometry.setCorner1A(p);
  line41Corner1B.findIntersection(circleOuterOuter, intersections);
  m_SectorSupportGeometry.setCorner1B(intersections[1]);
  m_SectorSupportGeometry.setCornerAngle(
    atan2(m_SectorSupportGeometry.getCorner1B().y() -
          m_SectorSupportGeometry.getCorner1A().y(),
          m_SectorSupportGeometry.getCorner1B().x() -
          m_SectorSupportGeometry.getCorner1A().x()) * CLHEP::rad);
  p.setX(m_SectorSupportPosition.getX() +
         m_SectorSupportGeometry.getThickness());
  p.setY(m_SectorSupportGeometry.getCorner1A().y() -
         m_SectorSupportGeometry.getThickness() *
         (1.0 / cos(m_SectorSupportGeometry.getCornerAngle()) -
          tan(m_SectorSupportGeometry.getCornerAngle())));
  p.setZ(0);
  m_SectorSupportGeometry.setCorner1AInner(p);
  Line2D lineCorner1(m_SectorSupportGeometry.getCorner1AInner().x(),
                     m_SectorSupportGeometry.getCorner1AInner().y(),
                     m_SectorSupportGeometry.getCorner1B().x() -
                     m_SectorSupportGeometry.getCorner1A().x(),
                     m_SectorSupportGeometry.getCorner1B().y() -
                     m_SectorSupportGeometry.getCorner1A().y());
  lineCorner1.findIntersection(circleOuterInner, intersections);
  m_SectorSupportGeometry.setCorner1BInner(intersections[1]);
  /* Corner 2. */
  line23Inner.findIntersection(circleOuterInner, intersections);
  m_SectorSupportGeometry.setCorner2Inner(intersections[1]);
  /* Corner 3. */
  line23Outer.findIntersection(circleInnerOuter, intersections);
  m_SectorSupportGeometry.setCorner3(intersections[1]);
  line23Inner.findIntersection(circleInnerInner, intersections);
  m_SectorSupportGeometry.setCorner3Inner(intersections[1]);
  line23Prism.findIntersection(circleInnerInner, intersections);
  p.setX(intersections[1].x());
  p.setY(m_SectorSupportPosition.getY() +
         m_SectorSupportGeometry.getThickness());
  p.setZ(0);
  m_SectorSupportGeometry.setCorner3Prism(p);
  /* Corner 4. */
  line41Outer.findIntersection(circleInnerOuter, intersections);
  m_SectorSupportGeometry.setCorner4(intersections[1]);
  line41Inner.findIntersection(circleInnerInner, intersections);
  m_SectorSupportGeometry.setCorner4Inner(intersections[1]);
  line41Prism.findIntersection(circleInnerInner, intersections);
  p.setX(m_SectorSupportPosition.getX() +
         m_SectorSupportGeometry.getThickness());
  p.setY(intersections[1].y());
  p.setZ(0);
  m_SectorSupportGeometry.setCorner4Prism(p);
}

calculateShieldGeometry()

void calculateShieldGeometry ( )

private

Calculate shield geometry data.

Definition at line 385 of file GeometryData.cc.

{
  int i;
  double r, l, dx, dy, xCenter, yCenter;
  const double asqrt2 = 1.0 / sqrt(2.0);
  HepGeom::Point3D<double> points[8];
  const ShieldDetailGeometry* detailA = m_ShieldGeometry.getDetailA();
  const ShieldDetailGeometry* detailB = m_ShieldGeometry.getDetailB();
  const ShieldDetailGeometry* detailC = m_ShieldGeometry.getDetailC();
  const ShieldDetailGeometry* detailD = m_ShieldGeometry.getDetailD();
  r = m_SectorSupportPosition.getInnerR() +
      m_SectorSupportGeometry.getThickness();
  /* Detail A. */
  EKLMPointToCLHEP(detailA->getPoint(0), points[0]);
  EKLMPointToCLHEP(detailA->getPoint(1), points[1]);
  EKLMPointToCLHEP(detailA->getPoint(2), points[2]);
  EKLMPointToCLHEP(detailA->getPoint(3), points[3]);
  EKLMPointToCLHEP(detailA->getPoint(4), points[4]);
  points[5].setX(detailA->getLengthX());
  points[5].setY(detailA->getLengthY());
  points[5].setZ(0);
  EKLMPointToCLHEP(detailA->getPoint(5), points[6]);
  EKLMPointToCLHEP(detailA->getPoint(6), points[7]);
  l = 0.5 * (detailA->getLengthX() + detailB->getLengthX());
  xCenter = -asqrt2 * l;
  yCenter = asqrt2 * l;
  for (i = 0; i < 8; i++)
    points[i] = HepGeom::Translate3D(xCenter, yCenter, 0) *
                HepGeom::RotateZ3D(-45.0 * CLHEP::deg) *
                HepGeom::Translate3D(-detailA->getLengthX() / 2,
                                     -detailA->getLengthY() / 2, 0) *
                points[i];
  getDetailDxDy(points, 8, r, 1, 1, dx, dy);
  m_ShieldGeometry.setDetailACenter(xCenter + dx, yCenter + dy);
  /* Details B, D, E. */
  points[0].setX(0);
  points[0].setY(-detailD->getLengthY());
  points[0].setZ(0);
  points[1].setX(detailD->getLengthX());
  points[1].setY(0);
  points[1].setZ(0);
  points[2].setX(detailB->getLengthX() - detailD->getLengthX());
  points[2].setY(0);
  points[2].setZ(0);
  points[3].setX(detailB->getLengthX());
  points[3].setY(-detailD->getLengthY());
  points[3].setZ(0);
  EKLMPointToCLHEP(detailB->getPoint(0), points[4]);
  EKLMPointToCLHEP(detailB->getPoint(1), points[5]);
  EKLMPointToCLHEP(detailB->getPoint(2), points[6]);
  EKLMPointToCLHEP(detailB->getPoint(3), points[7]);
  /* Detail B center coordinates before its shift are (0, 0). */
  for (i = 0; i < 8; i++)
    points[i] = HepGeom::RotateZ3D(-45.0 * CLHEP::deg) *
                HepGeom::Translate3D(-detailB->getLengthX() / 2,
                                     -detailB->getLengthY() / 2, 0) *
                points[i];
  getDetailDxDy(points, 8, r, 1, 1, dx, dy);
  m_ShieldGeometry.setDetailBCenter(dx, dy);
  /* Detail C. */
  EKLMPointToCLHEP(detailC->getPoint(0), points[0]);
  EKLMPointToCLHEP(detailC->getPoint(1), points[1]);
  EKLMPointToCLHEP(detailC->getPoint(2), points[2]);
  EKLMPointToCLHEP(detailC->getPoint(3), points[3]);
  EKLMPointToCLHEP(detailC->getPoint(4), points[4]);
  points[5].setX(detailC->getLengthX());
  points[5].setY(detailC->getLengthY());
  points[5].setZ(0);
  EKLMPointToCLHEP(detailC->getPoint(5), points[6]);
  EKLMPointToCLHEP(detailC->getPoint(6), points[7]);
  l = 0.5 * (detailB->getLengthX() + detailC->getLengthX());
  xCenter = asqrt2 * l;
  yCenter = -asqrt2 * l;
  for (i = 0; i < 8; i++)
    points[i] = HepGeom::Translate3D(xCenter, yCenter, 0) *
                HepGeom::RotateZ3D(-45.0 * CLHEP::deg) *
                HepGeom::RotateY3D(180.0 * CLHEP::deg) *
                HepGeom::Translate3D(-detailC->getLengthX() / 2,
                                     -detailC->getLengthY() / 2, 0) *
                points[i];
  getDetailDxDy(points, 8, r, 1, 1, dx, dy);
  m_ShieldGeometry.setDetailCCenter(xCenter + dx, yCenter + dy);
}

checkDetectorLayer()

void checkDetectorLayer ( int  section, int  layer  ) const

inherited

Check if detector layer number is correct (fatal error if not).

Section number must be checked separately.

Parameters

[in]	section	Section number.
[in]	layer	Layer number.

Definition at line 312 of file EKLMGeometry.cc.

{
  /* cppcheck-suppress variableScope */
  const char* sectionName[2] = {"backward", "forward"};
  if (layer < 0 || layer > m_NDetectorLayers[section - 1])
    B2FATAL("Number of layer must be less from 1 to the number of "
            "detector layers in the " << sectionName[section - 1] << " section ("
            << m_NDetectorLayers[section - 1] << ").");
}

checkDetectorLayerNumber()

void checkDetectorLayerNumber ( int  section, int  layer  ) const

inherited

Check if number of detector layers is correct (fatal error if not).

Section number must be checked separately.

Parameters

[in]	section	Section number.
[in]	layer	Layer number.

Definition at line 301 of file EKLMGeometry.cc.

{
  /* cppcheck-suppress variableScope */
  const char* sectionName[2] = {"backward", "forward"};
  if (layer < 0 || layer > m_NLayers ||
      layer > m_ElementNumbers->getMaximalDetectorLayerNumber(section))
    B2FATAL("Number of detector layers in the " << sectionName[section - 1] <<
            " section must be from 0 to the number of layers ( " <<
            m_NLayers << ").");
}

checkSegmentSupport()

void checkSegmentSupport ( int  support ) const

inherited

Check if segment support number is correct (fatal error if not).

Parameters

[in]

support

Segment support element number.

Definition at line 322 of file EKLMGeometry.cc.

{
  if (support <= 0 || support > m_ElementNumbers->getMaximalSegmentNumber() + 1)
    B2FATAL("Number of segment support element must be from 1 to " <<
            m_ElementNumbers->getMaximalSegmentNumber() + 1 << ".");
}

checkStripSegment()

void checkStripSegment ( int  strip ) const

inherited

Check if number of strip in a segment is correct (fatal error if not).

Parameters

[in]

strip

Strip number.

Definition at line 329 of file EKLMGeometry.cc.

{
  if (strip <= 0 || strip > m_ElementNumbers->getNStripsSegment())
    B2FATAL("Number of strip in a segment must be from 1 to " <<
            m_ElementNumbers->getNStripsSegment() << ".");
}

fillStripIndexArrays()

void fillStripIndexArrays ( )

private

Fill strip index arrays.

Definition at line 202 of file GeometryData.cc.

{
  const char err[] = "Strip sorting algorithm error.";
  int i;
  double l;
  std::vector<double> strips;
  std::vector<double>::iterator it;
  std::map<double, int> mapLengthStrip;
  std::map<double, int> mapLengthStrip2;
  std::map<double, int>::iterator itm;
  for (i = 0; i < m_NStrips; i++) {
    strips.push_back(m_StripPosition[i].getLength());
    mapLengthStrip.insert(
      std::pair<double, int>(m_StripPosition[i].getLength(), i));
  }
  sort(strips.begin(), strips.end(), compareLength);
  l = strips[0];
  m_nStripDifferent = 1;
  for (it = strips.begin(); it != strips.end(); ++it) {
    if ((*it) != l) {
      l = (*it);
      m_nStripDifferent++;
    }
  }
  m_StripLenToAll = (int*)malloc(m_nStripDifferent * sizeof(int));
  if (m_StripLenToAll == nullptr)
    B2FATAL(c_MemErr);
  i = 0;
  l = strips[0];
  itm = mapLengthStrip.find(l);
  if (itm == mapLengthStrip.end())
    B2FATAL(err);
  m_StripLenToAll[i] = itm->second;
  mapLengthStrip2.insert(std::pair<double, int>(l, i));
  for (it = strips.begin(); it != strips.end(); ++it) {
    if ((*it) != l) {
      l = (*it);
      i++;
      itm = mapLengthStrip.find(l);
      if (itm == mapLengthStrip.end())
        B2FATAL(err);
      m_StripLenToAll[i] = itm->second;
      mapLengthStrip2.insert(std::pair<double, int>(l, i));
    }
  }
  m_StripAllToLen = (int*)malloc(m_NStrips * sizeof(int));
  if (m_StripAllToLen == nullptr)
    B2FATAL(c_MemErr);
  for (i = 0; i < m_NStrips; i++) {
    itm = mapLengthStrip2.find(m_StripPosition[i].getLength());
    if (itm == mapLengthStrip2.end())
      B2FATAL(err);
    m_StripAllToLen[i] = itm->second;
  }
}

getEndcapStructureGeometry()

const EndcapStructureGeometry * getEndcapStructureGeometry ( ) const

inlineinherited

Get endcap structure geometry data.

Definition at line 1783 of file EKLMGeometry.h.

    {
      return &m_EndcapStructureGeometry;
    }

getLayerPosition()

const ElementPosition * getLayerPosition ( ) const

inlineinherited

Get position data for layers.

Definition at line 1799 of file EKLMGeometry.h.

    {
      return &m_LayerPosition;
    }

getLayerShiftZ()

double getLayerShiftZ ( ) const

inlineinherited

Get Z distance between two layers.

Definition at line 1807 of file EKLMGeometry.h.

    {
      return m_LayerShiftZ;
    }

getLayerTransform()

void getLayerTransform	(	HepGeom::Transform3D *	t,
		int	n
	)		const

Get layer transformation.

Parameters

[out]	t	Transformation.
[in]	n	Number of layer.

Numbers start from 0.

Definition at line 662 of file GeometryData.cc.

{
  *t = HepGeom::Translate3D(0.0, 0.0, m_SectionPosition.getLength() / 2.0 -
                            (n + 1) * m_LayerShiftZ +
                            0.5 * m_LayerPosition.getLength());
}

getMaximalStripLength()

double getMaximalStripLength ( ) const

inline

Get maximal strip length.

Definition at line 79 of file GeometryData.h.

      {
        return m_StripPosition[m_StripLenToAll[
                                 m_nStripDifferent - 1]].getLength();
      }

getNDetectorLayers()

int getNDetectorLayers ( int  section ) const

inherited

Get number of detector layers.

Parameters

[in]

section

Section number (1 - backward, 2 - forward).

Definition at line 293 of file EKLMGeometry.cc.

{
  m_ElementNumbers->checkSection(section);
  return m_NDetectorLayers[section - 1];
}

getNLayers()

int getNLayers ( ) const

inlineinherited

Get number of layers.

Definition at line 1693 of file EKLMGeometry.h.

    {
      return m_NLayers;
    }

getNPlanes()

int getNPlanes ( ) const

inlineinherited

Get number of planes.

Definition at line 1715 of file EKLMGeometry.h.

    {
      return m_NPlanes;
    }

getNSections()

int getNSections ( ) const

inlineinherited

Get number of sections.

Definition at line 1685 of file EKLMGeometry.h.

    {
      return m_NSections;
    }

getNSectors()

int getNSectors ( ) const

inlineinherited

Get number of sectors.

Definition at line 1707 of file EKLMGeometry.h.

    {
      return m_NSectors;
    }

getNSegments()

int getNSegments ( ) const

inlineinherited

Get number of segments.

Definition at line 1723 of file EKLMGeometry.h.

    {
      return m_NSegments;
    }

getNSegmentSupportElementsSector()

int getNSegmentSupportElementsSector ( ) const

inlineinherited

Get number of segment support elements *in one sector).

Definition at line 1731 of file EKLMGeometry.h.

    {
      return m_NSegmentSupportElementsSector;
    }

getNStrips()

int getNStrips ( ) const

inlineinherited

Get number of strips.

Definition at line 1739 of file EKLMGeometry.h.

    {
      return m_NStrips;
    }

getNStripsDifferentLength()

int getNStripsDifferentLength ( ) const

inline

Get number of strips with different lengths.

Definition at line 88 of file GeometryData.h.

      {
        return m_nStripDifferent;
      }

getPlanePosition()

const ElementPosition * getPlanePosition ( ) const

inlineinherited

Get position data for planes.

Definition at line 1839 of file EKLMGeometry.h.

    {
      return &m_PlanePosition;
    }

getPlaneTransform()

void getPlaneTransform	(	HepGeom::Transform3D *	t,
		int	n
	)		const

Get plane transformation.

Parameters

[out]	t	Transformation.
[in]	n	Number of plane.

Numbers start from 0.

Definition at line 690 of file GeometryData.cc.

{
  if (n == 0)
    *t = HepGeom::Translate3D(m_PlanePosition.getX(), m_PlanePosition.getY(),
                              m_PlanePosition.getZ()) *
         HepGeom::Rotate3D(180. * CLHEP::deg,
                           HepGeom::Vector3D<double>(1., 1., 0.));
  else
    *t = HepGeom::Translate3D(m_PlanePosition.getX(), m_PlanePosition.getY(),
                              -m_PlanePosition.getZ());
}

getPlasticSheetGeometry()

const PlasticSheetGeometry * getPlasticSheetGeometry ( ) const

inlineinherited

Get plastic sheet geometry data.

Definition at line 1847 of file EKLMGeometry.h.

    {
      return &m_PlasticSheetGeometry;
    }

getSectionPosition()

const ElementPosition * getSectionPosition ( ) const

inlineinherited

Get position data for sections.

Definition at line 1791 of file EKLMGeometry.h.

    {
      return &m_SectionPosition;
    }

getSectionTransform()

void getSectionTransform	(	HepGeom::Transform3D *	t,
		int	n
	)		const

Get section transformation.

Parameters

[out]	t	Transformation.
[in]	n	Number of section.

Numbers start from 0.

Definition at line 650 of file GeometryData.cc.

{
  if (n == 0)
    *t = HepGeom::Translate3D(m_SectionPosition.getX(), m_SectionPosition.getY(),
                              -m_SectionPosition.getZ() + m_SolenoidZ);
  else
    *t = HepGeom::Translate3D(m_SectionPosition.getX(), m_SectionPosition.getY(),
                              m_SectionPosition.getZ() + m_SolenoidZ) *
         HepGeom::RotateY3D(180.*CLHEP::deg);
}

getSectorPosition()

const ElementPosition * getSectorPosition ( ) const

inlineinherited

Get position data for sectors.

Definition at line 1815 of file EKLMGeometry.h.

    {
      return &m_SectorPosition;
    }

getSectorSupportGeometry()

const SectorSupportGeometry * getSectorSupportGeometry ( ) const

inlineinherited

Get sector support geometry data.

Definition at line 1831 of file EKLMGeometry.h.

    {
      return &m_SectorSupportGeometry;
    }

getSectorSupportPosition()

const ElementPosition * getSectorSupportPosition ( ) const

inlineinherited

Get position data for sector support structure.

Definition at line 1823 of file EKLMGeometry.h.

    {
      return &m_SectorSupportPosition;
    }

getSectorTransform()

void getSectorTransform	(	HepGeom::Transform3D *	t,
		int	n
	)		const

Get sector transformation.

Parameters

[out]	t	Transformation.
[in]	n	Number of sector.

Numbers start from 0.

Definition at line 670 of file GeometryData.cc.

{
  switch (n) {
    case 0:
      *t = HepGeom::Translate3D(0., 0., 0.);
      break;
    case 1:
      *t = HepGeom::RotateY3D(180.0 * CLHEP::deg);
      break;
    case 2:
      *t = HepGeom::RotateZ3D(90.0 * CLHEP::deg) *
           HepGeom::RotateY3D(180.0 * CLHEP::deg);
      break;
    case 3:
      *t = HepGeom::RotateZ3D(-90.0 * CLHEP::deg);
      break;
  }
}

getSegmentSupportGeometry()

const SegmentSupportGeometry * getSegmentSupportGeometry ( ) const

inlineinherited

Get segment support geometry data.

Definition at line 1855 of file EKLMGeometry.h.

    {
      return &m_SegmentSupportGeometry;
    }

getSegmentSupportPosition()

const EKLMGeometry::SegmentSupportPosition * getSegmentSupportPosition ( int  plane, int  support  ) const

inherited

Get position data for segment support structure.

Parameters

[in]	plane	Plane number.
[in]	support	Segment support element number.

Definition at line 339 of file EKLMGeometry.cc.

{
  m_ElementNumbers->checkPlane(plane);
  checkSegmentSupport(support);
  return &m_SegmentSupportPosition[(plane - 1) * (m_NSegments + 1) +
                                               support - 1];
}

getSheetTransform()

void getSheetTransform	(	HepGeom::Transform3D *	t,
		int	n
	)		const

Get plastic sheet element transformation.

Parameters

[out]	t	Transformation.
[in]	n	Number of list.

Numbers start from 0. This function is intended for construction of the list solids, so the transformation does not include z shift.

Definition at line 710 of file GeometryData.cc.

{
  double y;
  y = m_StripPosition[n].getY();
  if (n % m_ElementNumbers->getNStripsSegment() == 0)
    y = y + 0.5 * m_PlasticSheetGeometry.getDeltaL();
  else if (n % m_ElementNumbers->getNStripsSegment() ==
           m_ElementNumbers->getNStripsSegment() - 1)
    y = y - 0.5 * m_PlasticSheetGeometry.getDeltaL();
  *t = HepGeom::Translate3D(m_StripPosition[n].getX(), y, 0.0);
}

getShieldGeometry()

const ShieldGeometry * getShieldGeometry ( ) const

inlineinherited

Get shield layer details geometry data.

Definition at line 1885 of file EKLMGeometry.h.

    {
      return &m_ShieldGeometry;
    }

getSolenoidZ()

double getSolenoidZ ( ) const

inlineinherited

Get solenoid center Z coordinate.

Definition at line 1775 of file EKLMGeometry.h.

    {
      return m_SolenoidZ;
    }

getStripGeometry()

const StripGeometry * getStripGeometry ( ) const

inlineinherited

Get strip geometry data.

Definition at line 1871 of file EKLMGeometry.h.

    {
      return &m_StripGeometry;
    }

getStripLength()

double getStripLength ( int  strip ) const

inline

Get strip length.

Parameters

[in]

strip

Strip number.

Definition at line 71 of file GeometryData.h.

      {
        return m_StripPosition[strip - 1].getLength();
      }

getStripLengthIndex()

int getStripLengthIndex ( int  positionIndex ) const

inline

Get index in length-based array.

Parameters

[in]

positionIndex

index in position-based array.

Definition at line 97 of file GeometryData.h.

      {
        return m_StripAllToLen[positionIndex];
      }

getStripPosition()

const EKLMGeometry::ElementPosition * getStripPosition ( int  strip ) const

inherited

Get position data for strips.

Parameters

[in]

strip

Strip number.

Definition at line 348 of file EKLMGeometry.cc.

{
  m_ElementNumbers->checkStrip(strip);
  return &m_StripPosition[strip - 1];
}

getStripPositionIndex()

int getStripPositionIndex ( int  lengthIndex ) const

inline

Get index in position-based array.

Parameters

[in]

lengthIndex

index in length-based array.

Definition at line 106 of file GeometryData.h.

      {
        return m_StripLenToAll[lengthIndex];
      }

getStripTransform()

void getStripTransform	(	HepGeom::Transform3D *	t,
		int	n
	)		const

Get strip transformation.

Parameters

[out]	t	Transformation.
[in]	n	Number of strip.

Numbers start from 0.

Definition at line 703 of file GeometryData.cc.

{
  *t = HepGeom::Translate3D(m_StripPosition[n].getX(),
                            m_StripPosition[n].getY(), 0.0);
}

hitInEKLM()

bool hitInEKLM

(

double

z

)

const

Check if z coordinate may be in EKLM.

Parameters

[in]

z

Z coordinate (cm).

Warning: this is not a complete check!

Definition at line 639 of file GeometryData.cc.

{
  double zMm;
  zMm = z / Unit::cm * CLHEP::cm;
  return (zMm > m_MinZForward) || (zMm < m_MaxZBackward);
}

initializeFromDatabase()

void initializeFromDatabase ( )

private

Initialize from database.

Definition at line 596 of file GeometryData.cc.

{
  DBObjPtr<EKLMGeometry> eklmGeometry;
  if (!eklmGeometry.isValid())
    B2FATAL("No EKLM geometry data in the database.");
  EKLMGeometry::operator=(*eklmGeometry);
  m_Geometry = new EKLMGeometry(*this);
}

initializeFromGearbox()

void initializeFromGearbox ( const GearDir *  gearDir )

private

Initialize from Gearbox (XML).

Parameters

[in]

gearDir

XML data directory.

Definition at line 477 of file GeometryData.cc.

{
  int i, j, k;
  std::string name;
  ShieldDetailGeometry shieldDetailGeometry;
  GearDir gd(*gearDir);
  gd.append("/EKLM");
  /* Beam-background study. */
  m_BeamBackgroundStudy = gd.getBool("BeamBackgroundStudy");
  /* Numbers of elements. */
  m_NSections = gd.getInt("NSections");
  m_ElementNumbers->checkSection(m_NSections);
  m_NLayers = gd.getInt("NLayers");
  m_ElementNumbers->checkLayer(m_NLayers);
  m_NDetectorLayers = new int[m_NSections];
  m_NDetectorLayers[0] = gd.getInt("NDetectorLayersBackward");
  checkDetectorLayerNumber(1, m_NDetectorLayers[0]);
  if (m_NSections == 2) {
    m_NDetectorLayers[1] = gd.getInt("NDetectorLayersForward");
    checkDetectorLayerNumber(2, m_NDetectorLayers[1]);
  }
  m_NSectors = gd.getInt("NSectors");
  m_ElementNumbers->checkSector(m_NSectors);
  m_NPlanes = gd.getInt("NPlanes");
  m_ElementNumbers->checkPlane(m_NPlanes);
  m_NSegments = gd.getInt("NSegments");
  m_ElementNumbers->checkSegment(m_NSegments);
  m_NSegmentSupportElementsSector = (m_NSegments + 1) * m_NPlanes;
  m_NStrips = gd.getInt("NStrips");
  m_ElementNumbers->checkStrip(m_NStrips);
  /* Geometry parameters. */
  m_SolenoidZ = gd.getLength("SolenoidZ") * CLHEP::cm;
  readEndcapStructureGeometry(gd);
  GearDir section(gd);
  section.append("/Section");
  readPositionData(&m_SectionPosition, &section);
  readSizeData(&m_SectionPosition, &section);
  GearDir layer(gd);
  layer.append("/Layer");
  readSizeData(&m_LayerPosition, &layer);
  m_LayerShiftZ = layer.getLength("ShiftZ") * CLHEP::cm;
  GearDir sector(gd);
  sector.append("/Sector");
  readSizeData(&m_SectorPosition, &sector);
  GearDir sectorSupport(gd);
  sectorSupport.append("/SectorSupport");
  readPositionData(&m_SectorSupportPosition, &sectorSupport);
  readSizeData(&m_SectorSupportPosition, &sectorSupport);
  readSectorSupportGeometry(&m_SectorSupportGeometry, &sectorSupport);
  GearDir plane(gd);
  plane.append("/Plane");
  readPositionData(&m_PlanePosition, &plane);
  readSizeData(&m_PlanePosition, &plane);
  GearDir plasticSheet(gd);
  plasticSheet.append("/PlasticSheet");
  m_PlasticSheetGeometry.setWidth(plasticSheet.getLength("Width") * CLHEP::cm);
  m_PlasticSheetGeometry.setDeltaL(plasticSheet.getLength("DeltaL") *
                                   CLHEP::cm);
  GearDir segmentSupport(gd);
  segmentSupport.append("/SegmentSupport");
  m_SegmentSupportGeometry.setTopWidth(
    segmentSupport.getLength("TopWidth") * CLHEP::cm);
  m_SegmentSupportGeometry.setTopThickness(
    segmentSupport.getLength("TopThickness") * CLHEP::cm);
  m_SegmentSupportGeometry.setMiddleWidth(
    segmentSupport.getLength("MiddleWidth") * CLHEP::cm);
  m_SegmentSupportGeometry.setMiddleThickness(
    segmentSupport.getLength("MiddleThickness") * CLHEP::cm);
  try {
    m_SegmentSupportPosition =
      new EKLMGeometry::SegmentSupportPosition[m_NSegmentSupportElementsSector];
  } catch (std::bad_alloc& ba) {
    B2FATAL(c_MemErr);
  }
  for (j = 0; j < m_NPlanes; j++) {
    for (i = 0; i <= m_NSegments; i++) {
      k = j * (m_NSegments + 1) + i;
      GearDir segmentSupport2(segmentSupport);
      name = "/SegmentSupportPlane[" + std::to_string(j + 1) + "]";
      segmentSupport2.append(name);
      name = "/SegmentSupport[" + std::to_string(i + 1) + "]";
      segmentSupport2.append(name);
      m_SegmentSupportPosition[k].setLength(
        segmentSupport2.getLength("Length") * CLHEP::cm);
      m_SegmentSupportPosition[k].setX(
        segmentSupport2.getLength("X") * CLHEP::cm);
      m_SegmentSupportPosition[k].setY(
        segmentSupport2.getLength("Y") * CLHEP::cm);
      m_SegmentSupportPosition[k].setZ(
        segmentSupport2.getLength("Z") * CLHEP::cm);
      m_SegmentSupportPosition[k].setDeltaLRight(
        segmentSupport2.getLength("DeltaLRight") * CLHEP::cm);
      m_SegmentSupportPosition[k].setDeltaLLeft(
        segmentSupport2.getLength("DeltaLLeft") * CLHEP::cm);
    }
  }
  readXMLDataStrips(gd);
  GearDir shield(gd);
  shield.append("/Shield");
  m_ShieldGeometry.setThickness(shield.getLength("Thickness") * CLHEP::cm);
  GearDir shieldDetailA(shield);
  shieldDetailA.append("/Detail[@id=\"A\"]");
  readShieldDetailGeometry(&shieldDetailGeometry, &shieldDetailA);
  m_ShieldGeometry.setDetailA(shieldDetailGeometry);
  GearDir shieldDetailB(shield);
  shieldDetailB.append("/Detail[@id=\"B\"]");
  readShieldDetailGeometry(&shieldDetailGeometry, &shieldDetailB);
  m_ShieldGeometry.setDetailB(shieldDetailGeometry);
  GearDir shieldDetailC(shield);
  shieldDetailC.append("/Detail[@id=\"C\"]");
  readShieldDetailGeometry(&shieldDetailGeometry, &shieldDetailC);
  m_ShieldGeometry.setDetailC(shieldDetailGeometry);
  GearDir shieldDetailD(shield);
  shieldDetailD.append("/Detail[@id=\"D\"]");
  readShieldDetailGeometry(&shieldDetailGeometry, &shieldDetailD);
  m_ShieldGeometry.setDetailD(shieldDetailGeometry);
  m_Geometry = new EKLMGeometry(*this);
}

Instance()

const EKLM::GeometryData & Instance ( enum DataSource  dataSource = c_Database, const GearDir *  gearDir = nullptr  )

static

Instantiation.

Parameters

[in]	dataSource	Data source.
[in]	gearDir	XML data directory.

Definition at line 33 of file GeometryData.cc.

{
  static EKLM::GeometryData gd(dataSource, gearDir);
  return gd;
}

readEndcapStructureGeometry()

void readEndcapStructureGeometry ( const GearDir &  gd )

private

Read section structure geometry data.

Parameters

[in]

gd

Current GearDir.

Definition at line 469 of file GeometryData.cc.

{
  GearDir d(gd);
  d.append("/EndcapStructure");
  m_EndcapStructureGeometry.setPhi(d.getAngle("Phi") * CLHEP::rad);
  m_EndcapStructureGeometry.setNSides(d.getInt("NSides"));
}

readXMLDataStrips()

void readXMLDataStrips ( const GearDir &  gd )

private

Read strip parameters from XML database.

Definition at line 258 of file GeometryData.cc.

{
  int i;
  std::string name;
  GearDir Strips(gd);
  Strips.append("/Strip");
  m_StripGeometry.setWidth(Strips.getLength("Width") * CLHEP::cm);
  m_StripGeometry.setThickness(Strips.getLength("Thickness") * CLHEP::cm);
  m_StripGeometry.setGrooveDepth(Strips.getLength("GrooveDepth") * CLHEP::cm);
  m_StripGeometry.setGrooveWidth(Strips.getLength("GrooveWidth") * CLHEP::cm);
  m_StripGeometry.setNoScintillationThickness(
    Strips.getLength("NoScintillationThickness") * CLHEP::cm);
  m_StripGeometry.setRSSSize(Strips.getLength("RSSSize") * CLHEP::cm);
  try {
    m_StripPosition = new EKLMGeometry::ElementPosition[m_NStrips];
  } catch (std::bad_alloc& ba) {
    B2FATAL(c_MemErr);
  }
  for (i = 0; i < m_NStrips; i++) {
    GearDir StripContent(Strips);
    name = "/Strip[" + std::to_string(i + 1) + "]";
    StripContent.append(name);
    m_StripPosition[i].setLength(StripContent.getLength("Length") * CLHEP::cm);
    m_StripPosition[i].setX(StripContent.getLength("X") * CLHEP::cm);
    m_StripPosition[i].setY(StripContent.getLength("Y") * CLHEP::cm);
    m_StripPosition[i].setZ(StripContent.getLength("Z") * CLHEP::cm);
  }
}

saveToDatabase()

void saveToDatabase

(

const IntervalOfValidity &

iov

)

const

Save geometry data to database.

Definition at line 634 of file GeometryData.cc.

{
  Database::Instance().storeData("EKLMGeometry", m_Geometry, iov);
}

Member Data Documentation

m_BeamBackgroundStudy

bool m_BeamBackgroundStudy = false

protectedinherited

ROOT streamer.

Whether to perform beam-background study.

Definition at line 1896 of file EKLMGeometry.h.

m_ElementNumbers

const EKLMElementNumbers* m_ElementNumbers

protectedinherited

Element numbers.

Definition at line 1893 of file EKLMGeometry.h.

m_EndcapStructureGeometry

EndcapStructureGeometry m_EndcapStructureGeometry

protectedinherited

Section structure geometry data.

Definition at line 1926 of file EKLMGeometry.h.

m_Geometry

EKLMGeometry* m_Geometry

private

Copy of data in this class used to write it to database.

Definition at line 227 of file GeometryData.h.

m_LayerPosition

ElementPosition m_LayerPosition

protectedinherited

Position data for layers.

Definition at line 1932 of file EKLMGeometry.h.

m_LayerShiftZ

double m_LayerShiftZ

protectedinherited

Z distance between two layers.

Definition at line 1935 of file EKLMGeometry.h.

m_MaxZBackward

double m_MaxZBackward

private

Maximal z coordinate of the backward section.

Definition at line 242 of file GeometryData.h.

m_MinZForward

double m_MinZForward

private

Minimal z coordinate of the forward section.

Definition at line 239 of file GeometryData.h.

m_NDetectorLayers

int* m_NDetectorLayers

protectedinherited

Number of detector layers.

Definition at line 1905 of file EKLMGeometry.h.

m_NLayers

int m_NLayers

protectedinherited

Number of layers in one section.

Definition at line 1902 of file EKLMGeometry.h.

m_NPlanes

int m_NPlanes

protectedinherited

Number of planes in one sector.

Definition at line 1911 of file EKLMGeometry.h.

m_NSections

int m_NSections

protectedinherited

Number of sections.

Definition at line 1899 of file EKLMGeometry.h.

m_NSectors

int m_NSectors

protectedinherited

Number of sectors in one layer.

Definition at line 1908 of file EKLMGeometry.h.

m_NSegments

int m_NSegments

protectedinherited

Number of segments in one plane.

Definition at line 1914 of file EKLMGeometry.h.

m_NSegmentSupportElementsSector

int m_NSegmentSupportElementsSector

protectedinherited

Number of segment support elements in one sector.

Definition at line 1917 of file EKLMGeometry.h.

m_nStripDifferent

int m_nStripDifferent

private

Number of strips with different lengths in one plane.

Definition at line 230 of file GeometryData.h.

m_NStrips

int m_NStrips

protectedinherited

Number of strips in one plane.

Definition at line 1920 of file EKLMGeometry.h.

m_PlanePosition

ElementPosition m_PlanePosition

protectedinherited

Position data for planes.

Definition at line 1947 of file EKLMGeometry.h.

m_PlasticSheetGeometry

PlasticSheetGeometry m_PlasticSheetGeometry

protectedinherited

Plastic sheet geometry data.

Definition at line 1950 of file EKLMGeometry.h.

m_SectionPosition

ElementPosition m_SectionPosition

protectedinherited

Position data for sections.

Definition at line 1929 of file EKLMGeometry.h.

m_SectorPosition

ElementPosition m_SectorPosition

protectedinherited

Position data for sectors.

Definition at line 1938 of file EKLMGeometry.h.

m_SectorSupportGeometry

SectorSupportGeometry m_SectorSupportGeometry

protectedinherited

Sector support geometry data.

Definition at line 1944 of file EKLMGeometry.h.

m_SectorSupportPosition

ElementPosition m_SectorSupportPosition

protectedinherited

Position data for sector support structure.

Definition at line 1941 of file EKLMGeometry.h.

m_SegmentSupportGeometry

SegmentSupportGeometry m_SegmentSupportGeometry

protectedinherited

Segment support geometry data.

Definition at line 1953 of file EKLMGeometry.h.

m_SegmentSupportPosition

SegmentSupportPosition* m_SegmentSupportPosition

protectedinherited

Position data for segment support structure.

Definition at line 1956 of file EKLMGeometry.h.

m_ShieldGeometry

ShieldGeometry m_ShieldGeometry

protectedinherited

Shield layer details geometry data.

Definition at line 1965 of file EKLMGeometry.h.

m_SolenoidZ

double m_SolenoidZ

protectedinherited

Solenoid center Z coordinate.

Definition at line 1923 of file EKLMGeometry.h.

m_StripAllToLen

int* m_StripAllToLen

private

Number of strip in length-based array.

Definition at line 233 of file GeometryData.h.

m_StripGeometry

StripGeometry m_StripGeometry

protectedinherited

Strip geometry data.

Definition at line 1959 of file EKLMGeometry.h.

m_StripLenToAll

int* m_StripLenToAll

private

Number of strip in position-based array.

Definition at line 236 of file GeometryData.h.

m_StripPosition

ElementPosition* m_StripPosition

protectedinherited

Position data for strips.

Definition at line 1962 of file EKLMGeometry.h.

---

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

• klm/eklm/geometry/include/GeometryData.h
• klm/eklm/geometry/src/GeometryData.cc