photonIsolationHIProducer Class Reference

Inheritance diagram for photonIsolationHIProducer:

Public Member Functions
	photonIsolationHIProducer (const edm::ParameterSet &ps)
	~photonIsolationHIProducer () 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 Member Functions
void	produce (edm::Event &evt, const edm::EventSetup &es) override

Private Attributes
edm::EDGetTokenT< reco::BasicClusterCollection >	barrelClusters_
edm::EDGetTokenT< EcalRecHitCollection >	barrelEcalHits_
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<>
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 22 of file photonIsolationHIProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 45 of file photonIsolationHIProducer.cc.

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

photonIsolationHIProducer::~photonIsolationHIProducer ( )

override

Definition at line 59 of file photonIsolationHIProducer.cc.

{}

Member Function Documentation

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

static

Definition at line 137 of file photonIsolationHIProducer.cc.

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

                                                                                           {
  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);
}

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

overrideprivate

Definition at line 61 of file photonIsolationHIProducer.cc.

References HLT_2018_cff::barrelClusters, barrelClusters_, piZeroAnalyzer_cfi::barrelEcalHits, barrelEcalHits_, reco::HIPhotonIsolation::ecalClusterIsoR1(), reco::HIPhotonIsolation::ecalClusterIsoR2(), reco::HIPhotonIsolation::ecalClusterIsoR3(), reco::HIPhotonIsolation::ecalClusterIsoR4(), reco::HIPhotonIsolation::ecalClusterIsoR5(), cosmicPhotonAnalyzer_cfi::eMax, HLT_2018_cff::endcapClusters, endcapClusters_, piZeroAnalyzer_cfi::endcapEcalHits, endcapEcalHits_, trigObjTnPSource_cfi::filler, 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(), TrackInfoProducer_cfi::rechits, SurveyInfoScenario_cff::seed, reco::HIPhotonIsolation::seedTime(), reco::HIPhotonIsolation::swissCrx(), ntuplemaker::time, duplicaterechits_cfi::trackCollection, reco::HIPhotonIsolation::trackIsoR1PtCut20(), reco::HIPhotonIsolation::trackIsoR2PtCut20(), reco::HIPhotonIsolation::trackIsoR3PtCut20(), reco::HIPhotonIsolation::trackIsoR4PtCut20(), reco::HIPhotonIsolation::trackIsoR5PtCut20(), trackQuality_, and tracks_.

                                                                              {
  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(evt, es, trackCollection, trackQuality_);
  EcalClusterLazyTools lazyTool(evt, 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));
}

Member Data Documentation

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

private

Definition at line 38 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<EcalRecHitCollection> photonIsolationHIProducer::barrelEcalHits_

private

Definition at line 33 of file photonIsolationHIProducer.cc.

Referenced by produce().

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

private

Definition at line 39 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<EcalRecHitCollection> photonIsolationHIProducer::endcapEcalHits_

private

Definition at line 34 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<HBHERecHitCollection> photonIsolationHIProducer::hbhe_

private

Definition at line 35 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<HFRecHitCollection> photonIsolationHIProducer::hf_

private

Definition at line 36 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<HORecHitCollection> photonIsolationHIProducer::ho_

private

Definition at line 37 of file photonIsolationHIProducer.cc.

Referenced by produce().

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

private

Definition at line 32 of file photonIsolationHIProducer.cc.

Referenced by produce().

std::string photonIsolationHIProducer::trackQuality_

private

Definition at line 42 of file photonIsolationHIProducer.cc.

Referenced by produce().

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

private

Definition at line 40 of file photonIsolationHIProducer.cc.

Referenced by produce().