#include <EleIsoDetIdCollectionProducer.h>

Inheritance diagram for EleIsoDetIdCollectionProducer:

Public Member Functions
void	beginJob ()

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

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

	~EleIsoDetIdCollectionProducer ()

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 Attributes
edm::InputTag	emObjectLabel_

edm::EDGetToken	emObjectToken_

double	energyCut_

double	etCandCut_

double	etCut_

std::vector< int >	flagsexclEB_

std::vector< int >	flagsexclEE_

double	innerRadius_

std::string	interestingDetIdCollection_

double	outerRadius_

edm::InputTag	recHitsLabel_

edm::EDGetToken	recHitsToken_

std::vector< int >	severitiesexclEB_

std::vector< int >	severitiesexclEE_

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: Matthew LeBourgeois PH/CMG Modified from : RecoEcal/EgammaClusterProducers/{src,interface}/InterestingDetIdCollectionProducer.{h,cc} by Paolo Meridiani PH/CMG

Implementation: <Notes on="" implementation>="">

Definition at line 34 of file EleIsoDetIdCollectionProducer.h.

Constructor & Destructor Documentation

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

explicit

ctor

Definition at line 31 of file EleIsoDetIdCollectionProducer.cc.

References flagsexclEB_, flagsexclEE_, edm::ParameterSet::getParameter(), interestingDetIdCollection_, severitiesexclEB_, and severitiesexclEE_.

                                                                                            :
             recHitsToken_(consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag > ("recHitsLabel"))),
         emObjectToken_(consumes<reco::GsfElectronCollection>(iConfig.getParameter< edm::InputTag > ("emObjectLabel"))),
         recHitsLabel_(iConfig.getParameter< edm::InputTag > ("recHitsLabel")),
         emObjectLabel_(iConfig.getParameter< edm::InputTag > ("emObjectLabel")),
             energyCut_(iConfig.getParameter<double>("energyCut")),
             etCut_(iConfig.getParameter<double>("etCut")),
             etCandCut_(iConfig.getParameter<double> ("etCandCut")),
             outerRadius_(iConfig.getParameter<double>("outerRadius")),
             innerRadius_(iConfig.getParameter<double>("innerRadius")),
             interestingDetIdCollection_(iConfig.getParameter<std::string>("interestingDetIdCollection"))
  {
 
    const std::vector<std::string> flagnamesEB = 
      iConfig.getParameter<std::vector<std::string> >("RecHitFlagToBeExcludedEB");
    
    const std::vector<std::string> flagnamesEE =
      iConfig.getParameter<std::vector<std::string> >("RecHitFlagToBeExcludedEE");
    
    flagsexclEB_= 
      StringToEnumValue<EcalRecHit::Flags>(flagnamesEB);
    
    flagsexclEE_=
      StringToEnumValue<EcalRecHit::Flags>(flagnamesEE);
    
    const std::vector<std::string> severitynamesEB = 
      iConfig.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcludedEB");
    
    severitiesexclEB_= 
      StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEB);
 
    const std::vector<std::string> severitynamesEE = 
      iConfig.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcludedEE");
    
    severitiesexclEE_= 
      StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEE);
    
    //register your products
    produces< DetIdCollection > (interestingDetIdCollection_) ;
 }

EleIsoDetIdCollectionProducer::~EleIsoDetIdCollectionProducer ( )

Definition at line 72 of file EleIsoDetIdCollectionProducer.cc.

73 {}

Member Function Documentation

void EleIsoDetIdCollectionProducer::beginJob ( void )

Definition at line 75 of file EleIsoDetIdCollectionProducer.cc.

76 {}

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

virtual

producer

Implements edm::stream::EDProducerBase.

Definition at line 80 of file EleIsoDetIdCollectionProducer.cc.

References EcalRecHit::checkFlag(), EcalRecHit::checkFlags(), EcalRecHit::detid(), DetId::Ecal, EcalBarrel, EcalEndcap, emObjectToken_, EcalRecHit::energy(), energyCut_, eta, etCandCut_, etCut_, spr::find(), flagsexclEB_, flagsexclEE_, edm::EventSetup::get(), edm::Event::getByToken(), innerRadius_, edm::InputTag::instance(), interestingDetIdCollection_, GeomDetEnumerators::isBarrel(), EcalRecHit::kGood, mag(), outerRadius_, perp(), edm::Event::put(), recHitsLabel_, recHitsToken_, severitiesexclEB_, and severitiesexclEE_.

                                                                                      {
     using namespace edm;
     using namespace std;
 
     //Get EM Object
     Handle<reco::GsfElectronCollection> emObjectH;
     iEvent.getByToken(emObjectToken_,emObjectH);
 
     // take EcalRecHits
     Handle<EcalRecHitCollection> recHitsH;
     iEvent.getByToken(recHitsToken_,recHitsH);
 
     edm::ESHandle<CaloGeometry> pG;
     iSetup.get<CaloGeometryRecord>().get(pG);    
     const CaloGeometry* caloGeom = pG.product();
 
     //Get the channel status from the db
     //edm::ESHandle<EcalChannelStatus> chStatus;
     //iSetup.get<EcalChannelStatusRcd>().get(chStatus);
 
     edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
     iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
     const EcalSeverityLevelAlgo* sevLevel = sevlv.product();
 
     CaloDualConeSelector<EcalRecHit> *doubleConeSel_ = 0;
     if(recHitsLabel_.instance() == "EcalRecHitsEB")
         doubleConeSel_= new CaloDualConeSelector<EcalRecHit>(innerRadius_,outerRadius_, &*pG, DetId::Ecal, EcalBarrel);
     else if(recHitsLabel_.instance() == "EcalRecHitsEE")
         doubleConeSel_= new CaloDualConeSelector<EcalRecHit>(innerRadius_,outerRadius_, &*pG, DetId::Ecal, EcalEndcap);
 
     //Create empty output collections
     std::auto_ptr< DetIdCollection > detIdCollection (new DetIdCollection() ) ;
 
     reco::GsfElectronCollection::const_iterator emObj;
     if(doubleConeSel_) { //if cone selector was created
         for (emObj = emObjectH->begin(); emObj != emObjectH->end();  emObj++) { //Loop over candidates
 
             if(emObj->et() < etCandCut_) continue; //don't calculate if object hasn't enough energy
             
             GlobalPoint pclu (emObj->caloPosition().x(),emObj->caloPosition().y(),emObj->caloPosition().z());
             doubleConeSel_->selectCallback(pclu, *recHitsH, [&](const EcalRecHit& recIt) {
               if (recIt.energy() < energyCut_) 
                 return;  //dont fill if below E noise value
 
               double et = recIt.energy() * 
                 caloGeom->getPosition(recIt.detid()).perp() / 
                 caloGeom->getPosition(recIt.detid()).mag();
 
               bool isBarrel = false;
               if (fabs(caloGeom->getPosition(recIt.detid()).eta() < 1.479)) 
                 isBarrel = true;
               
               if (et < etCut_) 
                 return;  //dont fill if below ET noise value
 
               std::vector<int>::const_iterator sit;
               int severityFlag = sevLevel->severityLevel(recIt.detid(), *recHitsH);
               if (isBarrel) {
                 sit = std::find(severitiesexclEB_.begin(), severitiesexclEB_.end(), severityFlag);
                 if (sit != severitiesexclEB_.end())
                   return;
               } else {
                 sit = std::find(severitiesexclEE_.begin(), severitiesexclEE_.end(), severityFlag);
                 if (sit != severitiesexclEE_.end())
                   return;
               }
 
               std::vector<int>::const_iterator vit;
               if (isBarrel) {
                 // new rechit flag checks
                 //vit = std::find(flagsexclEB_.begin(), flagsexclEB_.end(), ((EcalRecHit*)(&*recIt))->recoFlag());
                 //if (vit != flagsexclEB_.end())
                 //  continue;
                 if (!(recIt.checkFlag(EcalRecHit::kGood))) {
                   if (recIt.checkFlags(flagsexclEB_)) {                
                     return;
                   }
                 }
               } else {
                 // new rechit flag checks
                 //vit = std::find(flagsexclEE_.begin(), flagsexclEE_.end(), ((EcalRecHit*)(&*recIt))->recoFlag());
                 //if (vit != flagsexclEE_.end())
                 //  continue;
                 if (!(recIt.checkFlag(EcalRecHit::kGood))) {
                   if (recIt.checkFlags(flagsexclEE_)) {                
                     return;
                   }
                 }
               }
 
               if (std::find(detIdCollection->begin(),detIdCollection->end(),recIt.detid()) == detIdCollection->end()) 
                 detIdCollection->push_back(recIt.detid()); 
             }); //end rechits
             
         } //end candidates
         
         delete doubleConeSel_;
     } //end if cone selector was created
     
     iEvent.put( detIdCollection, interestingDetIdCollection_ );
 }