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_
const edm::ESGetToken < CaloGeometry, CaloGeometryRecord >	geometryToken_
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 23 of file photonIsolationHIProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 48 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()},
      geometryToken_{esConsumes()},
      trackQuality_(config.getParameter<std::string>("trackQuality")) {
  produces<reco::HIPhotonIsolationMap>();
}

Member Function Documentation

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

static

Definition at line 145 of file photonIsolationHIProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_FULL_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 64 of file photonIsolationHIProducer.cc.

References HLT_FULL_cff::barrelClusters, barrelClusters_, ootPhotons_cff::barrelEcalHits, barrelEcalHits_, reco::HIPhotonIsolation::ecalClusterIsoR1(), reco::HIPhotonIsolation::ecalClusterIsoR2(), reco::HIPhotonIsolation::ecalClusterIsoR3(), reco::HIPhotonIsolation::ecalClusterIsoR4(), reco::HIPhotonIsolation::ecalClusterIsoR5(), ecalClusterToolsESGetTokens_, HLT_FULL_cff::endcapClusters, endcapClusters_, ootPhotons_cff::endcapEcalHits, endcapEcalHits_, relativeConstraints::geometry, geometryToken_, EcalClusterLazyToolsBase::ESGetTokens::get(), EcalClusterIsoCalculator::getBkgSubEcalClusterIso(), HcalRechitIsoCalculator::getBkgSubHcalRechitIso(), TrackIsoCalculator::getBkgSubTrackIso(), edm::Event::getByToken(), edm::EventSetup::getHandle(), hbhe_, reco::HIPhotonIsolation::hcalRechitIsoR1(), reco::HIPhotonIsolation::hcalRechitIsoR2(), reco::HIPhotonIsolation::hcalRechitIsoR3(), reco::HIPhotonIsolation::hcalRechitIsoR4(), reco::HIPhotonIsolation::hcalRechitIsoR5(), hf_, hcalSimParameters_cfi::ho, ho_, edm::ESHandleBase::isValid(), eostools::move(), photonProducer_, interactiveExample::photons, edm::ESHandle< class >::product(), edm::Event::put(), HI_PhotonSkim_cff::rechits, fileCollector::seed, reco::HIPhotonIsolation::seedTime(), reco::HIPhotonIsolation::swissCrx(), HLT_FULL_cff::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(barrelClusters, endcapClusters);
  edm::ESHandle<CaloGeometry> geometryHandle = es.getHandle(geometryToken_);
  const CaloGeometry* geometry = nullptr;
  if (geometryHandle.isValid()) {
    geometry = geometryHandle.product();
  }
  HcalRechitIsoCalculator RxC(geometry, 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));
}

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().

const EcalClusterLazyTools::ESGetTokens photonIsolationHIProducer::ecalClusterToolsESGetTokens_

private

Definition at line 42 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().

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> photonIsolationHIProducer::geometryToken_

private

Definition at line 43 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 45 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().