Implementation of HCal strip readout. More...

#include <HcalGeometry.h>

Public Types
enum class	ScintillatorOrientation { horizontal = 0 , vertical = 1 , depth = 2 }
	Encodes the orientation of a bar. More...

Public Member Functions
	~HcalGeometry ()=default
	Class destructor.

ScintillatorOrientation	getScintillatorOrientation (const ldmx::HcalID id) const

TVector3	getStripCenterPosition (ldmx::HcalID id) const
	Get a strip center position from a combined hcal ID.

std::map< ldmx::HcalID, TVector3 >	getStripPositionMap () const
	Get the strip position map.

bool	backLayerIsHorizontal (const int layer) const
	Check whether a given layer corresponds to a horizontal (scintillator length along the x-axis) or vertical layer in the back HCal.

double	getHalfTotalWidth (int isection, int layer=1) const
	Get the half total width of a layer for a given section(strip) for back(side) Hcal.

double	getScintillatorLength (ldmx::HcalID id) const
	Get the length of a scintillator bar.

double	getScintillatorWidth () const
	Get the scitillator width.

double	getScintillatorThickness () const
	Get the scitillator thickness.

int	getNumSections () const
	Get the number of sections.

int	getNumLayers (int isection) const
	Get the number of layers for that section.

int	getNumStrips (int isection, int layer=1) const
	Get the number of strips per layer for that section and layer.

int	getZeroStrip (int isection, int layer=1) const
	Get the location of the zeroStrip in a given section and layer.

double	getEcalDx () const
	Get the length of the Ecal in (x) for the side Hcal.

double	getEcalDy () const
	Get the length of the Ecal in (y) for the side Hcal.

bool	hasSide3DReadout () const
	Does the Side Hcal have 3D readout?

bool	isPrototype () const

std::vector< double >	rotateGlobalToLocalBarPosition (const std::vector< double > &globalPosition, const ldmx::HcalID &id) const
	Coordinates that are given by Geant4 are typically global.

Public Member Functions inherited from framework::ConditionsObject
	ConditionsObject (const std::string &name)
	Class constructor.

virtual	~ConditionsObject ()
	Destructor.

std::string	getName () const
	Get the name of this object.

Static Public Attributes
static constexpr const char *	CONDITIONS_OBJECT_NAME {"HcalGeometry"}
	Conditions object: The name of the python configuration calling this class (Hcal/python/HcalGeometry.py) needs to match the CONDITIONS_OBJECT_NAME exactly.

Private Member Functions
	HcalGeometry (const framework::config::Parameters &ps)
	Class constructor, for use only by the provider.

void	buildStripPositionMap ()
	Map builder of HcalID and position.

void	printPositionMap (int section) const
	Debugging utility, prints out the HcalID and corresponding value of all entries in the strip_position_map_ for a given section.

void	printPositionMap () const
	Debugging utility, prints out the HcalID and corresponding value of all entries in the strip_position_map_.

Private Attributes
int	verbose_ {0}
	Parameters that apply to all types of geometries Verbosity, not configurable but helpful if developing.

double	scint_thickness_
	Thickness of scintillator.

double	scint_width_
	Width of Scintillator Strip [mm].

std::vector< double >	zero_layer_
	Front of HCal relative to world geometry for each section [mm].

std::vector< double >	layer_thickness_
	Thickness of the layers in each section [mm].

std::vector< int >	num_layers_
	Number of layers in each section.

int	num_sections_
	Number of sections.

double	ecal_dx_
	Lenght of the Ecal (in x and y)

double	ecal_dy_

double	y_offset_

int	back_horizontal_parity_ {}

int	side_3d_readout_ {}

std::vector< std::vector< int > >	num_strips_
	Number of strips per layer in each section and each layer.

std::vector< std::vector< double > >	zero_strip_
	The plane of the zero'th strip of each section [mm].

std::vector< std::vector< double > >	half_total_width_
	Half Total Width of Strips [mm].

std::vector< std::vector< double > >	scint_length_
	Length of strips [mm].

bool	is_prototype_ {}

std::map< ldmx::HcalID, TVector3 >	strip_position_map_
	Map of the HcalID position of strip centers relative to world geometry.

Friends
class	hcal::HcalGeometryProvider

Detailed Description

Implementation of HCal strip readout.

Definition at line 34 of file HcalGeometry.h.

Member Enumeration Documentation

◆ ScintillatorOrientation

enum class ldmx::HcalGeometry::ScintillatorOrientation

strong

Encodes the orientation of a bar.

horizontal : The length of the bar is along the x-axis vertical : The length of the bar is along the y-axis depth : The length of the bar is along the z-axis

Definition at line 50 of file HcalGeometry.h.

                                     {
    horizontal = 0,
    vertical = 1,
    depth = 2
  };

Constructor & Destructor Documentation

◆ ~HcalGeometry()

ldmx::HcalGeometry::~HcalGeometry ( )

default

Class destructor.

Does nothing because the stl containers clean up automatically.

◆ HcalGeometry()

ldmx::HcalGeometry::HcalGeometry ( const framework::config::Parameters & ps )

private

Class constructor, for use only by the provider.

Parameters

ps	Parameters to configure the HcalGeometry

Definition at line 11 of file HcalGeometry.cxx.

    : framework::ConditionsObject(HcalGeometry::CONDITIONS_OBJECT_NAME) {
  scint_thickness_ = ps.getParameter<double>("scint_thickness");
  scint_width_ = ps.getParameter<double>("scint_width");
  zero_layer_ = ps.getParameter<std::vector<double>>("zero_layer");
  layer_thickness_ = ps.getParameter<std::vector<double>>("layer_thickness");
  num_layers_ = ps.getParameter<std::vector<int>>("num_layers");
  num_sections_ = ps.getParameter<int>("num_sections");
  ecal_dx_ = ps.getParameter<double>("ecal_dx");
  ecal_dy_ = ps.getParameter<double>("ecal_dy");
  verbose_ = ps.getParameter<int>("verbose");
  back_horizontal_parity_ = ps.getParameter<int>("back_horizontal_parity");
  side_3d_readout_ = ps.getParameter<int>("side_3d_readout");
  y_offset_ = ps.getParameter<double>("y_offset");
 
  auto detectors_valid =
      ps.getParameter<std::vector<std::string>>("detectors_valid");
  // If one of the strings in detectors_valid is "ldmx-hcal-prototype", we
  // will use prototype geometry initialization
  is_prototype_ = std::find_if(detectors_valid.cbegin(), detectors_valid.cend(),
                               [](const auto detector) {
                                 return detector.find("ldmx-hcal-prototype") !=
                                        std::string::npos;
                               }) != detectors_valid.cend();
 
  num_strips_ = ps.getParameter<std::vector<std::vector<int>>>("num_strips");
  half_total_width_ =
      ps.getParameter<std::vector<std::vector<double>>>("half_total_width");
  zero_strip_ = ps.getParameter<std::vector<std::vector<double>>>("zero_strip");
  scint_length_ =
      ps.getParameter<std::vector<std::vector<double>>>("scint_length");
 
  buildStripPositionMap();
 
  if (verbose_ > 0) {
    printPositionMap();
  }
}

References buildStripPositionMap(), ecal_dx_, framework::config::Parameters::getParameter(), half_total_width_, layer_thickness_, num_layers_, num_sections_, num_strips_, printPositionMap(), scint_length_, scint_thickness_, scint_width_, verbose_, zero_layer_, and zero_strip_.

Member Function Documentation

◆ backLayerIsHorizontal()

bool ldmx::HcalGeometry::backLayerIsHorizontal ( const int layer ) const

inline

Check whether a given layer corresponds to a horizontal (scintillator length along the x-axis) or vertical layer in the back HCal.

See the back_horizontal_parity_ member for details.

Definition at line 89 of file HcalGeometry.h.

                                                    {
    return layer % 2 == back_horizontal_parity_;
  }

◆ buildStripPositionMap()

void ldmx::HcalGeometry::buildStripPositionMap ( )

private

Map builder of HcalID and position.

To build the map we loop over the number of Hcal sections, layers and strips. The Hcal sections range from 0 to 4. (We hard-code the number of sections as seen in HcalID) The Hcal layers range from 1 to NumLayers_[section]. The Hcal strips range from 0 to NumStrips_[section].

Odd layers have horizontal strips. Even layers have vertical strips.

For back Hcal:

layers in z
strips occupy thickness of scintillator in z (e.g. 20mm)
strips orientation is in x(y) depending on back_horizontal parity

Now compute, y(x) position for horizontal(vertical) layers, relative to the center of detector. Strips enumeration starts from -y(-x) stripcenter will be large for +y(+x) and the half width of the strip needs to be subtracted The halfwidth of the scintillator is given by half_total_width_. The x(y) position is set to the center of the strip (0).

For side Hcal before 3D readout

layers in y(x)
all layers have strips in x(y) for top-bottom (left-right) sections
all layers have strips occupying width of scintillator in z (e.g. 50mm)

Definition at line 171 of file HcalGeometry.cxx.

                                         {
  // We hard-code the number of sections as seen in HcalID
  for (unsigned int section = 0; section < num_sections_; section++) {
    for (unsigned int layer = 1; layer <= num_layers_[section]; layer++) {
      for (unsigned int strip = 0; strip < getNumStrips(section, layer);
           strip++) {
        // initialize values
        double x{-99999}, y{-99999}, z{-99999};
 
        // get hcal section
        ldmx::HcalID::HcalSection hcalsection =
            (ldmx::HcalID::HcalSection)section;
 
        const ldmx::HcalID id{section, layer, strip};
        const auto orientation{getScintillatorOrientation(id)};
        // the center of a layer: (layer-1) * (layer_thickness) +
        // scint_thickness/2
        double layercenter =
            (layer - 1) * layer_thickness_.at(section) + 0.5 * scint_thickness_;
 
        // the center of a strip: (strip + 0.5) * (strip_dx)
        double stripcenter = (strip + 0.5) * scint_width_;
 
        if (hcalsection == ldmx::HcalID::HcalSection::BACK) {
          // z position: zero-layer(z) + layer_z + scint_thickness / 2
          z = zero_layer_.at(section) + layercenter;
 
          if (orientation == ScintillatorOrientation::horizontal) {
            y = stripcenter - getZeroStrip(section, layer);
            x = 0;
          } else {
            x = stripcenter - getZeroStrip(section, layer);
            y = 0;
          }
        } else {
          if (side_3d_readout_) {
            /*
             *
             * For 3D readout:
             * - odd layers have strips in z
             * - even layers have strips in x(y) for top-bottom (left-right)
             * sections
             * - odd layers have strips occupying width of scintillator in x(y)
             * - even layers have strips occupying width of scintillator in z
             *
             */
            switch (hcalsection) {
              case ldmx::HcalID::HcalSection::BACK:
                // Handled earlier in the code!
              case ldmx::HcalID::HcalSection::LEFT:
              case ldmx::HcalID::HcalSection::RIGHT:
                if (orientation == ScintillatorOrientation::vertical) {
                  x = zero_layer_[section] + 0.5 * scint_thickness_ +
                      (layer - 1) * layer_thickness_[section];
                  y = ecal_dy_ -
                      (getScintillatorLength({id.section(), 2, id.strip()}) -
                       getScintillatorLength(id)) /
                          2;
                  z = getZeroStrip(section, layer) +
                      (strip + 0.5) * getScintillatorWidth();
                } else if (orientation == ScintillatorOrientation::depth) {
                  x = zero_layer_[section] + 0.5 * scint_thickness_ +
                      layer_thickness_[section] * (layer - 1);
                  y = -ecal_dy_ / 2 + (strip + 0.5) * getScintillatorWidth();
                  z = getZeroStrip(section, layer + 1) +
                      getScintillatorLength(id) / 2;
                }
                if (section == ldmx::HcalID::HcalSection::LEFT) {
                  y *= -1;
                  x *= -1;
                }
                break;
 
              case ldmx::HcalID::HcalSection::BOTTOM:
              case ldmx::HcalID::HcalSection::TOP:
                if (orientation == ScintillatorOrientation::horizontal) {
                  //
                  // Second half of the expression is the difference between the
                  // longest strips (first module) and the current module.
                  //
                  // 22 mm extra for space for 1 absorber and one air box
                  x = -ecal_dx_ / 2 - 2 - 20 +
                      (getScintillatorLength({id.section(), 2, id.strip()}) -
                       getScintillatorLength(id)) /
                          2;
                  y = zero_layer_[section] + 0.5 * scint_thickness_ +
                      (layer - 1) * layer_thickness_[section];
                  z = getZeroStrip(section, layer) +
                      (strip + 0.5) * getScintillatorWidth();
                }
                if (orientation == ScintillatorOrientation::depth) {
                  x = (ecal_dx_ / 2) - (strip + 0.5) * getScintillatorWidth();
                  y = zero_layer_[section] + 0.5 * scint_thickness_ +
                      layer_thickness_[section] * (layer - 1);
                  z = getZeroStrip(section, layer + 1) +
                      getScintillatorLength(id) / 2;
                }
                if (section == ldmx::HcalID::HcalSection::BOTTOM) {
                  y *= -1;
                  x *= -1;
                }
                break;
            }
 
          } else {
            // z position: zero-strip(z) + strip_center(z)
            z = getZeroStrip(section, layer) + stripcenter;
            if (hcalsection == ldmx::HcalID::HcalSection::TOP or
                hcalsection == ldmx::HcalID::HcalSection::BOTTOM) {
              y = zero_layer_.at(section) + layercenter;
              x = getHalfTotalWidth(section, layer);
              if (hcalsection == ldmx::HcalID::HcalSection::BOTTOM) {
                y *= -1;
                x *= -1;
              }
 
            } else {
              x = zero_layer_.at(section) + layercenter;
              y = getHalfTotalWidth(section, layer);
              if (hcalsection == ldmx::HcalID::HcalSection::RIGHT) {
                x *= -1;
                y *= -1;
              }
            }
          }
        }
 
        y += y_offset_;
        TVector3 pos;
        pos.SetXYZ(x, y, z);
        strip_position_map_[ldmx::HcalID(section, layer, strip)] = pos;
      }  // loop over strips
    }    // loop over layers
  }      // loop over sections
}  // strip position map

References ecal_dx_, getHalfTotalWidth(), getNumStrips(), getScintillatorLength(), getScintillatorWidth(), getZeroStrip(), layer_thickness_, num_layers_, num_sections_, scint_thickness_, scint_width_, strip_position_map_, and zero_layer_.

Referenced by HcalGeometry().

◆ getEcalDx()

double ldmx::HcalGeometry::getEcalDx ( ) const

inline

Get the length of the Ecal in (x) for the side Hcal.

Definition at line 152 of file HcalGeometry.h.

152{ return ecal_dx_; }

References ecal_dx_.

◆ getEcalDy()

double ldmx::HcalGeometry::getEcalDy ( ) const

inline

Get the length of the Ecal in (y) for the side Hcal.

Definition at line 157 of file HcalGeometry.h.

157{ return ecal_dy_; }

◆ getHalfTotalWidth()

double ldmx::HcalGeometry::getHalfTotalWidth	(	int	isection,
		int	layer = `1`
	)		const

inline

Get the half total width of a layer for a given section(strip) for back(side) Hcal.

Parameters

section
layer

Returns: half total width [mm]

Definition at line 99 of file HcalGeometry.h.

                                                              {
    // Layer numbering starts at 1, but a vector is zero-indexed
    auto layer_index = layer - 1;
    return half_total_width_.at(isection).at(layer_index);
  }

References half_total_width_.

Referenced by buildStripPositionMap().

◆ getNumLayers()

int ldmx::HcalGeometry::getNumLayers ( int isection ) const

inline

Get the number of layers for that section.

Definition at line 131 of file HcalGeometry.h.

131{ return num_layers_.at(isection); }

References num_layers_.

◆ getNumSections()

int ldmx::HcalGeometry::getNumSections ( ) const

inline

Get the number of sections.

Definition at line 126 of file HcalGeometry.h.

126{ return num_sections_; }

References num_sections_.

◆ getNumStrips()

int ldmx::HcalGeometry::getNumStrips	(	int	isection,
		int	layer = `1`
	)		const

inline

Get the number of strips per layer for that section and layer.

Definition at line 136 of file HcalGeometry.h.

                                                      {
    auto layer_index = layer - 1;
    return num_strips_.at(isection).at(layer_index);
  }

References num_strips_.

Referenced by buildStripPositionMap(), and printPositionMap().

◆ getScintillatorLength()

double ldmx::HcalGeometry::getScintillatorLength ( ldmx::HcalID id ) const

inline

Get the length of a scintillator bar.

Parameters

id	The HcalID of the bar

Returns: The length of the bar with ID id [mm]

Definition at line 110 of file HcalGeometry.h.

                                                    {
    return scint_length_.at(id.section()).at(id.layer() - 1);
  }

References scint_length_.

Referenced by buildStripPositionMap().

◆ getScintillatorOrientation()

HcalGeometry::ScintillatorOrientation ldmx::HcalGeometry::getScintillatorOrientation ( const ldmx::HcalID id ) const

Definition at line 100 of file HcalGeometry.cxx.

                               {
  if (hasSide3DReadout()) {
    // v14 or later detector
    switch (id.section()) {
      case ldmx::HcalID::HcalSection::TOP:
      case ldmx::HcalID::HcalSection::BOTTOM:
        // Odd layers are in z/depth direction, even are in the x/horizontal
        // direction
        return id.layer() % 2 == 0 ? ScintillatorOrientation::horizontal
                                   : ScintillatorOrientation::depth;
 
      case ldmx::HcalID::HcalSection::LEFT:
      case ldmx::HcalID::HcalSection::RIGHT:
        // Odd layers are in the z/depth direction, even are in the y/vertical
        // direction
        return id.layer() % 2 == 0 ? ScintillatorOrientation::vertical
                                   : ScintillatorOrientation::depth;
      case ldmx::HcalID::HcalSection::BACK:
        // Configurable
        return id.layer() % 2 == back_horizontal_parity_
                   ? ScintillatorOrientation::horizontal
                   : ScintillatorOrientation::vertical;
    }  // V14 or later detector
  }
  if (isPrototype()) {
    // The prototype only has the back section. However, the orientation
    // depends on the configuration so we delegate to the
    // back_horizontal_parity parameter
    return id.layer() % 2 == back_horizontal_parity_
               ? ScintillatorOrientation::horizontal
               : ScintillatorOrientation::vertical;
  }  // Prototype detector
  // v13/v12
  switch (id.section()) {
    // For the v13 side hcal, the bars in each section have the same
    // orientation
    case ldmx::HcalID::HcalSection::TOP:
    case ldmx::HcalID::HcalSection::BOTTOM:
      return ScintillatorOrientation::horizontal;
    case ldmx::HcalID::HcalSection::LEFT:
    case ldmx::HcalID::HcalSection::RIGHT:
      return ScintillatorOrientation::vertical;
    case ldmx::HcalID::HcalSection::BACK:
      // Configurable
      return id.layer() % 2 == back_horizontal_parity_
                 ? ScintillatorOrientation::horizontal
                 : ScintillatorOrientation::vertical;
  }  // v13/v12 detector
  // Can only reach this part if we somehow didn't match any of the options
  // above. This could happen if someone introduces a new geometry but doesn't
  // patch this part.
  EXCEPTION_RAISE("InvalidRotation",
                  "Attempted to rotate into an invalid "
                  "orientation for a scintillator bar!");
}

◆ getScintillatorThickness()

double ldmx::HcalGeometry::getScintillatorThickness ( ) const

inline

Get the scitillator thickness.

Definition at line 121 of file HcalGeometry.h.

121{ return scint_thickness_; }

References scint_thickness_.

◆ getScintillatorWidth()

double ldmx::HcalGeometry::getScintillatorWidth ( ) const

inline

Get the scitillator width.

Definition at line 116 of file HcalGeometry.h.

116{ return scint_width_; }

References scint_width_.

Referenced by buildStripPositionMap().

◆ getStripCenterPosition()

TVector3 ldmx::HcalGeometry::getStripCenterPosition ( ldmx::HcalID id ) const

inline

Get a strip center position from a combined hcal ID.

Exceptions

std::out_of_range if HcalID is not on map.

Parameters

HcalID

Returns: A TVector3 with the X, Y and Z position of the center of the bar.

Definition at line 74 of file HcalGeometry.h.

                                                       {
    return strip_position_map_.at(id);
  }

References strip_position_map_.

Referenced by printPositionMap().

◆ getStripPositionMap()

std::map< ldmx::HcalID, TVector3 > ldmx::HcalGeometry::getStripPositionMap ( ) const

inline

Get the strip position map.

Definition at line 81 of file HcalGeometry.h.

                                                           {
    return strip_position_map_;
  }

References strip_position_map_.

◆ getZeroStrip()

int ldmx::HcalGeometry::getZeroStrip	(	int	isection,
		int	layer = `1`
	)		const

inline

Get the location of the zeroStrip in a given section and layer.

Definition at line 144 of file HcalGeometry.h.

                                                      {
    auto layer_index = layer - 1;
    return zero_strip_.at(isection).at(layer_index);
  }

References zero_strip_.

Referenced by buildStripPositionMap().

◆ hasSide3DReadout()

bool ldmx::HcalGeometry::hasSide3DReadout ( ) const

inline

Does the Side Hcal have 3D readout?

In other words, does the side hcal layers alter in scintillator direction (z vs x/y).

Definition at line 165 of file HcalGeometry.h.

165{ return side_3d_readout_; }

◆ isPrototype()

bool ldmx::HcalGeometry::isPrototype ( ) const

inline

Definition at line 171 of file HcalGeometry.h.

171{ return is_prototype_; }

◆ printPositionMap() [1/2]

void ldmx::HcalGeometry::printPositionMap ( ) const

inlineprivate

Debugging utility, prints out the HcalID and corresponding value of all entries in the strip_position_map_.

For printing only one of the sections, see the overloaded version of this function taking a section parameter.

Definition at line 231 of file HcalGeometry.h.

                                {
    for (int section = 0; section < num_sections_; ++section) {
      printPositionMap(section);
    }
  }

References num_sections_, and printPositionMap().

Referenced by HcalGeometry(), and printPositionMap().

◆ printPositionMap() [2/2]

void ldmx::HcalGeometry::printPositionMap ( int section ) const

private

Debugging utility, prints out the HcalID and corresponding value of all entries in the strip_position_map_ for a given section.

Parameters

section The section number to print, see HcalID for details.

Definition at line 156 of file HcalGeometry.cxx.

                                                     {
  // Note that layer numbering starts at 1 rather than 0
  for (int layer = 1; layer <= num_layers_[section]; ++layer) {
    for (int strip = 0; strip < getNumStrips(section, layer); ++strip) {
      HcalID id(section, layer, strip);
      auto centerPosition = getStripCenterPosition(id);
      auto x = centerPosition[0];
      auto y = centerPosition[1];
      auto z = centerPosition[2];
      std::cout << id << ": Center position: (" << x << ", " << y << ", " << z
                << ")\n";
    }
  }
}

References getNumStrips(), getStripCenterPosition(), and num_layers_.

◆ rotateGlobalToLocalBarPosition()

std::vector< double > ldmx::HcalGeometry::rotateGlobalToLocalBarPosition	(	const std::vector< double > &	globalPosition,
		const ldmx::HcalID &	id
	)		const

Coordinates that are given by Geant4 are typically global.

These can be transformed into corresponding local coordinates of a volume with a TopTransform. However, a TopTransform translates to the local bar but does not do the rotation (this is because we don't do a rotation when placing the bars in the GDML). This is used primarily for recording pre and post step positions in local coordinates of the volume in HcalSD.

the logic below does the rotation to the local coordiates where x : short transverse side of bar y : long transverse side of bar z : along length of bar

Note: This logic only applies to the v14 and prototype detector; however, support for v12 is not broken because no studies using these pre/post step positions have been (or should be) done with the v12 detector. /; : The native position type in Geant4 is typically a G4ThreeVector, which is a typedef for a CLHEP::Hep3Vector. However, DetDescr currently does not have a dependency on Geant4/CLHEP so we are taking the position as a vector of floats (which is what is used by SimCalorimeterHit)

Definition at line 49 of file HcalGeometry.cxx.

                                                                         {
  const auto orientation{getScintillatorOrientation(id)};
  switch (id.section()) {
    case ldmx::HcalID::HcalSection::BACK:
      switch (orientation) {
        case ScintillatorOrientation::horizontal:
          return {globalPosition[2], globalPosition[1], globalPosition[0]};
        case ScintillatorOrientation::vertical:
          return {globalPosition[2], globalPosition[0], globalPosition[1]};
        default:  // Should not be possible with current geometries
          EXCEPTION_RAISE("InvalidRotation",
                          "Attempted to rotate into an invalid "
                          "orientation for a scintillator bar!");
      }
    case ldmx::HcalID::HcalSection::TOP:
      [[fallthrough]];
    case ldmx::HcalID::HcalSection::BOTTOM:
      switch (orientation) {
        case ScintillatorOrientation::horizontal:
          return {globalPosition[1], globalPosition[2], globalPosition[0]};
        case ScintillatorOrientation::depth:
          return {globalPosition[1], globalPosition[0], globalPosition[2]};
        default:  // Should not be possible with current geometries
          EXCEPTION_RAISE("InvalidRotation",
                          "Attempted to rotate into an invalid "
                          "orientation for a scintillator bar!");
      }
    case ldmx::HcalID::HcalSection::LEFT:
      [[fallthrough]];
    case ldmx::HcalID::HcalSection::RIGHT:
      switch (orientation) {
        case ScintillatorOrientation::vertical:
          return {globalPosition[0], globalPosition[2], globalPosition[1]};
        case ScintillatorOrientation::depth:
          return globalPosition;
        default:  // Should not be possible with current geometries
          EXCEPTION_RAISE("InvalidRotation",
                          "Attempted to rotate into an invalid "
                          "orientation for a scintillator bar!");
      }
    default:
      // Can only reach this part if we somehow didn't match any of the options
      // above. This could happen if someone introduces a new geometry but
      // doesn't patch this part.
      EXCEPTION_RAISE("InvalidRotation",
                      "Attempted to rotate into an invalid "
                      "orientation for a scintillator bar!");
  }
}

Friends And Related Symbol Documentation

◆ hcal::HcalGeometryProvider

friend class hcal::HcalGeometryProvider

friend

Definition at line 205 of file HcalGeometry.h.

Member Data Documentation

◆ back_horizontal_parity_

int ldmx::HcalGeometry::back_horizontal_parity_ {}

private

Definition at line 270 of file HcalGeometry.h.

270{};

◆ CONDITIONS_OBJECT_NAME

constexpr const char* ldmx::HcalGeometry::CONDITIONS_OBJECT_NAME {"HcalGeometry"}

staticconstexpr

Conditions object: The name of the python configuration calling this class (Hcal/python/HcalGeometry.py) needs to match the CONDITIONS_OBJECT_NAME exactly.

Definition at line 42 of file HcalGeometry.h.

42{"HcalGeometry"};

Referenced by simcore::HcalSD::decodeCopyNumber(), hcal::HcalTriggerGeometryProvider::getCondition(), simcore::HcalSD::ProcessHits(), hcal::HcalClusterProducer::produce(), hcal::HcalDigiProducer::produce(), hcal::HcalDoubleEndRecProducer::produce(), hcal::HcalRecProducer::produce(), and hcal::HcalSingleEndRecProducer::produce().

◆ ecal_dx_

double ldmx::HcalGeometry::ecal_dx_

private

Lenght of the Ecal (in x and y)

Definition at line 261 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), getEcalDx(), and HcalGeometry().

◆ ecal_dy_

double ldmx::HcalGeometry::ecal_dy_

private

Definition at line 262 of file HcalGeometry.h.

◆ half_total_width_

std::vector<std::vector<double> > ldmx::HcalGeometry::half_total_width_

private

Half Total Width of Strips [mm].

Definition at line 280 of file HcalGeometry.h.

Referenced by getHalfTotalWidth(), and HcalGeometry().

◆ is_prototype_

bool ldmx::HcalGeometry::is_prototype_ {}

private

Definition at line 284 of file HcalGeometry.h.

284{};

◆ layer_thickness_

std::vector<double> ldmx::HcalGeometry::layer_thickness_

private

Thickness of the layers in each section [mm].

Definition at line 252 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), and HcalGeometry().

◆ num_layers_

std::vector<int> ldmx::HcalGeometry::num_layers_

private

Number of layers in each section.

Definition at line 255 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), getNumLayers(), HcalGeometry(), and printPositionMap().

◆ num_sections_

int ldmx::HcalGeometry::num_sections_

private

Number of sections.

Definition at line 258 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), getNumSections(), HcalGeometry(), and printPositionMap().

◆ num_strips_

std::vector<std::vector<int> > ldmx::HcalGeometry::num_strips_

private

Number of strips per layer in each section and each layer.

Definition at line 276 of file HcalGeometry.h.

Referenced by getNumStrips(), and HcalGeometry().

◆ scint_length_

std::vector<std::vector<double> > ldmx::HcalGeometry::scint_length_

private

Length of strips [mm].

Definition at line 282 of file HcalGeometry.h.

Referenced by getScintillatorLength(), and HcalGeometry().

◆ scint_thickness_

double ldmx::HcalGeometry::scint_thickness_

private

Thickness of scintillator.

Definition at line 243 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), getScintillatorThickness(), and HcalGeometry().

◆ scint_width_

double ldmx::HcalGeometry::scint_width_

private

Width of Scintillator Strip [mm].

Definition at line 246 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), getScintillatorWidth(), and HcalGeometry().

◆ side_3d_readout_

int ldmx::HcalGeometry::side_3d_readout_ {}

private

Definition at line 273 of file HcalGeometry.h.

273{};

◆ strip_position_map_

std::map<ldmx::HcalID, TVector3> ldmx::HcalGeometry::strip_position_map_

private

Map of the HcalID position of strip centers relative to world geometry.

The map is not configurable and is calculated by buildStripPositionMap().

Definition at line 290 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), getStripCenterPosition(), and getStripPositionMap().

◆ verbose_

int ldmx::HcalGeometry::verbose_ {0}

private

Parameters that apply to all types of geometries Verbosity, not configurable but helpful if developing.

Definition at line 240 of file HcalGeometry.h.

240{0};

Referenced by HcalGeometry().

◆ y_offset_

double ldmx::HcalGeometry::y_offset_

private

Definition at line 265 of file HcalGeometry.h.

◆ zero_layer_

std::vector<double> ldmx::HcalGeometry::zero_layer_

private

Front of HCal relative to world geometry for each section [mm].

Definition at line 249 of file HcalGeometry.h.

Referenced by buildStripPositionMap(), and HcalGeometry().

◆ zero_strip_

std::vector<std::vector<double> > ldmx::HcalGeometry::zero_strip_

private

The plane of the zero'th strip of each section [mm].

Definition at line 278 of file HcalGeometry.h.

Referenced by getZeroStrip(), and HcalGeometry().

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

DetDescr/include/DetDescr/HcalGeometry.h
DetDescr/src/DetDescr/HcalGeometry.cxx

Public Types

Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Friends

Detailed Description

Member Enumeration Documentation

◆ ScintillatorOrientation

Constructor & Destructor Documentation

◆ ~HcalGeometry()

◆ HcalGeometry()

Member Function Documentation

◆ backLayerIsHorizontal()

◆ buildStripPositionMap()

◆ getEcalDx()

◆ getEcalDy()

◆ getHalfTotalWidth()

◆ getNumLayers()

◆ getNumSections()

◆ getNumStrips()

◆ getScintillatorLength()

◆ getScintillatorOrientation()

◆ getScintillatorThickness()

◆ getScintillatorWidth()

◆ getStripCenterPosition()

◆ getStripPositionMap()

◆ getZeroStrip()

◆ hasSide3DReadout()

◆ isPrototype()

◆ printPositionMap() [1/2]

◆ printPositionMap() [2/2]

◆ rotateGlobalToLocalBarPosition()

Friends And Related Symbol Documentation

◆ hcal::HcalGeometryProvider

Member Data Documentation

◆ back_horizontal_parity_

◆ CONDITIONS_OBJECT_NAME

◆ ecal_dx_

◆ ecal_dy_

◆ half_total_width_

◆ is_prototype_

◆ layer_thickness_

◆ num_layers_

◆ num_sections_

◆ num_strips_

◆ scint_length_

◆ scint_thickness_

◆ scint_width_

◆ side_3d_readout_

◆ strip_position_map_

◆ verbose_

◆ y_offset_

◆ zero_layer_

◆ zero_strip_