HcalRecHitsMaker.h

Go to the documentation of this file.

#ifndef FastSimulation__HcalRecHitsMaker__h
#define FastSimulation__HcalRecHitsMaker__h

#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
#include "FastSimulation/Utilities/interface/GaussianTail.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include <map>
#include <vector>

//#include <boost/cstdint.hpp>

class CaloGeometry;
class RandomEngine;
// class HcalTPGCoder;
class HcalSimParameterMap;
class HcalDbService;
class HcalRespCorrs;

namespace edm { 
  class ParameterSet;
  class Event;
  class EventSetup;
}

class HcalRecHitsMaker
{
 public:
  HcalRecHitsMaker(edm::ParameterSet const & p,int,const RandomEngine* random);
  ~HcalRecHitsMaker();

  void loadHcalRecHits(edm::Event &iEvent, HBHERecHitCollection& hbheHits, HBHEDigiCollection& hbheDigis);
  void loadHcalRecHits(edm::Event &iEvent, HORecHitCollection &ho, HODigiCollection & hoDigis);
  void loadHcalRecHits(edm::Event &iEvent, HFRecHitCollection &hfHits, HFDigiCollection& hfDigis);
  void init(const edm::EventSetup &es,bool dodigis,bool domiscalib);

 private:
  unsigned createVectorsOfCells(const edm::EventSetup &es);
  unsigned createVectorOfSubdetectorCells( const CaloGeometry&,int subdetn,std::vector<int>&);
  unsigned noisifySubdet(std::vector<float >& theMap, std::vector<int>& theHits,const std::vector<int>& thecells, unsigned ncells, double  hcalHotFraction_, const GaussianTail *,double sigma,double threshold,double correctionfactor); 
  // Not currently used. Will probably be removed soon.
  //  void noisifySignal(std::map<uint32_t,std::pair<float,bool> >& theMap); 
  void noisify();
  double noiseInfCfromDB(const HcalDbService * conditions,const HcalDetId & detId);
  void Fill(int id,float energy, std::vector<int> & myHits,float noise,float correctionfactor);
  void loadPCaloHits(const edm::Event & iEvent);
  
  void clean();
  void cleanSubDet(std::vector<float>& hits,std::vector<int>& cells);
  // conversion for digitization
  int fCtoAdc(double fc) const;
  double fractionOOT(int time_slice);

 private:
  unsigned det_;
  std::vector<double> threshold_;
  std::vector<double> noise_;
  std::vector<double> corrfac_;
  std::vector<double> hcalHotFraction_;
  unsigned nnoise_;

  //  edm::ESHandle<CaloTowerConstituentsMap> calotowerMap_;
  edm::InputTag inputCol_;
  static bool initialized_;
  bool doDigis_;
  bool doMiscalib_;
  bool doSaturation_;
  bool noiseFromDb_;
  double refactor_;
  double refactor_mean_;
  std::string hcalfileinpath_;

  std::vector<float> hcalRecHits_;

  std::vector<int> firedCells_;

  static std::vector<HcalDetId> theDetIds_;
  static std::vector<float> miscalib_;

  // coefficients for fC to ADC conversion
  static std::vector<int> fctoadc_;

  static std::vector<float> peds_;
  static std::vector<float> gains_;
  static std::vector<float> sat_;
  static std::vector<float> noisesigma_;
  static std::vector<float> TPGFactor_;
 
  // the hashed indices
  static unsigned maxIndex_;
  static std::vector<int> hbhi_;
  static std::vector<int> hehi_;
  static std::vector<int> hohi_;
  static std::vector<int> hfhi_;
  unsigned nhbcells_,nhecells_,nhocells_,nhfcells_;

  const RandomEngine* random_;
  std::vector<GaussianTail*> myGaussianTailGenerators_;

  //  const HcalTPGCoder * myCoder_;
  //  HcalSimParameterMap * myHcalSimParameterMap_;

  // the access to the response corection factors
  const HcalRespCorrs* myRespCorr;
};

#endif