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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 49 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 68 of file photonIsolationHIProducer.cc.

{}

Member Function Documentation

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

static

Definition at line 150 of file photonIsolationHIProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, 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 71 of file photonIsolationHIProducer.cc.

References barrelClusters_, eleEcalExtractorBlocks_cff::barrelEcalHits, barrelEcalHits_, reco::HIPhotonIsolation::ecalClusterIsoR1(), reco::HIPhotonIsolation::ecalClusterIsoR2(), reco::HIPhotonIsolation::ecalClusterIsoR3(), reco::HIPhotonIsolation::ecalClusterIsoR4(), reco::HIPhotonIsolation::ecalClusterIsoR5(), photonPostprocessing_cfi::eMax, endcapClusters_, eleEcalExtractorBlocks_cff::endcapEcalHits, endcapEcalHits_, objects.autophobj::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_, muons_cff::photons, edm::Event::put(), TrackInfoProducer_cfi::rechits, SurveyInfoScenario_cff::seed, reco::HIPhotonIsolation::seedTime(), reco::HIPhotonIsolation::swissCrx(), ntuplemaker::time, findElectronsInSiStrips_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 41 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<EcalRecHitCollection> photonIsolationHIProducer::barrelEcalHits_

private

Definition at line 36 of file photonIsolationHIProducer.cc.

Referenced by produce().

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

private

Definition at line 42 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<EcalRecHitCollection> photonIsolationHIProducer::endcapEcalHits_

private

Definition at line 37 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<HBHERecHitCollection> photonIsolationHIProducer::hbhe_

private

Definition at line 38 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<HFRecHitCollection> photonIsolationHIProducer::hf_

private

Definition at line 39 of file photonIsolationHIProducer.cc.

Referenced by produce().

edm::EDGetTokenT<HORecHitCollection> photonIsolationHIProducer::ho_

private

Definition at line 40 of file photonIsolationHIProducer.cc.

Referenced by produce().

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

private

Definition at line 35 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 43 of file photonIsolationHIProducer.cc.

Referenced by produce().