Inheritance diagram for FixedGridRhoProducerFastjetFromRecHit:

Public Member Functions
	FixedGridRhoProducerFastjetFromRecHit (const edm::ParameterSet &iConfig)

	~FixedGridRhoProducerFastjetFromRecHit () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
std::array< double, 4 >	getHitP4 (const DetId &detId, const double hitE, const CaloGeometry &caloGeometry) const

bool	passedEcalNoiseCut (const EcalRecHit &hit, const EcalPFRecHitThresholds &thresholds) const

bool	passedHcalNoiseCut (const HBHERecHit &hit, const HcalPFCuts *) const

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

Private Attributes
fastjet::GridMedianBackgroundEstimator	bge_

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

const bool	cutsFromDB_

const edm::EDGetTokenT< EcalRecHitCollection >	ebRecHitsTag_

const edm::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd >	ecalPFRecHitThresholdsToken_

const edm::EDGetTokenT< EcalRecHitCollection >	eeRecHitsTag_

const EgammaHcalIsolation::arrayHB	eThresHB_

const EgammaHcalIsolation::arrayHE	eThresHE_

const edm::EDGetTokenT< HBHERecHitCollection >	hbheRecHitsTag_

edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd >	hcalCutsToken_

HcalPFCuts const *	paramPF_ = nullptr

const bool	skipECAL_

const bool	skipHCAL_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 33 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Constructor & Destructor Documentation

◆ FixedGridRhoProducerFastjetFromRecHit()

FixedGridRhoProducerFastjetFromRecHit::FixedGridRhoProducerFastjetFromRecHit ( const edm::ParameterSet & iConfig )

explicit

Definition at line 67 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References deDxTools::esConsumes().

     : bge_(iConfig.getParameter<double>("maxRapidity"), iConfig.getParameter<double>("gridSpacing")),
       hbheRecHitsTag_(consumes(iConfig.getParameter<edm::InputTag>("hbheRecHitsTag"))),
       ebRecHitsTag_(consumes(iConfig.getParameter<edm::InputTag>("ebRecHitsTag"))),
       eeRecHitsTag_(consumes(iConfig.getParameter<edm::InputTag>("eeRecHitsTag"))),
       eThresHB_(iConfig.getParameter<EgammaHcalIsolation::arrayHB>("eThresHB")),
       eThresHE_(iConfig.getParameter<EgammaHcalIsolation::arrayHE>("eThresHE")),
       skipHCAL_(iConfig.getParameter<bool>("skipHCAL")),
       skipECAL_(iConfig.getParameter<bool>("skipECAL")),
       ecalPFRecHitThresholdsToken_{esConsumes()},
       caloGeometryToken_{esConsumes()},
       cutsFromDB_(iConfig.getParameter<bool>("usePFThresholdsFromDB")) {
   if (skipHCAL_ && skipECAL_) {
     throw cms::Exception("FixedGridRhoProducerFastjetFromRecHit")
         << "skipHCAL and skipECAL both can't be True. Please make at least one of them False.";
   }
   if (cutsFromDB_) {
     hcalCutsToken_ = esConsumes<HcalPFCuts, HcalPFCutsRcd>(edm::ESInputTag("", "withTopo"));
   }
   produces<double>();
 }

◆ ~FixedGridRhoProducerFastjetFromRecHit()

FixedGridRhoProducerFastjetFromRecHit::~FixedGridRhoProducerFastjetFromRecHit ( )

overridedefault

Member Function Documentation

◆ fillDescriptions()

void FixedGridRhoProducerFastjetFromRecHit::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 89 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                                        {
   edm::ParameterSetDescription desc;
   //We use raw recHits and not the PF recHits, because in Phase1 muon and egamma HLT,
   //PF recHit collection is regional, while the full raw recHit collection is available.
   desc.add<edm::InputTag>("hbheRecHitsTag", edm::InputTag("hltHbhereco"));
   desc.add<edm::InputTag>("ebRecHitsTag", edm::InputTag("hltEcalRecHit", "EcalRecHitsEB"));
   desc.add<edm::InputTag>("eeRecHitsTag", edm::InputTag("hltEcalRecHit", "EcalRecHitsEE"));
   desc.add<bool>("skipHCAL", false);
   desc.add<bool>("skipECAL", false);
   //eThresHB/HE are from RecoParticleFlow/PFClusterProducer/python/particleFlowRecHitHBHE_cfi.py
   desc.add<std::vector<double> >("eThresHB", {0.1, 0.2, 0.3, 0.3});
   desc.add<std::vector<double> >("eThresHE", {0.1, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2});
   desc.add<double>("maxRapidity", 2.5);
   desc.add<double>("gridSpacing", 0.55);
   desc.add<bool>("usePFThresholdsFromDB", true);
   descriptions.addWithDefaultLabel(desc);
 }

◆ getHitP4()

std::array< double, 4 > FixedGridRhoProducerFastjetFromRecHit::getHitP4	(	const DetId &	detId,
		const double	hitE,
		const CaloGeometry &	caloGeometry
	)		const

private

Definition at line 153 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References funct::cos(), hcalRecHitTable_cff::detId, CaloGeometry::getSubdetectorGeometry(), funct::sin(), and GeomDetEnumerators::subDetGeom.

Referenced by produce().

                                                                                                               {
   const CaloSubdetectorGeometry *subDetGeom = caloGeometry.getSubdetectorGeometry(detId);
   const auto &gpPos = subDetGeom->getGeometry(detId)->repPos();
   const double thispt = hitE / cosh(gpPos.eta());
   const double thispx = thispt * cos(gpPos.phi());
   const double thispy = thispt * sin(gpPos.phi());
   const double thispz = thispt * sinh(gpPos.eta());
   return std::array<double, 4>{{thispx, thispy, thispz, hitE}};
 }

◆ passedEcalNoiseCut()

bool FixedGridRhoProducerFastjetFromRecHit::passedEcalNoiseCut	(	const EcalRecHit &	hit,
		const EcalPFRecHitThresholds &	thresholds
	)		const

private

Definition at line 180 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References particleFlowZeroSuppressionECAL_cff::thresholds.

Referenced by produce().

                                                                                                                {
   return (hit.energy() > thresholds[hit.detid()]);
 }

◆ passedHcalNoiseCut()

bool FixedGridRhoProducerFastjetFromRecHit::passedHcalNoiseCut	(	const HBHERecHit &	hit,
		const HcalPFCuts *	hcalcuts
	)		const

private

Definition at line 166 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References eThresHB_, eThresHE_, HcalCondObjectContainer< Item >::getValues(), HcalBarrel, HcalEndcap, hit::id, and B2GTnPMonitor_cfi::item.

Referenced by produce().

                                                                                                  {
   if (hcalcuts != nullptr) {  // using hcal cuts from DB
     const HcalPFCut *item = hcalcuts->getValues(hit.id().rawId());
     return (hit.energy() > item->noiseThreshold());
   } else {  // using hcal cuts from config file
     const auto thisDetId = hit.id();
     const auto thisDepth = thisDetId.depth();
     return (thisDetId.subdet() == HcalBarrel && hit.energy() > eThresHB_[thisDepth - 1]) ||
            (thisDetId.subdet() == HcalEndcap && hit.energy() > eThresHE_[thisDepth - 1]);
   }
 }

◆ produce()

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

overrideprivate

Definition at line 109 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References bge_, caloGeometryToken_, cutsFromDB_, ebRecHitsTag_, ecalPFRecHitThresholdsToken_, eeRecHitsTag_, edm::EventSetup::getData(), getHitP4(), ecalDrivenElectronSeeds_cfi::hbheRecHits, hbheRecHitsTag_, hcalCutsToken_, hit::id, iEvent, PixelMapPlotter::inputs, paramPF_, passedEcalNoiseCut(), passedHcalNoiseCut(), skipECAL_, skipHCAL_, and particleFlowZeroSuppressionECAL_cff::thresholds.

                                                                                                  {
   if (cutsFromDB_) {
     paramPF_ = &iSetup.getData(hcalCutsToken_);
   }
 
   std::vector<fastjet::PseudoJet> inputs;
   auto const &thresholds = iSetup.getData(ecalPFRecHitThresholdsToken_);
   auto const &caloGeometry = iSetup.getData(caloGeometryToken_);
 
   if (!skipHCAL_) {
     auto const &hbheRecHits = iEvent.get(hbheRecHitsTag_);
     inputs.reserve(inputs.size() + hbheRecHits.size());
     for (const auto &hit : hbheRecHits) {
       if (passedHcalNoiseCut(hit, paramPF_)) {
         const auto &hitp4 = getHitP4(hit.id(), hit.energy(), caloGeometry);
         inputs.emplace_back(fastjet::PseudoJet(hitp4[0], hitp4[1], hitp4[2], hitp4[3]));
       }
     }
   }
 
   if (!skipECAL_) {
     auto const &ebRecHits = iEvent.get(ebRecHitsTag_);
     inputs.reserve(inputs.size() + ebRecHits.size());
     for (const auto &hit : ebRecHits) {
       if (passedEcalNoiseCut(hit, thresholds)) {
         const auto &hitp4 = getHitP4(hit.id(), hit.energy(), caloGeometry);
         inputs.emplace_back(fastjet::PseudoJet(hitp4[0], hitp4[1], hitp4[2], hitp4[3]));
       }
     }
 
     auto const &eeRecHits = iEvent.get(eeRecHitsTag_);
     inputs.reserve(inputs.size() + eeRecHits.size());
     for (const auto &hit : eeRecHits) {
       if (passedEcalNoiseCut(hit, thresholds)) {
         const auto &hitp4 = getHitP4(hit.id(), hit.energy(), caloGeometry);
         inputs.emplace_back(fastjet::PseudoJet(hitp4[0], hitp4[1], hitp4[2], hitp4[3]));
       }
     }
   }
 
   bge_.set_particles(inputs);
   iEvent.put(std::make_unique<double>(bge_.rho()));
 }