CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
ElectronIdentifier Class Reference

#include <ElectronIdentifier.h>

Public Member Functions

float dEtaInSeed (const reco::GsfElectronPtr &ele)
 
 ElectronIdentifier (const edm::ParameterSet &c)
 
float isolation (const reco::GsfElectronPtr &ele)
 
bool passID (const reco::GsfElectronPtr &ele, edm::Handle< reco::BeamSpot > beamspot, edm::Handle< reco::ConversionCollection > conversions)
 
void setID (std::string ID)
 
void setRho (double rho)
 

Private Attributes

EffectiveAreas _effectiveAreas
 
std::array< std::array< std::array< double, 2 >, 4 >, 8 > cuts_
 
int ID_
 
double rho_
 

Detailed Description

Definition at line 40 of file ElectronIdentifier.h.

Constructor & Destructor Documentation

◆ ElectronIdentifier()

ElectronIdentifier::ElectronIdentifier ( const edm::ParameterSet &  c)

Definition at line 24 of file ElectronIdentifier.cc.

25  : _effectiveAreas((c.getParameter<edm::FileInPath>("effAreasConfigFile")).fullPath())
26 
27 {
28  rho_ = -1;
29  ID_ = -1;
30  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 0.0115;
31  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 0.011;
32  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.00998;
33  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 0.00998;
34 
35  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 0.037;
36  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 0.0314;
37  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.0298;
38  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 0.0292;
39 
40  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 0.00749;
41  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 0.00477;
42  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.00311;
43  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 0.00308;
44 
45  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 0.00895;
46  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 0.00868;
47  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.00609;
48  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 0.00605;
49 
50  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 0.228;
51  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 0.222;
52  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.103;
53  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 0.0816;
54 
55  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 0.213;
56  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 0.213;
57  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.045;
58  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 0.0394;
59 
60  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 0.356;
61  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 0.298;
62  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.253;
63  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 0.0414;
64 
65  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 0.211;
66  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 0.101;
67  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.0878;
68  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 0.0641;
69 
70  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 0.175;
71  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 0.0994;
72  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.0695;
73  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 0.0588;
74 
75  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 0.159;
76  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 0.107;
77  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.0821;
78  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 0.0571;
79 
80  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 0.299;
81  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 0.241;
82  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.134;
83  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 0.0129;
84 
85  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 0.15;
86  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 0.14;
87  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.13;
88  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 0.0129;
89 
90  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 2;
91  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 1;
92  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 1;
93  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 1;
94 
95  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 3;
96  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 1;
97  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 1;
98  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 1;
99 
100  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::VETO][EleIDEtaBins::BARREL] = 1;
101  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::LOOSE][EleIDEtaBins::BARREL] = 1;
102  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 1;
103  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::TIGHT][EleIDEtaBins::BARREL] = 1;
104 
105  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::VETO][EleIDEtaBins::ENDCAP] = 1;
106  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::LOOSE][EleIDEtaBins::ENDCAP] = 1;
107  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 1;
108  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::TIGHT][EleIDEtaBins::ENDCAP] = 1;
109 }

References BARREL, CONVERSION, cuts_, DETAINSEED, DPHIIN, ENDCAP, HOVERE, ID_, ISO, LOOSE, MEDIUM, MISSINGHITS, ONEOVERE, rho_, SIGMAIETA, TIGHT, and VETO.

Member Function Documentation

◆ dEtaInSeed()

float ElectronIdentifier::dEtaInSeed ( const reco::GsfElectronPtr &  ele)

Definition at line 132 of file ElectronIdentifier.cc.

132  {
133  return ele->superCluster().isNonnull() && ele->superCluster()->seed().isNonnull()
134  ? ele->deltaEtaSuperClusterTrackAtVtx() - ele->superCluster()->eta() + ele->superCluster()->seed()->eta()
135  : std::numeric_limits<float>::max();
136 }

References reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), edm::Ref< C, T, F >::isNonnull(), SiStripPI::max, and reco::GsfElectron::superCluster().

Referenced by passID().

◆ isolation()

float ElectronIdentifier::isolation ( const reco::GsfElectronPtr &  ele)

Definition at line 137 of file ElectronIdentifier.cc.

137  {
138  if (rho_ < 0) {
139  throw;
140  }
141  const reco::GsfElectron::PflowIsolationVariables& pfIso = ele->pfIsolationVariables();
142  const float chad = pfIso.sumChargedHadronPt;
143  const float nhad = pfIso.sumNeutralHadronEt;
144  const float pho = pfIso.sumPhotonEt;
145  const float eA = _effectiveAreas.getEffectiveArea(fabs(ele->superCluster()->eta()));
146  const float iso = chad + std::max(0.0, nhad + pho - rho_ * eA);
147 
148  // Apply the cut and return the result
149  // Scale by pT if the relative isolation is requested but avoid division by 0
150  return iso;
151 }

References _effectiveAreas, EffectiveAreas::getEffectiveArea(), SiStripPI::max, reco::GsfElectron::pfIsolationVariables(), rho_, reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, and reco::GsfElectron::superCluster().

Referenced by passID().

◆ passID()

bool ElectronIdentifier::passID ( const reco::GsfElectronPtr &  ele,
edm::Handle< reco::BeamSpot >  beamspot,
edm::Handle< reco::ConversionCollection >  conversions 
)

Definition at line 153 of file ElectronIdentifier.cc.

155  {
156  if (ID_ == -1)
157  throw;
158  unsigned int region = fabs(ele->superCluster()->eta()) < 1.479 ? EleIDEtaBins::BARREL : EleIDEtaBins::BARREL;
159 
160  if (ele->full5x5_sigmaIetaIeta() > cuts_[EleIDCutNames::SIGMAIETA][ID_][region])
161  return false;
162  if (dEtaInSeed(ele) > cuts_[EleIDCutNames::DETAINSEED][ID_][region])
163  return false;
164  if (std::abs(ele->deltaPhiSuperClusterTrackAtVtx()) > cuts_[EleIDCutNames::DPHIIN][ID_][region])
165  return false;
166  if (ele->hadronicOverEm() > cuts_[EleIDCutNames::HOVERE][ID_][region])
167  return false;
168  if (isolation(ele) / ele->pt() > cuts_[EleIDCutNames::ISO][ID_][region])
169  return false;
170  if (std::abs(1.0 - ele->eSuperClusterOverP()) / ele->ecalEnergy() > cuts_[EleIDCutNames::ONEOVERE][ID_][region])
171  return false;
172  if ((ele->gsfTrack()->hitPattern().numberOfAllHits(reco::HitPattern::MISSING_INNER_HITS)) >
173  cuts_[EleIDCutNames::MISSINGHITS][ID_][region])
174  return false;
175  if (ConversionTools::hasMatchedConversion(*ele, *conversions, beamspot->position()))
176  return false;
177 
178  return true;
179 }

References funct::abs(), BARREL, pwdgSkimBPark_cfi::conversions, cuts_, reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), DETAINSEED, dEtaInSeed(), DPHIIN, reco::GsfElectron::ecalEnergy(), reco::GsfElectron::eSuperClusterOverP(), reco::GsfElectron::full5x5_sigmaIetaIeta(), reco::GsfElectron::gsfTrack(), reco::GsfElectron::hadronicOverEm(), ConversionTools::hasMatchedConversion(), HOVERE, ID_, ISO, isolation(), reco::HitPattern::MISSING_INNER_HITS, MISSINGHITS, ONEOVERE, reco::LeafCandidate::pt(), HLT_FULL_cff::region, SIGMAIETA, and reco::GsfElectron::superCluster().

Referenced by ZCounting::analyzeElectrons().

◆ setID()

void ElectronIdentifier::setID ( std::string  ID)

Definition at line 120 of file ElectronIdentifier.cc.

120  {
121  if (ID == "TIGHT")
122  ID_ = EleIDWorkingPoints::TIGHT;
123  else if (ID == "MEDIUM")
124  ID_ = EleIDWorkingPoints::MEDIUM;
125  else if (ID == "LOOSE")
126  ID_ = EleIDWorkingPoints::LOOSE;
127  else if (ID == "VETO")
128  ID_ = EleIDWorkingPoints::VETO;
129  else
130  throw;
131 }

References ID_, LOOSE, MEDIUM, TIGHT, and VETO.

Referenced by ZCounting::ZCounting().

◆ setRho()

void ElectronIdentifier::setRho ( double  rho)

Definition at line 111 of file ElectronIdentifier.cc.

111  {
112  if (rho >= 0) {
113  rho_ = rho;
114  } else {
115  throw cms::Exception("ValueError") << "Encountered invalid value for energy density rho.\n"
116  << "Value: " << rho << "\n"
117  << "Rho should be a real, positive number.\n";
118  }
119 }

References Exception, rho, and rho_.

Referenced by ZCounting::analyzeElectrons().

Member Data Documentation

◆ _effectiveAreas

EffectiveAreas ElectronIdentifier::_effectiveAreas
private

Definition at line 57 of file ElectronIdentifier.h.

Referenced by isolation().

◆ cuts_

std::array<std::array<std::array<double, 2>, 4>, 8> ElectronIdentifier::cuts_
private

Definition at line 55 of file ElectronIdentifier.h.

Referenced by ElectronIdentifier(), and passID().

◆ ID_

int ElectronIdentifier::ID_
private

Definition at line 54 of file ElectronIdentifier.h.

Referenced by ElectronIdentifier(), passID(), and setID().

◆ rho_

double ElectronIdentifier::rho_
private

Definition at line 53 of file ElectronIdentifier.h.

Referenced by ElectronIdentifier(), isolation(), and setRho().

reco::GsfElectron::gsfTrack
GsfTrackRef gsfTrack() const override
reference to a GsfTrack
Definition: GsfElectron.h:164
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:236
TIGHT
Definition: ElectronIdentifier.h:38
reco::GsfElectron::PflowIsolationVariables
Definition: GsfElectron.h:669
ONEOVERE
Definition: ElectronIdentifier.h:34
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:508
BARREL
Definition: ElectronIdentifier.h:39
ElectronIdentifier::cuts_
std::array< std::array< std::array< double, 2 >, 4 >, 8 > cuts_
Definition: ElectronIdentifier.h:55
edm::FileInPath
Definition: FileInPath.h:61
reco::GsfElectron::PflowIsolationVariables::sumPhotonEt
float sumPhotonEt
sum pt of PF photons // old float photonIso ;
Definition: GsfElectron.h:673
MISSINGHITS
Definition: ElectronIdentifier.h:35
ElectronIdentifier::isolation
float isolation(const reco::GsfElectronPtr &ele)
Definition: ElectronIdentifier.cc:137
ConversionTools::hasMatchedConversion
static bool hasMatchedConversion(const reco::GsfElectron &ele, const reco::ConversionCollection &convCol, const math::XYZPoint &beamspot, bool allowCkfMatch=true, float lxyMin=2.0, float probMin=1e-6, unsigned int nHitsBeforeVtxMax=0)
Definition: ConversionTools.cc:181
DETAINSEED
Definition: ElectronIdentifier.h:30
beamspot
Definition: BeamSpotWrite2Txt.h:8
CONVERSION
Definition: ElectronIdentifier.h:36
ElectronIdentifier::dEtaInSeed
float dEtaInSeed(const reco::GsfElectronPtr &ele)
Definition: ElectronIdentifier.cc:132
DDAxes::rho
align::ID
uint32_t ID
Definition: Definitions.h:24
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:233
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
reco::GsfElectron::full5x5_sigmaIetaIeta
float full5x5_sigmaIetaIeta() const
Definition: GsfElectron.h:479
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:238
ElectronIdentifier::rho_
double rho_
Definition: ElectronIdentifier.h:53
DPHIIN
Definition: ElectronIdentifier.h:31
LOOSE
Definition: ElectronIdentifier.h:38
HLT_FULL_cff.region
region
Definition: HLT_FULL_cff.py:88286
HOVERE
Definition: ElectronIdentifier.h:32
VETO
Definition: ElectronIdentifier.h:38
ISO
Definition: ElectronIdentifier.h:33
reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt
float sumNeutralHadronEt
sum pt of neutral hadrons // old float neutralHadronIso ;
Definition: GsfElectron.h:672
SIGMAIETA
Definition: ElectronIdentifier.h:29
ElectronIdentifier::ID_
int ID_
Definition: ElectronIdentifier.h:54
Exception
Definition: hltDiff.cc:245
reco::HitPattern::MISSING_INNER_HITS
Definition: HitPattern.h:155
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:229
EffectiveAreas::getEffectiveArea
const float getEffectiveArea(float eta) const
Definition: EffectiveAreas.cc:44
ElectronIdentifier::_effectiveAreas
EffectiveAreas _effectiveAreas
Definition: ElectronIdentifier.h:57
reco::GsfElectron::superCluster
SuperClusterRef superCluster() const override
reference to a SuperCluster
Definition: GsfElectron.h:163
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::GsfElectron::pfIsolationVariables
const PflowIsolationVariables & pfIsolationVariables() const
Definition: GsfElectron.h:721
c
auto & c
Definition: CAHitNtupletGeneratorKernelsImpl.h:56
MEDIUM
Definition: ElectronIdentifier.h:38
pwdgSkimBPark_cfi.conversions
conversions
Definition: pwdgSkimBPark_cfi.py:10
reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt
float sumChargedHadronPt
sum-pt of charged Hadron // old float chargedHadronIso ;
Definition: GsfElectron.h:671
ENDCAP
Definition: ElectronIdentifier.h:39
reco::GsfElectron::ecalEnergy
float ecalEnergy() const
Definition: GsfElectron.h:883
Generated for CMSSW Reference Manual by doxygen 1.8.16