#include <EgammaIsoHcalDetIdCollectionProducer.h>

Inheritance diagram for EgammaIsoHcalDetIdCollectionProducer:

Public Member Functions
void	beginRun (edm::Run const &, const edm::EventSetup &) final

	EgammaIsoHcalDetIdCollectionProducer (const edm::ParameterSet &)

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

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Private Attributes
edm::EDGetTokenT< reco::GsfElectronCollection >	elesToken_

EGHcalRecHitSelector	hcalHitSelector_

std::string	interestingDetIdCollection_

float	minEleEt_

float	minPhoEt_

float	minSCEt_

edm::EDGetTokenT< reco::PhotonCollection >	phosToken_

edm::EDGetTokenT< HBHERecHitCollection >	recHitsToken_

edm::EDGetTokenT< reco::SuperClusterCollection >	superClustersToken_

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

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Original author: Sam Harper (RAL)

Make a collection of detids to be kept tipically in a AOD rechit collection Modified from the ECAL version "InterestingDetIdCollectionProducer" to be HCAL

Definition at line 39 of file EgammaIsoHcalDetIdCollectionProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file EgammaIsoHcalDetIdCollectionProducer.cc.

References elesToken_, edm::ParameterSet::getParameter(), interestingDetIdCollection_, minEleEt_, minPhoEt_, minSCEt_, phosToken_, recHitsToken_, AlCaHLTBitMon_QueryRunRegistry::string, and superClustersToken_.

     : hcalHitSelector_(iConfig.getParameter<edm::ParameterSet>("hitSelection")) {
   recHitsToken_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("recHitsLabel"));
   elesToken_ = consumes<reco::GsfElectronCollection>(iConfig.getParameter<edm::InputTag>("elesLabel"));
 
   phosToken_ = consumes<reco::PhotonCollection>(iConfig.getParameter<edm::InputTag>("phosLabel"));
 
   superClustersToken_ =
       consumes<reco::SuperClusterCollection>(iConfig.getParameter<edm::InputTag>("superClustersLabel"));
 
   minSCEt_ = iConfig.getParameter<double>("minSCEt");
   minEleEt_ = iConfig.getParameter<double>("minEleEt");
   minPhoEt_ = iConfig.getParameter<double>("minPhoEt");
 
   interestingDetIdCollection_ = iConfig.getParameter<std::string>("interestingDetIdCollection");
 
   //register your products
   produces<DetIdCollection>(interestingDetIdCollection_);
 }

Member Function Documentation

void EgammaIsoHcalDetIdCollectionProducer::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

final

Definition at line 43 of file EgammaIsoHcalDetIdCollectionProducer.cc.

References hcalHitSelector_, and EGHcalRecHitSelector::setup().

                                                                                                   {
   hcalHitSelector_.setup(iSetup);
 }

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

static

Definition at line 29 of file EgammaIsoHcalDetIdCollectionProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_2018_cff::InputTag, EGHcalRecHitSelector::makePSetDescription(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("recHitsLabel", edm::InputTag("hbhereco"));
   desc.add<edm::InputTag>("elesLabel", edm::InputTag("gedGsfElectrons"));
   desc.add<edm::InputTag>("phosLabel", edm::InputTag("gedPhotons"));
   desc.add<edm::InputTag>("superClustersLabel", edm::InputTag("particleFlowEGamma"));
   desc.add<double>("minSCEt", 20);
   desc.add<double>("minEleEt", 20);
   desc.add<double>("minPhoEt", 20);
   desc.add<std::string>("interestingDetIdCollection", "");
   desc.add<edm::ParameterSetDescription>("hitSelection", EGHcalRecHitSelector::makePSetDescription());
   descriptions.add(("interestingGedEgammaIsoHCALDetId"), desc);
 }

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

override

Definition at line 48 of file EgammaIsoHcalDetIdCollectionProducer.cc.

References EGHcalRecHitSelector::addDetIds(), elesToken_, SimCluster::energy(), edm::Event::getByToken(), hcalHitSelector_, interestingDetIdCollection_, edm::HandleBase::isValid(), minEleEt_, minPhoEt_, minSCEt_, eostools::move(), phosToken_, edm::Event::put(), FastTrackerRecHitMaskProducer_cfi::recHits, recHitsToken_, SimDataFormats::CaloAnalysis::sc, funct::sin(), HLT_2018_cff::superClusters, superClustersToken_, SimCluster::theta(), and tier0::unique().

                                                                                                 {
   edm::Handle<reco::SuperClusterCollection> superClusters;
   iEvent.getByToken(superClustersToken_, superClusters);
 
   edm::Handle<reco::GsfElectronCollection> eles;
   iEvent.getByToken(elesToken_, eles);
 
   edm::Handle<reco::PhotonCollection> phos;
   iEvent.getByToken(phosToken_, phos);
 
   edm::Handle<HBHERecHitCollection> recHits;
   iEvent.getByToken(recHitsToken_, recHits);
 
   std::vector<DetId> indexToStore;
   indexToStore.reserve(100);
 
   if (eles.isValid() && recHits.isValid()) {
     for (auto& ele : *eles) {
       float scEt = ele.superCluster()->energy() * std::sin(ele.superCluster()->position().theta());
       if (scEt > minEleEt_ || ele.et() > minEleEt_) {
         hcalHitSelector_.addDetIds(*ele.superCluster(), *recHits, indexToStore);
       }
     }
   }
   if (phos.isValid() && recHits.isValid()) {
     for (auto& pho : *phos) {
       float scEt = pho.superCluster()->energy() * std::sin(pho.superCluster()->position().theta());
       if (scEt > minPhoEt_ || pho.et() > minPhoEt_) {
         hcalHitSelector_.addDetIds(*pho.superCluster(), *recHits, indexToStore);
       }
     }
   }
   if (superClusters.isValid() && recHits.isValid()) {
     for (auto& sc : *superClusters) {
       float scEt = sc.energy() * std::sin(sc.position().theta());
       if (scEt > minSCEt_) {
         hcalHitSelector_.addDetIds(sc, *recHits, indexToStore);
       }
     }
   }
 
   //unify the vector
   std::sort(indexToStore.begin(), indexToStore.end());
   std::unique(indexToStore.begin(), indexToStore.end());
 
   auto detIdCollection = std::make_unique<DetIdCollection>(indexToStore);
 
   iEvent.put(std::move(detIdCollection), interestingDetIdCollection_);
 }