CaloSegment.h

Go to the documentation of this file.

#ifndef CaloSegment_h
#define CaloSegment_h

//FAMOS headers
#include "FastSimulation/CaloGeometryTools/interface/CaloPoint.h"

#include <string>

class CaloGeometryHelper;
class CaloSegment {
public:
  typedef math::XYZVector XYZVector;
  typedef math::XYZVector XYZPoint;

  enum Material { PbWO4 = 0, CRACK = 1, GAP = 2, PS = 3, HCAL = 4, ECALHCALGAP = 5, PSEEGAP = 6 };

  CaloSegment(const CaloPoint& in,
              const CaloPoint& out,
              double si,
              double siX0,
              double liX0,
              Material mat,
              const CaloGeometryHelper*);
  ~CaloSegment() { ; }
  inline double sEntrance() const { return sentrance_; };
  inline double sExit() const { return sexit_; };
  inline double sX0Entrance() const { return sX0entrance_; };
  inline double sX0Exit() const { return sX0exit_; };
  inline double sL0Entrance() const { return sL0entrance_; };
  inline double sL0Exit() const { return sL0exit_; };
  inline double length() const { return length_; };
  inline double X0length() const { return X0length_; };
  inline double L0length() const { return L0length_; };
  inline const CaloPoint& entrance() const { return entrance_; };
  inline const CaloPoint& exit() const { return exit_; };

  inline bool operator<(const CaloSegment& s) const { return sentrance_ < s.sEntrance() && sexit_ < sExit(); }
  inline Material material() const { return material_; };
  inline DetId::Detector whichDetector() const { return detector_; };
  XYZPoint positionAtDepthincm(double depth) const;
  XYZPoint positionAtDepthinX0(double depth) const;
  XYZPoint positionAtDepthinL0(double depth) const;

  double x0FromCm(double cm) const;

private:
  //  static ECALProperties myCaloProperties;
  CaloPoint entrance_;
  CaloPoint exit_;
  double sentrance_;
  double sexit_;
  double sX0entrance_;
  double sX0exit_;
  double length_;
  double X0length_;
  double sL0entrance_;
  double sL0exit_;
  double L0length_;
  Material material_;
  DetId::Detector detector_;

public:
  class inX0Segment {
  public:
    //      inSegment(const CaloSegment & ref):segment_(ref){;};
    inX0Segment(double depth)
        : ref_(depth) {
            //std::cout << "inSegment " << std::endl;
          };
    ~inX0Segment() { ; };
    // in X0 !!!
    //      bool operator() (double value) const
    //  {
    //    return (value>segment_.sX0Entrance()&&value<segment_.sX0Exit());
    //  }
    bool operator()(const CaloSegment& segment) const {
      return (ref_ > segment.sX0Entrance() && ref_ < segment.sX0Exit());
    }

  private:
    //      const CaloSegment & segment_;
    double ref_;
  };

  class inL0Segment {
  public:
    //      inSegment(const CaloSegment & ref):segment_(ref){;};
    inL0Segment(double depth)
        : ref_(depth) {
            //std::cout << "inSegment " << std::endl;
          };
    ~inL0Segment() { ; };
    // in X0 !!!
    //      bool operator() (double value) const
    //  {
    //    return (value>segment_.sX0Entrance()&&value<segment_.sX0Exit());
    //  }
    bool operator()(const CaloSegment& segment) const {
      return (ref_ > segment.sL0Entrance() && ref_ < segment.sL0Exit());
    }

  private:
    //      const CaloSegment & segment_;
    double ref_;
  };

  class inSegment {
  public:
    //      inSegment(const CaloSegment & ref):segment_(ref){;};
    inSegment(double depth)
        : ref_(depth) {
            //std::cout << "inSegment " << std::endl;
          };
    ~inSegment() { ; };
    // in X0 !!!
    //      bool operator() (double value) const
    //  {
    //    return (value>segment_.sX0Entrance()&&value<segment_.sX0Exit());
    //  }
    bool operator()(const CaloSegment& segment) const {
      //      std::cout << " Entrance " << segment.sEntrance() << " Exit " << segment.sExit() << " " << ref_ << " " << segment.whichDetector() << std::endl;
      return (ref_ > segment.sEntrance() && ref_ < segment.sExit());
    }

  private:
    //      const CaloSegment & segment_;
    double ref_;
  };
};
#include <iosfwd>
std::ostream& operator<<(std::ostream& o, const CaloSegment& cid);

#endif