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

#include <EGammaMvaEleEstimatorCSA14.h>

Public Types

enum  MVAType { kTrig = 0, kNonTrig = 1, kNonTrigPhys14 = 2 }
 

Public Member Functions

void bindVariables ()
 
 EGammaMvaEleEstimatorCSA14 ()
 
UInt_t GetMVABin (double eta, double pt) const
 
void initialize (std::string methodName, std::string weightsfile, EGammaMvaEleEstimatorCSA14::MVAType type)
 
void initialize (std::string methodName, EGammaMvaEleEstimatorCSA14::MVAType type, Bool_t useBinnedVersion, std::vector< std::string > weightsfiles)
 
Bool_t isInitialized () const
 
Double_t mvaValue (const reco::GsfElectron &ele, const reco::Vertex &vertex, const TransientTrackBuilder &transientTrackBuilder, noZS::EcalClusterLazyTools myEcalCluster, bool printDebug=kFALSE)
 
Double_t mvaValue (const pat::Electron &ele, bool printDebug)
 
 ~EGammaMvaEleEstimatorCSA14 ()
 

Private Attributes

Bool_t fisInitialized
 
std::string fMethodname
 
MVAType fMVAType
 
Float_t fMVAVar_abseta
 
Float_t fMVAVar_ChargedIso_DR0p0To0p1
 
Float_t fMVAVar_ChargedIso_DR0p1To0p2
 
Float_t fMVAVar_ChargedIso_DR0p2To0p3
 
Float_t fMVAVar_ChargedIso_DR0p3To0p4
 
Float_t fMVAVar_ChargedIso_DR0p4To0p5
 
Float_t fMVAVar_d0
 
Float_t fMVAVar_deta
 
Float_t fMVAVar_detacalo
 
Float_t fMVAVar_dphi
 
Float_t fMVAVar_eleEoPout
 
Float_t fMVAVar_EoP
 
Float_t fMVAVar_EoPout
 
Float_t fMVAVar_eta
 
Float_t fMVAVar_etawidth
 
Float_t fMVAVar_fbrem
 
Float_t fMVAVar_GammaIso_DR0p0To0p1
 
Float_t fMVAVar_GammaIso_DR0p1To0p2
 
Float_t fMVAVar_GammaIso_DR0p2To0p3
 
Float_t fMVAVar_GammaIso_DR0p3To0p4
 
Float_t fMVAVar_GammaIso_DR0p4To0p5
 
Float_t fMVAVar_gsfchi2
 
Float_t fMVAVar_HoE
 
Float_t fMVAVar_IoEmIoP
 
Float_t fMVAVar_ip3d
 
Float_t fMVAVar_ip3dSig
 
Float_t fMVAVar_isBarrel
 
Float_t fMVAVar_isEndcap
 
Float_t fMVAVar_kfchi2
 
Float_t fMVAVar_kfhits
 
Float_t fMVAVar_kfhitsall
 
Float_t fMVAVar_NeutralHadronIso_DR0p0To0p1
 
Float_t fMVAVar_NeutralHadronIso_DR0p1To0p2
 
Float_t fMVAVar_NeutralHadronIso_DR0p2To0p3
 
Float_t fMVAVar_NeutralHadronIso_DR0p3To0p4
 
Float_t fMVAVar_NeutralHadronIso_DR0p4To0p5
 
Float_t fMVAVar_OneMinusE1x5E5x5
 
Float_t fMVAVar_phiwidth
 
Float_t fMVAVar_PreShowerOverRaw
 
Float_t fMVAVar_pt
 
Float_t fMVAVar_R9
 
Float_t fMVAVar_rho
 
Float_t fMVAVar_SCeta
 
Float_t fMVAVar_see
 
Float_t fMVAVar_spp
 
UInt_t fNMVABins
 
std::vector< TMVA::MethodBase * > fTMVAMethod
 
std::vector< TMVA::Reader * > fTMVAReader
 
Bool_t fUseBinnedVersion
 

Detailed Description

Definition at line 41 of file EGammaMvaEleEstimatorCSA14.h.

Member Enumeration Documentation

enum EGammaMvaEleEstimatorCSA14::MVAType
Enumerator
kTrig 
kNonTrig 
kNonTrigPhys14 

Definition at line 46 of file EGammaMvaEleEstimatorCSA14.h.

46  {
47  kTrig = 0, // MVA for triggering electrons
48  kNonTrig = 1, // MVA for non-triggering electrons
49  kNonTrigPhys14 = 2, // MVA for non-triggering electrons in Phys14
50  };

Constructor & Destructor Documentation

EGammaMvaEleEstimatorCSA14::EGammaMvaEleEstimatorCSA14 ( )

Definition at line 11 of file EGammaMvaEleEstimatorCSA14.cc.

11  :
12 fMethodname("BDTG method"),
13 fisInitialized(kFALSE),
14 fMVAType(kTrig),
15 fUseBinnedVersion(kTRUE),
16 fNMVABins(0)
17 {
18  // Constructor.
19 }
EGammaMvaEleEstimatorCSA14::~EGammaMvaEleEstimatorCSA14 ( )

Definition at line 22 of file EGammaMvaEleEstimatorCSA14.cc.

References fTMVAReader, and mps_fire::i.

23 {
24  for (unsigned int i=0;i<fTMVAReader.size(); ++i) {
25  if (fTMVAReader[i]) delete fTMVAReader[i];
26  }
27 }
mps_fire.i
i
Definition: mps_fire.py:156
EGammaMvaEleEstimatorCSA14::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimatorCSA14.h:84

Member Function Documentation

void EGammaMvaEleEstimatorCSA14::bindVariables ( )

Definition at line 487 of file EGammaMvaEleEstimatorCSA14.cc.

References fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_kfchi2, fMVAVar_OneMinusE1x5E5x5, fMVAVar_R9, fMVAVar_spp, and edm::detail::isnan().

Referenced by isInitialized(), and mvaValue().

487  {
488 
489  // this binding is needed for variables that sometime diverge.
490 
491 
492  if(fMVAVar_fbrem < -1.)
493  fMVAVar_fbrem = -1.;
494 
495  fMVAVar_deta = fabs(fMVAVar_deta);
496  if(fMVAVar_deta > 0.06)
497  fMVAVar_deta = 0.06;
498 
499 
500  fMVAVar_dphi = fabs(fMVAVar_dphi);
501  if(fMVAVar_dphi > 0.6)
502  fMVAVar_dphi = 0.6;
503 
504 
505  if(fMVAVar_EoP > 20.)
506  fMVAVar_EoP = 20.;
507 
508  if(fMVAVar_eleEoPout > 20.)
509  fMVAVar_eleEoPout = 20.;
510 
511 
512  fMVAVar_detacalo = fabs(fMVAVar_detacalo);
513  if(fMVAVar_detacalo > 0.2)
514  fMVAVar_detacalo = 0.2;
515 
516  if(fMVAVar_OneMinusE1x5E5x5 < -1.)
517  fMVAVar_OneMinusE1x5E5x5 = -1;
518 
519  if(fMVAVar_OneMinusE1x5E5x5 > 2.)
520  fMVAVar_OneMinusE1x5E5x5 = 2.;
521 
522 
523 
524  if(fMVAVar_R9 > 5)
525  fMVAVar_R9 = 5;
526 
527  if(fMVAVar_gsfchi2 > 200.)
528  fMVAVar_gsfchi2 = 200;
529 
530 
531  if(fMVAVar_kfchi2 > 10.)
532  fMVAVar_kfchi2 = 10.;
533 
534 
535  // Needed for a bug in CMSSW_420, fixed in more recent CMSSW versions
536  if(std::isnan(fMVAVar_spp))
537  fMVAVar_spp = 0.;
538 
539 
540  return;
541 }
edm::detail::isnan
bool isnan(float x)
Definition: math.h:13
UInt_t EGammaMvaEleEstimatorCSA14::GetMVABin ( double  eta,
double  pt 
) const

Definition at line 206 of file EGammaMvaEleEstimatorCSA14.cc.

References stringResolutionProvider_cfi::bin, fMVAType, kNonTrig, kNonTrigPhys14, and kTrig.

Referenced by isInitialized(), and mvaValue().

206  {
207 
208  //Default is to return the first bin
209  unsigned int bin = 0;
210 
211 
212  if (fMVAType == EGammaMvaEleEstimatorCSA14::kNonTrig ) {
213  bin = 0;
214  if (pt < 10 && fabs(eta) < 1.479) bin = 0;
215  if (pt < 10 && fabs(eta) >= 1.479) bin = 1;
216  if (pt >= 10 && fabs(eta) < 1.479) bin = 2;
217  if (pt >= 10 && fabs(eta) >= 1.479) bin = 3;
218  }
219 
220 
221  if (fMVAType == EGammaMvaEleEstimatorCSA14::kTrig
222  ) {
223  bin = 0;
224  if (pt >= 10 && fabs(eta) < 1.479) bin = 0;
225  if (pt >= 10 && fabs(eta) >= 1.479) bin = 1;
226  }
227 
228  if (fMVAType == EGammaMvaEleEstimatorCSA14::kNonTrigPhys14 ){
229  bin = 0;
230  if (pt < 10 && fabs(eta) < 0.8) bin = 0;
231  if (pt < 10 && fabs(eta) >= 0.8 && fabs(eta) < 1.479) bin = 1;
232  if (pt < 10 && fabs(eta) >= 1.479) bin = 2;
233  if (pt >= 10 && fabs(eta) < 0.8) bin = 3;
234  if (pt >= 10 && fabs(eta) >= 0.8 && fabs(eta) < 1.479) bin = 4;
235  if (pt >= 10 && fabs(eta) >= 1.479) bin = 5;
236  }
237 
238  return bin;
239 }
stringResolutionProvider_cfi.eta
eta
Definition: stringResolutionProvider_cfi.py:14
stringResolutionProvider_cfi.bin
bin
set the eta bin as selection string.
Definition: stringResolutionProvider_cfi.py:11
void EGammaMvaEleEstimatorCSA14::initialize ( std::string  methodName,
std::string  weightsfile,
EGammaMvaEleEstimatorCSA14::MVAType  type 
)

Definition at line 30 of file EGammaMvaEleEstimatorCSA14.cc.

Referenced by heppy::EGammaMvaEleEstimatorFWLite::initialize().

33 {
34 
35  std::vector<std::string> tempWeightFileVector;
36  tempWeightFileVector.push_back(weightsfile);
37  initialize(methodName,type,kFALSE,tempWeightFileVector);
38 }
type
type
Definition: HCALResponse.h:21
EGammaMvaEleEstimatorCSA14::initialize
void initialize(std::string methodName, std::string weightsfile, EGammaMvaEleEstimatorCSA14::MVAType type)
Definition: EGammaMvaEleEstimatorCSA14.cc:30
void EGammaMvaEleEstimatorCSA14::initialize ( std::string  methodName,
EGammaMvaEleEstimatorCSA14::MVAType  type,
Bool_t  useBinnedVersion,
std::vector< std::string >  weightsfiles 
)

Definition at line 42 of file EGammaMvaEleEstimatorCSA14.cc.

References gather_cfg::cout, data, Exception, FrontierConditions_GlobalTag_cff::file, fisInitialized, fMethodname, fMVAType, fMVAVar_abseta, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_isBarrel, fMVAVar_isEndcap, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_SCeta, fMVAVar_see, fMVAVar_spp, fNMVABins, fTMVAMethod, fTMVAReader, fUseBinnedVersion, mps_fire::i, kNonTrig, kNonTrigPhys14, and kTrig.

46  {
47 
48  //clean up first
49  for (unsigned int i=0;i<fTMVAReader.size(); ++i) {
50  if (fTMVAReader[i]) delete fTMVAReader[i];
51  }
52  fTMVAReader.clear();
53  fTMVAMethod.clear();
54  //initialize
55  fisInitialized = kTRUE;
56  fMVAType = type;
57  fMethodname = methodName;
58  fUseBinnedVersion = useBinnedVersion;
59 
60  //Define expected number of bins
61  UInt_t ExpectedNBins = 0;
62  if (type == kTrig) {
63  ExpectedNBins = 2;
64  }
65  else if (type == kNonTrig) {
66  ExpectedNBins = 4;
67  }
68  else if (type ==kNonTrigPhys14) {
69  ExpectedNBins = 6;
70  }
71 
72  fNMVABins = ExpectedNBins;
73 
74  //Check number of weight files given
75  if (fNMVABins != weightsfiles.size() ) {
76  std::cout << "Error: Expected Number of bins = " << fNMVABins << " does not equal to weightsfiles.size() = "
77  << weightsfiles.size() << std::endl;
78 
79  #ifndef STANDALONE
80  assert(fNMVABins == weightsfiles.size());
81  #endif
82  }
83 
84  //Loop over all bins
85  for (unsigned int i=0;i<fNMVABins; ++i) {
86 
87  TMVA::Reader *tmpTMVAReader = new TMVA::Reader( "!Color:!Silent:Error" );
88  tmpTMVAReader->SetVerbose(kTRUE);
89 
90  if (type == kTrig) {
91  // Pure tracking variables
92  // Pure tracking variables
93  tmpTMVAReader->AddVariable("fBrem", &fMVAVar_fbrem);
94  tmpTMVAReader->AddVariable("kfchi2", &fMVAVar_kfchi2);
95  tmpTMVAReader->AddVariable("kfhits", &fMVAVar_kfhits);
96  tmpTMVAReader->AddVariable("gsfChi2", &fMVAVar_gsfchi2);
97 
98  // Geometrical matchings
99  tmpTMVAReader->AddVariable("eledeta", &fMVAVar_deta);
100  tmpTMVAReader->AddVariable("eledphi", &fMVAVar_dphi);
101  tmpTMVAReader->AddVariable("detacalo", &fMVAVar_detacalo);
102 
103  // Pure ECAL -> shower shapes
104  tmpTMVAReader->AddVariable("noZSsee", &fMVAVar_see);
105  tmpTMVAReader->AddVariable("noZSspp", &fMVAVar_spp);
106  tmpTMVAReader->AddVariable("etawidth", &fMVAVar_etawidth);
107  tmpTMVAReader->AddVariable("phiwidth", &fMVAVar_phiwidth);
108  tmpTMVAReader->AddVariable("noZSe1x5e5x5", &fMVAVar_OneMinusE1x5E5x5);
109  tmpTMVAReader->AddVariable("noZSr9", &fMVAVar_R9);
110 
111  // Energy matching
112  tmpTMVAReader->AddVariable("HtoE", &fMVAVar_HoE);
113  tmpTMVAReader->AddVariable("EoP", &fMVAVar_EoP);
114  tmpTMVAReader->AddVariable("IoEmIoP", &fMVAVar_IoEmIoP);
115  tmpTMVAReader->AddVariable("EEleoPout", &fMVAVar_eleEoPout);
116  if(i == 1 ) tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
117 
118  tmpTMVAReader->AddSpectator("pt", &fMVAVar_pt);
119  tmpTMVAReader->AddSpectator("absEta", &fMVAVar_abseta);
120  }
121 
122 
123  if ((type == kNonTrig)||(type == kNonTrigPhys14)) {
124 
125  tmpTMVAReader->AddVariable("ele_kfhits", &fMVAVar_kfhits);
126  // Pure ECAL -> shower shapes
127  tmpTMVAReader->AddVariable("ele_oldsigmaietaieta", &fMVAVar_see);
128  tmpTMVAReader->AddVariable("ele_oldsigmaiphiiphi", &fMVAVar_spp);
129  tmpTMVAReader->AddVariable("ele_oldcircularity", &fMVAVar_OneMinusE1x5E5x5);
130  tmpTMVAReader->AddVariable("ele_oldr9", &fMVAVar_R9);
131  tmpTMVAReader->AddVariable("ele_scletawidth", &fMVAVar_etawidth);
132  tmpTMVAReader->AddVariable("ele_sclphiwidth", &fMVAVar_phiwidth);
133  tmpTMVAReader->AddVariable("ele_he", &fMVAVar_HoE);
134  if ((type == kNonTrig)&&(i == 1 || i == 3)) tmpTMVAReader->AddVariable("ele_psEoverEraw",&fMVAVar_PreShowerOverRaw);
135  if ((type == kNonTrigPhys14)&&(i == 2 || i == 5)) tmpTMVAReader->AddVariable("ele_psEoverEraw",&fMVAVar_PreShowerOverRaw);
136 
137 
138  //Pure tracking variables
139  tmpTMVAReader->AddVariable("ele_kfchi2", &fMVAVar_kfchi2);
140  tmpTMVAReader->AddVariable("ele_chi2_hits", &fMVAVar_gsfchi2);
141  // Energy matching
142  tmpTMVAReader->AddVariable("ele_fbrem", &fMVAVar_fbrem);
143  tmpTMVAReader->AddVariable("ele_ep", &fMVAVar_EoP);
144  tmpTMVAReader->AddVariable("ele_eelepout", &fMVAVar_eleEoPout);
145  tmpTMVAReader->AddVariable("ele_IoEmIop", &fMVAVar_IoEmIoP);
146 
147  // Geometrical matchings
148  tmpTMVAReader->AddVariable("ele_deltaetain", &fMVAVar_deta);
149  tmpTMVAReader->AddVariable("ele_deltaphiin", &fMVAVar_dphi);
150  tmpTMVAReader->AddVariable("ele_deltaetaseed", &fMVAVar_detacalo);
151 
152 
153 
154  tmpTMVAReader->AddSpectator("ele_pT", &fMVAVar_pt);
155  tmpTMVAReader->AddSpectator("ele_isbarrel", &fMVAVar_isBarrel);
156  tmpTMVAReader->AddSpectator("ele_isendcap", &fMVAVar_isEndcap);
157  if (type == kNonTrigPhys14) tmpTMVAReader->AddSpectator("scl_eta", &fMVAVar_SCeta);
158 
159  }
160 
161 
162 
163 
164 #ifndef STANDALONE
165  if ((fMethodname.find("BDT") == 0) && (weightsfiles[i].rfind(".xml.gz") == weightsfiles[i].length()-strlen(".xml.gz"))) {
166  gzFile file = gzopen(weightsfiles[i].c_str(), "rb");
167  if (file == nullptr) { std::cout << "Error opening gzip file associated to " << weightsfiles[i] << std::endl; throw cms::Exception("Configuration","Error reading zipped XML file"); }
168  std::vector<char> data;
169  data.reserve(1024*1024*10);
170  unsigned int bufflen = 32*1024;
171  char *buff = reinterpret_cast<char *>(malloc(bufflen));
172  if (buff == nullptr) { std::cout << "Error creating buffer for " << weightsfiles[i] << std::endl; gzclose(file); throw cms::Exception("Configuration","Error reading zipped XML file"); }
173  int read;
174  while ((read = gzread(file, buff, bufflen)) != 0) {
175  if (read == -1) { std::cout << "Error reading gzip file associated to " << weightsfiles[i] << ": " << gzerror(file,&read) << std::endl; gzclose(file); free(buff); throw cms::Exception("Configuration","Error reading zipped XML file"); }
176  data.insert(data.end(), buff, buff+read);
177  }
178  if (gzclose(file) != Z_OK) { std::cout << "Error closing gzip file associated to " << weightsfiles[i] << std::endl; }
179  free(buff);
180  data.push_back('\0'); // IMPORTANT
181  fTMVAMethod.push_back(dynamic_cast<TMVA::MethodBase*>(tmpTMVAReader->BookMVA(TMVA::Types::kBDT, &data[0])));
182  } else {
183  if (weightsfiles[i].rfind(".xml.gz") == weightsfiles[i].length()-strlen(".xml.gz")) {
184  std::cout << "Error: xml.gz unsupported for method " << fMethodname << ", weight file " << weightsfiles[i] << std::endl; throw cms::Exception("Configuration","Error reading zipped XML file");
185  }
186  fTMVAMethod.push_back(dynamic_cast<TMVA::MethodBase*>(tmpTMVAReader->BookMVA(fMethodname , weightsfiles[i])));
187  }
188 #else
189  if (weightsfiles[i].rfind(".xml.gz") == weightsfiles[i].length()-strlen(".xml.gz")) {
190  std::cout << "Error: xml.gz unsupported for method " << fMethodname << ", weight file " << weightsfiles[i] << std::endl; abort();
191  }
192  fTMVAMethod.push_back(dynamic_cast<TMVA::MethodBase*>(tmpTMVAReader->BookMVA(fMethodname , weightsfiles[i])));
193 #endif
194  std::cout << "MVABin " << i << " : MethodName = " << fMethodname
195  << " , type == " << type << " , "
196  << "Load weights file : " << weightsfiles[i]
197  << std::endl;
198  fTMVAReader.push_back(tmpTMVAReader);
199  }
200  std::cout << "Electron ID MVA Completed\n";
201 
202 }
type
type
Definition: HCALResponse.h:21
mps_fire.i
i
Definition: mps_fire.py:156
EGammaMvaEleEstimatorCSA14::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimatorCSA14.h:85
EGammaMvaEleEstimatorCSA14::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimatorCSA14.h:84
FrontierConditions_GlobalTag_cff.file
file
Definition: FrontierConditions_GlobalTag_cff.py:13
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
gather_cfg.cout
cout
Definition: gather_cfg.py:145
Bool_t EGammaMvaEleEstimatorCSA14::isInitialized ( ) const
inline

Definition at line 60 of file EGammaMvaEleEstimatorCSA14.h.

References bindVariables(), stringResolutionProvider_cfi::eta, fisInitialized, GetMVABin(), mvaValue(), and EnergyCorrector::pt.

60 { return fisInitialized; }
Double_t EGammaMvaEleEstimatorCSA14::mvaValue ( const reco::GsfElectron &  ele,
const reco::Vertex &  vertex,
const TransientTrackBuilder &  transientTrackBuilder,
noZS::EcalClusterLazyTools  myEcalCluster,
bool  printDebug = kFALSE 
)

Definition at line 248 of file EGammaMvaEleEstimatorCSA14.cc.

References IPTools::absoluteImpactParameter3D(), stringResolutionProvider_cfi::bin, bindVariables(), TransientTrackBuilder::build(), reco::GsfElectron::closestCtfTrackRef(), gather_cfg::cout, reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e1x5(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e3x3(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e5x5(), reco::GsfElectron::e5x5(), reco::GsfElectron::ecalEnergy(), reco::GsfElectron::eEleClusterOverPout(), reco::GsfElectron::eSuperClusterOverP(), reco::GsfElectron::fbrem(), fisInitialized, fMethodname, fMVAType, fMVAVar_abseta, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_ip3d, fMVAVar_ip3dSig, fMVAVar_isBarrel, fMVAVar_isEndcap, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_kfhitsall, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_SCeta, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, GetMVABin(), reco::GsfElectron::gsfTrack(), reco::GsfElectron::hadronicOverEm(), edm::Ref< C, T, F >::isAvailable(), edm::detail::isnan(), edm::Ref< C, T, F >::isNonnull(), kNonTrig, kNonTrigPhys14, kTrig, EcalClusterLazyToolsT< EcalClusterToolsImpl >::localCovariances(), reco::LeafCandidate::p(), reco::Vertex::position(), reco::LeafCandidate::pt(), mathSSE::sqrt(), reco::GsfElectron::superCluster(), and groupFilesInBlocks::tt.

Referenced by isInitialized(), and heppy::EGammaMvaEleEstimatorFWLite::mvaValue().

252  {
253 
254  if (!fisInitialized) {
255  std::cout << "Error: EGammaMvaEleEstimatorCSA14 not properly initialized.\n";
256  return -9999;
257  }
258 
259  if ( (fMVAType != EGammaMvaEleEstimatorCSA14::kTrig) && (fMVAType != EGammaMvaEleEstimatorCSA14::kNonTrig) && (fMVAType != EGammaMvaEleEstimatorCSA14::kNonTrigPhys14) ) {
260  std::cout << "Error: This method should be called for kTrig or kNonTrig or kNonTrigPhys14 MVA only" <<std::endl;
261  return -9999;
262  }
263 
264  bool validKF= false;
265  reco::TrackRef myTrackRef = ele.closestCtfTrackRef();
266  validKF = (myTrackRef.isAvailable());
267  validKF = (myTrackRef.isNonnull());
268 
269  // Pure tracking variables
270  fMVAVar_fbrem = ele.fbrem();
271  fMVAVar_kfchi2 = (validKF) ? myTrackRef->normalizedChi2() : 0 ;
272  fMVAVar_kfhits = (validKF) ? myTrackRef->hitPattern().trackerLayersWithMeasurement() : -1. ;
273  fMVAVar_kfhitsall = (validKF) ? myTrackRef->numberOfValidHits() : -1. ; // save also this in your ntuple as possible alternative
274  fMVAVar_gsfchi2 = ele.gsfTrack()->normalizedChi2();
275 
276 
277  // Geometrical matchings
278  fMVAVar_deta = ele.deltaEtaSuperClusterTrackAtVtx();
279  fMVAVar_dphi = ele.deltaPhiSuperClusterTrackAtVtx();
280  fMVAVar_detacalo = ele.deltaEtaSeedClusterTrackAtCalo();
281 
282 
283  // Pure ECAL -> shower shapes
284  std::vector<float> vCov = myEcalCluster.localCovariances(*(ele.superCluster()->seed())) ;
285  if (!isnan(vCov[0])) fMVAVar_see = sqrt (vCov[0]); //EleSigmaIEtaIEta
286  else fMVAVar_see = 0.;
287  if (!isnan(vCov[2])) fMVAVar_spp = sqrt (vCov[2]); //EleSigmaIPhiIPhi
288  else fMVAVar_spp = 0.;
289 
290  fMVAVar_etawidth = ele.superCluster()->etaWidth();
291  fMVAVar_phiwidth = ele.superCluster()->phiWidth();
292  fMVAVar_OneMinusE1x5E5x5 = (ele.e5x5()) !=0. ? 1.-(myEcalCluster.e1x5(*(ele.superCluster()->seed()))/myEcalCluster.e5x5(*(ele.superCluster()->seed()))) : -1. ;
293  fMVAVar_R9 = myEcalCluster.e3x3(*(ele.superCluster()->seed())) / ele.superCluster()->rawEnergy();
294 
295  // Energy matching
296  fMVAVar_HoE = ele.hadronicOverEm();
297  fMVAVar_EoP = ele.eSuperClusterOverP();
298  fMVAVar_IoEmIoP = (1.0/ele.ecalEnergy()) - (1.0 / ele.p()); // in the future to be changed with ele.gsfTrack()->p()
299  fMVAVar_eleEoPout = ele.eEleClusterOverPout();
300  fMVAVar_PreShowerOverRaw= ele.superCluster()->preshowerEnergy() / ele.superCluster()->rawEnergy();
301 
302 
303  // Spectators
304  fMVAVar_eta = ele.superCluster()->eta();
305  fMVAVar_abseta = fabs(ele.superCluster()->eta());
306  fMVAVar_pt = ele.pt();
307  fMVAVar_isBarrel = (ele.superCluster()->eta()<1.479);
308  fMVAVar_isEndcap = (ele.superCluster()->eta()>1.479);
309  fMVAVar_SCeta = ele.superCluster()->eta();
310 
311 
312  // for triggering electrons get the impact parameteres
313  if(fMVAType == kTrig) {
314  //d0
315  if (ele.gsfTrack().isNonnull()) {
316  fMVAVar_d0 = (-1.0)*ele.gsfTrack()->dxy(vertex.position());
317  } else if (ele.closestCtfTrackRef().isNonnull()) {
318  fMVAVar_d0 = (-1.0)*ele.closestCtfTrackRef()->dxy(vertex.position());
319  } else {
320  fMVAVar_d0 = -9999.0;
321  }
322 
323  //default values for IP3D
324  fMVAVar_ip3d = -999.0;
325  fMVAVar_ip3dSig = 0.0;
326  if (ele.gsfTrack().isNonnull()) {
327  const double gsfsign = ( (-ele.gsfTrack()->dxy(vertex.position())) >=0 ) ? 1. : -1.;
328 
329  const reco::TransientTrack &tt = transientTrackBuilder.build(ele.gsfTrack());
330  const std::pair<bool,Measurement1D> &ip3dpv = IPTools::absoluteImpactParameter3D(tt,vertex);
331  if (ip3dpv.first) {
332  double ip3d = gsfsign*ip3dpv.second.value();
333  double ip3derr = ip3dpv.second.error();
334  fMVAVar_ip3d = ip3d;
335  fMVAVar_ip3dSig = ip3d/ip3derr;
336  }
337  }
338  }
339 
340  // evaluate
341  bindVariables();
342  Double_t mva = -9999;
343  if (fUseBinnedVersion) {
344  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
345  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
346  } else {
347  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
348  }
349 
350 
351 
352  if(printDebug) {
353  std::cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType <<std::endl;
354  std::cout << " fbrem " << fMVAVar_fbrem
355  << " kfchi2 " << fMVAVar_kfchi2
356  << " mykfhits " << fMVAVar_kfhits
357  << " gsfchi2 " << fMVAVar_gsfchi2
358  << " deta " << fMVAVar_deta
359  << " dphi " << fMVAVar_dphi
360  << " detacalo " << fMVAVar_detacalo
361  << " see " << fMVAVar_see
362  << " spp " << fMVAVar_spp
363  << " etawidth " << fMVAVar_etawidth
364  << " phiwidth " << fMVAVar_phiwidth
365  << " OneMinusE1x5E5x5 " << fMVAVar_OneMinusE1x5E5x5
366  << " R9 " << fMVAVar_R9
367  << " HoE " << fMVAVar_HoE
368  << " EoP " << fMVAVar_EoP
369  << " IoEmIoP " << fMVAVar_IoEmIoP
370  << " eleEoPout " << fMVAVar_eleEoPout
371  << " d0 " << fMVAVar_d0
372  << " ip3d " << fMVAVar_ip3d
373  << " eta " << fMVAVar_eta
374  << " pt " << fMVAVar_pt <<std::endl;
375  std::cout << " ### MVA " << mva <<std::endl;
376  }
377 
378 
379 
380  return mva;
381 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:576
reco::LeafCandidate::pt
virtual double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
EcalClusterLazyToolsT::e5x5
float e5x5(const reco::BasicCluster &cluster)
Definition: EcalClusterLazyTools.h:255
EGammaMvaEleEstimatorCSA14::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimatorCSA14.h:85
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:245
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:16
IPTools::absoluteImpactParameter3D
std::pair< bool, Measurement1D > absoluteImpactParameter3D(const reco::TransientTrack &transientTrack, const reco::Vertex &vertex)
Definition: IPTools.cc:37
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:750
EcalClusterLazyToolsT::e3x3
float e3x3(const reco::BasicCluster &cluster)
Definition: EcalClusterLazyTools.h:243
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:109
EGammaMvaEleEstimatorCSA14::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimatorCSA14.h:84
reco::GsfElectron::closestCtfTrackRef
TrackRef closestCtfTrackRef() const
Definition: GsfElectron.h:201
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:249
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:487
edm::detail::isnan
bool isnan(float x)
Definition: math.h:13
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
reco::GsfElectron::superCluster
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
Definition: GsfElectron.h:184
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:252
reco::LeafCandidate::p
virtual double p() const final
magnitude of momentum vector
Definition: LeafCandidate.h:108
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:248
stringResolutionProvider_cfi.bin
bin
set the eta bin as selection string.
Definition: stringResolutionProvider_cfi.py:11
EcalClusterLazyToolsT::localCovariances
std::vector< float > localCovariances(const reco::BasicCluster &cluster, float w0=4.7)
Definition: EcalClusterLazyTools.h:363
reco::GsfElectron::ecalEnergy
float ecalEnergy() const
Definition: GsfElectron.h:837
reco::GsfElectron::e5x5
float e5x5() const
Definition: GsfElectron.h:438
EGammaMvaEleEstimatorCSA14::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimatorCSA14.cc:206
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:250
EcalClusterLazyToolsT::e1x5
float e1x5(const reco::BasicCluster &cluster)
Definition: EcalClusterLazyTools.h:218
gather_cfg.cout
cout
Definition: gather_cfg.py:145
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:185
Double_t EGammaMvaEleEstimatorCSA14::mvaValue ( const pat::Electron &  ele,
bool  printDebug 
)

Definition at line 385 of file EGammaMvaEleEstimatorCSA14.cc.

References stringResolutionProvider_cfi::bin, bindVariables(), pat::Electron::closestCtfTrackRef(), gather_cfg::cout, reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), reco::GsfElectron::ecalEnergy(), reco::GsfElectron::eEleClusterOverPout(), reco::GsfElectron::eSuperClusterOverP(), reco::GsfElectron::fbrem(), fisInitialized, fMethodname, fMVAType, fMVAVar_abseta, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_isBarrel, fMVAVar_isEndcap, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_kfhitsall, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_SCeta, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, reco::GsfElectron::full5x5_e1x5(), reco::GsfElectron::full5x5_e5x5(), reco::GsfElectron::full5x5_r9(), reco::GsfElectron::full5x5_sigmaIetaIeta(), reco::GsfElectron::full5x5_sigmaIphiIphi(), fUseBinnedVersion, GetMVABin(), pat::Electron::gsfTrack(), reco::GsfElectron::hadronicOverEm(), edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), kNonTrig, kNonTrigPhys14, kTrig, reco::LeafCandidate::p(), reco::LeafCandidate::pt(), and pat::Electron::superCluster().

386  {
387 
388  if (!fisInitialized) {
389  std::cout << "Error: EGammaMvaEleEstimatorCSA14 not properly initialized.\n";
390  return -9999;
391  }
392 
393  if ( (fMVAType != EGammaMvaEleEstimatorCSA14::kTrig) && (fMVAType != EGammaMvaEleEstimatorCSA14::kNonTrig) && (fMVAType != EGammaMvaEleEstimatorCSA14::kNonTrigPhys14) ) {
394  std::cout << "Error: This method should be called for kTrig or kNonTrig or kNonTrigPhys14 MVA only" <<std::endl;
395  return -9999;
396  }
397 
398  bool validKF= false;
399  reco::TrackRef myTrackRef = ele.closestCtfTrackRef();
400  validKF = (myTrackRef.isAvailable());
401  validKF = (myTrackRef.isNonnull());
402 
403  // Pure tracking variables
404  fMVAVar_fbrem = ele.fbrem();
405  fMVAVar_kfchi2 = (validKF) ? myTrackRef->normalizedChi2() : 0 ;
406  fMVAVar_kfhits = (validKF) ? myTrackRef->hitPattern().trackerLayersWithMeasurement() : -1. ;
407  fMVAVar_kfhitsall = (validKF) ? myTrackRef->numberOfValidHits() : -1. ; // save also this in your ntuple as possible alternative
408  fMVAVar_gsfchi2 = ele.gsfTrack()->normalizedChi2();
409 
410 
411  // Geometrical matchings
412  fMVAVar_deta = ele.deltaEtaSuperClusterTrackAtVtx();
413  fMVAVar_dphi = ele.deltaPhiSuperClusterTrackAtVtx();
414  fMVAVar_detacalo = ele.deltaEtaSeedClusterTrackAtCalo();
415 
416 
417  // Pure ECAL -> shower shapes
418  fMVAVar_see = ele.full5x5_sigmaIetaIeta(); //EleSigmaIEtaIEta
419  fMVAVar_spp = ele.full5x5_sigmaIphiIphi(); //EleSigmaIPhiIPhi
420 
421  fMVAVar_etawidth = ele.superCluster()->etaWidth();
422  fMVAVar_phiwidth = ele.superCluster()->phiWidth();
423  fMVAVar_OneMinusE1x5E5x5 = (ele.full5x5_e5x5()) !=0. ? 1.-(ele.full5x5_e1x5()/ele.full5x5_e5x5()) : -1. ;
424  fMVAVar_R9 = ele.full5x5_r9();
425 
426  // Energy matching
427  fMVAVar_HoE = ele.hadronicOverEm();
428  fMVAVar_EoP = ele.eSuperClusterOverP();
429  fMVAVar_IoEmIoP = (1.0/ele.ecalEnergy()) - (1.0 / ele.p()); // in the future to be changed with ele.gsfTrack()->p()
430  fMVAVar_eleEoPout = ele.eEleClusterOverPout();
431  fMVAVar_PreShowerOverRaw= ele.superCluster()->preshowerEnergy() / ele.superCluster()->rawEnergy();
432 
433 
434  // Spectators
435  fMVAVar_eta = ele.superCluster()->eta();
436  fMVAVar_abseta = fabs(ele.superCluster()->eta());
437  fMVAVar_pt = ele.pt();
438  fMVAVar_isBarrel = (ele.superCluster()->eta()<1.479);
439  fMVAVar_isEndcap = (ele.superCluster()->eta()>1.479);
440  fMVAVar_SCeta = ele.superCluster()->eta();
441 
442 
443 
444  // evaluate
445  bindVariables();
446  Double_t mva = -9999;
447  if (fUseBinnedVersion) {
448  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
449  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
450  } else {
451  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
452  }
453 
454 
455 
456  if(printDebug) {
457  std::cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType <<std::endl;
458  std::cout << " fbrem " << fMVAVar_fbrem
459  << " kfchi2 " << fMVAVar_kfchi2
460  << " mykfhits " << fMVAVar_kfhits
461  << " gsfchi2 " << fMVAVar_gsfchi2
462  << " deta " << fMVAVar_deta
463  << " dphi " << fMVAVar_dphi
464  << " detacalo " << fMVAVar_detacalo
465  << " see " << fMVAVar_see
466  << " spp " << fMVAVar_spp
467  << " etawidth " << fMVAVar_etawidth
468  << " phiwidth " << fMVAVar_phiwidth
469  << " OneMinusE1x5E5x5 " << fMVAVar_OneMinusE1x5E5x5
470  << " R9 " << fMVAVar_R9
471  << " HoE " << fMVAVar_HoE
472  << " EoP " << fMVAVar_EoP
473  << " IoEmIoP " << fMVAVar_IoEmIoP
474  << " eleEoPout " << fMVAVar_eleEoPout
475  << " eta " << fMVAVar_eta
476  << " pt " << fMVAVar_pt <<std::endl;
477  std::cout << " ### MVA " << mva <<std::endl;
478  }
479 
480 
481 
482  return mva;
483 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:576
reco::LeafCandidate::pt
virtual double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
EGammaMvaEleEstimatorCSA14::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimatorCSA14.h:85
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:245
reco::GsfElectron::full5x5_e5x5
float full5x5_e5x5() const
Definition: GsfElectron.h:459
reco::GsfElectron::full5x5_e1x5
float full5x5_e1x5() const
Definition: GsfElectron.h:457
reco::GsfElectron::full5x5_sigmaIphiIphi
float full5x5_sigmaIphiIphi() const
Definition: GsfElectron.h:456
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:750
EGammaMvaEleEstimatorCSA14::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimatorCSA14.h:84
reco::GsfElectron::full5x5_sigmaIetaIeta
float full5x5_sigmaIetaIeta() const
Definition: GsfElectron.h:455
pat::Electron::superCluster
reco::SuperClusterRef superCluster() const
override the reco::GsfElectron::superCluster method, to access the internal storage of the superclust...
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:249
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:487
pat::Electron::gsfTrack
reco::GsfTrackRef gsfTrack() const
override the reco::GsfElectron::gsfTrack method, to access the internal storage of the supercluster ...
pat::Electron::closestCtfTrackRef
reco::TrackRef closestCtfTrackRef() const
override the reco::GsfElectron::closestCtfTrackRef method, to access the internal storage of the trac...
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:252
reco::LeafCandidate::p
virtual double p() const final
magnitude of momentum vector
Definition: LeafCandidate.h:108
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:248
stringResolutionProvider_cfi.bin
bin
set the eta bin as selection string.
Definition: stringResolutionProvider_cfi.py:11
reco::GsfElectron::ecalEnergy
float ecalEnergy() const
Definition: GsfElectron.h:837
reco::GsfElectron::full5x5_r9
float full5x5_r9() const
Definition: GsfElectron.h:460
EGammaMvaEleEstimatorCSA14::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimatorCSA14.cc:206
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:250
gather_cfg.cout
cout
Definition: gather_cfg.py:145

Member Data Documentation

Bool_t EGammaMvaEleEstimatorCSA14::fisInitialized
private

Definition at line 87 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), isInitialized(), and mvaValue().

std::string EGammaMvaEleEstimatorCSA14::fMethodname
private

Definition at line 86 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

MVAType EGammaMvaEleEstimatorCSA14::fMVAType
private

Definition at line 88 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by GetMVABin(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_abseta
private

Definition at line 121 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_ChargedIso_DR0p0To0p1
private

Definition at line 128 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_ChargedIso_DR0p1To0p2
private

Definition at line 129 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_ChargedIso_DR0p2To0p3
private

Definition at line 130 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_ChargedIso_DR0p3To0p4
private

Definition at line 131 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_ChargedIso_DR0p4To0p5
private

Definition at line 132 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_d0
private

Definition at line 116 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_deta
private

Definition at line 98 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_detacalo
private

Definition at line 100 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_dphi
private

Definition at line 99 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_eleEoPout
private

Definition at line 112 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_EoP
private

Definition at line 110 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_EoPout
private

Definition at line 113 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_eta
private

Definition at line 120 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_etawidth
private

Definition at line 104 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_fbrem
private

Definition at line 92 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_GammaIso_DR0p0To0p1
private

Definition at line 133 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_GammaIso_DR0p1To0p2
private

Definition at line 134 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_GammaIso_DR0p2To0p3
private

Definition at line 135 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_GammaIso_DR0p3To0p4
private

Definition at line 136 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_GammaIso_DR0p4To0p5
private

Definition at line 137 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_gsfchi2
private

Definition at line 96 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_HoE
private

Definition at line 109 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_IoEmIoP
private

Definition at line 111 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_ip3d
private

Definition at line 117 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_ip3dSig
private

Definition at line 118 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_isBarrel
private

Definition at line 124 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_isEndcap
private

Definition at line 125 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_kfchi2
private

Definition at line 93 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_kfhits
private

Definition at line 94 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_kfhitsall
private

Definition at line 95 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_NeutralHadronIso_DR0p0To0p1
private

Definition at line 138 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_NeutralHadronIso_DR0p1To0p2
private

Definition at line 139 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_NeutralHadronIso_DR0p2To0p3
private

Definition at line 140 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_NeutralHadronIso_DR0p3To0p4
private

Definition at line 141 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_NeutralHadronIso_DR0p4To0p5
private

Definition at line 142 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_OneMinusE1x5E5x5
private

Definition at line 106 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_phiwidth
private

Definition at line 105 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_PreShowerOverRaw
private

Definition at line 114 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_pt
private

Definition at line 122 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_R9
private

Definition at line 107 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_rho
private

Definition at line 123 of file EGammaMvaEleEstimatorCSA14.h.

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_SCeta
private

Definition at line 126 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_see
private

Definition at line 102 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Float_t EGammaMvaEleEstimatorCSA14::fMVAVar_spp
private

Definition at line 103 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by bindVariables(), initialize(), and mvaValue().

UInt_t EGammaMvaEleEstimatorCSA14::fNMVABins
private

Definition at line 90 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize().

std::vector<TMVA::MethodBase*> EGammaMvaEleEstimatorCSA14::fTMVAMethod
private

Definition at line 85 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

std::vector<TMVA::Reader*> EGammaMvaEleEstimatorCSA14::fTMVAReader
private

Definition at line 84 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), mvaValue(), and ~EGammaMvaEleEstimatorCSA14().

Bool_t EGammaMvaEleEstimatorCSA14::fUseBinnedVersion
private

Definition at line 89 of file EGammaMvaEleEstimatorCSA14.h.

Referenced by initialize(), and mvaValue().

Generated for CMSSW Reference Manual by doxygen 1.8.11