CrystalPad.h

Go to the documentation of this file.

#ifndef FastSimulation_CaloGeometryTools_CrystalPad
#define FastSimulation_CaloGeometryTools_CrystalPad

#include "Geometry/CaloTopology/interface/CaloDirection.h"

#include "CLHEP/Vector/TwoVector.h"
#include "DataFormats/Math/interface/Vector3D.h"
//#include "Math/GenVector/Transform3D.h"
#include "FastSimulation/CaloGeometryTools/interface/Transform3DPJ.h"
//C++ headers
#include <vector>


class CrystalPad
{
 public:

  typedef math::XYZVector XYZVector;
  typedef math::XYZVector XYZPoint;
  typedef ROOT::Math::Transform3DPJ Transform3D;
  typedef ROOT::Math::Transform3DPJ::Point Point;


  CrystalPad() { dummy_ = true;};
  CrystalPad(unsigned number, const std::vector<CLHEP::Hep2Vector>& corners);
  CrystalPad(unsigned number, int onEcal, 
         const std::vector<XYZPoint>& corners, 
         const XYZPoint& origin, 
         const XYZVector& vec1, 
         const XYZVector& vec2);

  CrystalPad(unsigned number, 
         const std::vector<XYZPoint>& corners,
         const Transform3D&,
         double scaf=1.,bool bothdirections=false);

  CrystalPad(const CrystalPad& right);

  CrystalPad& operator = (const CrystalPad& rhs );

  ~CrystalPad(){;};
  
  bool inside(const CLHEP::Hep2Vector & point,bool debug=false) const;
  //  bool globalinside(XYZPoint) const;

  //  CLHEP::Hep2Vector localPoint(XYZPoint point) const;

  inline const std::vector<CLHEP::Hep2Vector> & getCorners() const {return corners_;}

  void resetCorners();

  void print() const;

  inline double survivalProbability() const { return survivalProbability_;};
  inline void setSurvivalProbability(double val) {survivalProbability_=val;};

  CLHEP::Hep2Vector& edge(unsigned iside,int n) ;

  CLHEP::Hep2Vector & edge(CaloDirection);

  inline unsigned getNumber() const{return number_;};

  /*
  inline XYZPoint originalCoordinates(CLHEP::Hep2Vector point) const
    {
      XYZPoint p(point.x(),point.y(),0.);
      return trans_.Inverse() * p;
    }
  */
  
  inline bool operator==(const CrystalPad& quad) const
    {
      //      std::cout << " First " << quad.getCellID() << " Second " << this->getCellID() << std::endl;
      return quad.getNumber()==this->getNumber();
    }

  inline bool operator<(const CrystalPad& quad) const
    {
      return (center_.mag()<quad.center().mag());
    }

  void extrems(double &xmin,double& xmax,double &ymin, double& ymax) const;

  inline const CLHEP::Hep2Vector& center() const {return center_;}

  void getDrawingCoordinates(std::vector<float>&x, std::vector<float>&y) const;   

  
 private:

  static std::vector<CLHEP::Hep2Vector> aVector;

  std::vector<CLHEP::Hep2Vector> corners_;
  std::vector<CLHEP::Hep2Vector> dir_;
  unsigned number_; 
  Transform3D trans_;
  ROOT::Math::Rotation3D rotation_;
  XYZVector translation_;
  double survivalProbability_;
  CLHEP::Hep2Vector center_;
  double epsilon_;
  bool dummy_;
  double yscalefactor_;

 public:
    class padEqual
      {
      public:
    padEqual(unsigned cell):ref_(cell) 
      {
        //      std::cout << " quadEqual " << ref_ << std::endl;
      };
    ~padEqual(){;};
    inline bool operator() (const CrystalPad & quad) const
      {
        return (ref_==quad.getNumber());
      }
      private:
    unsigned ref_;
      };
  


};

#include<iosfwd>
std::ostream& operator <<(std::ostream& o ,  CrystalPad & quad);

#endif