CMS 3D CMS Logo
ElectronIdentifier.cc
Go to the documentation of this file.
1 #include "DQMOffline/Lumi/interface/ElectronIdentifier.h"
2 #include "FWCore/Framework/interface/Event.h"
3 #include "FWCore/Common/interface/TriggerNames.h"
4 #include "FWCore/Utilities/interface/RegexMatch.h"
5 #include "DataFormats/Common/interface/Handle.h"
6 #include "DataFormats/Common/interface/TriggerResults.h"
7 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
8 #include "DataFormats/VertexReco/interface/Vertex.h"
9 #include "DataFormats/MuonReco/interface/Muon.h"
10 #include "DataFormats/MuonReco/interface/MuonSelectors.h"
11 
12 #include "DataFormats/EgammaCandidates/interface/Electron.h"
13 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
14 #include "DataFormats/EgammaReco/interface/SuperCluster.h"
15 #include "RecoEgamma/EgammaTools/interface/ConversionTools.h"
16 
17 #include "DQMOffline/Lumi/interface/TTrigger.h"
18 #include "DQMOffline/Lumi/interface/TriggerTools.h"
19 
20 #include <boost/foreach.hpp>
21 #include <TLorentzVector.h>
22 #include <TMath.h>
23 #include <algorithm>
24 
25 ElectronIdentifier::ElectronIdentifier (const edm::ParameterSet& c) :
26  _effectiveAreas( (c.getParameter<edm::FileInPath>("effAreasConfigFile")).fullPath())
27 
28 {
29  rho_ = -1;
30  ID_ = -1;
31  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 0.0115;
32  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 0.011;
33  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.00998;
34  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 0.00998;
35 
36  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 0.037;
37  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 0.0314;
38  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.0298;
39  cuts_[EleIDCutNames::SIGMAIETA][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 0.0292;
40 
41  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 0.00749;
42  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 0.00477;
43  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.00311;
44  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 0.00308;
45 
46  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 0.00895;
47  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 0.00868;
48  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.00609;
49  cuts_[EleIDCutNames::DETAINSEED][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 0.00605;
50 
51  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 0.228;
52  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 0.222;
53  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.103;
54  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 0.0816;
55 
56  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 0.213;
57  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 0.213;
58  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.045;
59  cuts_[EleIDCutNames::DPHIIN][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 0.0394;
60 
61  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 0.356;
62  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 0.298;
63  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.253;
64  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 0.0414;
65 
66  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 0.211;
67  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 0.101;
68  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.0878;
69  cuts_[EleIDCutNames::HOVERE][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 0.0641;
70 
71  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 0.175;
72  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 0.0994;
73  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.0695;
74  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 0.0588;
75 
76  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 0.159;
77  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 0.107;
78  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.0821;
79  cuts_[EleIDCutNames::ISO][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 0.0571;
80 
81  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 0.299;
82  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 0.241;
83  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 0.134;
84  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 0.0129;
85 
86  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 0.15;
87  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 0.14;
88  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 0.13;
89  cuts_[EleIDCutNames::ONEOVERE][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 0.0129;
90 
91  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 2;
92  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 1;
93  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 1;
94  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 1;
95 
96  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 3;
97  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 1;
98  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 1;
99  cuts_[EleIDCutNames::MISSINGHITS][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 1;
100 
101  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::VETO] [EleIDEtaBins::BARREL] = 1;
102  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::BARREL] = 1;
103  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::BARREL] = 1;
104  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::BARREL] = 1;
105 
106  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::VETO] [EleIDEtaBins::ENDCAP] = 1;
107  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::LOOSE] [EleIDEtaBins::ENDCAP] = 1;
108  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::MEDIUM][EleIDEtaBins::ENDCAP] = 1;
109  cuts_[EleIDCutNames::CONVERSION][EleIDWorkingPoints::TIGHT] [EleIDEtaBins::ENDCAP] = 1;
110 }
111 
112 void ElectronIdentifier::setRho(double rho) {
113  if(rho >= 0) {
114  rho_ = rho;
115  } else {
116  throw cms::Exception("ValueError")
117  << "Encountered invalid value for energy density rho.\n"
118  << "Value: " << rho << "\n"
119  << "Rho should be a real, positive number.\n";
120  }
121 }
122 void ElectronIdentifier::setID(std::string ID) {
123  if(ID=="TIGHT") ID_ = EleIDWorkingPoints::TIGHT;
124  else if(ID=="MEDIUM") ID_ = EleIDWorkingPoints::MEDIUM;
125  else if(ID=="LOOSE") ID_ = EleIDWorkingPoints::LOOSE;
126  else if(ID=="VETO") ID_ = EleIDWorkingPoints::VETO;
127  else throw;
128 }
129 float ElectronIdentifier::dEtaInSeed(const reco::GsfElectronPtr& ele){
130  return ele->superCluster().isNonnull() && ele->superCluster()->seed().isNonnull() ?
131  ele->deltaEtaSuperClusterTrackAtVtx() - ele->superCluster()->eta() + ele->superCluster()->seed()->eta() : std::numeric_limits<float>::max();
132  }
133 float ElectronIdentifier::isolation(const reco::GsfElectronPtr& ele) {
134  if(rho_ < 0 ) {
135  throw;
136  }
137  const reco::GsfElectron::PflowIsolationVariables& pfIso = ele->pfIsolationVariables();
138  const float chad = pfIso.sumChargedHadronPt;
139  const float nhad = pfIso.sumNeutralHadronEt;
140  const float pho = pfIso.sumPhotonEt;
141  const float eA = _effectiveAreas.getEffectiveArea( fabs(ele->superCluster()->eta()) );
142  const float iso = chad + std::max(0.0, nhad + pho - rho_*eA);
143 
144  // Apply the cut and return the result
145  // Scale by pT if the relative isolation is requested but avoid division by 0
146  return iso;
147 }
148 
149 
150 bool ElectronIdentifier::passID(const reco::GsfElectronPtr& ele,edm::Handle<reco::BeamSpot> beamspot,edm::Handle<reco::ConversionCollection> conversions) {
151  if(ID_ == -1) throw;
152  unsigned int region = fabs(ele->superCluster()->eta()) < 1.479 ? EleIDEtaBins::BARREL : EleIDEtaBins::BARREL;
153 
154  if( ele->full5x5_sigmaIetaIeta() > cuts_[EleIDCutNames::SIGMAIETA][ID_][region]) return false;
155  if( dEtaInSeed(ele) > cuts_[EleIDCutNames::DETAINSEED][ID_][region]) return false;
156  if( std::abs(ele->deltaPhiSuperClusterTrackAtVtx()) > cuts_[EleIDCutNames::DPHIIN][ID_][region]) return false;
157  if( ele->hadronicOverEm() > cuts_[EleIDCutNames::HOVERE][ID_][region]) return false;
158  if( isolation(ele)/ele->pt() > cuts_[EleIDCutNames::ISO][ID_][region]) return false;
159  if( std::abs(1.0 - ele->eSuperClusterOverP())/ele->ecalEnergy() > cuts_[EleIDCutNames::ONEOVERE][ID_][region]) return false;
160  if( (ele->gsfTrack()->hitPattern().numberOfAllHits(reco::HitPattern::MISSING_INNER_HITS)) > cuts_[EleIDCutNames::MISSINGHITS][ID_][region]) return false;
161  if( ConversionTools::hasMatchedConversion(*ele,conversions,beamspot->position())) return false;
162 
163  return true;
164 }
165 
166 
167 
168 
reco::GsfElectron::pfIsolationVariables
const PflowIsolationVariables & pfIsolationVariables() const
Definition: GsfElectron.h:673
reco::GsfElectron::gsfTrack
GsfTrackRef gsfTrack() const override
reference to a GsfTrack
Definition: GsfElectron.h:185
ElectronIdentifier::dEtaInSeed
float dEtaInSeed(const reco::GsfElectronPtr &ele)
Definition: ElectronIdentifier.cc:129
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:253
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:245
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
align::ID
uint32_t ID
Definition: Definitions.h:26
ElectronIdentifier::cuts_
std::array< std::array< std::array< double, 2 >, 4 >, 8 > cuts_
Definition: ElectronIdentifier.h:59
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
EffectiveAreas::getEffectiveArea
const float getEffectiveArea(float eta) const
Definition: EffectiveAreas.cc:50
conversions_cfi.conversions
conversions
Definition: conversions_cfi.py:8
reco::GsfElectron::full5x5_sigmaIetaIeta
float full5x5_sigmaIetaIeta() const
Definition: GsfElectron.h:455
ElectronIdentifier::_effectiveAreas
EffectiveAreas _effectiveAreas
Definition: ElectronIdentifier.h:61
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:249
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:487
reco::GsfElectron::PflowIsolationVariables::sumPhotonEt
float sumPhotonEt
sum pt of PF photons // old float photonIso ;
Definition: GsfElectron.h:627
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:252
ElectronIdentifier::passID
bool passID(const reco::GsfElectronPtr &ele, edm::Handle< reco::BeamSpot > beamspot, edm::Handle< reco::ConversionCollection > conversions)
Definition: ElectronIdentifier.cc:150
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt
float sumNeutralHadronEt
sum pt of neutral hadrons // old float neutralHadronIso ;
Definition: GsfElectron.h:626
ConversionTools::hasMatchedConversion
static bool hasMatchedConversion(const reco::GsfElectron &ele, const edm::Handle< 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:149
ElectronIdentifier::ElectronIdentifier
ElectronIdentifier(const edm::ParameterSet &c)
Definition: ElectronIdentifier.cc:25
ElectronIdentifier::setRho
void setRho(double rho)
Definition: ElectronIdentifier.cc:112
reco::GsfElectron::ecalEnergy
float ecalEnergy() const
Definition: GsfElectron.h:844
edm
HLT enums.
Definition: AlignableModifier.h:17
ElectronIdentifier::setID
void setID(std::string ID)
Definition: ElectronIdentifier.cc:122
reco::GsfElectron::superCluster
SuperClusterRef superCluster() const override
reference to a SuperCluster
Definition: GsfElectron.h:184
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
ElectronIdentifier::isolation
float isolation(const reco::GsfElectronPtr &ele)
Definition: ElectronIdentifier.cc:133
reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt
float sumChargedHadronPt
sum-pt of charged Hadron // old float chargedHadronIso ;
Definition: GsfElectron.h:625
Generated for CMSSW Reference Manual by doxygen 1.8.11