HGCDigitizer.h

Go to the documentation of this file.

#ifndef HGCalSimProducers_HGCDigitizer_h
#define HGCalSimProducers_HGCDigitizer_h
 
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/ForwardDetId/interface/ForwardSubdetector.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
 
#include "SimDataFormats/CaloHit/interface/PCaloHit.h"
#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"
 
#include "SimCalorimetry/HGCalSimProducers/interface/HGCEEDigitizer.h"
#include "SimCalorimetry/HGCalSimProducers/interface/HGCHEfrontDigitizer.h"
#include "SimCalorimetry/HGCalSimProducers/interface/HGCHEbackDigitizer.h"
#include "SimCalorimetry/HGCalSimProducers/interface/HFNoseDigitizer.h"
#include "DataFormats/HGCDigi/interface/HGCDigiCollections.h"
#include "DataFormats/HGCDigi/interface/PHGCSimAccumulator.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "Geometry/HGCalGeometry/interface/HGCalGeometry.h"
#include "Geometry/HcalTowerAlgo/interface/HcalGeometry.h"
 
#include <vector>
#include <map>
#include <unordered_set>
#include <memory>
#include <tuple>
 
class PCaloHit;
class PileUpEventPrincipal;
 
class HGCDigitizer {
public:
  HGCDigitizer(const edm::ParameterSet& ps, edm::ConsumesCollector& iC);
  ~HGCDigitizer() {}
 
  // index , det id, time
  typedef std::tuple<int, uint32_t, float> HGCCaloHitTuple_t;
  static bool orderByDetIdThenTime(const HGCCaloHitTuple_t& a, const HGCCaloHitTuple_t& b) {
    unsigned int detId_a(std::get<1>(a)), detId_b(std::get<1>(b));
 
    if (detId_a < detId_b)
      return true;
    if (detId_a > detId_b)
      return false;
 
    double time_a(std::get<2>(a)), time_b(std::get<2>(b));
    if (time_a < time_b)
      return true;
 
    return false;
  }
 
  void accumulate(edm::Event const& e, edm::EventSetup const& c, CLHEP::HepRandomEngine* hre);
  void accumulate_forPreMix(edm::Event const& e, edm::EventSetup const& c, CLHEP::HepRandomEngine* hre);
 
  void accumulate(PileUpEventPrincipal const& e, edm::EventSetup const& c, CLHEP::HepRandomEngine* hre);
  void accumulate_forPreMix(PileUpEventPrincipal const& e, edm::EventSetup const& c, CLHEP::HepRandomEngine* hre);
 
  template <typename GEOM>
  void accumulate(edm::Handle<edm::PCaloHitContainer> const& hits,
                  int bxCrossing,
                  const GEOM* geom,
                  CLHEP::HepRandomEngine* hre);
  template <typename GEOM>
  void accumulate_forPreMix(edm::Handle<edm::PCaloHitContainer> const& hits,
                            int bxCrossing,
                            const GEOM* geom,
                            CLHEP::HepRandomEngine* hre);
 
  void accumulate_forPreMix(const PHGCSimAccumulator& simAccumulator, const bool minbiasFlag);
  void initializeEvent(edm::Event const& e, edm::EventSetup const& c);
  void finalizeEvent(edm::Event& e, edm::EventSetup const& c, CLHEP::HepRandomEngine* hre);
 
  bool producesEEDigis() { return ((mySubDet_ == ForwardSubdetector::HGCEE) || (myDet_ == DetId::HGCalEE)); }
  bool producesHEfrontDigis() { return ((mySubDet_ == ForwardSubdetector::HGCHEF) || (myDet_ == DetId::HGCalHSi)); }
  bool producesHEbackDigis() { return ((mySubDet_ == ForwardSubdetector::HGCHEB) || (myDet_ == DetId::HGCalHSc)); }
  bool producesHFNoseDigis() { return ((mySubDet_ == ForwardSubdetector::HFNose) && (myDet_ == DetId::Forward)); }
  std::string digiCollection() { return digiCollection_; }
  int geometryType() { return geometryType_; }
 
  void beginRun(const edm::EventSetup& es);
  void endRun();
 
private:
  uint32_t getType() const;
  bool getWeight(std::array<float, 3>& tdcForToAOnset, float& keV2fC) const;
  std::string hitCollection_, digiCollection_;
 
  //geometry type (0 pre-TDR; 1 TDR)
  int geometryType_;
 
  //digitization type (it's up to the specializations to decide it's meaning)
  int digitizationType_;
 
  // if true, we're running mixing in premixing stage1 and have to produce the output differently
  bool premixStage1_;
 
  // Minimum charge threshold for premixing stage1
  double premixStage1MinCharge_;
  // Maximum charge for packing in premixing stage1
  double premixStage1MaxCharge_;
 
  //handle sim hits
  int maxSimHitsAccTime_;
  double bxTime_, ev_per_eh_pair_;
  std::unique_ptr<hgc::HGCSimHitDataAccumulator> simHitAccumulator_;
  std::unique_ptr<hgc::HGCPUSimHitDataAccumulator> pusimHitAccumulator_;
  void resetSimHitDataAccumulator();
  void resetPUSimHitDataAccumulator();
  //debug position
  void checkPosition(const HGCalDigiCollection* digis) const;
 
  //digitizers
  std::unique_ptr<HGCEEDigitizer> theHGCEEDigitizer_;
  std::unique_ptr<HGCHEbackDigitizer> theHGCHEbackDigitizer_;
  std::unique_ptr<HGCHEfrontDigitizer> theHGCHEfrontDigitizer_;
  std::unique_ptr<HFNoseDigitizer> theHFNoseDigitizer_;
 
  //geometries
  std::unordered_set<DetId> validIds_;
  const HGCalGeometry* gHGCal_;
  const HcalGeometry* gHcal_;
 
  //detector and subdetector id
  DetId::Detector myDet_;
  ForwardSubdetector mySubDet_;
 
  //misc switches
  uint32_t verbosity_;
 
  //reference speed to evaluate time of arrival at the sensititive detector, assuming the center of CMS
  float refSpeed_;
 
  //delay to apply after evaluating time of arrival at the sensitive detector
  float tofDelay_;
 
  //average occupancies
  std::array<double, 4> averageOccupancies_;
  uint32_t nEvents_;
 
  //maxBx limit beyond which the Digitizer should filter out all hits
  static const unsigned int maxBx_ = 14;
  static const unsigned int thisBx_ = 9;
  std::vector<float> cce_;
  std::unordered_map<uint32_t, std::vector<std::pair<float, float> > > hitRefs_bx0;
  std::unordered_map<uint32_t, std::vector<std::tuple<float, float, float> > > PhitRefs_bx0;
  std::unordered_map<uint32_t, bool> hitOrder_monitor;
};
 
#endif