PhotonIsolationCalculator Class Reference

#include <PhotonIsolationCalculator.h>

Public Member Functions
void	calculate (const reco::Photon *, const edm::Event &, const edm::EventSetup &es, reco::Photon::FiducialFlags &phofid, reco::Photon::IsolationVariables &phoisolR03, reco::Photon::IsolationVariables &phoisolR04) const
	PhotonIsolationCalculator ()
void	setup (const edm::ParameterSet &conf, std::vector< int > const &flagsEB_, std::vector< int > const &flagsEE_, std::vector< int > const &severitiesEB_, std::vector< int > const &severitiesEE_, edm::ConsumesCollector &&iC)
	~PhotonIsolationCalculator ()

Private Member Functions
double	calculateEcalRecHitIso (const reco::Photon *photon, const edm::Event &iEvent, const edm::EventSetup &iSetup, double RCone, double RConeInner, double etaSlice, double eMin, double etMin, bool vetoClusteredHits, bool useNumCrystals) const
template<bool isoBC>
double	calculateHcalRecHitIso (const reco::Photon *photon, const CaloGeometry &geometry, const HcalTopology &hcalTopology, const HcalChannelQuality &hcalChStatus, const HcalSeverityLevelComputer &hcalSevLvlComputer, const CaloTowerConstituentsMap &towerMap, const HBHERecHitCollection &hbheRecHits, double RCone, double RConeInner, int depth) const
void	calculateTrackIso (const reco::Photon *photon, const edm::Event &e, double &trkCone, int &ntrkCone, double pTThresh=0, double RCone=.4, double RinnerCone=.1, double etaSlice=0.015, double lip=0.2, double d0=0.1) const

Static Private Member Functions
static void	classify (const reco::Photon *photon, bool &isEBPho, bool &isEEPho, bool &isEBEtaGap, bool &isEBPhiGap, bool &isEERingGap, bool &isEEDeeGap, bool &isEBEEGap)

Private Attributes
edm::EDGetToken	barrelecalCollection_
edm::EDGetToken	beamSpotProducerTag_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
std::array< double, 5 >	ecalIsoBarrelRadiusA_
std::array< double, 5 >	ecalIsoBarrelRadiusB_
std::array< double, 5 >	ecalIsoEndcapRadiusA_
std::array< double, 5 >	ecalIsoEndcapRadiusB_
edm::ESGetToken < EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	ecalSevLvlToken_
edm::EDGetToken	endcapecalCollection_
std::vector< int >	flagsEB_
std::vector< int >	flagsEE_
edm::EDGetTokenT < HBHERecHitCollection >	hbheRecHitsTag_
edm::ESGetToken < HcalChannelQuality, HcalChannelQualityRcd >	hcalChannelQualityToken_
EgammaHcalIsolation::arrayHB	hcalIsoEThresHB_
EgammaHcalIsolation::arrayHE	hcalIsoEThresHE_
std::array< double, 7 >	hcalIsoInnerRadAEB_
std::array< double, 7 >	hcalIsoInnerRadAEE_
std::array< double, 7 >	hcalIsoInnerRadBEB_
std::array< double, 7 >	hcalIsoInnerRadBEE_
std::array< double, 7 >	hcalIsoOuterRadAEB_
std::array< double, 7 >	hcalIsoOuterRadAEE_
std::array< double, 7 >	hcalIsoOuterRadBEB_
std::array< double, 7 >	hcalIsoOuterRadBEE_
edm::ESGetToken < HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd >	hcalSevLvlComputerToken_
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	hcalTopologyToken_
int	maxHcalSeverity_
std::vector< double >	moduleEtaBoundary_
double	modulePhiBoundary_
std::vector< int >	severityExclEB_
std::vector< int >	severityExclEE_
edm::ESGetToken < CaloTowerConstituentsMap, CaloGeometryRecord >	towerMapToken_
edm::EDGetToken	trackInputTag_
std::array< double, 6 >	trkIsoBarrelRadiusA_
std::array< double, 6 >	trkIsoBarrelRadiusB_
std::array< double, 6 >	trkIsoEndcapRadiusA_
std::array< double, 6 >	trkIsoEndcapRadiusB_
bool	useNumCrystals_
bool	vetoClusteredEcalHits_

Detailed Description

Determine and Set quality information on Photon Objects

Author

A. Askew, N. Marinelli, M.B. Anderson

Definition at line 31 of file PhotonIsolationCalculator.h.

Constructor & Destructor Documentation

PhotonIsolationCalculator::PhotonIsolationCalculator ( )

inline

Definition at line 33 of file PhotonIsolationCalculator.h.

{}

PhotonIsolationCalculator::~PhotonIsolationCalculator ( )

inline

Definition at line 35 of file PhotonIsolationCalculator.h.

{}

Member Function Documentation

void PhotonIsolationCalculator::calculate	(	const reco::Photon *	pho,
		const edm::Event &	e,
		const edm::EventSetup &	es,
		reco::Photon::FiducialFlags &	phofid,
		reco::Photon::IsolationVariables &	phoisolR03,
		reco::Photon::IsolationVariables &	phoisolR04
	)		const

Definition at line 145 of file PhotonIsolationCalculator.cc.

References calculateEcalRecHitIso(), calculateTrackIso(), caloGeometryToken_, classify(), HLT_FULL_cff::depth, EcalBarrel, ecalIsoBarrelRadiusA_, ecalIsoBarrelRadiusB_, ecalIsoEndcapRadiusA_, ecalIsoEndcapRadiusB_, reco::Photon::IsolationVariables::ecalRecHitSumEt, fbc, edm::Event::get(), edm::EventSetup::getData(), ootPhotons_cff::hbheRecHits, hbheRecHitsTag_, hcalChannelQualityToken_, hcalIsoInnerRadAEB_, hcalIsoInnerRadAEE_, hcalIsoInnerRadBEB_, hcalIsoInnerRadBEE_, hcalIsoOuterRadAEB_, hcalIsoOuterRadAEE_, hcalIsoOuterRadBEB_, hcalIsoOuterRadBEE_, reco::Photon::IsolationVariables::hcalRecHitSumEt, reco::Photon::IsolationVariables::hcalRecHitSumEtBc, hcalSevLvlComputerToken_, hcalTopologyToken_, reco::CaloCluster::hitsAndFractions(), gpuClustering::id, SurfaceOrientation::inner, reco::Photon::FiducialFlags::isEB, reco::Photon::FiducialFlags::isEBEEGap, reco::Photon::FiducialFlags::isEBEtaGap, reco::Photon::FiducialFlags::isEBPhiGap, reco::Photon::FiducialFlags::isEE, reco::Photon::FiducialFlags::isEEDeeGap, reco::Photon::FiducialFlags::isEERingGap, reco::Photon::IsolationVariables::nTrkHollowCone, reco::Photon::IsolationVariables::nTrkSolidCone, SurfaceOrientation::outer, reco::Photon::IsolationVariables::pre7DepthHcal, DetId::subdetId(), reco::Photon::superCluster(), towerMapToken_, trkIsoBarrelRadiusA_, trkIsoBarrelRadiusB_, trkIsoEndcapRadiusA_, trkIsoEndcapRadiusB_, reco::Photon::IsolationVariables::trkSumPtHollowCone, reco::Photon::IsolationVariables::trkSumPtSolidCone, useNumCrystals_, and vetoClusteredEcalHits_.

Referenced by PhotonProducer::fillPhotonCollection().

                                                                                         {
  //Get fiducial flags. This does not really belong here
  bool isEBPho = false;
  bool isEEPho = false;
  bool isEBEtaGap = false;
  bool isEBPhiGap = false;
  bool isEERingGap = false;
  bool isEEDeeGap = false;
  bool isEBEEGap = false;
  classify(pho, isEBPho, isEEPho, isEBEtaGap, isEBPhiGap, isEERingGap, isEEDeeGap, isEBEEGap);
  phofid.isEB = isEBPho;
  phofid.isEE = isEEPho;
  phofid.isEBEtaGap = isEBEtaGap;
  phofid.isEBPhiGap = isEBPhiGap;
  phofid.isEERingGap = isEERingGap;
  phofid.isEEDeeGap = isEEDeeGap;
  phofid.isEBEEGap = isEBEEGap;

  auto const& caloGeometry = es.getData(caloGeometryToken_);
  auto const& hcalTopology = &es.getData(hcalTopologyToken_);
  auto const& hcalChannelQuality = &es.getData(hcalChannelQualityToken_);
  auto const& hcalSevLvlComputer = &es.getData(hcalSevLvlComputerToken_);
  auto const& towerMap = es.getData(towerMapToken_);

  // Calculate isolation variables. cone sizes and thresholds
  // are set for Barrel and endcap separately

  reco::SuperClusterRef scRef = pho->superCluster();
  const reco::BasicCluster& seedCluster = *(scRef->seed());
  DetId seedXtalId = seedCluster.hitsAndFractions()[0].first;
  int detector = seedXtalId.subdetId();

  //Isolation parameters variables
  double photonEcalRecHitConeInnerRadiusA_;
  double photonEcalRecHitConeOuterRadiusA_;
  double photonEcalRecHitEtaSliceA_;
  double photonEcalRecHitThreshEA_;
  double photonEcalRecHitThreshEtA_;
  double trackConeOuterRadiusA_;
  double trackConeInnerRadiusA_;
  double isolationtrackThresholdA_;
  double isolationtrackEtaSliceA_;
  double trackLipRadiusA_;
  double trackD0RadiusA_;
  double photonEcalRecHitConeInnerRadiusB_;
  double photonEcalRecHitConeOuterRadiusB_;
  double photonEcalRecHitEtaSliceB_;
  double photonEcalRecHitThreshEB_;
  double photonEcalRecHitThreshEtB_;
  double trackConeOuterRadiusB_;
  double trackConeInnerRadiusB_;
  double isolationtrackThresholdB_;
  double isolationtrackEtaSliceB_;
  double trackLipRadiusB_;
  double trackD0RadiusB_;

  if (detector == EcalBarrel) {
    trackConeOuterRadiusA_ = trkIsoBarrelRadiusA_[0];
    trackConeInnerRadiusA_ = trkIsoBarrelRadiusA_[1];
    isolationtrackThresholdA_ = trkIsoBarrelRadiusA_[2];
    trackLipRadiusA_ = trkIsoBarrelRadiusA_[3];
    trackD0RadiusA_ = trkIsoBarrelRadiusA_[4];
    isolationtrackEtaSliceA_ = trkIsoBarrelRadiusA_[5];

    photonEcalRecHitConeInnerRadiusA_ = ecalIsoBarrelRadiusA_[0];
    photonEcalRecHitConeOuterRadiusA_ = ecalIsoBarrelRadiusA_[1];
    photonEcalRecHitEtaSliceA_ = ecalIsoBarrelRadiusA_[2];
    photonEcalRecHitThreshEA_ = ecalIsoBarrelRadiusA_[3];
    photonEcalRecHitThreshEtA_ = ecalIsoBarrelRadiusA_[4];

    trackConeOuterRadiusB_ = trkIsoBarrelRadiusB_[0];
    trackConeInnerRadiusB_ = trkIsoBarrelRadiusB_[1];
    isolationtrackThresholdB_ = trkIsoBarrelRadiusB_[2];
    trackLipRadiusB_ = trkIsoBarrelRadiusB_[3];
    trackD0RadiusB_ = trkIsoBarrelRadiusB_[4];
    isolationtrackEtaSliceB_ = trkIsoBarrelRadiusB_[5];

    photonEcalRecHitConeInnerRadiusB_ = ecalIsoBarrelRadiusB_[0];
    photonEcalRecHitConeOuterRadiusB_ = ecalIsoBarrelRadiusB_[1];
    photonEcalRecHitEtaSliceB_ = ecalIsoBarrelRadiusB_[2];
    photonEcalRecHitThreshEB_ = ecalIsoBarrelRadiusB_[3];
    photonEcalRecHitThreshEtB_ = ecalIsoBarrelRadiusB_[4];
  } else {
    // detector==EcalEndcap

    trackConeOuterRadiusA_ = trkIsoEndcapRadiusA_[0];
    trackConeInnerRadiusA_ = trkIsoEndcapRadiusA_[1];
    isolationtrackThresholdA_ = trkIsoEndcapRadiusA_[2];
    trackLipRadiusA_ = trkIsoEndcapRadiusA_[3];
    trackD0RadiusA_ = trkIsoEndcapRadiusA_[4];
    isolationtrackEtaSliceA_ = trkIsoEndcapRadiusA_[5];

    photonEcalRecHitConeInnerRadiusA_ = ecalIsoEndcapRadiusA_[0];
    photonEcalRecHitConeOuterRadiusA_ = ecalIsoEndcapRadiusA_[1];
    photonEcalRecHitEtaSliceA_ = ecalIsoEndcapRadiusA_[2];
    photonEcalRecHitThreshEA_ = ecalIsoEndcapRadiusA_[3];
    photonEcalRecHitThreshEtA_ = ecalIsoEndcapRadiusA_[4];

    trackConeOuterRadiusB_ = trkIsoEndcapRadiusB_[0];
    trackConeInnerRadiusB_ = trkIsoEndcapRadiusB_[1];
    isolationtrackThresholdB_ = trkIsoEndcapRadiusB_[2];
    trackLipRadiusB_ = trkIsoEndcapRadiusB_[3];
    trackD0RadiusB_ = trkIsoEndcapRadiusB_[4];
    isolationtrackEtaSliceB_ = trkIsoEndcapRadiusB_[5];

    photonEcalRecHitConeInnerRadiusB_ = ecalIsoEndcapRadiusB_[0];
    photonEcalRecHitConeOuterRadiusB_ = ecalIsoEndcapRadiusB_[1];
    photonEcalRecHitEtaSliceB_ = ecalIsoEndcapRadiusB_[2];
    photonEcalRecHitThreshEB_ = ecalIsoEndcapRadiusB_[3];
    photonEcalRecHitThreshEtB_ = ecalIsoEndcapRadiusB_[4];
  }

  //Calculate hollow cone track isolation, CONE A
  int ntrkA = 0;
  double trkisoA = 0;
  calculateTrackIso(pho,
                    e,
                    trkisoA,
                    ntrkA,
                    isolationtrackThresholdA_,
                    trackConeOuterRadiusA_,
                    trackConeInnerRadiusA_,
                    isolationtrackEtaSliceA_,
                    trackLipRadiusA_,
                    trackD0RadiusA_);

  //Calculate solid cone track isolation, CONE A
  int sntrkA = 0;
  double strkisoA = 0;
  calculateTrackIso(pho,
                    e,
                    strkisoA,
                    sntrkA,
                    isolationtrackThresholdA_,
                    trackConeOuterRadiusA_,
                    0.,
                    isolationtrackEtaSliceA_,
                    trackLipRadiusA_,
                    trackD0RadiusA_);

  phoisolR1.nTrkHollowCone = ntrkA;
  phoisolR1.trkSumPtHollowCone = trkisoA;
  phoisolR1.nTrkSolidCone = sntrkA;
  phoisolR1.trkSumPtSolidCone = strkisoA;

  //Calculate hollow cone track isolation, CONE B
  int ntrkB = 0;
  double trkisoB = 0;
  calculateTrackIso(pho,
                    e,
                    trkisoB,
                    ntrkB,
                    isolationtrackThresholdB_,
                    trackConeOuterRadiusB_,
                    trackConeInnerRadiusB_,
                    isolationtrackEtaSliceB_,
                    trackLipRadiusB_,
                    trackD0RadiusB_);

  //Calculate solid cone track isolation, CONE B
  int sntrkB = 0;
  double strkisoB = 0;
  calculateTrackIso(pho,
                    e,
                    strkisoB,
                    sntrkB,
                    isolationtrackThresholdB_,
                    trackConeOuterRadiusB_,
                    0.,
                    isolationtrackEtaSliceB_,
                    trackLipRadiusB_,
                    trackD0RadiusB_);

  phoisolR2.nTrkHollowCone = ntrkB;
  phoisolR2.trkSumPtHollowCone = trkisoB;
  phoisolR2.nTrkSolidCone = sntrkB;
  phoisolR2.trkSumPtSolidCone = strkisoB;

  //   std::cout << "Output from solid cone track isolation: ";
  //   std::cout << " Sum pT: " << strkiso << " ntrk: " << sntrk << std::endl;

  double EcalRecHitIsoA = calculateEcalRecHitIso(pho,
                                                 e,
                                                 es,
                                                 photonEcalRecHitConeOuterRadiusA_,
                                                 photonEcalRecHitConeInnerRadiusA_,
                                                 photonEcalRecHitEtaSliceA_,
                                                 photonEcalRecHitThreshEA_,
                                                 photonEcalRecHitThreshEtA_,
                                                 vetoClusteredEcalHits_,
                                                 useNumCrystals_);
  phoisolR1.ecalRecHitSumEt = EcalRecHitIsoA;

  double EcalRecHitIsoB = calculateEcalRecHitIso(pho,
                                                 e,
                                                 es,
                                                 photonEcalRecHitConeOuterRadiusB_,
                                                 photonEcalRecHitConeInnerRadiusB_,
                                                 photonEcalRecHitEtaSliceB_,
                                                 photonEcalRecHitThreshEB_,
                                                 photonEcalRecHitThreshEtB_,
                                                 vetoClusteredEcalHits_,
                                                 useNumCrystals_);
  phoisolR2.ecalRecHitSumEt = EcalRecHitIsoB;

  if (not hbheRecHitsTag_.isUninitialized()) {
    auto const& hbheRecHits = e.get(hbheRecHitsTag_);

    auto fcone = [this,
                  pho,
                  &caloGeometry,
                  &hcalTopo = *hcalTopology,
                  &hcalQual = *hcalChannelQuality,
                  &hcalSev = *hcalSevLvlComputer,
                  &towerMap,
                  &hbheRecHits](double outer, double inner, int depth) {
      return calculateHcalRecHitIso<false>(
          pho, caloGeometry, hcalTopo, hcalQual, hcalSev, towerMap, hbheRecHits, outer, inner, depth);
    };

    auto fbc = [this,
                pho,
                &caloGeometry,
                &hcalTopo = *hcalTopology,
                &hcalQual = *hcalChannelQuality,
                &hcalSev = *hcalSevLvlComputer,
                &towerMap,
                &hbheRecHits](double outer, int depth) {
      return calculateHcalRecHitIso<true>(
          pho, caloGeometry, hcalTopo, hcalQual, hcalSev, towerMap, hbheRecHits, outer, 0., depth);
    };

    for (size_t id = 0; id < phoisolR1.hcalRecHitSumEt.size(); ++id) {
      const auto& outerA = detector == EcalBarrel ? hcalIsoOuterRadAEB_[id] : hcalIsoOuterRadAEE_[id];
      const auto& outerB = detector == EcalBarrel ? hcalIsoOuterRadBEB_[id] : hcalIsoOuterRadBEE_[id];
      const auto& innerA = detector == EcalBarrel ? hcalIsoInnerRadAEB_[id] : hcalIsoInnerRadAEE_[id];
      const auto& innerB = detector == EcalBarrel ? hcalIsoInnerRadBEB_[id] : hcalIsoInnerRadBEE_[id];

      phoisolR1.hcalRecHitSumEt[id] = fcone(outerA, innerA, id + 1);
      phoisolR2.hcalRecHitSumEt[id] = fcone(outerB, innerB, id + 1);

      phoisolR1.hcalRecHitSumEtBc[id] = fbc(outerA, id + 1);
      phoisolR2.hcalRecHitSumEtBc[id] = fbc(outerB, id + 1);
    }
  }

  phoisolR1.pre7DepthHcal = false;
  phoisolR2.pre7DepthHcal = false;
}

double PhotonIsolationCalculator::calculateEcalRecHitIso ( const reco::Photon *  photon, const edm::Event &  iEvent, const edm::EventSetup &  iSetup, double  RCone, double  RConeInner, double  etaSlice, double  eMin, double  etMin, bool  vetoClusteredHits, bool  useNumCrystals  ) const

private

Definition at line 484 of file PhotonIsolationCalculator.cc.

References barrelecalCollection_, caloGeometryToken_, EgammaRecHitIsolation::doFlagChecks(), EgammaRecHitIsolation::doSeverityChecks(), DetId::Ecal, ecalSevLvlToken_, endcapecalCollection_, flagsEB_, flagsEE_, edm::Event::getByToken(), edm::EventSetup::getData(), EgammaRecHitIsolation::getEtSum(), edm::EventSetup::getHandle(), edm::Handle< T >::product(), EgammaRecHitIsolation::setUseNumCrystals(), EgammaRecHitIsolation::setVetoClustered(), severityExclEB_, and severityExclEE_.

Referenced by calculate().

                                                                                 {
  edm::Handle<EcalRecHitCollection> ecalhitsCollEB;
  edm::Handle<EcalRecHitCollection> ecalhitsCollEE;
  iEvent.getByToken(endcapecalCollection_, ecalhitsCollEE);

  iEvent.getByToken(barrelecalCollection_, ecalhitsCollEB);

  const EcalRecHitCollection* rechitsCollectionEE_ = ecalhitsCollEE.product();
  const EcalRecHitCollection* rechitsCollectionEB_ = ecalhitsCollEB.product();

  const EcalSeverityLevelAlgo* sevLevel = &iSetup.getData(ecalSevLvlToken_);

  edm::ESHandle<CaloGeometry> geoHandle = iSetup.getHandle(caloGeometryToken_);

  EgammaRecHitIsolation phoIsoEB(
      RCone, RConeInner, etaSlice, etMin, eMin, geoHandle, *rechitsCollectionEB_, sevLevel, DetId::Ecal);

  phoIsoEB.setVetoClustered(vetoClusteredHits);
  phoIsoEB.setUseNumCrystals(useNumXtals);
  phoIsoEB.doSeverityChecks(ecalhitsCollEB.product(), severityExclEB_);
  phoIsoEB.doFlagChecks(flagsEB_);
  double ecalIsolEB = phoIsoEB.getEtSum(photon);

  EgammaRecHitIsolation phoIsoEE(
      RCone, RConeInner, etaSlice, etMin, eMin, geoHandle, *rechitsCollectionEE_, sevLevel, DetId::Ecal);

  phoIsoEE.setVetoClustered(vetoClusteredHits);
  phoIsoEE.setUseNumCrystals(useNumXtals);
  phoIsoEE.doSeverityChecks(ecalhitsCollEE.product(), severityExclEE_);
  phoIsoEE.doFlagChecks(flagsEE_);

  double ecalIsolEE = phoIsoEE.getEtSum(photon);
  //  delete phoIso;
  double ecalIsol = ecalIsolEB + ecalIsolEE;

  return ecalIsol;
}

template<bool isoBC>

double PhotonIsolationCalculator::calculateHcalRecHitIso ( const reco::Photon *  photon, const CaloGeometry &  geometry, const HcalTopology &  hcalTopology, const HcalChannelQuality &  hcalChStatus, const HcalSeverityLevelComputer &  hcalSevLvlComputer, const CaloTowerConstituentsMap &  towerMap, const HBHERecHitCollection &  hbheRecHits, double  RCone, double  RConeInner, int  depth  ) const

private

Definition at line 532 of file PhotonIsolationCalculator.cc.

References hcalIsoEThresHB_, hcalIsoEThresHE_, EgammaHcalIsolation::isBehindClusterSeed, maxHcalSeverity_, and EgammaHcalIsolation::withinConeAroundCluster.

                                                                          {
  const EgammaHcalIsolation::arrayHB e04{{0., 0., 0., 0.}};
  const EgammaHcalIsolation::arrayHE e07{{0., 0., 0., 0., 0., 0., 0.}};

  if constexpr (isoBC) {
    auto hcaliso = EgammaHcalIsolation(EgammaHcalIsolation::InclusionRule::withinConeAroundCluster,
                                       RCone,
                                       EgammaHcalIsolation::InclusionRule::isBehindClusterSeed,
                                       RConeInner,
                                       hcalIsoEThresHB_,
                                       e04,
                                       maxHcalSeverity_,
                                       hcalIsoEThresHE_,
                                       e07,
                                       maxHcalSeverity_,
                                       hbheRecHits,
                                       geometry,
                                       hcalTopology,
                                       hcalChStatus,
                                       hcalSevLvlComputer,
                                       towerMap);

    return hcaliso.getHcalEtSumBc(photon, depth);
  } else {
    auto hcaliso = EgammaHcalIsolation(EgammaHcalIsolation::InclusionRule::withinConeAroundCluster,
                                       RCone,
                                       EgammaHcalIsolation::InclusionRule::withinConeAroundCluster,
                                       RConeInner,
                                       hcalIsoEThresHB_,
                                       e04,
                                       maxHcalSeverity_,
                                       hcalIsoEThresHE_,
                                       e07,
                                       maxHcalSeverity_,
                                       hbheRecHits,
                                       geometry,
                                       hcalTopology,
                                       hcalChStatus,
                                       hcalSevLvlComputer,
                                       towerMap);

    return hcaliso.getHcalEtSum(photon, depth);
  }
}

void PhotonIsolationCalculator::calculateTrackIso ( const reco::Photon *  photon, const edm::Event &  e, double &  trkCone, int &  ntrkCone, double  pTThresh = 0, double  RCone = .4, double  RinnerCone = .1, double  etaSlice = 0.015, double  lip = 0.2, double  d0 = 0.1  ) const

private

Definition at line 445 of file PhotonIsolationCalculator.cc.

References beamSpotProducerTag_, edm::Event::getByToken(), PhotonTkIsolation::getIso(), edm::HandleBase::isValid(), edm::Handle< T >::product(), HLT_FULL_cff::trackCollection, trackInputTag_, tracks, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by calculate().

                                                                   {
  ntrkCone = 0;
  trkCone = 0;
  //get the tracks
  edm::Handle<reco::TrackCollection> tracks;
  e.getByToken(trackInputTag_, tracks);
  if (!tracks.isValid()) {
    return;
  }
  const reco::TrackCollection* trackCollection = tracks.product();
  //Photon Eta and Phi.  Hope these are correct.
  reco::BeamSpot vertexBeamSpot;
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  e.getByToken(beamSpotProducerTag_, recoBeamSpotHandle);
  vertexBeamSpot = *recoBeamSpotHandle;

  PhotonTkIsolation phoIso(RCone,
                           RinnerCone,
                           etaSlice,
                           pTThresh,
                           lip,
                           d0,
                           trackCollection,
                           math::XYZPoint(vertexBeamSpot.x0(), vertexBeamSpot.y0(), vertexBeamSpot.z0()));

  std::pair<int, double> res = phoIso.getIso(photon);
  ntrkCone = res.first;
  trkCone = res.second;
}

void PhotonIsolationCalculator::classify ( const reco::Photon *  photon, bool &  isEBPho, bool &  isEEPho, bool &  isEBEtaGap, bool &  isEBPhiGap, bool &  isEERingGap, bool &  isEEDeeGap, bool &  isEBEEGap  )

staticprivate

Definition at line 400 of file PhotonIsolationCalculator.cc.

References EcalBarrel, EcalEndcap, PVValHelper::eta, reco::CaloCluster::hitsAndFractions(), EEDetId::isNextToDBoundary(), EBDetId::isNextToEtaBoundary(), EBDetId::isNextToPhiBoundary(), EEDetId::isNextToRingBoundary(), DetId::subdetId(), and reco::Photon::superCluster().

Referenced by calculate().

                                                          {
  const reco::CaloCluster& seedCluster = *(photon->superCluster()->seed());
  // following line is temporary until the D. Evans or F. Ferri submit the new tag for ClusterAlgos
  // DEfinitive will be
  // DetId seedXtalId = scRef->seed()->seed();
  DetId seedXtalId = seedCluster.hitsAndFractions()[0].first;
  int detector = seedXtalId.subdetId();

  //Set fiducial flags for this photon.
  double eta = photon->superCluster()->position().eta();
  double feta = fabs(eta);

  if (detector == EcalBarrel) {
    isEBPho = true;
    if (EBDetId::isNextToEtaBoundary(EBDetId(seedXtalId))) {
      if (fabs(feta - 1.479) < .1)
        isEBEEGap = true;
      else
        isEBEtaGap = true;
    }

    if (EBDetId::isNextToPhiBoundary(EBDetId(seedXtalId)))
      isEBPhiGap = true;

  } else if (detector == EcalEndcap) {
    isEEPho = true;
    if (EEDetId::isNextToRingBoundary(EEDetId(seedXtalId))) {
      if (fabs(feta - 1.479) < .1)
        isEBEEGap = true;
      else
        isEERingGap = true;
    }

    if (EEDetId::isNextToDBoundary(EEDetId(seedXtalId)))
      isEEDeeGap = true;
  }
}

void PhotonIsolationCalculator::setup	(	const edm::ParameterSet &	conf,
		std::vector< int > const &	flagsEB_,
		std::vector< int > const &	flagsEE_,
		std::vector< int > const &	severitiesEB_,
		std::vector< int > const &	severitiesEE_,
		edm::ConsumesCollector &&	iC
	)

Isolation parameters for barrel and for two different cone sizes

Isolation parameters for Endcap and for two different cone sizes

Definition at line 29 of file PhotonIsolationCalculator.cc.

References barrelecalCollection_, beamSpotProducerTag_, caloGeometryToken_, ecalIsoBarrelRadiusA_, ecalIsoBarrelRadiusB_, ecalIsoEndcapRadiusA_, ecalIsoEndcapRadiusB_, ecalSevLvlToken_, endcapecalCollection_, flagsEB_, flagsEE_, edm::ParameterSet::getParameter(), hbheRecHitsTag_, hcalChannelQualityToken_, hcalIsoEThresHB_, hcalIsoEThresHE_, hcalIsoInnerRadAEB_, hcalIsoInnerRadAEE_, hcalIsoInnerRadBEB_, hcalIsoInnerRadBEE_, hcalIsoOuterRadAEB_, hcalIsoOuterRadAEE_, hcalIsoOuterRadBEB_, hcalIsoOuterRadBEE_, hcalSevLvlComputerToken_, hcalTopologyToken_, mps_fire::i, maxHcalSeverity_, moduleEtaBoundary_, modulePhiBoundary_, severityExclEB_, severityExclEE_, towerMapToken_, trackInputTag_, trkIsoBarrelRadiusA_, trkIsoBarrelRadiusB_, trkIsoEndcapRadiusA_, trkIsoEndcapRadiusB_, useNumCrystals_, and vetoClusteredEcalHits_.

Referenced by PhotonProducer::PhotonProducer().

                                                                 {
  trackInputTag_ = iC.consumes<reco::TrackCollection>(conf.getParameter<edm::InputTag>("trackProducer"));
  beamSpotProducerTag_ = iC.consumes<reco::BeamSpot>(conf.getParameter<edm::InputTag>("beamSpotProducer"));
  barrelecalCollection_ =
      iC.consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("barrelEcalRecHitCollection"));
  endcapecalCollection_ =
      iC.consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("endcapEcalRecHitCollection"));

  auto hbhetag = conf.getParameter<edm::InputTag>("HBHERecHitCollection");
  if (not hbhetag.label().empty())
    hbheRecHitsTag_ = iC.consumes<HBHERecHitCollection>(hbhetag);

  caloGeometryToken_ = decltype(caloGeometryToken_){iC.esConsumes()};
  hcalTopologyToken_ = decltype(hcalTopologyToken_){iC.esConsumes()};
  hcalChannelQualityToken_ = decltype(hcalChannelQualityToken_){iC.esConsumes(edm::ESInputTag("", "withTopo"))};
  hcalSevLvlComputerToken_ = decltype(hcalSevLvlComputerToken_){iC.esConsumes()};
  towerMapToken_ = decltype(towerMapToken_){iC.esConsumes()};
  ecalSevLvlToken_ = iC.esConsumes();

  //  gsfRecoInputTag_ = conf.getParameter<edm::InputTag>("GsfRecoCollection");
  modulePhiBoundary_ = conf.getParameter<double>("modulePhiBoundary");
  moduleEtaBoundary_ = conf.getParameter<std::vector<double>>("moduleEtaBoundary");
  //
  vetoClusteredEcalHits_ = conf.getParameter<bool>("vetoClustered");
  useNumCrystals_ = conf.getParameter<bool>("useNumCrystals");

  int i = 0;
  trkIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("TrackConeOuterRadiusA_Barrel"));
  trkIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("TrackConeInnerRadiusA_Barrel"));
  trkIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("isolationtrackThresholdA_Barrel"));
  trkIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("longImpactParameterA_Barrel"));
  trkIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("transImpactParameterA_Barrel"));
  trkIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("isolationtrackEtaSliceA_Barrel"));

  i = 0;
  ecalIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitInnerRadiusA_Barrel"));
  ecalIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitOuterRadiusA_Barrel"));
  ecalIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitEtaSliceA_Barrel"));
  ecalIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitThreshEA_Barrel"));
  ecalIsoBarrelRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitThreshEtA_Barrel"));

  hcalIsoInnerRadAEB_ = conf.getParameter<std::array<double, 7>>("HcalRecHitInnerRadiusA_Barrel");
  hcalIsoOuterRadAEB_ = conf.getParameter<std::array<double, 7>>("HcalRecHitOuterRadiusA_Barrel");

  i = 0;
  trkIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("TrackConeOuterRadiusB_Barrel"));
  trkIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("TrackConeInnerRadiusB_Barrel"));
  trkIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("isolationtrackThresholdB_Barrel"));
  trkIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("longImpactParameterB_Barrel"));
  trkIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("transImpactParameterB_Barrel"));
  trkIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("isolationtrackEtaSliceB_Barrel"));

  i = 0;
  ecalIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitInnerRadiusB_Barrel"));
  ecalIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitOuterRadiusB_Barrel"));
  ecalIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitEtaSliceB_Barrel"));
  ecalIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitThreshEB_Barrel"));
  ecalIsoBarrelRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitThreshEtB_Barrel"));

  hcalIsoInnerRadBEB_ = conf.getParameter<std::array<double, 7>>("HcalRecHitInnerRadiusB_Barrel");
  hcalIsoOuterRadBEB_ = conf.getParameter<std::array<double, 7>>("HcalRecHitOuterRadiusB_Barrel");

  i = 0;
  trkIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("TrackConeOuterRadiusA_Endcap"));
  trkIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("TrackConeInnerRadiusA_Endcap"));
  trkIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("isolationtrackThresholdA_Endcap"));
  trkIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("longImpactParameterA_Endcap"));
  trkIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("transImpactParameterA_Endcap"));
  trkIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("isolationtrackEtaSliceA_Endcap"));

  i = 0;
  ecalIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitInnerRadiusA_Endcap"));
  ecalIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitOuterRadiusA_Endcap"));
  ecalIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitEtaSliceA_Endcap"));
  ecalIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitThreshEA_Endcap"));
  ecalIsoEndcapRadiusA_[i++] = (conf.getParameter<double>("EcalRecHitThreshEtA_Endcap"));

  hcalIsoInnerRadAEE_ = conf.getParameter<std::array<double, 7>>("HcalRecHitInnerRadiusA_Endcap");
  hcalIsoOuterRadAEE_ = conf.getParameter<std::array<double, 7>>("HcalRecHitOuterRadiusA_Endcap");

  i = 0;
  trkIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("TrackConeOuterRadiusB_Endcap"));
  trkIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("TrackConeInnerRadiusB_Endcap"));
  trkIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("isolationtrackThresholdB_Endcap"));
  trkIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("longImpactParameterB_Endcap"));
  trkIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("transImpactParameterB_Endcap"));
  trkIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("isolationtrackEtaSliceB_Endcap"));

  i = 0;
  ecalIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitInnerRadiusB_Endcap"));
  ecalIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitOuterRadiusB_Endcap"));
  ecalIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitEtaSliceB_Endcap"));
  ecalIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitThreshEB_Endcap"));
  ecalIsoEndcapRadiusB_[i++] = (conf.getParameter<double>("EcalRecHitThreshEtB_Endcap"));

  hcalIsoInnerRadBEE_ = conf.getParameter<std::array<double, 7>>("HcalRecHitInnerRadiusB_Endcap");
  hcalIsoOuterRadBEE_ = conf.getParameter<std::array<double, 7>>("HcalRecHitOuterRadiusB_Endcap");

  //Pick up the variables for the spike removal
  flagsEB_ = flagsEB;
  flagsEE_ = flagsEE;
  severityExclEB_ = severitiesEB;
  severityExclEE_ = severitiesEE;

  hcalIsoEThresHB_ = conf.getParameter<EgammaHcalIsolation::arrayHB>("recHitEThresholdHB");
  hcalIsoEThresHE_ = conf.getParameter<EgammaHcalIsolation::arrayHE>("recHitEThresholdHE");
  maxHcalSeverity_ = conf.getParameter<int>("maxHcalRecHitSeverity");
}

Member Data Documentation

edm::EDGetToken PhotonIsolationCalculator::barrelecalCollection_

private

Definition at line 96 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

edm::EDGetToken PhotonIsolationCalculator::beamSpotProducerTag_

private

Definition at line 108 of file PhotonIsolationCalculator.h.

Referenced by calculateTrackIso(), and setup().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> PhotonIsolationCalculator::caloGeometryToken_

private

Definition at line 100 of file PhotonIsolationCalculator.h.

Referenced by calculate(), calculateEcalRecHitIso(), and setup().

std::array<double, 5> PhotonIsolationCalculator::ecalIsoBarrelRadiusA_

private

Definition at line 115 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 5> PhotonIsolationCalculator::ecalIsoBarrelRadiusB_

private

Definition at line 117 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 5> PhotonIsolationCalculator::ecalIsoEndcapRadiusA_

private

Definition at line 120 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 5> PhotonIsolationCalculator::ecalIsoEndcapRadiusB_

private

Definition at line 122 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> PhotonIsolationCalculator::ecalSevLvlToken_

private

Definition at line 105 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

edm::EDGetToken PhotonIsolationCalculator::endcapecalCollection_

private

Definition at line 97 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

std::vector<int> PhotonIsolationCalculator::flagsEB_

private

Definition at line 140 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

std::vector<int> PhotonIsolationCalculator::flagsEE_

private

Definition at line 141 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

edm::EDGetTokenT<HBHERecHitCollection> PhotonIsolationCalculator::hbheRecHitsTag_

private

Definition at line 98 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

edm::ESGetToken<HcalChannelQuality, HcalChannelQualityRcd> PhotonIsolationCalculator::hcalChannelQualityToken_

private

Definition at line 102 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

EgammaHcalIsolation::arrayHB PhotonIsolationCalculator::hcalIsoEThresHB_

private

Definition at line 136 of file PhotonIsolationCalculator.h.

Referenced by calculateHcalRecHitIso(), and setup().

EgammaHcalIsolation::arrayHE PhotonIsolationCalculator::hcalIsoEThresHE_

private

Definition at line 137 of file PhotonIsolationCalculator.h.

Referenced by calculateHcalRecHitIso(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoInnerRadAEB_

private

Definition at line 124 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoInnerRadAEE_

private

Definition at line 130 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoInnerRadBEB_

private

Definition at line 127 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoInnerRadBEE_

private

Definition at line 133 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoOuterRadAEB_

private

Definition at line 125 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoOuterRadAEE_

private

Definition at line 131 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoOuterRadBEB_

private

Definition at line 128 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 7> PhotonIsolationCalculator::hcalIsoOuterRadBEE_

private

Definition at line 134 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

edm::ESGetToken<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd> PhotonIsolationCalculator::hcalSevLvlComputerToken_

private

Definition at line 103 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> PhotonIsolationCalculator::hcalTopologyToken_

private

Definition at line 101 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

int PhotonIsolationCalculator::maxHcalSeverity_

private

Definition at line 138 of file PhotonIsolationCalculator.h.

Referenced by calculateHcalRecHitIso(), and setup().

std::vector<double> PhotonIsolationCalculator::moduleEtaBoundary_

private

Definition at line 110 of file PhotonIsolationCalculator.h.

Referenced by setup().

double PhotonIsolationCalculator::modulePhiBoundary_

private

Definition at line 109 of file PhotonIsolationCalculator.h.

Referenced by setup().

std::vector<int> PhotonIsolationCalculator::severityExclEB_

private

Definition at line 142 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

std::vector<int> PhotonIsolationCalculator::severityExclEE_

private

Definition at line 143 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

edm::ESGetToken<CaloTowerConstituentsMap, CaloGeometryRecord> PhotonIsolationCalculator::towerMapToken_

private

Definition at line 104 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

edm::EDGetToken PhotonIsolationCalculator::trackInputTag_

private

Definition at line 107 of file PhotonIsolationCalculator.h.

Referenced by calculateTrackIso(), and setup().

std::array<double, 6> PhotonIsolationCalculator::trkIsoBarrelRadiusA_

private

Definition at line 114 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 6> PhotonIsolationCalculator::trkIsoBarrelRadiusB_

private

Definition at line 116 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 6> PhotonIsolationCalculator::trkIsoEndcapRadiusA_

private

Definition at line 119 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::array<double, 6> PhotonIsolationCalculator::trkIsoEndcapRadiusB_

private

Definition at line 121 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

bool PhotonIsolationCalculator::useNumCrystals_

private

Definition at line 112 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

bool PhotonIsolationCalculator::vetoClusteredEcalHits_

private

Definition at line 111 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().