Inheritance diagram for HLTRechitInRegionsProducer< T1 >:

Public Member Functions
	HLTRechitInRegionsProducer (const edm::ParameterSet &ps)

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

	~HLTRechitInRegionsProducer () override=default

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 Types
typedef std::vector< T1 >	T1Collection

typedef T1::const_iterator	T1iterator

Private Member Functions
void	eventSetupConsumes ()

template<>
void	eventSetupConsumes ()

void	getEtaPhiRegions (std::vector< RectangularEtaPhiRegion > *, T1Collection const &, edm::EventSetup const &, bool const)

template<>
void	getEtaPhiRegions (std::vector< RectangularEtaPhiRegion > *theRegions, T1Collection const &theCandidateCollection, edm::EventSetup const &eventSetup, bool const isolatedCase)

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

const bool	doIsolated_

const std::vector< edm::InputTag >	hitLabels

std::vector< edm::EDGetTokenT< EcalRecHitCollection > >	hitTokens

edm::ESGetToken< L1CaloGeometry, L1CaloGeometryRecord >	l1CaloGeometryToken_

const double	l1LowerThr_

const double	l1LowerThrIgnoreIsolation_

const edm::EDGetTokenT< T1Collection >	l1TokenIsolated_

const edm::EDGetTokenT< T1Collection >	l1TokenNonIsolated_

const double	l1UpperThr_

const std::vector< std::string >	productLabels

const double	regionEtaMargin_

const double	regionPhiMargin_

std::vector< edm::EDGetTokenT< EcalUncalibratedRecHitCollection > >	uncalibHitTokens

const bool	useUncalib_

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

template<typename T1>
class HLTRechitInRegionsProducer< T1 >

Definition at line 70 of file HLTRechitInRegionsProducer.cc.

Member Typedef Documentation

◆ T1Collection

template<typename T1 >

typedef std::vector<T1> HLTRechitInRegionsProducer< T1 >::T1Collection

private

Definition at line 71 of file HLTRechitInRegionsProducer.cc.

◆ T1iterator

template<typename T1 >

typedef T1::const_iterator HLTRechitInRegionsProducer< T1 >::T1iterator

private

Definition at line 72 of file HLTRechitInRegionsProducer.cc.

Constructor & Destructor Documentation

◆ HLTRechitInRegionsProducer()

template<typename T1 >

HLTRechitInRegionsProducer< T1 >::HLTRechitInRegionsProducer ( const edm::ParameterSet & ps )

Definition at line 110 of file HLTRechitInRegionsProducer.cc.

References HLTRechitInRegionsProducer< T1 >::eventSetupConsumes(), HLTRechitInRegionsProducer< T1 >::hitLabels, HLTRechitInRegionsProducer< T1 >::hitTokens, mps_fire::i, HLTRechitInRegionsProducer< T1 >::productLabels, HLTRechitInRegionsProducer< T1 >::uncalibHitTokens, and HLTRechitInRegionsProducer< T1 >::useUncalib_.

     : useUncalib_(ps.getParameter<bool>("useUncalib")),
       doIsolated_(ps.getParameter<bool>("doIsolated")),
       l1TokenIsolated_(doIsolated_ ? consumes<T1Collection>(ps.getParameter<edm::InputTag>("l1TagIsolated"))
                                    : edm::EDGetTokenT<T1Collection>()),
       l1TokenNonIsolated_(consumes<T1Collection>(ps.getParameter<edm::InputTag>("l1TagNonIsolated"))),
       l1LowerThr_(ps.getParameter<double>("l1LowerThr")),
       l1UpperThr_(ps.getParameter<double>("l1UpperThr")),
       l1LowerThrIgnoreIsolation_(ps.getParameter<double>("l1LowerThrIgnoreIsolation")),
       regionEtaMargin_(ps.getParameter<double>("regionEtaMargin")),
       regionPhiMargin_(ps.getParameter<double>("regionPhiMargin")),
       hitLabels(ps.getParameter<std::vector<edm::InputTag>>("ecalhitLabels")),
       productLabels(ps.getParameter<std::vector<std::string>>("productLabels")) {
   eventSetupConsumes();
 
   if (useUncalib_) {
     for (unsigned int i = 0; i < hitLabels.size(); i++) {
       uncalibHitTokens.push_back(consumes<EcalUncalibratedRecHitCollection>(hitLabels[i]));
       produces<EcalUncalibratedRecHitCollection>(productLabels[i]);
     }
   } else {
     for (unsigned int i = 0; i < hitLabels.size(); i++) {
       hitTokens.push_back(consumes<EcalRecHitCollection>(hitLabels[i]));
       produces<EcalRecHitCollection>(productLabels[i]);
     }
   }
 }

◆ ~HLTRechitInRegionsProducer()

template<typename T1 >

HLTRechitInRegionsProducer< T1 >::~HLTRechitInRegionsProducer ( )

overridedefault

Member Function Documentation

◆ eventSetupConsumes() [1/2]

template<typename T1 >

void HLTRechitInRegionsProducer< T1 >::eventSetupConsumes ( )

private

Definition at line 145 of file HLTRechitInRegionsProducer.cc.

References deDxTools::esConsumes().

Referenced by HLTRechitInRegionsProducer< T1 >::HLTRechitInRegionsProducer().

                                                         {
   caloGeometryToken_ = esConsumes();
 }

◆ eventSetupConsumes() [2/2]

template<>

void HLTRechitInRegionsProducer< l1extra::L1EmParticle >::eventSetupConsumes ( )

private

Definition at line 139 of file HLTRechitInRegionsProducer.cc.

References deDxTools::esConsumes().

                                                                          {
   caloGeometryToken_ = esConsumes();
   l1CaloGeometryToken_ = esConsumes();
 }

◆ fillDescriptions()

template<typename T1 >

void HLTRechitInRegionsProducer< T1 >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 150 of file HLTRechitInRegionsProducer.cc.

References edm::ConfigurationDescriptions::add(), defaultModuleLabel(), submitPVResolutionJobs::desc, HLT_2022v12_cff::InputTag, and HLTMuonOfflineAnalyzer_cfi::inputTags.

                                                                                                 {
   edm::ParameterSetDescription desc;
   std::vector<std::string> productTags;
   productTags.push_back("EcalRegionalRecHitsEB");
   productTags.push_back("EcalRegionalRecHitsEE");
   desc.add<std::vector<std::string>>("productLabels", productTags);
   std::vector<edm::InputTag> inputTags;
   inputTags.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEB"));
   inputTags.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEE"));
   inputTags.push_back(edm::InputTag("hltESRegionalEgammaRecHit:EcalRecHitsES"));
   desc.add<std::vector<edm::InputTag>>("ecalhitLabels", inputTags);
   desc.add<edm::InputTag>("l1TagIsolated", edm::InputTag("l1extraParticles", "Isolated"));
   desc.add<edm::InputTag>("l1TagNonIsolated", edm::InputTag("l1extraParticles", "NonIsolated"));
   desc.add<bool>("useUncalib", true);
   desc.add<bool>("doIsolated", true);
   desc.add<double>("l1LowerThr", 5.0);
   desc.add<double>("l1UpperThr", 999.);
   desc.add<double>("l1LowerThrIgnoreIsolation", 0.0);
   desc.add<double>("regionEtaMargin", 0.14);
   desc.add<double>("regionPhiMargin", 0.4);
   descriptions.add(defaultModuleLabel<HLTRechitInRegionsProducer<T1>>(), desc);
 }

◆ getEtaPhiRegions() [1/2]

template<typename T1 >

void HLTRechitInRegionsProducer< T1 >::getEtaPhiRegions	(	std::vector< RectangularEtaPhiRegion > *	theRegions,
		T1Collection const &	theCandidateCollection,
		edm::EventSetup const &	,
		bool const
	)

private

Definition at line 322 of file HLTRechitInRegionsProducer.cc.

References l1tSlwPFPuppiJets_cfi::phiLow.

                                                                   {
   for (unsigned int candItr = 0; candItr < theCandidateCollection.size(); candItr++) {
     T1 emItr = theCandidateCollection.at(candItr);
     if ((emItr.et() > l1LowerThr_) and (emItr.et() < l1UpperThr_)) {
       double etaLow = emItr.eta() - regionEtaMargin_;
       double etaHigh = emItr.eta() + regionEtaMargin_;
       double phiLow = emItr.phi() - regionPhiMargin_;
       double phiHigh = emItr.phi() + regionPhiMargin_;
 
       theRegions->push_back(RectangularEtaPhiRegion(etaLow, etaHigh, phiLow, phiHigh));
     }
   }
 }

◆ getEtaPhiRegions() [2/2]

template<>

void HLTRechitInRegionsProducer< l1extra::L1EmParticle >::getEtaPhiRegions	(	std::vector< RectangularEtaPhiRegion > *	theRegions,
		T1Collection const &	theCandidateCollection,
		edm::EventSetup const &	eventSetup,
		bool const	isolatedCase
	)

private

Definition at line 285 of file HLTRechitInRegionsProducer.cc.

References reco::LeafCandidate::et(), L1GctEmCand::etaIndex(), options_cfi::eventSetup, l1extra::L1EmParticle::gctEmCand(), L1GctEmCand::phiIndex(), and l1tSlwPFPuppiJets_cfi::phiLow.

                              {
   auto const& l1CaloGeom = eventSetup.getData(l1CaloGeometryToken_);
 
   for (unsigned int candItr = 0; candItr < theCandidateCollection.size(); candItr++) {
     l1extra::L1EmParticle emItr = theCandidateCollection.at(candItr);
 
     if (!isolatedCase) {
       if (doIsolated_ and (emItr.et() < l1LowerThrIgnoreIsolation_))
         continue;
     }
 
     if ((emItr.et() > l1LowerThr_) and (emItr.et() < l1UpperThr_)) {
       // Access the GCT hardware object corresponding to the L1Extra EM object.
       int etaIndex = emItr.gctEmCand()->etaIndex();
       int phiIndex = emItr.gctEmCand()->phiIndex();
 
       // Use the L1CaloGeometry to find the eta, phi bin boundaries.
       double etaLow = l1CaloGeom.etaBinLowEdge(etaIndex);
       double etaHigh = l1CaloGeom.etaBinHighEdge(etaIndex);
       double phiLow = l1CaloGeom.emJetPhiBinLowEdge(phiIndex);
       double phiHigh = l1CaloGeom.emJetPhiBinHighEdge(phiIndex);
 
       etaLow -= regionEtaMargin_;
       etaHigh += regionEtaMargin_;
       phiLow -= regionPhiMargin_;
       phiHigh += regionPhiMargin_;
 
       theRegions->push_back(RectangularEtaPhiRegion(etaLow, etaHigh, phiLow, phiHigh));
     }
   }
 }

◆ produce()

template<typename T1 >

void HLTRechitInRegionsProducer< T1 >::produce	(	edm::Event &	evt,
		edm::EventSetup const &	eventSetup
	)

override

Definition at line 174 of file HLTRechitInRegionsProducer.cc.

References edm::SortedCollection< T, SORT >::begin(), DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), options_cfi::eventSetup, relativeConstraints::geometry, edm::Event::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getGeometry(), hfClusterShapes_cfi::hits, mps_fire::i, sistrip::SpyUtilities::isValid(), eostools::move(), or, edm::Handle< T >::product(), HLT_2022v12_cff::productLabels, edm::Event::put(), FastTrackerRecHitMaskProducer_cfi::recHits, and HLT_2022v12_cff::region.

                                                                                            {
   // get the collection geometry:
   auto const& geometry = eventSetup.getData(caloGeometryToken_);
   const CaloSubdetectorGeometry* geometry_p;
   std::unique_ptr<const CaloSubdetectorTopology> topology;
 
   //Get the L1 EM Particle Collection
   edm::Handle<T1Collection> emIsolColl;
   if (doIsolated_) {
     evt.getByToken(l1TokenIsolated_, emIsolColl);
   }
 
   std::vector<RectangularEtaPhiRegion> regions;
   if (doIsolated_)
     getEtaPhiRegions(&regions, *emIsolColl, eventSetup, true);
 
   if (!doIsolated_ or (l1LowerThrIgnoreIsolation_ < 64))
     getEtaPhiRegions(&regions, evt.get(l1TokenNonIsolated_), eventSetup, false);
 
   if (useUncalib_) {
     edm::Handle<EcalUncalibratedRecHitCollection> urhcH[3];
     for (unsigned int i = 0; i < hitLabels.size(); i++) {
       auto uhits = std::make_unique<EcalUncalibratedRecHitCollection>();
 
       evt.getByToken(uncalibHitTokens[i], urhcH[i]);
       if (!(urhcH[i].isValid())) {
         edm::LogError("ProductNotFound") << "could not get a handle on the EcalRecHitCollection! ("
                                          << hitLabels[i].encode() << ")" << std::endl;
         return;
       }
       const EcalUncalibratedRecHitCollection* uncalibRecHits = urhcH[i].product();
 
       if (!uncalibRecHits->empty()) {
         if ((*uncalibRecHits)[0].id().subdetId() == EcalBarrel) {
           geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
           topology = std::make_unique<EcalBarrelTopology>(geometry);
         } else if ((*uncalibRecHits)[0].id().subdetId() == EcalEndcap) {
           geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
           topology = std::make_unique<EcalEndcapTopology>(geometry);
         } else if ((*uncalibRecHits)[0].id().subdetId() == EcalPreshower) {
           geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
           topology = std::make_unique<EcalPreshowerTopology>();
         } else
           throw(std::runtime_error("\n\nProducer encountered invalied ecalhitcollection type.\n\n"));
 
         if (!regions.empty()) {
           EcalUncalibratedRecHitCollection::const_iterator it;
 
           for (it = uncalibRecHits->begin(); it != uncalibRecHits->end(); it++) {
             auto this_cell = geometry_p->getGeometry(it->id());
 
             std::vector<RectangularEtaPhiRegion>::const_iterator region;
             for (region = regions.begin(); region != regions.end(); region++) {
               if (region->inRegion(this_cell->etaPos(), this_cell->phiPos())) {
                 uhits->push_back(*it);
                 break;
               }
             }
           }
         }
       }
       evt.put(std::move(uhits), productLabels[i]);
     }
 
   } else {
     edm::Handle<EcalRecHitCollection> rhcH[3];
     for (unsigned int i = 0; i < hitLabels.size(); i++) {
       auto hits = std::make_unique<EcalRecHitCollection>();
 
       evt.getByToken(hitTokens[i], rhcH[i]);
       if (!(rhcH[i].isValid())) {
         edm::LogError("ProductNotFound") << "could not get a handle on the EcalRecHitCollection! ("
                                          << hitLabels[i].encode() << ")" << std::endl;
         return;
       }
       const EcalRecHitCollection* recHits = rhcH[i].product();
 
       if (!recHits->empty()) {
         if ((*recHits)[0].id().subdetId() == EcalBarrel) {
           geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
           topology = std::make_unique<EcalBarrelTopology>(geometry);
         } else if ((*recHits)[0].id().subdetId() == EcalEndcap) {
           geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
           topology = std::make_unique<EcalEndcapTopology>(geometry);
         } else if ((*recHits)[0].id().subdetId() == EcalPreshower) {
           geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
           topology = std::make_unique<EcalPreshowerTopology>();
         } else
           throw(std::runtime_error("\n\nProducer encountered invalied ecalhitcollection type.\n\n"));
 
         if (!regions.empty()) {
           EcalRecHitCollection::const_iterator it;
           for (it = recHits->begin(); it != recHits->end(); it++) {
             auto this_cell = geometry_p->getGeometry(it->id());
 
             std::vector<RectangularEtaPhiRegion>::const_iterator region;
             for (region = regions.begin(); region != regions.end(); region++) {
               if (region->inRegion(this_cell->etaPos(), this_cell->phiPos())) {
                 hits->push_back(*it);
                 break;
               }
             }
           }
         }
       }
       evt.put(std::move(hits), productLabels[i]);
     }
   }
 }