#include <HLTHcalPFClusterIsolationProducer.h>

Inheritance diagram for HLTHcalPFClusterIsolationProducer< T1 >:

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

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

	~HLTHcalPFClusterIsolationProducer ()

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 Types
typedef std::vector< T1 >	T1Collection

typedef edm::AssociationMap < edm::OneToValue< std::vector < T1 >, float > >	T1IsolationMap

typedef edm::Ref< T1Collection >	T1Ref

Private Attributes
bool	doRhoCorrection_

double	drMax_

double	drVetoBarrel_

double	drVetoEndcap_

float	effectiveAreaBarrel_

float	effectiveAreaEndcap_

double	energyBarrel_

double	energyEndcap_

double	etaStripBarrel_

double	etaStripEndcap_

edm::EDGetTokenT < reco::PFClusterCollection >	pfClusterProducerHCAL_

edm::EDGetTokenT < reco::PFClusterCollection >	pfClusterProducerHFEM_

edm::EDGetTokenT < reco::PFClusterCollection >	pfClusterProducerHFHAD_

edm::EDGetTokenT< T1Collection >	recoCandidateProducer_

float	rhoMax_

edm::EDGetTokenT< double >	rhoProducer_

float	rhoScale_

bool	useHF_

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

template<typename T1>
class HLTHcalPFClusterIsolationProducer< T1 >

Definition at line 29 of file HLTHcalPFClusterIsolationProducer.h.

Member Typedef Documentation

template<typename T1 >

typedef std::vector<T1> HLTHcalPFClusterIsolationProducer< T1 >::T1Collection

private

Definition at line 31 of file HLTHcalPFClusterIsolationProducer.h.

template<typename T1 >

typedef edm::AssociationMap<edm::OneToValue<std::vector<T1>, float > > HLTHcalPFClusterIsolationProducer< T1 >::T1IsolationMap

private

Definition at line 33 of file HLTHcalPFClusterIsolationProducer.h.

template<typename T1 >

typedef edm::Ref<T1Collection> HLTHcalPFClusterIsolationProducer< T1 >::T1Ref

private

Definition at line 32 of file HLTHcalPFClusterIsolationProducer.h.

Constructor & Destructor Documentation

template<typename T1 >

HLTHcalPFClusterIsolationProducer< T1 >::HLTHcalPFClusterIsolationProducer ( const edm::ParameterSet & config )

explicit

Definition at line 20 of file HLTHcalPFClusterIsolationProducer.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                       {
     
   std::string recoCandidateProducerName = "recoCandidateProducer";
   if ((typeid(HLTHcalPFClusterIsolationProducer<T1>) == typeid(HLTHcalPFClusterIsolationProducer<reco::RecoEcalCandidate>))) recoCandidateProducerName = "recoEcalCandidateProducer";
   
   recoCandidateProducer_ = consumes<T1Collection>(config.getParameter<edm::InputTag>(recoCandidateProducerName));
   pfClusterProducerHCAL_     = consumes<reco::PFClusterCollection>(config.getParameter<edm::InputTag>("pfClusterProducerHCAL"));
   useHF_                     = config.getParameter<bool>("useHF");  
   if (useHF_) {
     pfClusterProducerHFEM_     = consumes<reco::PFClusterCollection>(config.getParameter<edm::InputTag>("pfClusterProducerHFEM"));
     pfClusterProducerHFHAD_    = consumes<reco::PFClusterCollection>(config.getParameter<edm::InputTag>("pfClusterProducerHFHAD"));
   }
 
   drMax_          = config.getParameter<double>("drMax");
   drVetoBarrel_   = config.getParameter<double>("drVetoBarrel");
   drVetoEndcap_   = config.getParameter<double>("drVetoEndcap");
   etaStripBarrel_ = config.getParameter<double>("etaStripBarrel");
   etaStripEndcap_ = config.getParameter<double>("etaStripEndcap");
   energyBarrel_   = config.getParameter<double>("energyBarrel");
   energyEndcap_   = config.getParameter<double>("energyEndcap");
 
   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"); 
   effectiveAreaBarrel_            = config.getParameter<double>("effectiveAreaBarrel");
   effectiveAreaEndcap_            = config.getParameter<double>("effectiveAreaEndcap");
 
   produces <T1IsolationMap >();
 }

template<typename T1 >

HLTHcalPFClusterIsolationProducer< T1 >::~HLTHcalPFClusterIsolationProducer ( )

Definition at line 54 of file HLTHcalPFClusterIsolationProducer.cc.

55 {}

Member Function Documentation

template<typename T1 >

void HLTHcalPFClusterIsolationProducer< T1 >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 58 of file HLTHcalPFClusterIsolationProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::ifValue(), mergeVDriftHistosByStation::name, or, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                        {
  
   std::string recoCandidateProducerName = "recoCandidateProducer";
   if ((typeid(HLTHcalPFClusterIsolationProducer<T1>) == typeid(HLTHcalPFClusterIsolationProducer<reco::RecoEcalCandidate>))) recoCandidateProducerName = "recoEcalCandidateProducer";
 
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>(recoCandidateProducerName, edm::InputTag("hltL1SeededRecoEcalCandidatePF"));
   desc.add<edm::InputTag>("pfClusterProducerHCAL", edm::InputTag("hltParticleFlowClusterHCAL"));
   desc.ifValue(edm::ParameterDescription<bool>("useHF", false, true),
            true >> (edm::ParameterDescription<edm::InputTag>("pfClusterProducerHFEM", edm::InputTag("hltParticleFlowClusterHFEM"), true) and
             edm::ParameterDescription<edm::InputTag>("pfClusterProducerHFHAD", edm::InputTag("hltParticleFlowClusterHFHAD"), true)) or
            false >> (edm::ParameterDescription<edm::InputTag>("pfClusterProducerHFEM", edm::InputTag(""), true) and
              edm::ParameterDescription<edm::InputTag>("pfClusterProducerHFHAD", edm::InputTag(""), true)));
   desc.add<edm::InputTag>("rhoProducer", edm::InputTag("fixedGridRhoFastjetAllCalo"));
   desc.add<bool>("doRhoCorrection", false);
   desc.add<double>("rhoMax", 9.9999999E7); 
   desc.add<double>("rhoScale", 1.0); 
   desc.add<double>("effectiveAreaBarrel", 0.101);
   desc.add<double>("effectiveAreaEndcap", 0.046);
   desc.add<double>("drMax", 0.3);
   desc.add<double>("drVetoBarrel", 0.0);
   desc.add<double>("drVetoEndcap", 0.0);
   desc.add<double>("etaStripBarrel", 0.0);
   desc.add<double>("etaStripEndcap", 0.0);
   desc.add<double>("energyBarrel", 0.0);
   desc.add<double>("energyEndcap", 0.0);
   descriptions.add(std::string("hlt")+std::string(typeid(HLTHcalPFClusterIsolationProducer<T1>).name()), desc);
 }

template<typename T1 >

void HLTHcalPFClusterIsolationProducer< T1 >::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::EDProducer.

Definition at line 88 of file HLTHcalPFClusterIsolationProducer.cc.

References deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::PFCluster::energy(), reco::CaloCluster::eta(), edm::Event::getByToken(), i, edm::AssociationMap< Tag >::insert(), reco::CaloCluster::phi(), edm::Handle< T >::product(), reco::PFCluster::pt(), edm::Event::put(), and rho.

                                                                                                 {
   
   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_;
 
   edm::Handle<T1Collection> recoCandHandle;
   edm::Handle<reco::PFClusterCollection> clusterHcalHandle;
   edm::Handle<reco::PFClusterCollection> clusterHfemHandle;
   edm::Handle<reco::PFClusterCollection> clusterHfhadHandle;
 
   iEvent.getByToken(recoCandidateProducer_,recoCandHandle);
   iEvent.getByToken(pfClusterProducerHCAL_, clusterHcalHandle);
   const reco::PFClusterCollection* forIsolationHcal = clusterHcalHandle.product();
 
   if (useHF_) {
     iEvent.getByToken(pfClusterProducerHFEM_, clusterHfemHandle);
     iEvent.getByToken(pfClusterProducerHFHAD_, clusterHfhadHandle);
   }
 
   T1IsolationMap recoCandMap;
   
   float dRVeto = -1.;
   float etaStrip = -1;
   
   for (unsigned int iReco = 0; iReco < recoCandHandle->size(); iReco++) {
     T1Ref candRef(recoCandHandle, iReco);
     
     if (fabs(candRef->eta()) < 1.479) {
       dRVeto = drVetoBarrel_;
       etaStrip = etaStripBarrel_;
     } else {
       dRVeto = drVetoEndcap_;
       etaStrip = etaStripEndcap_;
     }
     
     float sum = 0;
     
     // Loop over the 3 types of PFClusters
 
     for(unsigned i=0; i<forIsolationHcal->size(); i++) {
       const reco::PFCluster& pfclu = (*forIsolationHcal)[i];
       
       if (fabs(candRef->eta()) < 1.479) {
     if (fabs(pfclu.pt()) < energyBarrel_)
       continue;
       } else {
     if (fabs(pfclu.energy()) < energyEndcap_)
       continue;
       }
 
       float dEta = fabs(candRef->eta() - pfclu.eta());
       if(dEta < etaStrip) continue;
       
       float dR = deltaR(candRef->eta(), candRef->phi(), pfclu.eta(), pfclu.phi());
       if(dR > drMax_ || dR < dRVeto) continue;
       
       sum += pfclu.pt();
     }
 
     if (useHF_) {
       const reco::PFClusterCollection* forIsolationHfem = clusterHfemHandle.product();
       const reco::PFClusterCollection* forIsolationHfhad = clusterHfhadHandle.product();
       
       for(unsigned i=0; i<forIsolationHfem->size(); i++) {
     const reco::PFCluster& pfclu = (*forIsolationHfem)[i];
     
     if (fabs(candRef->eta()) < 1.479) {
       if (fabs(pfclu.pt()) < energyBarrel_)
         continue;
     } else {
       if (fabs(pfclu.energy()) < energyEndcap_)
         continue;
     }
     
     float dEta = fabs(candRef->eta() - pfclu.eta());
     if(dEta < etaStrip) continue;
     
     float dR = deltaR(candRef->eta(), candRef->phi(), pfclu.eta(), pfclu.phi());
     if(dR > drMax_ || dR < dRVeto) continue;
     
     sum += pfclu.pt();
       }
       
       for(unsigned i=0; i<forIsolationHfhad->size(); i++) {
     const reco::PFCluster& pfclu = (*forIsolationHfhad)[i];
     
     if (fabs(candRef->eta()) < 1.479) {
       if (fabs(pfclu.pt()) < energyBarrel_)
         continue;
     } else {
       if (fabs(pfclu.energy()) < energyEndcap_)
         continue;
     }
     
     float dEta = fabs(candRef->eta() - pfclu.eta());
     if(dEta < etaStrip) continue;
     
     float dR = deltaR(candRef->eta(), candRef->phi(), pfclu.eta(), pfclu.phi());
     if(dR > drMax_ || dR < dRVeto) continue;
     
     sum += pfclu.pt();
       }
     }
  
     if (doRhoCorrection_) {
       if (fabs(candRef->eta()) < 1.479) 
     sum = sum - rho*effectiveAreaBarrel_;
       else
     sum = sum - rho*effectiveAreaEndcap_;
     }
 
     recoCandMap.insert(candRef, sum);
   }
   
   std::auto_ptr<T1IsolationMap> mapForEvent(new T1IsolationMap(recoCandMap));
   iEvent.put(mapForEvent);
 }