Crystal.h

Go to the documentation of this file.

#ifndef FastSimulation_GeometryTool_Crystal_h
#define FastSimulation_GeometryTool_Crystal_h

// Data Formats
#include "Math/GenVector/Plane3D.h"
#include "DataFormats/Math/interface/Vector3D.h"

// Unfortunately, GlobalPoints are also needed
#include "Geometry/CaloTopology/interface/CaloDirection.h"
#include "DataFormats/DetId/interface/DetId.h"

//FAMOS
#include "FastSimulation/CaloGeometryTools/interface/CrystalNeighbour.h"
#include "FastSimulation/CaloGeometryTools/interface/BaseCrystal.h"

#include <vector>

class DetId;

class Crystal
{

 public:

  typedef math::XYZVector XYZVector;
  typedef math::XYZVector XYZPoint;
  typedef ROOT::Math::Plane3D Plane3D;

  // side numbering
  //  enum CrystalSide{EAST=0,NORTH=1,WEST=2,SOUTH=3,FRONT=4,BACK=5};
  Crystal(){number_ = 0;};
  Crystal(const DetId&  cell,const BaseCrystal* bc=nullptr);

  inline const XYZPoint& getCorner(unsigned i) const { return myCrystal_->getCorner(i);};
  inline const XYZPoint& getCenter() const {return myCrystal_->getCenter();};
  inline const XYZPoint& getFrontCenter() const {return myCrystal_->getFrontCenter();};
  inline const XYZPoint & getBackCenter() const {return myCrystal_->getBackCenter();}
  inline const XYZVector& getFirstEdge() const {return  myCrystal_->getFirstEdge();}
  inline const XYZVector& getFifthEdge() const {return myCrystal_->getFifthEdge();}
  inline const DetId & getDetId() const {return cellid_;};
  inline const int getSubdetNumber() const {return myCrystal_->getSubdetNumber();}
  void print() const {return myCrystal_->print();}
  void getLateralEdges(unsigned i,XYZPoint& a,XYZPoint& b) const {myCrystal_->getLateralEdges(i,a,b); };
  void getFrontSide(XYZPoint &a,XYZPoint &b,XYZPoint &c,XYZPoint &d) const {myCrystal_->getFrontSide(a,b,c,d); }
  void getFrontSide(std::vector<XYZPoint>& corners) const {myCrystal_->getFrontSide(corners);}
  void getBackSide(XYZPoint &a,XYZPoint &b,XYZPoint &c,XYZPoint &d) const {myCrystal_->getBackSide(a,b,c,d);}
  void getBackSide(std::vector<XYZPoint>& corners) const {myCrystal_->getBackSide(corners);}
  void getLateralSide(unsigned i,XYZPoint &a,XYZPoint &b,XYZPoint &c,XYZPoint &d) const {myCrystal_->getLateralSide(i,a,b,c,d);}
  void getLateralSide(unsigned i,std::vector<XYZPoint>& corners) const {myCrystal_->getLateralSide(i,corners);}
  void getSide(const CaloDirection& side, XYZPoint &a,XYZPoint &b,XYZPoint &c,XYZPoint &d) const
  {myCrystal_->getSide(side,a,b,c,d);}
  void getSide(const CaloDirection& side, std::vector<XYZPoint>& corners) const
  {myCrystal_->getSide(side,corners);}


  const Plane3D& getFrontPlane() const {return myCrystal_->getFrontPlane();}
  const Plane3D& getBackPlane() const {return myCrystal_->getBackPlane();}
  const Plane3D& getLateralPlane(unsigned i) const {return myCrystal_->getLateralPlane(i);}
  const Plane3D& getPlane(const CaloDirection& side) const { return myCrystal_->getPlane(side);}

 inline const XYZVector& getLateralEdge(unsigned i) const {return myCrystal_->getLateralEdge(i);}

 inline const XYZVector& exitingNormal(const CaloDirection& side) const {return myCrystal_->exitingNormal(side);}

  static unsigned oppositeDirection(unsigned iside);

  void getDrawingCoordinates(std::vector<float> &x,std::vector<float> &y,std::vector<float> &z) const
  {
    myCrystal_->getDrawingCoordinates(x,y,z);
  }

  inline void setNumber(unsigned n) {number_=n;};

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

  inline CrystalNeighbour & crystalNeighbour(unsigned iq) { return  neighbours_[iq];}
  
  inline const XYZVector & getAxis() const { return myCrystal_->getAxis();}

  inline void setX0Back(double val) {X0back_=val;}

  inline double getX0Back() const { return X0back_;}

  ~Crystal(){;};

 private:
  unsigned number_;
  DetId cellid_ ;
  std::vector<CrystalNeighbour> neighbours_;
  double X0back_;
  const BaseCrystal * myCrystal_;

 public:
  class crystalEqual
    {
    public:
      crystalEqual(const DetId & cell):ref_(cell)
    {;};
      ~crystalEqual(){;};
      inline bool operator() (const Crystal& xtal) const
    {
      return (ref_==xtal.getDetId());
    }
    private:
      const DetId& ref_;
    };
};
#endif