CaloHitResponse.h

Go to the documentation of this file.

#ifndef CaloSimAlgos_CaloHitResponse_h
#define CaloSimAlgos_CaloHitResponse_h

#include "CalibFormats/CaloObjects/interface/CaloSamples.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "SimDataFormats/CrossingFrame/interface/MixCollection.h"
#include "SimDataFormats/CaloHit/interface/PCaloHit.h"
#include "CLHEP/Random/RandPoissonQ.h"
#include "SimCalorimetry/CaloSimAlgos/interface/CaloVPECorrection.h"

#include<map>
#include<vector>

#define ChangeHcalEnergyScale

class CaloVShape;
class CaloShapes;
class CaloVSimParameterMap;
class CaloVHitCorrection;
class CaloVHitFilter;
class CaloSimParameters;

class CaloHitResponse 
{
public:
  typedef std::map<DetId, CaloSamples> AnalogSignalMap;
  // get this from somewhere external
  enum {BUNCHSPACE=25};

  CaloHitResponse(const CaloVSimParameterMap * parameterMap, const CaloVShape * shape);
  CaloHitResponse(const CaloVSimParameterMap * parameterMap, const CaloShapes * shapes);

  virtual ~CaloHitResponse();

  // change HBHE scale
  void initHBHEScale();
  void setHBHEScale(std::string &); //GMA

  void setBunchRange(int minBunch, int maxBunch);

  void setGeometry(const CaloGeometry * geometry) { theGeometry = geometry; }

  virtual bool keepBlank() const { return true ; }

  virtual void initializeHits() {}

  virtual void finalizeHits() {}

  virtual void run(MixCollection<PCaloHit> & hits);

  virtual void add(const PCaloHit & hit);

  void add(const CaloSamples & signal);

  void setHitFilter(const CaloVHitFilter * filter) {
    theHitFilter = filter;
  }

  void setHitCorrection(const CaloVHitCorrection * hitCorrection) {
    theHitCorrection = hitCorrection;
  }

  void setPECorrection(const CaloVPECorrection * peCorrection) {
    thePECorrection = peCorrection;
  }

  virtual void setRandomEngine(CLHEP::HepRandomEngine & engine);

  void clear() {theAnalogSignalMap.clear();}
 
  void addHit(const PCaloHit * hit, CaloSamples & frame) const;

  virtual CaloSamples makeAnalogSignal(const PCaloHit & inputHit) const;

  double analogSignalAmplitude(const DetId & id, float energy, const CaloSimParameters & parameters) const;

  CaloSamples * findSignal(const DetId & detId);

  int nSignals() const {return theAnalogSignalMap.size();}

  CaloSamples makeBlankSignal(const DetId & detId) const;


  double timeOfFlight(const DetId & detId) const;

  void setPhaseShift(const double & thePhaseShift) { thePhaseShift_ = thePhaseShift; }


  bool withinBunchRange(int bunchCrossing) const {
    return(bunchCrossing >= theMinBunch && bunchCrossing <= theMaxBunch);
  }

protected:

  AnalogSignalMap theAnalogSignalMap;

  const CaloVSimParameterMap * theParameterMap;
  const CaloShapes * theShapes;
  const CaloVShape * theShape;
  const CaloVHitCorrection * theHitCorrection;
  const CaloVPECorrection * thePECorrection;
  const CaloVHitFilter * theHitFilter;

  const CaloGeometry * theGeometry;

  mutable CLHEP::RandPoissonQ * theRandPoisson;

  int theMinBunch;
  int theMaxBunch;

  double thePhaseShift_;

  // private : 
  bool  changeScale;
#ifdef ChangeHcalEnergyScale
  float hcal_en_scale[100][72][4];
#endif  
};

#endif