HcalHitMaker.h

Go to the documentation of this file.

#ifndef HcalHitMaker_h
#define HcalHitMaker_h

#include "FastSimulation/Event/interface/FSimTrack.h"
#include "FastSimulation/CaloHitMakers/interface/CaloHitMaker.h"
#include "FastSimulation/CaloHitMakers/interface/EcalHitMaker.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"

//#include "Math/GenVector/Transform3D.h"
#include "FastSimulation/CaloGeometryTools/interface/Transform3DPJ.h"

//#include <boost/cstdint.hpp>

class CaloGeometryHelper;

class HcalHitMaker : public CaloHitMaker
{
 public:
 
  typedef math::XYZVector XYZVector;
  typedef math::XYZVector XYZPoint;
  typedef ROOT::Math::Transform3DPJ Transform3D;

  HcalHitMaker(EcalHitMaker &, unsigned );
  ~HcalHitMaker() {;}
  
  inline void setSpotEnergy(double e) { spotEnergy=e;} 
  
  bool addHit(double r,double phi,unsigned layer=0);
  
  bool addHit(const XYZPoint & point ,unsigned layer=0);

  // get the hits
  const std::map<uint32_t,float>& getHits() { return hitMap_ ;} ;

  bool setDepth(double,bool inCm=false); 

 private:
    EcalHitMaker& myGrid;
    
    const FSimTrack * myTrack;
    XYZPoint ecalEntrance_;
    XYZVector particleDirection;
    int onHcal;
    
    double currentDepth_;
    Transform3D locToGlobal_;
    double radiusFactor_;
    bool mapCalculated_;
    
 public:
    static int getSubHcalDet(const FSimTrack* t)
      {
    //  std::cout << " getSubHcalDet " << std::endl;
    // According to  DataFormats/ HcalDetId/ interface/ HcalSubdetector.h
    //  std::cout << " onHcal " << t->onHcal() << " onVFcal " << t->onVFcal() << std::endl;
    if(t->onHcal()==1) return HcalBarrel;
    if(t->onHcal()==2) return HcalEndcap;
    if(t->onVFcal()==2) return HcalForward;
    return -1;
      } 

};
#endif