HLTEcalPFClusterIsolationProducer< T1 > Class Template Reference

#include <HLTEcalPFClusterIsolationProducer.h>

Inheritance diagram for HLTEcalPFClusterIsolationProducer< T1 >:

Public Member Functions
	HLTEcalPFClusterIsolationProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~HLTEcalPFClusterIsolationProducer ()
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)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::stream::EDProducerBase
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 Member Functions
bool	computedRVeto (T1Ref candRef, reco::PFClusterRef pfclu)

Private Attributes
const std::vector< double >	absEtaLowEdges_
const bool	doRhoCorrection_
const double	drMax_
double	drVeto2_
const double	drVetoBarrel_
const double	drVetoEndcap_
const std::vector< double >	effectiveAreas_
const double	energyBarrel_
const double	energyEndcap_
const double	etaStripBarrel_
const double	etaStripEndcap_
const edm::EDGetTokenT< reco::PFClusterCollection >	pfClusterProducer_
edm::EDGetTokenT< T1Collection >	recoCandidateProducer_
const double	rhoMax_
const edm::EDGetTokenT< double >	rhoProducer_
const double	rhoScale_

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
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 HLTEcalPFClusterIsolationProducer< T1 >

Definition at line 28 of file HLTEcalPFClusterIsolationProducer.h.

Member Typedef Documentation

template<typename T1 >

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

private

Definition at line 30 of file HLTEcalPFClusterIsolationProducer.h.

template<typename T1 >

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

private

Definition at line 32 of file HLTEcalPFClusterIsolationProducer.h.

template<typename T1 >

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

private

Definition at line 31 of file HLTEcalPFClusterIsolationProducer.h.

Constructor & Destructor Documentation

template<typename T1 >

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

explicit

Definition at line 29 of file HLTEcalPFClusterIsolationProducer.cc.

References HLTEcalPFClusterIsolationProducer< T1 >::absEtaLowEdges_, HLTEcalPFClusterIsolationProducer< T1 >::doRhoCorrection_, HLTEcalPFClusterIsolationProducer< T1 >::effectiveAreas_, Exception, edm::ParameterSet::getParameter(), HLTEcalPFClusterIsolationProducer< T1 >::recoCandidateProducer_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                     :
  pfClusterProducer_  (consumes<reco::PFClusterCollection>(config.getParameter<edm::InputTag>("pfClusterProducer"))),
  rhoProducer_        (consumes<double>(config.getParameter<edm::InputTag>("rhoProducer"))),
  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")),
  rhoMax_             (config.getParameter<double>("rhoMax")),
  rhoScale_           (config.getParameter<double>("rhoScale")),
  effectiveAreas_     (config.getParameter<std::vector<double> >("effectiveAreas")),
  absEtaLowEdges_     (config.getParameter<std::vector<double> >("absEtaLowEdges")) {

  if (doRhoCorrection_) {
    if (absEtaLowEdges_.size() != effectiveAreas_.size())
      throw cms::Exception("IncompatibleVects") << "absEtaLowEdges and effectiveAreas should be of the same size. \n";

    if (absEtaLowEdges_.at(0) != 0.0)
      throw cms::Exception("IncompleteCoverage") << "absEtaLowEdges should start from 0. \n";

    for (unsigned int aIt = 0; aIt < absEtaLowEdges_.size() - 1; aIt++) {
      if ( !(absEtaLowEdges_.at( aIt ) < absEtaLowEdges_.at( aIt + 1 )) )
        throw cms::Exception("ImproperBinning") << "absEtaLowEdges entries should be in increasing order. \n";
    }
  }

  std::string recoCandidateProducerName = "recoCandidateProducer";
  if ((typeid(HLTEcalPFClusterIsolationProducer<T1>) == typeid(HLTEcalPFClusterIsolationProducer<reco::RecoEcalCandidate>))) recoCandidateProducerName = "recoEcalCandidateProducer";
    
  recoCandidateProducer_ = consumes<T1Collection>(config.getParameter<edm::InputTag>(recoCandidateProducerName));
  produces <T1IsolationMap>();

}

template<typename T1 >

HLTEcalPFClusterIsolationProducer< T1 >::~HLTEcalPFClusterIsolationProducer

(

)

Definition at line 67 of file HLTEcalPFClusterIsolationProducer.cc.

{}

Member Function Documentation

template<typename T1 >

bool HLTEcalPFClusterIsolationProducer< T1 >::computedRVeto ( T1Ref  candRef, reco::PFClusterRef  pfclu  )

private

template<typename T1 >

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

static

Definition at line 71 of file HLTEcalPFClusterIsolationProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), defaultModuleLabel(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                       {

  std::string recoCandidateProducerName = "recoCandidateProducer";
  if ((typeid(HLTEcalPFClusterIsolationProducer<T1>) == typeid(HLTEcalPFClusterIsolationProducer<reco::RecoEcalCandidate>))) recoCandidateProducerName = "recoEcalCandidateProducer";
  
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>(recoCandidateProducerName, edm::InputTag("hltL1SeededRecoEcalCandidatePF"));
  desc.add<edm::InputTag>("pfClusterProducer", edm::InputTag("hltParticleFlowClusterECAL")); 
  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>("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);
  desc.add<std::vector<double> >("effectiveAreas", {0.29, 0.21}); // 2016 post-ichep sinEle default
  desc.add<std::vector<double> >("absEtaLowEdges", {0.0, 1.479}); // Barrel, Endcap
  descriptions.add(defaultModuleLabel<HLTEcalPFClusterIsolationProducer<T1>>(), desc);
}

template<typename T1 >

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

override

Definition at line 96 of file HLTEcalPFClusterIsolationProducer.cc.

References funct::abs(), HLTEcalPFClusterIsolationProducer< T1 >::absEtaLowEdges_, HLTEcalPFClusterIsolationProducer< T1 >::doRhoCorrection_, HLTEcalPFClusterIsolationProducer< T1 >::drMax_, HLTEcalPFClusterIsolationProducer< T1 >::drVetoBarrel_, HLTEcalPFClusterIsolationProducer< T1 >::drVetoEndcap_, HLTEcalPFClusterIsolationProducer< T1 >::effectiveAreas_, HLTEcalPFClusterIsolationProducer< T1 >::energyBarrel_, HLTEcalPFClusterIsolationProducer< T1 >::energyEndcap_, HLTEcalPFClusterIsolationProducer< T1 >::etaStripBarrel_, HLTEcalPFClusterIsolationProducer< T1 >::etaStripEndcap_, edm::Event::getByToken(), EcalPFClusterIsolation< T1 >::getSum(), edm::AssociationMap< Tag >::insert(), HLTEcalPFClusterIsolationProducer< T1 >::pfClusterProducer_, edm::Handle< T >::product(), edm::Event::put(), HLTEcalPFClusterIsolationProducer< T1 >::recoCandidateProducer_, rho, HLTEcalPFClusterIsolationProducer< T1 >::rhoMax_, HLTEcalPFClusterIsolationProducer< T1 >::rhoProducer_, and HLTEcalPFClusterIsolationProducer< T1 >::rhoScale_.

                                                                                                 {

  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> clusterHandle;

  iEvent.getByToken(recoCandidateProducer_,recoCandHandle);
  iEvent.getByToken(pfClusterProducer_, clusterHandle);

  EcalPFClusterIsolation<T1> isoAlgo(drMax_, drVetoBarrel_, drVetoEndcap_, etaStripBarrel_, etaStripEndcap_, energyBarrel_, energyEndcap_);
  T1IsolationMap recoCandMap(recoCandHandle);
    
  for (unsigned int iReco = 0; iReco < recoCandHandle->size(); iReco++) {
    T1Ref candRef(recoCandHandle, iReco);
    
    float sum = isoAlgo.getSum(candRef, clusterHandle);

    if (doRhoCorrection_) {
      int iEA = -1;
      auto cEta = std::abs(candRef->eta());
      for (int bIt = absEtaLowEdges_.size() - 1; bIt > -1; bIt--) {
        if ( cEta > absEtaLowEdges_.at(bIt) ) {
          iEA = bIt;
          break;
        }
      }

      sum = sum - rho*effectiveAreas_.at(iEA);
    }

    recoCandMap.insert(candRef, sum);
  }
  
  iEvent.put(std::make_unique<T1IsolationMap>(recoCandMap));
}

Member Data Documentation

template<typename T1 >

const std::vector<double> HLTEcalPFClusterIsolationProducer< T1 >::absEtaLowEdges_

private

Definition at line 62 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::HLTEcalPFClusterIsolationProducer(), and HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const bool HLTEcalPFClusterIsolationProducer< T1 >::doRhoCorrection_

private

Definition at line 58 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::HLTEcalPFClusterIsolationProducer(), and HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::drMax_

private

Definition at line 50 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

double HLTEcalPFClusterIsolationProducer< T1 >::drVeto2_

private

Definition at line 49 of file HLTEcalPFClusterIsolationProducer.h.

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::drVetoBarrel_

private

Definition at line 51 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::drVetoEndcap_

private

Definition at line 52 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const std::vector<double> HLTEcalPFClusterIsolationProducer< T1 >::effectiveAreas_

private

Definition at line 61 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::HLTEcalPFClusterIsolationProducer(), and HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::energyBarrel_

private

Definition at line 55 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::energyEndcap_

private

Definition at line 56 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::etaStripBarrel_

private

Definition at line 53 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::etaStripEndcap_

private

Definition at line 54 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const edm::EDGetTokenT<reco::PFClusterCollection> HLTEcalPFClusterIsolationProducer< T1 >::pfClusterProducer_

private

Definition at line 46 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

edm::EDGetTokenT<T1Collection> HLTEcalPFClusterIsolationProducer< T1 >::recoCandidateProducer_

private

Definition at line 45 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::HLTEcalPFClusterIsolationProducer(), and HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::rhoMax_

private

Definition at line 59 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const edm::EDGetTokenT<double> HLTEcalPFClusterIsolationProducer< T1 >::rhoProducer_

private

Definition at line 47 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().

template<typename T1 >

const double HLTEcalPFClusterIsolationProducer< T1 >::rhoScale_

private

Definition at line 60 of file HLTEcalPFClusterIsolationProducer.h.

Referenced by HLTEcalPFClusterIsolationProducer< T1 >::produce().