#include <EgammaHLTBcHcalIsolationProducersRegional.h>

Inheritance diagram for EgammaHLTBcHcalIsolationProducersRegional:

Public Member Functions
	EgammaHLTBcHcalIsolationProducersRegional (const edm::ParameterSet &)

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

	~EgammaHLTBcHcalIsolationProducersRegional ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
	EgammaHLTBcHcalIsolationProducersRegional (const EgammaHLTBcHcalIsolationProducersRegional &rhs)

EgammaHLTBcHcalIsolationProducersRegional &	operator= (const EgammaHLTBcHcalIsolationProducersRegional &rhs)

Private Attributes
edm::EDGetTokenT < CaloTowerCollection >	caloTowerProducer_

int	depth_

bool	doEtSum_

bool	doRhoCorrection_

float	effectiveAreaBarrel_

float	effectiveAreaEndcap_

float	etMin_

ElectronHcalHelper::Configuration	hcalCfg_

ElectronHcalHelper *	hcalHelper_

float	innerCone_

float	outerCone_

edm::EDGetTokenT < reco::RecoEcalCandidateCollection >	recoEcalCandidateProducer_

float	rhoMax_

edm::EDGetTokenT< double >	rhoProducer_

float	rhoScale_

bool	useSingleTower_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

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

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

Definition at line 34 of file EgammaHLTBcHcalIsolationProducersRegional.h.

Constructor & Destructor Documentation

EgammaHLTBcHcalIsolationProducersRegional::EgammaHLTBcHcalIsolationProducersRegional ( const edm::ParameterSet & config )

explicit

Definition at line 26 of file EgammaHLTBcHcalIsolationProducersRegional.cc.

References caloTowerProducer_, depth_, doEtSum_, doRhoCorrection_, effectiveAreaBarrel_, effectiveAreaEndcap_, etMin_, edm::ParameterSet::getParameter(), hcalCfg_, hcalHelper_, ElectronHcalHelper::Configuration::hcalTowers, ElectronHcalHelper::Configuration::hOverEConeSize, ElectronHcalHelper::Configuration::hOverEPtMin, innerCone_, outerCone_, recoEcalCandidateProducer_, rhoMax_, rhoProducer_, rhoScale_, useSingleTower_, and ElectronHcalHelper::Configuration::useTowers.

                                                                                                                   {
 
   recoEcalCandidateProducer_ = consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"));
   caloTowerProducer_         = consumes<CaloTowerCollection>(config.getParameter<edm::InputTag>("caloTowerProducer"));
 
   doRhoCorrection_           = config.getParameter<bool>("doRhoCorrection");
   if (doRhoCorrection_)
     rhoProducer_               = consumes<double>(config.getParameter<edm::InputTag>("rhoProducer"));
 
   rhoMax_                    = config.getParameter<double>("rhoMax"); 
   rhoScale_                  = config.getParameter<double>("rhoScale"); 
 
   etMin_                     = config.getParameter<double>("etMin");  
   innerCone_                 = config.getParameter<double>("innerCone");
   outerCone_                 = config.getParameter<double>("outerCone");
   depth_                     = config.getParameter<int>("depth");
   doEtSum_                   = config.getParameter<bool>("doEtSum"); 
   effectiveAreaBarrel_       = config.getParameter<double>("effectiveAreaBarrel");
   effectiveAreaEndcap_       = config.getParameter<double>("effectiveAreaEndcap");
 
   useSingleTower_            = config.getParameter<bool>("useSingleTower");
   
   hcalCfg_.hOverEConeSize = outerCone_;
   hcalCfg_.useTowers = true;
   hcalCfg_.hcalTowers = caloTowerProducer_;
   hcalCfg_.hOverEPtMin = etMin_;
 
   hcalHelper_ = new ElectronHcalHelper(hcalCfg_);
 
   produces <reco::RecoEcalCandidateIsolationMap>(); 
 }

EgammaHLTBcHcalIsolationProducersRegional::~EgammaHLTBcHcalIsolationProducersRegional ( )

Definition at line 58 of file EgammaHLTBcHcalIsolationProducersRegional.cc.

References hcalHelper_.

                                                                                       {
   delete hcalHelper_;
 }

EgammaHLTBcHcalIsolationProducersRegional::EgammaHLTBcHcalIsolationProducersRegional ( const EgammaHLTBcHcalIsolationProducersRegional & rhs )

inlineprivate

Definition at line 40 of file EgammaHLTBcHcalIsolationProducersRegional.h.

40 {}

Member Function Documentation

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

static

Definition at line 62 of file EgammaHLTBcHcalIsolationProducersRegional.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                                            {
 
   edm::ParameterSetDescription desc;
   
   desc.add<edm::InputTag>(("recoEcalCandidateProducer"), edm::InputTag("hltRecoEcalCandidate"));
   desc.add<edm::InputTag>(("caloTowerProducer"), edm::InputTag("hltTowerMakerForAll"));
   desc.add<edm::InputTag>(("rhoProducer"), edm::InputTag("fixedGridRhoFastjetAllCalo"));
   desc.add<bool>(("doRhoCorrection"), false);
   desc.add<double>(("rhoMax"), 999999.); 
   desc.add<double>(("rhoScale"), 1.0); 
   desc.add<double>(("etMin"), -1.0);  
   desc.add<double>(("innerCone"), 0);
   desc.add<double>(("outerCone"), 0.15);
   desc.add<int>(("depth"), -1);
   desc.add<bool>(("doEtSum"), false); 
   desc.add<double>(("effectiveAreaBarrel"), 0.021);
   desc.add<double>(("effectiveAreaEndcap"), 0.040);
   desc.add<bool>(("useSingleTower"), false);
   descriptions.add(("hltEgammaHLTBcHcalIsolationProducersRegional"), desc);  
 }

EgammaHLTBcHcalIsolationProducersRegional& EgammaHLTBcHcalIsolationProducersRegional::operator= ( const EgammaHLTBcHcalIsolationProducersRegional & rhs )

inlineprivate

Definition at line 41 of file EgammaHLTBcHcalIsolationProducersRegional.h.

41 { return *this; }

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

virtual

Implements edm::EDProducer.

Definition at line 83 of file EgammaHLTBcHcalIsolationProducersRegional.cc.

References caloTowerProducer_, ElectronHcalHelper::checkSetup(), depth_, doEtSum_, doRhoCorrection_, effectiveAreaBarrel_, effectiveAreaEndcap_, etMin_, edm::Event::getByToken(), EgammaTowerIsolation::getTowerEtSum(), ElectronHcalHelper::hcalESum(), ElectronHcalHelper::hcalESumDepth1BehindClusters(), ElectronHcalHelper::hcalESumDepth2BehindClusters(), hcalHelper_, ElectronHcalHelper::hcalTowersBehindClusters(), innerCone_, edm::AssociationMap< Tag >::insert(), outerCone_, edm::Handle< T >::product(), edm::Event::put(), ElectronHcalHelper::readEvent(), recoEcalCandidateProducer_, rho, rhoMax_, rhoProducer_, rhoScale_, and useSingleTower_.

                                                                                                      {
 
   // Get the HLT filtered objects
   edm::Handle<reco::RecoEcalCandidateCollection> recoEcalCandHandle;
   iEvent.getByToken(recoEcalCandidateProducer_, recoEcalCandHandle);
 
   edm::Handle<CaloTowerCollection> caloTowersHandle;
   iEvent.getByToken(caloTowerProducer_, caloTowersHandle);
 
   edm::Handle<double> rhoHandle;
   double rho = 0.0;
 
   if (doRhoCorrection_) {
     iEvent.getByToken(rhoProducer_, rhoHandle);
     rho = *(rhoHandle.product());
   }
 
   if (rho > rhoMax_)
     rho = rhoMax_;
   
   rho = rho*rhoScale_;
   
   hcalHelper_->checkSetup(iSetup);
   hcalHelper_->readEvent(iEvent);
 
   reco::RecoEcalCandidateIsolationMap isoMap;
   
   for(unsigned int iRecoEcalCand=0; iRecoEcalCand <recoEcalCandHandle->size(); iRecoEcalCand++) {
     
     reco::RecoEcalCandidateRef recoEcalCandRef(recoEcalCandHandle, iRecoEcalCand);
     
     float isol = 0;
     
     std::vector<CaloTowerDetId> towersBehindCluster;
     
     if (useSingleTower_)
       towersBehindCluster = hcalHelper_->hcalTowersBehindClusters(*(recoEcalCandRef->superCluster()));
     
     if (doEtSum_) { //calculate hcal isolation excluding the towers behind the cluster which will be used for H for H/E
       EgammaTowerIsolation isolAlgo(outerCone_, innerCone_, etMin_, depth_, caloTowersHandle.product());
       if (useSingleTower_) 
     isol = isolAlgo.getTowerEtSum(&(*recoEcalCandRef), &(towersBehindCluster)); // towersBehindCluster are excluded from the isolation sum
       else
     isol = isolAlgo.getTowerEtSum(&(*recoEcalCandRef));
       
       if (doRhoCorrection_) {
     if (fabs(recoEcalCandRef->superCluster()->eta()) < 1.442) 
       isol = isol - rho*effectiveAreaBarrel_;
     else
       isol = isol - rho*effectiveAreaEndcap_;
       }
     } else { //calcuate H for H/E
       if (useSingleTower_) 
     isol = hcalHelper_->hcalESumDepth1BehindClusters(towersBehindCluster) + hcalHelper_->hcalESumDepth2BehindClusters(towersBehindCluster);
       else
     isol = hcalHelper_->hcalESum(*(recoEcalCandRef->superCluster()));
     }
 
     isoMap.insert(recoEcalCandRef, isol);
   }
   
   std::auto_ptr<reco::RecoEcalCandidateIsolationMap> isolMap(new reco::RecoEcalCandidateIsolationMap(isoMap));
   iEvent.put(isolMap);
 }