CaloHitMaker.h

Go to the documentation of this file.

#ifndef CaloHitMaker_h
#define CaloHitMaker_h

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/Math/interface/Vector3D.h"
#include "Math/GenVector/Plane3D.h"
#include "SimG4CMS/Calo/interface/CaloHitID.h"

//CLHEP headers
//#include "CLHEP/Geometry/Point3D.h"
//#include "CLHEP/Geometry/Plane3D.h"


//STL headers
#include <string>
#include <map>

class CaloGeometryHelper;
class CalorimeterProperties;

class CaloHitMaker
{
 public:

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

  CaloHitMaker(const CaloGeometryHelper * calo,DetId::Detector det,int subdetn,int cal,unsigned sht=0);
  virtual ~CaloHitMaker(){;}
  
  virtual bool addHit(double r,double phi,unsigned layer=0)=0;
  virtual void setSpotEnergy(double e)=0;
  virtual const std::map<CaloHitID,float>& getHits()=0; 

  const CaloGeometryHelper * getCalorimeter() const 
    {
      //      std::cout << "CaloHitMaker is returning myCalorimeter " << myCalorimeter << std::endl;
      return myCalorimeter;
    }

 protected:
  static XYZPoint intersect(const Plane3D& p,const XYZPoint& a,const XYZPoint& b,double& t,bool segment,bool debug=false);

  const CaloGeometryHelper * myCalorimeter;    
  const CalorimeterProperties * theCaloProperties;
  double moliereRadius;
  double interactionLength;
  double spotEnergy;

  bool EMSHOWER;
  bool HADSHOWER;
  bool MIP;

 private:
  DetId::Detector base_;
  int subdetn_;
  int onCal_;


 protected:
  unsigned showerType_;
  std::map<CaloHitID,float> hitMap_;
  
};

#endif