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

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

Private Attributes
fastjet::GridMedianBackgroundEstimator	bge_

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

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_

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 31 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Constructor & Destructor Documentation

◆ FixedGridRhoProducerFastjetFromRecHit()

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

explicit

Definition at line 60 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()} {
   if (skipHCAL_ && skipECAL_) {
     throw cms::Exception("FixedGridRhoProducerFastjetFromRecHit")
         << "skipHCAL and skipECAL both can't be True. Please make at least one of them False.";
   }
   produces<double>();
 }

◆ ~FixedGridRhoProducerFastjetFromRecHit()

FixedGridRhoProducerFastjetFromRecHit::~FixedGridRhoProducerFastjetFromRecHit ( )

overridedefault

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 78 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);
   descriptions.addWithDefaultLabel(desc);
 }

◆ getHitP4()

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

private

Definition at line 137 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 158 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

private

Definition at line 150 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References eThresHB_, eThresHE_, HcalBarrel, HcalEndcap, and hit::id.

Referenced by produce().

                                                                                           {
   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 97 of file FixedGridRhoProducerFastjetFromRecHit.cc.

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

                                                                                                  {
   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)) {
         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()));
 }