EgammaHLTHGCalIDVarProducer Class Reference

Inheritance diagram for EgammaHLTHGCalIDVarProducer:

Classes
class	PCAAssocMap

Public Member Functions
	EgammaHLTHGCalIDVarProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~EgammaHLTHGCalIDVarProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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

Private Attributes
const edm::EDGetTokenT< reco::PFRecHitCollection >	hgcalRecHitToken_
float	hOverECone_
const edm::EDGetTokenT< reco::CaloClusterCollection >	layerClusterToken_
std::vector< PCAAssocMap >	pcaAssocMaps_
float	rCylinder_
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection >	recoEcalCandidateToken_
HGCalShowerShapeHelper	ssCalc_

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

Detailed Description

Definition at line 20 of file EgammaHLTHGCalIDVarProducer.cc.

Constructor & Destructor Documentation

EgammaHLTHGCalIDVarProducer()

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

explicit

Definition at line 60 of file EgammaHLTHGCalIDVarProducer.cc.

    : rCylinder_(config.getParameter<double>("rCylinder")),
      hOverECone_(config.getParameter<double>("hOverECone")),
      recoEcalCandidateToken_(
          consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"))),
      hgcalRecHitToken_(consumes<reco::PFRecHitCollection>(config.getParameter<edm::InputTag>("hgcalRecHits"))),
      layerClusterToken_(consumes<reco::CaloClusterCollection>(config.getParameter<edm::InputTag>("layerClusters"))),
      ssCalc_(consumesCollector()) {
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2xx, "sigma2xx"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2yy, "sigma2yy"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2zz, "sigma2zz"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2xy, "sigma2xy"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2yz, "sigma2yz"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2zx, "sigma2zx"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2uu, "sigma2uu"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2vv, "sigma2vv"));
  pcaAssocMaps_.emplace_back(PCAAssocMap(&HGCalShowerShapeHelper::ShowerWidths::sigma2ww, "sigma2ww"));
 
  produces<reco::RecoEcalCandidateIsolationMap>("rVar");
  produces<reco::RecoEcalCandidateIsolationMap>("hForHOverE");
  for (auto& var : pcaAssocMaps_) {
    produces<reco::RecoEcalCandidateIsolationMap>(var.name());
  }
}

References pcaAssocMaps_, HGCalShowerShapeHelper::ShowerWidths::sigma2uu, HGCalShowerShapeHelper::ShowerWidths::sigma2vv, HGCalShowerShapeHelper::ShowerWidths::sigma2ww, HGCalShowerShapeHelper::ShowerWidths::sigma2xx, HGCalShowerShapeHelper::ShowerWidths::sigma2xy, HGCalShowerShapeHelper::ShowerWidths::sigma2yy, HGCalShowerShapeHelper::ShowerWidths::sigma2yz, HGCalShowerShapeHelper::ShowerWidths::sigma2zx, HGCalShowerShapeHelper::ShowerWidths::sigma2zz, and trigObjTnPSource_cfi::var.

~EgammaHLTHGCalIDVarProducer()

EgammaHLTHGCalIDVarProducer::~EgammaHLTHGCalIDVarProducer ( )

override

Definition at line 85 of file EgammaHLTHGCalIDVarProducer.cc.

{}

Member Function Documentation

fillDescriptions()

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

static

Definition at line 87 of file EgammaHLTHGCalIDVarProducer.cc.

                                                                                             {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("recoEcalCandidateProducer", edm::InputTag("hltL1SeededRecoEcalCandidate"));
  desc.add<edm::InputTag>("hgcalRecHits", edm::InputTag("hgcalRecHits"));
  desc.add<edm::InputTag>("layerClusters", edm::InputTag("layerClusters"));
  desc.add<double>("rCylinder", 2.8);
  desc.add<double>("hOverECone", 0.15);
  descriptions.add(("hltEgammaHLTHGCalIDVarProducer"), desc);
}

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

produce()

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

override

Definition at line 97 of file EgammaHLTHGCalIDVarProducer.cc.

                                                                                       {
  auto recoEcalCandHandle = iEvent.getHandle(recoEcalCandidateToken_);
  const auto& hgcalRecHits = iEvent.get(hgcalRecHitToken_);
  const auto& layerClusters = iEvent.get(layerClusterToken_);
 
  ssCalc_.initPerEvent(iSetup, hgcalRecHits);
 
  auto rVarMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
  auto hForHoverEMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
  for (auto& pcaMap : pcaAssocMaps_) {
    pcaMap.initMap(recoEcalCandHandle);
  }
 
  for (size_t candNr = 0; candNr < recoEcalCandHandle->size(); candNr++) {
    reco::RecoEcalCandidateRef candRef(recoEcalCandHandle, candNr);
    ssCalc_.initPerObject(candRef->superCluster()->hitsAndFractions());
    rVarMap->insert(candRef, ssCalc_.getRvar(rCylinder_, candRef->superCluster()->energy()));
 
    float hForHoverE = HGCalClusterTools::hadEnergyInCone(
        candRef->superCluster()->eta(), candRef->superCluster()->phi(), layerClusters, 0., hOverECone_, 0., 0.);
    hForHoverEMap->insert(candRef, hForHoverE);
    auto pcaWidths = ssCalc_.getPCAWidths(rCylinder_);
    for (auto& pcaMap : pcaAssocMaps_) {
      pcaMap.insert(candRef, pcaWidths);
    }
  }
  iEvent.put(std::move(rVarMap), "rVar");
  iEvent.put(std::move(hForHoverEMap), "hForHOverE");
  for (auto& pcaMap : pcaAssocMaps_) {
    iEvent.put(pcaMap.releaseMap(), pcaMap.name());
  }
}

References HGCalShowerShapeHelper::getPCAWidths(), HGCalShowerShapeHelper::getRvar(), HGCalClusterTools::hadEnergyInCone(), HLTEgPhaseIITestSequence_cff::hgcalRecHits, hgcalRecHitToken_, hOverECone_, iEvent, HGCalShowerShapeHelper::initPerEvent(), HGCalShowerShapeHelper::initPerObject(), HLTEgPhaseIITestSequence_cff::layerClusters, layerClusterToken_, eostools::move(), pcaAssocMaps_, rCylinder_, recoEcalCandidateToken_, and ssCalc_.

Member Data Documentation

hgcalRecHitToken_

const edm::EDGetTokenT<reco::PFRecHitCollection> EgammaHLTHGCalIDVarProducer::hgcalRecHitToken_

private

Definition at line 55 of file EgammaHLTHGCalIDVarProducer.cc.

Referenced by produce().

hOverECone_

float EgammaHLTHGCalIDVarProducer::hOverECone_

private

Definition at line 52 of file EgammaHLTHGCalIDVarProducer.cc.

Referenced by produce().

layerClusterToken_

const edm::EDGetTokenT<reco::CaloClusterCollection> EgammaHLTHGCalIDVarProducer::layerClusterToken_

private

Definition at line 56 of file EgammaHLTHGCalIDVarProducer.cc.

Referenced by produce().

pcaAssocMaps_

std::vector<PCAAssocMap> EgammaHLTHGCalIDVarProducer::pcaAssocMaps_

private

Definition at line 53 of file EgammaHLTHGCalIDVarProducer.cc.

Referenced by EgammaHLTHGCalIDVarProducer(), and produce().

rCylinder_

float EgammaHLTHGCalIDVarProducer::rCylinder_

private

Definition at line 51 of file EgammaHLTHGCalIDVarProducer.cc.

Referenced by produce().

recoEcalCandidateToken_

const edm::EDGetTokenT<reco::RecoEcalCandidateCollection> EgammaHLTHGCalIDVarProducer::recoEcalCandidateToken_

private

Definition at line 54 of file EgammaHLTHGCalIDVarProducer.cc.

Referenced by produce().

ssCalc_

HGCalShowerShapeHelper EgammaHLTHGCalIDVarProducer::ssCalc_

private

Definition at line 57 of file EgammaHLTHGCalIDVarProducer.cc.

Referenced by produce().