Inheritance diagram for HLTRecHitInAllL1RegionsProducer< RecHitType >:

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

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

	~HLTRecHitInAllL1RegionsProducer () 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 Types
using	RecHitCollectionType = edm::SortedCollection< RecHitType >

Private Member Functions
L1RegionDataBase *	createL1RegionData (const std::string &, const edm::ParameterSet &, edm::ConsumesCollector &&)

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

std::vector< std::unique_ptr< L1RegionDataBase > >	l1RegionData_

std::vector< std::string >	productLabels_

std::vector< edm::InputTag >	recHitLabels_

std::vector< edm::EDGetTokenT< RecHitCollectionType > >	recHitTokens_

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 RecHitType>
class HLTRecHitInAllL1RegionsProducer< RecHitType >

Definition at line 110 of file HLTRecHitInAllL1RegionsProducer.cc.

Member Typedef Documentation

◆ RecHitCollectionType

template<typename RecHitType >

using HLTRecHitInAllL1RegionsProducer< RecHitType >::RecHitCollectionType = edm::SortedCollection<RecHitType>

private

Definition at line 111 of file HLTRecHitInAllL1RegionsProducer.cc.

Constructor & Destructor Documentation

◆ HLTRecHitInAllL1RegionsProducer()

template<typename RecHitType >

HLTRecHitInAllL1RegionsProducer< RecHitType >::HLTRecHitInAllL1RegionsProducer ( const edm::ParameterSet & ps )

Definition at line 136 of file HLTRecHitInAllL1RegionsProducer.cc.

References deDxTools::esConsumes().

     : caloGeometryToken_{esConsumes()} {
   const std::vector<edm::ParameterSet> l1InputRegions =
       para.getParameter<std::vector<edm::ParameterSet>>("l1InputRegions");
   for (auto& pset : l1InputRegions) {
     const std::string type = pset.getParameter<std::string>("type");
     // meh I was going to use a factory but it was going to be overly complex for my needs
     l1RegionData_.emplace_back(createL1RegionData(type, pset, consumesCollector()));
   }
   recHitLabels_ = para.getParameter<std::vector<edm::InputTag>>("recHitLabels");
   productLabels_ = para.getParameter<std::vector<std::string>>("productLabels");
 
   for (unsigned int collNr = 0; collNr < recHitLabels_.size(); collNr++) {
     recHitTokens_.push_back(consumes<RecHitCollectionType>(recHitLabels_[collNr]));
     produces<RecHitCollectionType>(productLabels_[collNr]);
   }
 }

◆ ~HLTRecHitInAllL1RegionsProducer()

template<typename RecHitType >

HLTRecHitInAllL1RegionsProducer< RecHitType >::~HLTRecHitInAllL1RegionsProducer ( )

inlineoverride

Definition at line 115 of file HLTRecHitInAllL1RegionsProducer.cc.

115 {}

Member Function Documentation

◆ createL1RegionData()

template<typename RecHitType >

L1RegionDataBase * HLTRecHitInAllL1RegionsProducer< RecHitType >::createL1RegionData	(	const std::string &	type,
		const edm::ParameterSet &	para,
		edm::ConsumesCollector &&	consumesColl
	)

private

Definition at line 262 of file HLTRecHitInAllL1RegionsProducer.cc.

References Exception.

                                                                                              {
   if (type == "L1EmParticle") {
     return new L1RegionData<l1extra::L1EmParticleCollection>(para, consumesColl);
   } else if (type == "L1JetParticle") {
     return new L1RegionData<l1extra::L1JetParticleCollection>(para, consumesColl);
   } else if (type == "L1MuonParticle") {
     return new L1RegionData<l1extra::L1MuonParticleCollection>(para, consumesColl);
   } else if (type == "EGamma") {
     return new L1RegionData<l1t::EGammaBxCollection>(para, consumesColl);
   } else if (type == "Jet") {
     return new L1RegionData<l1t::JetBxCollection>(para, consumesColl);
   } else if (type == "Muon") {
     return new L1RegionData<l1t::MuonBxCollection>(para, consumesColl);
   } else if (type == "Tau") {
     return new L1RegionData<l1t::TauBxCollection>(para, consumesColl);
   } else {
     //this is a major issue and could lead to rather subtle efficiency losses, so if its incorrectly configured, we're aborting the job!
     throw cms::Exception("InvalidConfig")
         << " type " << type
         << " is not recognised, this means the rec-hit you think you are keeping may not be and you should fix this "
            "error as it can lead to hard to find efficiency loses"
         << std::endl;
   }
 }

◆ fillDescriptions()

template<typename RecHitType >

void HLTRecHitInAllL1RegionsProducer< RecHitType >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 155 of file HLTRecHitInAllL1RegionsProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSet::addParameter(), defaultModuleLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, JetMonitor_cfi::jetPSet, HLT_2023v12_cff::l1InputRegions, HLT_2023v12_cff::recHitLabels, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                              {
   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> recHitLabels;
   recHitLabels.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEB"));
   recHitLabels.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEE"));
   recHitLabels.push_back(edm::InputTag("hltESRegionalEgammaRecHit:EcalRecHitsES"));
   desc.add<std::vector<edm::InputTag>>("recHitLabels", recHitLabels);
   std::vector<edm::ParameterSet> l1InputRegions;
 
   edm::ParameterSet emIsoPSet;
   emIsoPSet.addParameter<std::string>("type", "L1EmParticle");
   emIsoPSet.addParameter<double>("minEt", 5);
   emIsoPSet.addParameter<double>("maxEt", 999);
   emIsoPSet.addParameter<double>("regionEtaMargin", 0.14);
   emIsoPSet.addParameter<double>("regionPhiMargin", 0.4);
   emIsoPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltL1extraParticles:NonIsolated"));
   l1InputRegions.push_back(emIsoPSet);
   edm::ParameterSet emNonIsoPSet;
   emNonIsoPSet.addParameter<std::string>("type", "L1EmParticle");
   emNonIsoPSet.addParameter<double>("minEt", 5);
   emNonIsoPSet.addParameter<double>("maxEt", 999);
   emNonIsoPSet.addParameter<double>("regionEtaMargin", 0.14);
   emNonIsoPSet.addParameter<double>("regionPhiMargin", 0.4);
   emNonIsoPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltL1extraParticles:Isolated"));
   l1InputRegions.push_back(emNonIsoPSet);
 
   // Why no Central Jets here? They are present in the python config, e.g. OnLine_HLT_GRun.py
   // SHarper: because these are the default parameters designed to reproduce the original (no jets) behaviour
   //
   edm::ParameterSet egPSet;
   egPSet.addParameter<std::string>("type", "EGamma");
   egPSet.addParameter<double>("minEt", 5);
   egPSet.addParameter<double>("maxEt", 999);
   egPSet.addParameter<double>("regionEtaMargin", 0.4);
   egPSet.addParameter<double>("regionPhiMargin", 0.5);
   egPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltCaloStage2Digis"));
   l1InputRegions.push_back(egPSet);
 
   edm::ParameterSet jetPSet;
   jetPSet.addParameter<std::string>("type", "EGamma");
   jetPSet.addParameter<double>("minEt", 200);
   jetPSet.addParameter<double>("maxEt", 999);
   jetPSet.addParameter<double>("regionEtaMargin", 0.4);
   jetPSet.addParameter<double>("regionPhiMargin", 0.5);
   jetPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltCaloStage2Digis"));
   l1InputRegions.push_back(jetPSet);
 
   edm::ParameterSetDescription l1InputRegionDesc;
   l1InputRegionDesc.add<std::string>("type");
   l1InputRegionDesc.add<double>("minEt");
   l1InputRegionDesc.add<double>("maxEt");
   l1InputRegionDesc.add<double>("regionEtaMargin");
   l1InputRegionDesc.add<double>("regionPhiMargin");
   l1InputRegionDesc.add<edm::InputTag>("inputColl");
   desc.addVPSet("l1InputRegions", l1InputRegionDesc, l1InputRegions);
 
   descriptions.add(defaultModuleLabel<HLTRecHitInAllL1RegionsProducer<RecHitType>>(), desc);
 }

◆ produce()

template<typename RecHitType >

void HLTRecHitInAllL1RegionsProducer< RecHitType >::produce	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

override

Definition at line 219 of file HLTRecHitInAllL1RegionsProducer.cc.

References edmPickEvents::event, input, eostools::move(), Skims_PA_cff::name, rpcPointValidation_cfi::recHit, FastTrackerRecHitMaskProducer_cfi::recHits, HLT_2023v12_cff::region, singleTopDQM_cfi::setup, and GeomDetEnumerators::subDetGeom.

                                                                                                      {
   // get the collection geometry:
   auto const& caloGeom = setup.getData(caloGeometryToken_);
 
   std::vector<RectangularEtaPhiRegion> regions;
   std::for_each(l1RegionData_.begin(),
                 l1RegionData_.end(),
                 [&event, &setup, &regions](const std::unique_ptr<L1RegionDataBase>& input) {
                   input->getEtaPhiRegions(event, setup, regions);
                 });
 
   for (size_t recHitCollNr = 0; recHitCollNr < recHitTokens_.size(); recHitCollNr++) {
     edm::Handle<RecHitCollectionType> recHits;
     event.getByToken(recHitTokens_[recHitCollNr], recHits);
 
     if (!(recHits.isValid())) {
       edm::LogError("ProductNotFound") << "could not get a handle on the " << typeid(RecHitCollectionType).name()
                                        << " named " << recHitLabels_[recHitCollNr].encode() << std::endl;
       continue;
     }
 
     auto filteredRecHits = std::make_unique<RecHitCollectionType>();
 
     if (!recHits->empty()) {
       const CaloSubdetectorGeometry* subDetGeom = caloGeom.getSubdetectorGeometry(recHits->front().id());
       if (!regions.empty()) {
         for (const RecHitType& recHit : *recHits) {
           auto this_cell = subDetGeom->getGeometry(recHit.id());
           for (const auto& region : regions) {
             if (region.inRegion(this_cell->etaPos(), this_cell->phiPos())) {
               filteredRecHits->push_back(recHit);
               break;
             }
           }
         }
       }  //end check of empty regions
     }    //end check of empty rec-hits
     //   std::cout <<"putting fileter coll in "<<filteredRecHits->size()<<std::endl;
     event.put(std::move(filteredRecHits), productLabels_[recHitCollNr]);
   }  //end loop over all rec hit collections
 }