photonIsolationHIProducer Class Reference

Inheritance diagram for photonIsolationHIProducer:

Public Member Functions
	photonIsolationHIProducer (const edm::ParameterSet &ps)
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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

Private Member Functions
void	produce (edm::Event &evt, const edm::EventSetup &es) override

Private Attributes
edm::EDGetTokenT< reco::BasicClusterCollection >	barrelClusters_
edm::EDGetTokenT< EcalRecHitCollection >	barrelEcalHits_
const EcalClusterLazyTools::ESGetTokens	ecalClusterToolsESGetTokens_
edm::EDGetTokenT< reco::BasicClusterCollection >	endcapClusters_
edm::EDGetTokenT< EcalRecHitCollection >	endcapEcalHits_
edm::EDGetTokenT< HBHERecHitCollection >	hbhe_
edm::EDGetTokenT< HFRecHitCollection >	hf_
edm::EDGetTokenT< HORecHitCollection >	ho_
edm::EDGetTokenT< reco::PhotonCollection >	photonProducer_
std::string	trackQuality_
edm::EDGetTokenT< reco::TrackCollection >	tracks_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 22 of file photonIsolationHIProducer.cc.

Constructor & Destructor Documentation

photonIsolationHIProducer()

photonIsolationHIProducer::photonIsolationHIProducer ( const edm::ParameterSet &  ps )

explicit

Definition at line 46 of file photonIsolationHIProducer.cc.

    : photonProducer_(consumes(config.getParameter<edm::InputTag>("photonProducer"))),
      barrelEcalHits_(consumes(config.getParameter<edm::InputTag>("ebRecHitCollection"))),
      endcapEcalHits_(consumes(config.getParameter<edm::InputTag>("eeRecHitCollection"))),
      hbhe_(consumes(config.getParameter<edm::InputTag>("hbhe"))),
      hf_(consumes(config.getParameter<edm::InputTag>("hf"))),
      ho_(consumes(config.getParameter<edm::InputTag>("ho"))),
      barrelClusters_(consumes(config.getParameter<edm::InputTag>("basicClusterBarrel"))),
      endcapClusters_(consumes(config.getParameter<edm::InputTag>("basicClusterEndcap"))),
      tracks_(consumes(config.getParameter<edm::InputTag>("trackCollection"))),
      ecalClusterToolsESGetTokens_{consumesCollector()},
      trackQuality_(config.getParameter<std::string>("trackQuality")) {
  produces<reco::HIPhotonIsolationMap>();
}

Member Function Documentation

fillDescriptions()

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

static

Definition at line 137 of file photonIsolationHIProducer.cc.

                                                                                           {
  edm::ParameterSetDescription desc;
 
  desc.add<edm::InputTag>("photonProducer", edm::InputTag("photons"));
  desc.add<edm::InputTag>("ebRecHitCollection", edm::InputTag("ecalRecHit:EcalRecHitsEB"));
  desc.add<edm::InputTag>("eeRecHitCollection", edm::InputTag("ecalRecHit:EcalRecHitsEE"));
  desc.add<edm::InputTag>("hbhe", edm::InputTag("hbhereco"));
  desc.add<edm::InputTag>("hf", edm::InputTag("hfreco"));
  desc.add<edm::InputTag>("ho", edm::InputTag("horeco"));
  desc.add<edm::InputTag>("basicClusterBarrel", edm::InputTag("islandBasicClusters:islandBarrelBasicClusters"));
  desc.add<edm::InputTag>("basicClusterEndcap", edm::InputTag("islandBasicClusters:islandEndcapBasicClusters"));
  desc.add<edm::InputTag>("trackCollection", edm::InputTag("hiGeneralTracks"));
  desc.add<std::string>("trackQuality", "highPurity");
 
  descriptions.add("photonIsolationHIProducer", desc);
}

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

produce()

void photonIsolationHIProducer::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

overrideprivate

Definition at line 61 of file photonIsolationHIProducer.cc.

                                                                              {
  edm::Handle<reco::PhotonCollection> photons;
  evt.getByToken(photonProducer_, photons);
  edm::Handle<EcalRecHitCollection> barrelEcalHits;
  evt.getByToken(barrelEcalHits_, barrelEcalHits);
  edm::Handle<EcalRecHitCollection> endcapEcalHits;
  evt.getByToken(endcapEcalHits_, endcapEcalHits);
  edm::Handle<HBHERecHitCollection> hbhe;
  evt.getByToken(hbhe_, hbhe);
  edm::Handle<HFRecHitCollection> hf;
  evt.getByToken(hf_, hf);
  edm::Handle<HORecHitCollection> ho;
  evt.getByToken(ho_, ho);
  edm::Handle<reco::BasicClusterCollection> barrelClusters;
  evt.getByToken(barrelClusters_, barrelClusters);
  edm::Handle<reco::BasicClusterCollection> endcapClusters;
  evt.getByToken(endcapClusters_, endcapClusters);
  edm::Handle<reco::TrackCollection> trackCollection;
  evt.getByToken(tracks_, trackCollection);
 
  auto outputMap = std::make_unique<reco::HIPhotonIsolationMap>();
  reco::HIPhotonIsolationMap::Filler filler(*outputMap);
  std::vector<reco::HIPhotonIsolation> isoVector;
 
  EcalClusterIsoCalculator CxC(evt, es, barrelClusters, endcapClusters);
  HcalRechitIsoCalculator RxC(evt, es, hbhe, hf, ho);
  TrackIsoCalculator TxC(*trackCollection, trackQuality_);
  EcalClusterLazyTools lazyTool(evt, ecalClusterToolsESGetTokens_.get(es), barrelEcalHits_, endcapEcalHits_);
 
  for (reco::PhotonCollection::const_iterator phoItr = photons->begin(); phoItr != photons->end(); ++phoItr) {
    reco::HIPhotonIsolation iso;
    // HI-style isolation info
    iso.ecalClusterIsoR1(CxC.getBkgSubEcalClusterIso(phoItr->superCluster(), 1, 0));
    iso.ecalClusterIsoR2(CxC.getBkgSubEcalClusterIso(phoItr->superCluster(), 2, 0));
    iso.ecalClusterIsoR3(CxC.getBkgSubEcalClusterIso(phoItr->superCluster(), 3, 0));
    iso.ecalClusterIsoR4(CxC.getBkgSubEcalClusterIso(phoItr->superCluster(), 4, 0));
    iso.ecalClusterIsoR5(CxC.getBkgSubEcalClusterIso(phoItr->superCluster(), 5, 0));
 
    iso.hcalRechitIsoR1(RxC.getBkgSubHcalRechitIso(phoItr->superCluster(), 1, 0));
    iso.hcalRechitIsoR2(RxC.getBkgSubHcalRechitIso(phoItr->superCluster(), 2, 0));
    iso.hcalRechitIsoR3(RxC.getBkgSubHcalRechitIso(phoItr->superCluster(), 3, 0));
    iso.hcalRechitIsoR4(RxC.getBkgSubHcalRechitIso(phoItr->superCluster(), 4, 0));
    iso.hcalRechitIsoR5(RxC.getBkgSubHcalRechitIso(phoItr->superCluster(), 5, 0));
 
    iso.trackIsoR1PtCut20(TxC.getBkgSubTrackIso(*phoItr, 1, 2));
    iso.trackIsoR2PtCut20(TxC.getBkgSubTrackIso(*phoItr, 2, 2));
    iso.trackIsoR3PtCut20(TxC.getBkgSubTrackIso(*phoItr, 3, 2));
    iso.trackIsoR4PtCut20(TxC.getBkgSubTrackIso(*phoItr, 4, 2));
    iso.trackIsoR5PtCut20(TxC.getBkgSubTrackIso(*phoItr, 5, 2));
 
    // ecal spike rejection info (seed timing)
    const reco::CaloClusterPtr seed = phoItr->superCluster()->seed();
    const DetId& id = lazyTool.getMaximum(*seed).first;
    float time = -999.;
    const EcalRecHitCollection& rechits = (phoItr->isEB() ? *barrelEcalHits : *endcapEcalHits);
    EcalRecHitCollection::const_iterator it = rechits.find(id);
    if (it != rechits.end()) {
      time = it->time();
    }
    iso.seedTime(time);
 
    // ecal spike rejectino info (swiss cross)
    float eMax = lazyTool.eMax(*seed);
    float eRight = lazyTool.eRight(*seed);
    float eLeft = lazyTool.eLeft(*seed);
    float eTop = lazyTool.eTop(*seed);
    float eBottom = lazyTool.eBottom(*seed);
    iso.swissCrx(1 - (eRight + eLeft + eTop + eBottom) / eMax);
 
    isoVector.push_back(iso);
  }
  filler.insert(photons, isoVector.begin(), isoVector.end());
  filler.fill();
  evt.put(std::move(outputMap));
}

References HLT_FULL_cff::barrelClusters, barrelClusters_, piZeroAnalyzer_cfi::barrelEcalHits, barrelEcalHits_, reco::HIPhotonIsolation::ecalClusterIsoR1(), reco::HIPhotonIsolation::ecalClusterIsoR2(), reco::HIPhotonIsolation::ecalClusterIsoR3(), reco::HIPhotonIsolation::ecalClusterIsoR4(), reco::HIPhotonIsolation::ecalClusterIsoR5(), ecalClusterToolsESGetTokens_, cosmicPhotonAnalyzer_cfi::eMax, HLT_FULL_cff::endcapClusters, endcapClusters_, piZeroAnalyzer_cfi::endcapEcalHits, endcapEcalHits_, trigObjTnPSource_cfi::filler, EcalClusterLazyToolsBase::ESGetTokens::get(), EcalClusterIsoCalculator::getBkgSubEcalClusterIso(), HcalRechitIsoCalculator::getBkgSubHcalRechitIso(), TrackIsoCalculator::getBkgSubTrackIso(), edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, hbhe_, reco::HIPhotonIsolation::hcalRechitIsoR1(), reco::HIPhotonIsolation::hcalRechitIsoR2(), reco::HIPhotonIsolation::hcalRechitIsoR3(), reco::HIPhotonIsolation::hcalRechitIsoR4(), reco::HIPhotonIsolation::hcalRechitIsoR5(), photonIsolationHIProducer_cfi::hf, hf_, photonIsolationHIProducer_cfi::ho, ho_, eostools::move(), photonProducer_, BPHMonitor_cfi::photons, edm::Event::put(), HI_PhotonSkim_cff::rechits, fileCollector::seed, reco::HIPhotonIsolation::seedTime(), reco::HIPhotonIsolation::swissCrx(), protons_cff::time, duplicaterechits_cfi::trackCollection, reco::HIPhotonIsolation::trackIsoR1PtCut20(), reco::HIPhotonIsolation::trackIsoR2PtCut20(), reco::HIPhotonIsolation::trackIsoR3PtCut20(), reco::HIPhotonIsolation::trackIsoR4PtCut20(), reco::HIPhotonIsolation::trackIsoR5PtCut20(), trackQuality_, and tracks_.

Member Data Documentation

barrelClusters_

edm::EDGetTokenT<reco::BasicClusterCollection> photonIsolationHIProducer::barrelClusters_

private

Definition at line 37 of file photonIsolationHIProducer.cc.

Referenced by produce().

barrelEcalHits_

edm::EDGetTokenT<EcalRecHitCollection> photonIsolationHIProducer::barrelEcalHits_

private

Definition at line 32 of file photonIsolationHIProducer.cc.

Referenced by produce().

ecalClusterToolsESGetTokens_

const EcalClusterLazyTools::ESGetTokens photonIsolationHIProducer::ecalClusterToolsESGetTokens_

private

Definition at line 41 of file photonIsolationHIProducer.cc.

Referenced by produce().

endcapClusters_

edm::EDGetTokenT<reco::BasicClusterCollection> photonIsolationHIProducer::endcapClusters_

private

Definition at line 38 of file photonIsolationHIProducer.cc.

Referenced by produce().

endcapEcalHits_

edm::EDGetTokenT<EcalRecHitCollection> photonIsolationHIProducer::endcapEcalHits_

private

Definition at line 33 of file photonIsolationHIProducer.cc.

Referenced by produce().

hbhe_

edm::EDGetTokenT<HBHERecHitCollection> photonIsolationHIProducer::hbhe_

private

Definition at line 34 of file photonIsolationHIProducer.cc.

Referenced by produce().

hf_

edm::EDGetTokenT<HFRecHitCollection> photonIsolationHIProducer::hf_

private

Definition at line 35 of file photonIsolationHIProducer.cc.

Referenced by produce().

ho_

edm::EDGetTokenT<HORecHitCollection> photonIsolationHIProducer::ho_

private

Definition at line 36 of file photonIsolationHIProducer.cc.

Referenced by produce().

photonProducer_

edm::EDGetTokenT<reco::PhotonCollection> photonIsolationHIProducer::photonProducer_

private

Definition at line 31 of file photonIsolationHIProducer.cc.

Referenced by produce().

trackQuality_

std::string photonIsolationHIProducer::trackQuality_

private

Definition at line 43 of file photonIsolationHIProducer.cc.

Referenced by produce().

tracks_

edm::EDGetTokenT<reco::TrackCollection> photonIsolationHIProducer::tracks_

private

Definition at line 39 of file photonIsolationHIProducer.cc.

Referenced by produce().