DetectorGeometry.h

Go to the documentation of this file.

 
#ifndef EVENTDISPLAY_DETECTORGEOMETRY_H
#define EVENTDISPLAY_DETECTORGEOMETRY_H
 
// STL
#include <cmath>     //sqrt
#include <iostream>  //cerr
#include <map>       //storage maps
#include <memory>    //unique_ptr
#include <utility>   //BoundingBox
#include <vector>    //BoundingBox
 
// LDMX Framework
#include "DetDescr/EcalHexReadout.h"
#include "DetDescr/EcalID.h"
#include "DetDescr/HcalID.h"
#include "Ecal/Event/EcalHit.h"
#include "Hcal/Event/HcalHit.h"
#include "SimCore/Event/SimTrackerHit.h"  //recoil hits
 
namespace eventdisplay {
 
typedef std::vector<std::pair<double, double> > BoundingBox;
 
struct HexPrism {
  double x;
  double y;
  double z;
  double height;
  double radius;
};
 
class DetectorGeometry {
 public:
  static const DetectorGeometry &getInstance() {
    static const DetectorGeometry DETECTOR_GEOMETRY;
    return DETECTOR_GEOMETRY;
  }
 
  BoundingBox getBoundingBox(const ldmx::HcalHit &hit) const;
 
  BoundingBox getBoundingBox(const std::vector<ldmx::HcalHit> &hitVec) const;
 
  BoundingBox getBoundingBox(ldmx::HcalID::HcalSection section) const;
 
  HexPrism getHexPrism(const ldmx::EcalID &id) const;
 
  HexPrism getHexTower(int towerIndex) const;
 
  double getRotAngle(int layerID, int moduleID) const;
 
  BoundingBox getBoundingBox(int layerID, int moduleID) const;
 
  BoundingBox getBoundingBox(const ldmx::SimTrackerHit &recoilHit) const;
 
 private:
  DetectorGeometry();
 
  // HCAL
 
  std::map<ldmx::HcalID::HcalSection, int> hcalNLayers_;
 
  std::map<ldmx::HcalID::HcalSection, int> hcalNStrips_;
 
  std::map<ldmx::HcalID::HcalSection, double> hcalLengthScint_;
 
  std::map<ldmx::HcalID::HcalSection, double> hcalZeroLayer_;
 
  std::map<ldmx::HcalID::HcalSection, double> hcalZeroStrip_;
 
  std::map<ldmx::HcalID::HcalSection, double> hcalLayerThickness_;
 
  int hcalParityVertical_;
 
  double hcalUncertaintyTimingPos_;
 
  double hcalThicknessScint_;
 
  double hcalWidthScint_;
 
  // ECAL
 
  double ecalSiThickness_;
 
  double ecalDepth_;
 
  double ecalZeroLayer_;
 
  std::unique_ptr<ldmx::EcalHexReadout> ecalHexReader_;
 
  // RECOIL TRACKER
 
  double recoilStereoStripLength_;
 
  double recoilStereoXWidth_;
 
  double recoilStereoYWidth_;
 
  double recoilStereoSeparation_;
 
  double recoilStereoAngle_;
 
  double recoilMonoStripLength_;
 
  double recoilMonoXWidth_;
 
  double recoilMonoYWidth_;
 
  double recoilMonoSeparation_;
 
  double recoilSensorThickness_;
 
  std::map<int, std::vector<double> > recoilModulePos_;
 
  std::map<int, double> recoilModuleAngle_;
};
 
}  // namespace eventdisplay
 
#endif /* EVENTDISPLAY_DETECTORGEOMETRY_H */