#include <EgammaIsoHcalDetIdCollectionProducer.h>

Inheritance diagram for EgammaIsoHcalDetIdCollectionProducer:

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

	EgammaIsoHcalDetIdCollectionProducer (const edm::ParameterSet &)
	ctor More...

virtual void	produce (edm::Event &, const edm::EventSetup &)
	producer More...

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	addDetIds (const reco::SuperCluster &superClus, const HBHERecHitCollection &recHits, std::vector< DetId > &detIdsToStore)

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

std::string	interestingDetIdCollection_

int	maxDIEta_

int	maxDIPhi_

float	minEleEt_

float	minEnergyHCAL_

float	minPhoEt_

float	minSCEt_

edm::EDGetTokenT < reco::PhotonCollection >	phosToken_

edm::EDGetTokenT < HBHERecHitCollection >	recHitsToken_

edm::EDGetTokenT < reco::SuperClusterCollection >	superClustersToken_

edm::ESHandle < CaloTowerConstituentsMap >	towerMap_

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 44 of file EgammaIsoHcalDetIdCollectionProducer.h.

Constructor & Destructor Documentation

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

explicit

ctor

Definition at line 13 of file EgammaIsoHcalDetIdCollectionProducer.cc.

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

 {
 
   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");
   
   maxDIEta_ = iConfig.getParameter<int>("maxDIEta");
   maxDIPhi_ = iConfig.getParameter<int>("maxDIPhi");
 
   minEnergyHCAL_= iConfig.getParameter<double>("minEnergyHCAL");
   
    //register your products
   produces< DetIdCollection > (interestingDetIdCollection_) ;
 
 
 }

Member Function Documentation

void EgammaIsoHcalDetIdCollectionProducer::addDetIds	(	const reco::SuperCluster &	superClus,
		const HBHERecHitCollection &	recHits,
		std::vector< DetId > &	detIdsToStore
	)

private

Definition at line 132 of file EgammaIsoHcalDetIdCollectionProducer.cc.

References funct::abs(), calDIEta(), calDIPhi(), DetId::det(), DetId::Ecal, maxDIEta_, maxDIPhi_, minEnergyHCAL_, reco::SuperCluster::seed(), ecaldqm::towerId(), and towerMap_.

Referenced by produce().

 {
   DetId seedId = superClus.seed()->seed();
   if(seedId.det() != DetId::Ecal) {
     edm::LogError("EgammaIsoHcalDetIdCollectionProducerError") << "Somehow the supercluster has a seed which is not ECAL, something is badly wrong";
   }
   //so we are using CaloTowers to get the iEta/iPhi of the HCAL rec hit behind the seed cluster, this might go funny on tower 28 but shouldnt matter there
  
   CaloTowerDetId towerId(towerMap_->towerOf(seedId)); 
   int seedHcalIEta = towerId.ieta();
   int seedHcalIPhi = towerId.iphi();
 
   for(auto& recHit : recHits){
     int dIEtaAbs = std::abs(calDIEta(seedHcalIEta,recHit.id().ieta()));
     int dIPhiAbs = std::abs(calDIPhi(seedHcalIPhi,recHit.id().iphi()));
   
     if(dIEtaAbs<=maxDIEta_ && dIPhiAbs<=maxDIPhi_ && recHit.energy()>minEnergyHCAL_) detIdsToStore.push_back(recHit.id().rawId());
   }
 
 }

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

finaloverridevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 46 of file EgammaIsoHcalDetIdCollectionProducer.cc.

References edm::EventSetup::get(), and towerMap_.

 {
    iSetup.get<IdealGeometryRecord>().get(towerMap_);
    //  std::cout <<" got geom "<<towerMap_.isValid()<<" "<<&(*towerMap_)<<std::endl;
 }

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

virtual

producer

Implements edm::stream::EDProducerBase.

Definition at line 54 of file EgammaIsoHcalDetIdCollectionProducer.cc.

References addDetIds(), elesToken_, edm::Event::getByToken(), interestingDetIdCollection_, edm::HandleBase::isValid(), minEleEt_, minPhoEt_, minSCEt_, phosToken_, edm::Event::put(), HLT_25ns14e33_v3_cff::recHits, recHitsToken_, funct::sin(), python.multivaluedict::sort(), and superClustersToken_.

 {
 
    // take BasicClusters
   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);
 
   //Create empty output collections
   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_) 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_) 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_) addDetIds(sc,*recHits,indexToStore);
     }
   }
   
   //unify the vector
   std::sort(indexToStore.begin(),indexToStore.end());
   std::unique(indexToStore.begin(),indexToStore.end());
   
   std::auto_ptr< DetIdCollection > detIdCollection (new DetIdCollection(indexToStore) ) ;
    
   iEvent.put( detIdCollection, interestingDetIdCollection_ );
 
 }