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