PFCTRecHitProducer.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFClusterProducer_PFCTRecHitProducer_h_
#define RecoParticleFlow_PFClusterProducer_PFCTRecHitProducer_h_

// system include files
#include <memory>
#include <vector>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/ParticleFlowReco/interface/PFRecHitFwd.h"

#include "CondFormats/HcalObjects/interface/HcalPFCorrs.h"
#include "CondFormats/HcalObjects/interface/HcalChannelQuality.h"
#include "CondFormats/EcalObjects/interface/EcalChannelStatus.h"
#include "Geometry/CaloTopology/interface/CaloTowerConstituentsMap.h"
#include "Geometry/CaloTopology/interface/HcalTopology.h"

#include "RecoParticleFlow/PFClusterProducer/interface/PFRecHitNavigatorBase.h"
#include "DataFormats/HcalRecHit/interface/HFRecHit.h"
#include "DataFormats/HcalRecHit/interface/HBHERecHit.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"

class CaloSubdetectorTopology;
class CaloSubdetectorGeometry;
class DetId;

class dso_hidden PFCTRecHitProducer final : public edm::stream::EDProducer<> {
public:
  explicit PFCTRecHitProducer(const edm::ParameterSet&);
  ~PFCTRecHitProducer() override;

  void beginLuminosityBlock(const edm::LuminosityBlock& lumi, const edm::EventSetup& es) override;

  void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override;

  reco::PFRecHit* createHcalRecHit(const DetId& detid,
                                   double energy,
                                   PFLayer::Layer layer,
                                   const CaloSubdetectorGeometry* geom,
                                   const CaloTowerDetId& newDetId);

protected:
  double thresh_Barrel_;
  double thresh_Endcap_;
  const HcalChannelQuality* theHcalChStatus;
  const EcalChannelStatus* theEcalChStatus;
  const CaloTowerConstituentsMap* theTowerConstituentsMap;
  const HcalTopology* theHcalTopology;

  // ----------access to event data
  edm::EDGetTokenT<HBHERecHitCollection> hcalToken_;
  edm::EDGetTokenT<HFRecHitCollection> hfToken_;
  edm::EDGetTokenT<CaloTowerCollection> towersToken_;

  double thresh_HF_;

  // Navigation in HF:  False = no real clustering in HF; True  = do clustering
  bool navigation_HF_;
  double weight_HFem_;
  double weight_HFhad_;

  // Apply HCAL DPG rechit calibration
  bool HCAL_Calib_;
  bool HF_Calib_;
  float HCAL_Calib_29;
  float HF_Calib_29;

  // Don't allow large energy in short fibres if there is no energy in long fibres
  double shortFibre_Cut;
  double longFibre_Fraction;

  // Don't allow large energy in long fibres if there is no energy in short fibres
  double longFibre_Cut;
  double shortFibre_Fraction;

  // Also apply HCAL DPG cleaning
  bool applyLongShortDPG_;

  // Don't allow too large timing excursion if energy in long/short fibres is large enough
  double longShortFibre_Cut;
  double minShortTiming_Cut;
  double maxShortTiming_Cut;
  double minLongTiming_Cut;
  double maxLongTiming_Cut;

  bool applyTimeDPG_;
  bool applyPulseDPG_;
  int HcalMaxAllowedHFLongShortSev_;
  int HcalMaxAllowedHFDigiTimeSev_;
  int HcalMaxAllowedHFInTimeWindowSev_;
  int HcalMaxAllowedChannelStatusSev_;

  int hcalHFLongShortFlagValue_;
  int hcalHFDigiTimeFlagValue_;
  int hcalHFInTimeWindowFlagValue_;

  // Compensate for dead ECAL channels
  bool ECAL_Compensate_;
  double ECAL_Threshold_;
  double ECAL_Compensation_;
  unsigned int ECAL_Dead_Code_;

  // Depth correction for EM and HAD rechits in the HF
  double EM_Depth_;
  double HAD_Depth_;

  int m_maxDepthHB;
  int m_maxDepthHE;

  std::unique_ptr<PFRecHitNavigatorBase> navigator_;
};

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(PFCTRecHitProducer);

#endif