CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
AntiElectronIDMVA3 Class Reference

#include <AntiElectronIDMVA3.h>

Public Member Functions

 AntiElectronIDMVA3 ()
 
void Initialize_from_file (const std::string &methodName, const std::string &oneProngNoEleMatch0Pi0woGSF_BL, const std::string &oneProngNoEleMatch0Pi0wGSF_BL, const std::string &oneProngNoEleMatch1Pi0woGSF_BL, const std::string &oneProngNoEleMatch1Pi0wGSF_BL, const std::string &oneProng0Pi0woGSF_BL, const std::string &oneProng0Pi0wGSF_BL, const std::string &oneProng1Pi0woGSF_BL, const std::string &oneProng1Pi0wGSF_BL, const std::string &oneProngNoEleMatch0Pi0woGSF_EC, const std::string &oneProngNoEleMatch0Pi0wGSF_EC, const std::string &oneProngNoEleMatch1Pi0woGSF_EC, const std::string &oneProngNoEleMatch1Pi0wGSF_EC, const std::string &oneProng0Pi0woGSF_EC, const std::string &oneProng0Pi0wGSF_EC, const std::string &oneProng1Pi0woGSF_EC, const std::string &oneProng1Pi0wGSF_EC)
 
void Initialize_from_string (const std::string &methodName, const std::string &oneProngNoEleMatch0Pi0woGSF_BL, const std::string &oneProngNoEleMatch0Pi0wGSF_BL, const std::string &oneProngNoEleMatch1Pi0woGSF_BL, const std::string &oneProngNoEleMatch1Pi0wGSF_BL, const std::string &oneProng0Pi0woGSF_BL, const std::string &oneProng0Pi0wGSF_BL, const std::string &oneProng1Pi0woGSF_BL, const std::string &oneProng1Pi0wGSF_BL, const std::string &oneProngNoEleMatch0Pi0woGSF_EC, const std::string &oneProngNoEleMatch0Pi0wGSF_EC, const std::string &oneProngNoEleMatch1Pi0woGSF_EC, const std::string &oneProngNoEleMatch1Pi0wGSF_EC, const std::string &oneProng0Pi0woGSF_EC, const std::string &oneProng0Pi0wGSF_EC, const std::string &oneProng1Pi0woGSF_EC, const std::string &oneProng1Pi0wGSF_EC)
 
double MVAValue (Float_t TauEtaAtEcalEntrance, Float_t TauPt, Float_t TaudCrackEta, Float_t TaudCrackPhi, Float_t TauEmFraction, Float_t TauSignalPFGammaCands, Float_t TauLeadPFChargedHadrHoP, Float_t TauLeadPFChargedHadrEoP, Float_t TauVisMass, Float_t TauHadrMva, const std::vector< Float_t > &GammasdEta, const std::vector< Float_t > &GammasdPhi, const std::vector< Float_t > &GammasPt, Float_t TauKFNumHits, Float_t TauGSFNumHits, Float_t TauGSFChi2, Float_t TauGSFTrackResol, Float_t TauGSFTracklnPt, Float_t TauGSFTrackEta, Float_t TauPhi, Float_t TauSignalPFChargedCands, Float_t TauHasGsf, Float_t ElecEta, Float_t ElecPhi, Float_t ElecPt, Float_t ElecEe, Float_t ElecEgamma, Float_t ElecPin, Float_t ElecPout, Float_t ElecFbrem, Float_t ElecChi2GSF, Float_t ElecGSFNumHits, Float_t ElecGSFTrackResol, Float_t ElecGSFTracklnPt, Float_t ElecGSFTrackEta)
 
double MVAValue (Float_t TauEtaAtEcalEntrance, Float_t TauPt, Float_t TaudCrackEta, Float_t TaudCrackPhi, Float_t TauEmFraction, Float_t TauSignalPFGammaCands, Float_t TauLeadPFChargedHadrHoP, Float_t TauLeadPFChargedHadrEoP, Float_t TauVisMass, Float_t TauHadrMva, Float_t TauGammaEtaMom, Float_t TauGammaPhiMom, Float_t TauGammaEnFrac, Float_t TauKFNumHits, Float_t TauGSFNumHits, Float_t TauGSFChi2, Float_t TauGSFTrackResol, Float_t TauGSFTracklnPt, Float_t TauGSFTrackEta, Float_t TauPhi, Float_t TauSignalPFChargedCands, Float_t TauHasGsf, Float_t ElecEta, Float_t ElecPhi, Float_t ElecPt, Float_t ElecEe, Float_t ElecEgamma, Float_t ElecPin, Float_t ElecPout, Float_t ElecFbrem, Float_t ElecChi2GSF, Float_t ElecGSFNumHits, Float_t ElecGSFTrackResol, Float_t ElecGSFTracklnPt, Float_t ElecGSFTrackEta)
 
double MVAValue (const reco::PFTau &thePFTau, const reco::GsfElectron &theGsfEle)
 
double MVAValue (const reco::PFTau &thePFTau)
 
 ~AntiElectronIDMVA3 ()
 

Private Member Functions

void bookMVAs ()
 
double dCrackEta (double eta)
 
double dCrackPhi (double phi, double eta)
 
double minimum (double a, double b)
 

Private Attributes

Float_t Elec_Chi2GSF_
 
Float_t Elec_EgammaOverPdif_
 
Float_t Elec_EtotOverPin_
 
Float_t Elec_Fbrem_
 
Float_t Elec_GSFNumHits_
 
Float_t Elec_GSFTrackEta_
 
Float_t Elec_GSFTracklnPt_
 
Float_t Elec_GSFTrackResol_
 
TMVA::Reader * fTMVAReader_ [16]
 
Float_t GammadEta_
 
Float_t GammadPhi_
 
Float_t GammadPt_
 
Bool_t isInitialized_
 
std::string methodName_
 
Float_t Tau_dCrackEta_
 
Float_t Tau_dCrackPhi_
 
Float_t Tau_EmFraction_
 
Float_t Tau_EtaAtEcalEntrance_
 
Float_t Tau_GammaEnFrac_
 
Float_t Tau_GammaEtaMom_
 
Float_t Tau_GammaPhiMom_
 
Float_t Tau_GSFChi2_
 
Float_t Tau_GSFTrackEta_
 
Float_t Tau_GSFTracklnPt_
 
Float_t Tau_GSFTrackResol_
 
Float_t Tau_HadrEoP_
 
Float_t Tau_HadrHoP_
 
Float_t Tau_HadrMva_
 
Float_t Tau_NumGammaCands_
 
Float_t Tau_NumHitsVariable_
 
Float_t Tau_Pt_
 
Float_t Tau_VisMass_
 
int verbosity_
 

Detailed Description

Definition at line 27 of file AntiElectronIDMVA3.h.

Constructor & Destructor Documentation

AntiElectronIDMVA3::AntiElectronIDMVA3 ( )

Definition at line 6 of file AntiElectronIDMVA3.cc.

References fTMVAReader_, i, and verbosity_.

7  : isInitialized_(kFALSE),
8  methodName_("BDTG")
9 {
10  for ( unsigned i = 0; i < 16; ++i ) {
11  fTMVAReader_[i] = 0;
12  }
13 
14  verbosity_ = 1;
15 }
i
int i
Definition: DBlmapReader.cc:9
AntiElectronIDMVA3::fTMVAReader_
TMVA::Reader * fTMVAReader_[16]
Definition: AntiElectronIDMVA3.h:156
AntiElectronIDMVA3::~AntiElectronIDMVA3 ( )

Definition at line 17 of file AntiElectronIDMVA3.cc.

References fTMVAReader_, and i.

18 {
19  for ( unsigned i = 0; i < 16; ++i ) {
20  if ( fTMVAReader_[i] ) delete fTMVAReader_[i];
21  }
22 }
i
int i
Definition: DBlmapReader.cc:9
AntiElectronIDMVA3::fTMVAReader_
TMVA::Reader * fTMVAReader_[16]
Definition: AntiElectronIDMVA3.h:156

Member Function Documentation

void AntiElectronIDMVA3::bookMVAs ( )
private

Definition at line 134 of file AntiElectronIDMVA3.cc.

References Elec_Chi2GSF_, Elec_EgammaOverPdif_, Elec_EtotOverPin_, Elec_Fbrem_, Elec_GSFNumHits_, Elec_GSFTrackEta_, Elec_GSFTracklnPt_, Elec_GSFTrackResol_, fTMVAReader_, GammadEta_, GammadPhi_, GammadPt_, k_NoEleMatchwGwGSF_BL, k_NoEleMatchwGwGSF_EC, k_NoEleMatchwGwoGSF_BL, k_NoEleMatchwGwoGSF_EC, k_NoEleMatchwoGwGSF_BL, k_NoEleMatchwoGwGSF_EC, k_NoEleMatchwoGwoGSF_BL, k_NoEleMatchwoGwoGSF_EC, k_wGwGSF_BL, k_wGwGSF_EC, k_wGwoGSF_BL, k_wGwoGSF_EC, k_woGwGSF_BL, k_woGwGSF_EC, k_woGwoGSF_BL, k_woGwoGSF_EC, Tau_dCrackEta_, Tau_dCrackPhi_, Tau_EmFraction_, Tau_EtaAtEcalEntrance_, Tau_GammaEnFrac_, Tau_GammaEtaMom_, Tau_GammaPhiMom_, Tau_GSFChi2_, Tau_GSFTrackEta_, Tau_GSFTracklnPt_, Tau_GSFTrackResol_, Tau_HadrEoP_, Tau_HadrHoP_, Tau_HadrMva_, Tau_NumGammaCands_, Tau_NumHitsVariable_, Tau_Pt_, Tau_VisMass_, and verbosity_.

Referenced by Initialize_from_file(), and Initialize_from_string().

135 {
136 
137  GammadEta_ = 0.;
138  GammadPhi_ = 0.;
139  GammadPt_ = 0.;
140 
141  Tau_EtaAtEcalEntrance_ = 0.;
142  Tau_Pt_ = 0.;
143  Tau_dCrackEta_ = 0.;
144  Tau_dCrackPhi_ = 0.;
145  Tau_EmFraction_ = 0.;
146  Tau_NumGammaCands_ = 0.;
147  Tau_HadrHoP_ = 0.;
148  Tau_HadrEoP_ = 0.;
149  Tau_VisMass_ = 0.;
150  Tau_HadrMva_ = 0.;
151  Tau_GammaEtaMom_ = 0.;
152  Tau_GammaPhiMom_ = 0.;
153  Tau_GammaEnFrac_ = 0.;
154  Tau_GSFChi2_ = 0.;
155  Tau_NumHitsVariable_ = 0.;
156  Tau_GSFTrackResol_ = 0.;
157  Tau_GSFTracklnPt_ = 0.;
158  Tau_GSFTrackEta_ = 0.;
159 
160  Elec_EtotOverPin_ = 0.;
161  Elec_EgammaOverPdif_ = 0.;
162  Elec_Fbrem_ = 0.;
163  Elec_Chi2GSF_ = 0.;
164  Elec_GSFNumHits_ = 0.;
165  Elec_GSFTrackResol_ = 0.;
166  Elec_GSFTracklnPt_ = 0.;
167  Elec_GSFTrackEta_ = 0.;
168 
169  //TMVA::Tools::Instance();
170 
171  TMVA::Reader *readerNoEleMatchwoGwoGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
172  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
173  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
174  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
175  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
176  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
177  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
178  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
179  readerNoEleMatchwoGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
180  readerNoEleMatchwoGwoGSF_BL->SetVerbose(verbosity_);
181 
182  TMVA::Reader *readerNoEleMatchwoGwGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
183  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
184  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
185  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
186  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
187  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
188  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
189  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
190  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
191  readerNoEleMatchwoGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
192  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
193  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
194  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
195  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
196  readerNoEleMatchwoGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
197  readerNoEleMatchwoGwGSF_BL->SetVerbose(verbosity_);
198 
199  TMVA::Reader *readerNoEleMatchwGwoGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
200  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
201  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
202  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
203  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
204  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
205  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
206  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
207  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
208  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
209  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
210  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
211  readerNoEleMatchwGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
212  readerNoEleMatchwGwoGSF_BL->SetVerbose(verbosity_);
213 
214  TMVA::Reader *readerNoEleMatchwGwGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
215  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
216  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
217  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
218  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
219  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
220  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
221  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
222  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
223  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
224  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
225  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
226  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
227  readerNoEleMatchwGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
228  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
229  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
230  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
231  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
232  readerNoEleMatchwGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
233  readerNoEleMatchwGwGSF_BL->SetVerbose(verbosity_);
234 
235  TMVA::Reader *readerwoGwoGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
236  readerwoGwoGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
237  readerwoGwoGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
238  readerwoGwoGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
239  readerwoGwoGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
240  readerwoGwoGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
241  readerwoGwoGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
242  readerwoGwoGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
243  readerwoGwoGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
244  readerwoGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
245  readerwoGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
246  readerwoGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
247  readerwoGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
248  readerwoGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
249  readerwoGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
250  readerwoGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
251  readerwoGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
252  readerwoGwoGSF_BL->SetVerbose(verbosity_);
253 
254  TMVA::Reader *readerwoGwGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
255  readerwoGwGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
256  readerwoGwGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
257  readerwoGwGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
258  readerwoGwGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
259  readerwoGwGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
260  readerwoGwGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
261  readerwoGwGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
262  readerwoGwGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
263  readerwoGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
264  readerwoGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
265  readerwoGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
266  readerwoGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
267  readerwoGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
268  readerwoGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
269  readerwoGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
270  readerwoGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
271  readerwoGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
272  readerwoGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
273  readerwoGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
274  readerwoGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
275  readerwoGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
276  readerwoGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
277  readerwoGwGSF_BL->SetVerbose(verbosity_);
278 
279  TMVA::Reader *readerwGwoGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
280  readerwGwoGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
281  readerwGwoGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
282  readerwGwoGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
283  readerwGwoGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
284  readerwGwoGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
285  readerwGwoGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
286  readerwGwoGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
287  readerwGwoGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
288  readerwGwoGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
289  readerwGwoGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
290  readerwGwoGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
291  readerwGwoGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
292  readerwGwoGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
293  readerwGwoGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
294  readerwGwoGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
295  readerwGwoGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
296  readerwGwoGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
297  readerwGwoGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
298  readerwGwoGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
299  readerwGwoGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
300  readerwGwoGSF_BL->SetVerbose(verbosity_);
301 
302  TMVA::Reader *readerwGwGSF_BL = new TMVA::Reader( "!Color:Silent:Error" );
303  readerwGwGSF_BL->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
304  readerwGwGSF_BL->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
305  readerwGwGSF_BL->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
306  readerwGwGSF_BL->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
307  readerwGwGSF_BL->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
308  readerwGwGSF_BL->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
309  readerwGwGSF_BL->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
310  readerwGwGSF_BL->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
311  readerwGwGSF_BL->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
312  readerwGwGSF_BL->AddVariable("Tau_Pt",&Tau_Pt_);
313  readerwGwGSF_BL->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
314  readerwGwGSF_BL->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
315  readerwGwGSF_BL->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
316  readerwGwGSF_BL->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
317  readerwGwGSF_BL->AddVariable("Tau_VisMass",&Tau_VisMass_);
318  readerwGwGSF_BL->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
319  readerwGwGSF_BL->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
320  readerwGwGSF_BL->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
321  readerwGwGSF_BL->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
322  readerwGwGSF_BL->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
323  readerwGwGSF_BL->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
324  readerwGwGSF_BL->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
325  readerwGwGSF_BL->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
326  readerwGwGSF_BL->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
327  readerwGwGSF_BL->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
328  readerwGwGSF_BL->AddVariable("Tau_dCrackPhi",&Tau_dCrackPhi_);
329  readerwGwGSF_BL->SetVerbose(verbosity_);
330 
332 
333  TMVA::Reader *readerNoEleMatchwoGwoGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
334  readerNoEleMatchwoGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
335  readerNoEleMatchwoGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
336  readerNoEleMatchwoGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
337  readerNoEleMatchwoGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
338  readerNoEleMatchwoGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
339  readerNoEleMatchwoGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
340  readerNoEleMatchwoGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
341  readerNoEleMatchwoGwoGSF_EC->SetVerbose(verbosity_);
342 
343  TMVA::Reader *readerNoEleMatchwoGwGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
344  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
345  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
346  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
347  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
348  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
349  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
350  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
351  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
352  readerNoEleMatchwoGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
353  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
354  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
355  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
356  readerNoEleMatchwoGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
357  readerNoEleMatchwoGwGSF_EC->SetVerbose(verbosity_);
358 
359  TMVA::Reader *readerNoEleMatchwGwoGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
360  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
361  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
362  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
363  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
364  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
365  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
366  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
367  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
368  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
369  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
370  readerNoEleMatchwGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
371  readerNoEleMatchwGwoGSF_EC->SetVerbose(verbosity_);
372 
373  TMVA::Reader *readerNoEleMatchwGwGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
374  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
375  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
376  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
377  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
378  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
379  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
380  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
381  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
382  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
383  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
384  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
385  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
386  readerNoEleMatchwGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
387  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
388  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
389  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
390  readerNoEleMatchwGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
391  readerNoEleMatchwGwGSF_EC->SetVerbose(verbosity_);
392 
393  TMVA::Reader *readerwoGwoGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
394  readerwoGwoGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
395  readerwoGwoGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
396  readerwoGwoGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
397  readerwoGwoGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
398  readerwoGwoGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
399  readerwoGwoGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
400  readerwoGwoGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
401  readerwoGwoGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
402  readerwoGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
403  readerwoGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
404  readerwoGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
405  readerwoGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
406  readerwoGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
407  readerwoGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
408  readerwoGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
409  readerwoGwoGSF_EC->SetVerbose(verbosity_);
410 
411  TMVA::Reader *readerwoGwGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
412  readerwoGwGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
413  readerwoGwGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
414  readerwoGwGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
415  readerwoGwGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
416  readerwoGwGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
417  readerwoGwGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
418  readerwoGwGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
419  readerwoGwGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
420  readerwoGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
421  readerwoGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
422  readerwoGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
423  readerwoGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
424  readerwoGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
425  readerwoGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
426  readerwoGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
427  readerwoGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
428  readerwoGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
429  readerwoGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
430  readerwoGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
431  readerwoGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
432  readerwoGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
433  readerwoGwGSF_EC->SetVerbose(verbosity_);
434 
435  TMVA::Reader *readerwGwoGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
436  readerwGwoGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
437  readerwGwoGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
438  readerwGwoGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
439  readerwGwoGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
440  readerwGwoGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
441  readerwGwoGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
442  readerwGwoGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
443  readerwGwoGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
444  readerwGwoGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
445  readerwGwoGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
446  readerwGwoGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
447  readerwGwoGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
448  readerwGwoGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
449  readerwGwoGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
450  readerwGwoGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
451  readerwGwoGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
452  readerwGwoGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
453  readerwGwoGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
454  readerwGwoGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
455  readerwGwoGSF_EC->SetVerbose(verbosity_);
456 
457  TMVA::Reader *readerwGwGSF_EC = new TMVA::Reader( "!Color:Silent:Error" );
458  readerwGwGSF_EC->AddVariable("Elec_EtotOverPin",&Elec_EtotOverPin_);
459  readerwGwGSF_EC->AddVariable("Elec_EgammaOverPdif",&Elec_EgammaOverPdif_);
460  readerwGwGSF_EC->AddVariable("Elec_Fbrem",&Elec_Fbrem_);
461  readerwGwGSF_EC->AddVariable("Elec_Chi2GSF",&Elec_Chi2GSF_);
462  readerwGwGSF_EC->AddVariable("Elec_GSFNumHits",&Elec_GSFNumHits_);
463  readerwGwGSF_EC->AddVariable("Elec_GSFTrackResol",&Elec_GSFTrackResol_);
464  readerwGwGSF_EC->AddVariable("Elec_GSFTracklnPt",&Elec_GSFTracklnPt_);
465  readerwGwGSF_EC->AddVariable("Elec_GSFTrackEta",&Elec_GSFTrackEta_);
466  readerwGwGSF_EC->AddVariable("Tau_EtaAtEcalEntrance",&Tau_EtaAtEcalEntrance_);
467  readerwGwGSF_EC->AddVariable("Tau_Pt",&Tau_Pt_);
468  readerwGwGSF_EC->AddVariable("Tau_EmFraction",&Tau_EmFraction_);
469  readerwGwGSF_EC->AddVariable("Tau_NumGammaCands",&Tau_NumGammaCands_);
470  readerwGwGSF_EC->AddVariable("Tau_HadrHoP",&Tau_HadrHoP_);
471  readerwGwGSF_EC->AddVariable("Tau_HadrEoP",&Tau_HadrEoP_);
472  readerwGwGSF_EC->AddVariable("Tau_VisMass",&Tau_VisMass_);
473  readerwGwGSF_EC->AddVariable("Tau_HadrMva",&Tau_HadrMva_);
474  readerwGwGSF_EC->AddVariable("Tau_GammaEtaMom",&Tau_GammaEtaMom_);
475  readerwGwGSF_EC->AddVariable("Tau_GammaPhiMom",&Tau_GammaPhiMom_);
476  readerwGwGSF_EC->AddVariable("Tau_GammaEnFrac",&Tau_GammaEnFrac_);
477  readerwGwGSF_EC->AddVariable("Tau_GSFChi2",&Tau_GSFChi2_);
478  readerwGwGSF_EC->AddVariable("(Tau_GSFNumHits - Tau_KFNumHits)/(Tau_GSFNumHits + Tau_KFNumHits)",&Tau_NumHitsVariable_);
479  readerwGwGSF_EC->AddVariable("Tau_GSFTrackResol",&Tau_GSFTrackResol_);
480  readerwGwGSF_EC->AddVariable("Tau_GSFTracklnPt",&Tau_GSFTracklnPt_);
481  readerwGwGSF_EC->AddVariable("Tau_GSFTrackEta",&Tau_GSFTrackEta_);
482  readerwGwGSF_EC->AddVariable("Tau_dCrackEta",&Tau_dCrackEta_);
483  readerwGwGSF_EC->SetVerbose(verbosity_);
484 
485  fTMVAReader_[k_NoEleMatchwoGwoGSF_BL] = readerNoEleMatchwoGwoGSF_BL;
486  fTMVAReader_[k_NoEleMatchwoGwGSF_BL] = readerNoEleMatchwoGwGSF_BL;
487  fTMVAReader_[k_NoEleMatchwGwoGSF_BL] = readerNoEleMatchwGwoGSF_BL;
488  fTMVAReader_[k_NoEleMatchwGwGSF_BL] = readerNoEleMatchwGwGSF_BL;
489  fTMVAReader_[k_woGwoGSF_BL] = readerwoGwoGSF_BL;
490  fTMVAReader_[k_woGwGSF_BL] = readerwoGwGSF_BL;
491  fTMVAReader_[k_wGwoGSF_BL] = readerwGwoGSF_BL;
492  fTMVAReader_[k_wGwGSF_BL] = readerwGwGSF_BL;
493  fTMVAReader_[k_NoEleMatchwoGwoGSF_EC] = readerNoEleMatchwoGwoGSF_EC;
494  fTMVAReader_[k_NoEleMatchwoGwGSF_EC] = readerNoEleMatchwoGwGSF_EC;
495  fTMVAReader_[k_NoEleMatchwGwoGSF_EC] = readerNoEleMatchwGwoGSF_EC;
496  fTMVAReader_[k_NoEleMatchwGwGSF_EC] = readerNoEleMatchwGwGSF_EC;
497  fTMVAReader_[k_woGwoGSF_EC] = readerwoGwoGSF_EC;
498  fTMVAReader_[k_woGwGSF_EC] = readerwoGwGSF_EC;
499  fTMVAReader_[k_wGwoGSF_EC] = readerwGwoGSF_EC;
500  fTMVAReader_[k_wGwGSF_EC] = readerwGwGSF_EC;
501 
502 }
AntiElectronIDMVA3::fTMVAReader_
TMVA::Reader * fTMVAReader_[16]
Definition: AntiElectronIDMVA3.h:156
double AntiElectronIDMVA3::dCrackEta ( double  eta)
private

Definition at line 1035 of file AntiElectronIDMVA3.cc.

References i, m, and minimum().

Referenced by MVAValue().

1035  {
1036 
1037  //Location of the eta-cracks
1038  double cracks[5] = {0, 4.44747e-01, 7.92824e-01, 1.14090e+00, 1.47464e+00};
1039 
1040  double m=99.; //the result
1041 
1042  for(int i=0;i<5;i++){
1043  double d = minimum(eta-cracks[i], eta+cracks[i]);
1044  if (TMath::Abs(d)<TMath::Abs(m)){
1045  m=d;
1046  }
1047  }
1048 
1049  return TMath::Abs(m);
1050 }
i
int i
Definition: DBlmapReader.cc:9
eta
T eta() const
Definition: Basic3DVectorLD.h:172
AntiElectronIDMVA3::minimum
double minimum(double a, double b)
Definition: AntiElectronIDMVA3.cc:969
double AntiElectronIDMVA3::dCrackPhi ( double  phi,
double  eta 
)
private

Definition at line 976 of file AntiElectronIDMVA3.cc.

References i, m, minimum(), Pi, and pi.

Referenced by MVAValue().

976  {
977 
978  double pi= TMath::Pi();// 3.14159265358979323846;
979 
980  //Location of the 18 phi-cracks
981  static std::vector<double> cPhi;
982  if(cPhi.size()==0)
983  {
984  cPhi.resize(18,0);
985  cPhi[0]=2.97025;
986  for(unsigned i=1;i<=17;++i) cPhi[i]=cPhi[0]-2*i*pi/18;
987  }
988 
989  //Shift of this location if eta<0
990  double delta_cPhi=0.00638;
991 
992  double m; //the result
993 
994  if (eta>=- 1.47464 && eta<= 1.47464){
995 
996  //the location is shifted
997  if(eta<0) phi +=delta_cPhi;
998 
999  if (phi>=-pi && phi<=pi){
1000 
1001  //the problem of the extrema
1002  if (phi<cPhi[17] || phi>=cPhi[0]){
1003  if (phi<0) phi+= 2*pi;
1004  m = minimum(phi -cPhi[0],phi-cPhi[17]-2*pi);
1005  }
1006 
1007  //between these extrema...
1008  else{
1009  bool OK = false;
1010  unsigned i=16;
1011  while(!OK){
1012  if (phi<cPhi[i]){
1013  m=minimum(phi-cPhi[i+1],phi-cPhi[i]);
1014  OK=true;
1015  }
1016  else i-=1;
1017  }
1018  }
1019  }
1020  else{
1021  m=0.; //if there is a problem, we assum that we are in a crack
1022  //std::cout<<"Problem in dminphi"<<std::endl;
1023  }
1024  }
1025  else{
1026  return -99.;
1027  //std::cout<<"Encap region"<<std::endl;
1028  }
1029 
1030  return TMath::Abs(m);
1031 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
i
int i
Definition: DBlmapReader.cc:9
eta
T eta() const
Definition: Basic3DVectorLD.h:172
pi
const Double_t pi
Definition: trackSplitPlot.h:35
AntiElectronIDMVA3::minimum
double minimum(double a, double b)
Definition: AntiElectronIDMVA3.cc:969
phi
Definition: DDAxes.h:10
void AntiElectronIDMVA3::Initialize_from_file ( const std::string &  methodName,
const std::string &  oneProngNoEleMatch0Pi0woGSF_BL,
const std::string &  oneProngNoEleMatch0Pi0wGSF_BL,
const std::string &  oneProngNoEleMatch1Pi0woGSF_BL,
const std::string &  oneProngNoEleMatch1Pi0wGSF_BL,
const std::string &  oneProng0Pi0woGSF_BL,
const std::string &  oneProng0Pi0wGSF_BL,
const std::string &  oneProng1Pi0woGSF_BL,
const std::string &  oneProng1Pi0wGSF_BL,
const std::string &  oneProngNoEleMatch0Pi0woGSF_EC,
const std::string &  oneProngNoEleMatch0Pi0wGSF_EC,
const std::string &  oneProngNoEleMatch1Pi0woGSF_EC,
const std::string &  oneProngNoEleMatch1Pi0wGSF_EC,
const std::string &  oneProng0Pi0woGSF_EC,
const std::string &  oneProng0Pi0wGSF_EC,
const std::string &  oneProng1Pi0woGSF_EC,
const std::string &  oneProng1Pi0wGSF_EC 
)

Definition at line 42 of file AntiElectronIDMVA3.cc.

References bookMVAs(), fTMVAReader_, i, isInitialized_, k_NoEleMatchwGwGSF_BL, k_NoEleMatchwGwGSF_EC, k_NoEleMatchwGwoGSF_BL, k_NoEleMatchwGwoGSF_EC, k_NoEleMatchwoGwGSF_BL, k_NoEleMatchwoGwGSF_EC, k_NoEleMatchwoGwoGSF_BL, k_NoEleMatchwoGwoGSF_EC, k_wGwGSF_BL, k_wGwGSF_EC, k_wGwoGSF_BL, k_wGwoGSF_EC, k_woGwGSF_BL, k_woGwGSF_EC, k_woGwoGSF_BL, k_woGwoGSF_EC, reco::details::loadTMVAWeights(), and methodName_.

59 {
60  for ( unsigned i = 0; i < 16; ++i ) {
61  if ( fTMVAReader_[i] ) delete fTMVAReader_[i];
62  }
63 
64  isInitialized_ = kTRUE;
65  methodName_ = methodName;
66 
67  bookMVAs();
68 
69  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwoGwoGSF_BL], "BDTG", oneProngNoEleMatch0Pi0woGSF_BL);
70  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwoGwGSF_BL], "BDTG", oneProngNoEleMatch0Pi0wGSF_BL);
71  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwGwoGSF_BL], "BDTG", oneProngNoEleMatch1Pi0woGSF_BL);
72  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwGwGSF_BL], "BDTG", oneProngNoEleMatch1Pi0wGSF_BL);
73  reco::details::loadTMVAWeights(fTMVAReader_[k_woGwoGSF_BL], "BDTG", oneProng0Pi0woGSF_BL);
74  reco::details::loadTMVAWeights(fTMVAReader_[k_woGwGSF_BL], "BDTG", oneProng0Pi0wGSF_BL);
75  reco::details::loadTMVAWeights(fTMVAReader_[k_wGwoGSF_BL], "BDTG", oneProng1Pi0woGSF_BL);
76  reco::details::loadTMVAWeights(fTMVAReader_[k_wGwGSF_BL], "BDTG", oneProng1Pi0wGSF_BL);
77  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwoGwoGSF_EC], "BDTG", oneProngNoEleMatch0Pi0woGSF_EC);
78  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwoGwGSF_EC], "BDTG", oneProngNoEleMatch0Pi0wGSF_EC);
79  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwGwoGSF_EC], "BDTG", oneProngNoEleMatch1Pi0woGSF_EC);
80  reco::details::loadTMVAWeights(fTMVAReader_[k_NoEleMatchwGwGSF_EC], "BDTG", oneProngNoEleMatch1Pi0wGSF_EC);
81  reco::details::loadTMVAWeights(fTMVAReader_[k_woGwoGSF_EC], "BDTG", oneProng0Pi0woGSF_EC);
82  reco::details::loadTMVAWeights(fTMVAReader_[k_woGwGSF_EC], "BDTG", oneProng0Pi0wGSF_EC);
83  reco::details::loadTMVAWeights(fTMVAReader_[k_wGwoGSF_EC], "BDTG", oneProng1Pi0woGSF_EC);
84  reco::details::loadTMVAWeights(fTMVAReader_[k_wGwGSF_EC], "BDTG", oneProng1Pi0wGSF_EC);
85 
86 }
i
int i
Definition: DBlmapReader.cc:9
AntiElectronIDMVA3::fTMVAReader_
TMVA::Reader * fTMVAReader_[16]
Definition: AntiElectronIDMVA3.h:156
reco::details::loadTMVAWeights
void loadTMVAWeights(TMVA::Reader *reader, const std::string &method, const std::string &weightFile, bool verbose=false)
Definition: TMVAZipReader.cc:50
void AntiElectronIDMVA3::Initialize_from_string ( const std::string &  methodName,
const std::string &  oneProngNoEleMatch0Pi0woGSF_BL,
const std::string &  oneProngNoEleMatch0Pi0wGSF_BL,
const std::string &  oneProngNoEleMatch1Pi0woGSF_BL,
const std::string &  oneProngNoEleMatch1Pi0wGSF_BL,
const std::string &  oneProng0Pi0woGSF_BL,
const std::string &  oneProng0Pi0wGSF_BL,
const std::string &  oneProng1Pi0woGSF_BL,
const std::string &  oneProng1Pi0wGSF_BL,
const std::string &  oneProngNoEleMatch0Pi0woGSF_EC,
const std::string &  oneProngNoEleMatch0Pi0wGSF_EC,
const std::string &  oneProngNoEleMatch1Pi0woGSF_EC,
const std::string &  oneProngNoEleMatch1Pi0wGSF_EC,
const std::string &  oneProng0Pi0woGSF_EC,
const std::string &  oneProng0Pi0wGSF_EC,
const std::string &  oneProng1Pi0woGSF_EC,
const std::string &  oneProng1Pi0wGSF_EC 
)

Definition at line 88 of file AntiElectronIDMVA3.cc.

References bookMVAs(), fTMVAReader_, i, isInitialized_, k_NoEleMatchwGwGSF_BL, k_NoEleMatchwGwGSF_EC, k_NoEleMatchwGwoGSF_BL, k_NoEleMatchwGwoGSF_EC, k_NoEleMatchwoGwGSF_BL, k_NoEleMatchwoGwGSF_EC, k_NoEleMatchwoGwoGSF_BL, k_NoEleMatchwoGwoGSF_EC, k_wGwGSF_BL, k_wGwGSF_EC, k_wGwoGSF_BL, k_wGwoGSF_EC, k_woGwGSF_BL, k_woGwGSF_EC, k_woGwoGSF_BL, k_woGwoGSF_EC, and methodName_.

105 {
106  for ( unsigned i = 0; i < 16; ++i ) {
107  if ( fTMVAReader_[i] ) delete fTMVAReader_[i];
108  }
109 
110  isInitialized_ = kTRUE;
111  methodName_ = methodName;
112 
113  bookMVAs();
114 
115  fTMVAReader_[k_NoEleMatchwoGwoGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch0Pi0woGSF_BL.data());
116  fTMVAReader_[k_NoEleMatchwoGwGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch0Pi0wGSF_BL.data());
117  fTMVAReader_[k_NoEleMatchwGwoGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch1Pi0woGSF_BL.data());
118  fTMVAReader_[k_NoEleMatchwGwGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch1Pi0wGSF_BL.data());
119  fTMVAReader_[k_woGwoGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProng0Pi0woGSF_BL.data());
120  fTMVAReader_[k_woGwGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProng0Pi0wGSF_BL.data());
121  fTMVAReader_[k_wGwoGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProng1Pi0woGSF_BL.data());
122  fTMVAReader_[k_wGwGSF_BL]->BookMVA(TMVA::Types::kBDT, oneProng1Pi0wGSF_BL.data());
123  fTMVAReader_[k_NoEleMatchwoGwoGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch0Pi0woGSF_EC.data());
124  fTMVAReader_[k_NoEleMatchwoGwGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch0Pi0wGSF_EC.data());
125  fTMVAReader_[k_NoEleMatchwGwoGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch1Pi0woGSF_EC.data());
126  fTMVAReader_[k_NoEleMatchwGwGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProngNoEleMatch1Pi0wGSF_EC.data());
127  fTMVAReader_[k_woGwoGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProng0Pi0woGSF_EC.data());
128  fTMVAReader_[k_woGwGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProng0Pi0wGSF_EC.data());
129  fTMVAReader_[k_wGwoGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProng1Pi0woGSF_EC.data());
130  fTMVAReader_[k_wGwGSF_EC]->BookMVA(TMVA::Types::kBDT, oneProng1Pi0wGSF_EC.data());
131 
132 }
i
int i
Definition: DBlmapReader.cc:9
AntiElectronIDMVA3::fTMVAReader_
TMVA::Reader * fTMVAReader_[16]
Definition: AntiElectronIDMVA3.h:156
double AntiElectronIDMVA3::minimum ( double  a,
double  b 
)
private

Definition at line 969 of file AntiElectronIDMVA3.cc.

References a, and b.

Referenced by dCrackEta(), and dCrackPhi().

969  {
970  if(TMath::Abs(b)<TMath::Abs(a)) return b;
971  else return a;
972 }
b
double b
Definition: hdecay.h:120
a
double a
Definition: hdecay.h:121
double AntiElectronIDMVA3::MVAValue ( Float_t  TauEtaAtEcalEntrance,
Float_t  TauPt,
Float_t  TaudCrackEta,
Float_t  TaudCrackPhi,
Float_t  TauEmFraction,
Float_t  TauSignalPFGammaCands,
Float_t  TauLeadPFChargedHadrHoP,
Float_t  TauLeadPFChargedHadrEoP,
Float_t  TauVisMass,
Float_t  TauHadrMva,
const std::vector< Float_t > &  GammasdEta,
const std::vector< Float_t > &  GammasdPhi,
const std::vector< Float_t > &  GammasPt,
Float_t  TauKFNumHits,
Float_t  TauGSFNumHits,
Float_t  TauGSFChi2,
Float_t  TauGSFTrackResol,
Float_t  TauGSFTracklnPt,
Float_t  TauGSFTrackEta,
Float_t  TauPhi,
Float_t  TauSignalPFChargedCands,
Float_t  TauHasGsf,
Float_t  ElecEta,
Float_t  ElecPhi,
Float_t  ElecPt,
Float_t  ElecEe,
Float_t  ElecEgamma,
Float_t  ElecPin,
Float_t  ElecPout,
Float_t  ElecFbrem,
Float_t  ElecChi2GSF,
Float_t  ElecGSFNumHits,
Float_t  ElecGSFTrackResol,
Float_t  ElecGSFTracklnPt,
Float_t  ElecGSFTrackEta 
)

Definition at line 504 of file AntiElectronIDMVA3.cc.

References dPhi(), i, and Pi.

Referenced by MVAValue().

539 {
540  double sumPt = 0.;
541  double dEta = 0.;
542  double dEta2 = 0.;
543  double dPhi = 0.;
544  double dPhi2 = 0.;
545  double sumPt2 = 0.;
546  for ( unsigned int i = 0 ; i < GammasPt.size() ; ++i ) {
547  double pt_i = GammasPt[i];
548  double phi_i = GammasdPhi[i];
549  if ( GammasdPhi[i] > TMath::Pi() ) phi_i = GammasdPhi[i] - 2*TMath::Pi();
550  else if ( GammasdPhi[i] < -TMath::Pi() ) phi_i = GammasdPhi[i] + 2*TMath::Pi();
551  double eta_i = GammasdEta[i];
552  sumPt += pt_i;
553  sumPt2 += (pt_i*pt_i);
554  dEta += (pt_i*eta_i);
555  dEta2 += (pt_i*eta_i*eta_i);
556  dPhi += (pt_i*phi_i);
557  dPhi2 += (pt_i*phi_i*phi_i);
558  }
559 
560  Float_t TauGammaEnFrac = sumPt/TauPt;
561 
562  if ( sumPt > 0. ) {
563  dEta /= sumPt;
564  dPhi /= sumPt;
565  dEta2 /= sumPt;
566  dPhi2 /= sumPt;
567  }
568 
569  Float_t TauGammaEtaMom = TMath::Sqrt(dEta2)*TMath::Sqrt(TauGammaEnFrac)*TauPt;
570  Float_t TauGammaPhiMom = TMath::Sqrt(dPhi2)*TMath::Sqrt(TauGammaEnFrac)*TauPt;
571 
572  return MVAValue(TauEtaAtEcalEntrance,
573  TauPt,
574  TaudCrackEta,
575  TaudCrackPhi,
576  TauEmFraction,
577  TauSignalPFGammaCands,
578  TauLeadPFChargedHadrHoP,
579  TauLeadPFChargedHadrEoP,
580  TauVisMass,
581  TauHadrMva,
582  TauGammaEtaMom,
583  TauGammaPhiMom,
584  TauGammaEnFrac,
585  TauKFNumHits,
586  TauGSFNumHits,
587  TauGSFChi2,
588  TauGSFTrackResol,
589  TauGSFTracklnPt,
590  TauGSFTrackEta,
591  TauPhi,
592  TauSignalPFChargedCands,
593  TauHasGsf,
594  ElecEta,
595  ElecPhi,
596  ElecPt,
597  ElecEe,
598  ElecEgamma,
599  ElecPin,
600  ElecPout,
601  ElecFbrem,
602  ElecChi2GSF,
603  ElecGSFNumHits,
604  ElecGSFTrackResol,
605  ElecGSFTracklnPt,
606  ElecGSFTrackEta);
607 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
i
int i
Definition: DBlmapReader.cc:9
dPhi
double dPhi(double phi1, double phi2)
Definition: JetUtil.h:30
AntiElectronIDMVA3::MVAValue
double MVAValue(Float_t TauEtaAtEcalEntrance, Float_t TauPt, Float_t TaudCrackEta, Float_t TaudCrackPhi, Float_t TauEmFraction, Float_t TauSignalPFGammaCands, Float_t TauLeadPFChargedHadrHoP, Float_t TauLeadPFChargedHadrEoP, Float_t TauVisMass, Float_t TauHadrMva, const std::vector< Float_t > &GammasdEta, const std::vector< Float_t > &GammasdPhi, const std::vector< Float_t > &GammasPt, Float_t TauKFNumHits, Float_t TauGSFNumHits, Float_t TauGSFChi2, Float_t TauGSFTrackResol, Float_t TauGSFTracklnPt, Float_t TauGSFTrackEta, Float_t TauPhi, Float_t TauSignalPFChargedCands, Float_t TauHasGsf, Float_t ElecEta, Float_t ElecPhi, Float_t ElecPt, Float_t ElecEe, Float_t ElecEgamma, Float_t ElecPin, Float_t ElecPout, Float_t ElecFbrem, Float_t ElecChi2GSF, Float_t ElecGSFNumHits, Float_t ElecGSFTrackResol, Float_t ElecGSFTracklnPt, Float_t ElecGSFTrackEta)
Definition: AntiElectronIDMVA3.cc:504
double AntiElectronIDMVA3::MVAValue ( Float_t  TauEtaAtEcalEntrance,
Float_t  TauPt,
Float_t  TaudCrackEta,
Float_t  TaudCrackPhi,
Float_t  TauEmFraction,
Float_t  TauSignalPFGammaCands,
Float_t  TauLeadPFChargedHadrHoP,
Float_t  TauLeadPFChargedHadrEoP,
Float_t  TauVisMass,
Float_t  TauHadrMva,
Float_t  TauGammaEtaMom,
Float_t  TauGammaPhiMom,
Float_t  TauGammaEnFrac,
Float_t  TauKFNumHits,
Float_t  TauGSFNumHits,
Float_t  TauGSFChi2,
Float_t  TauGSFTrackResol,
Float_t  TauGSFTracklnPt,
Float_t  TauGSFTrackEta,
Float_t  TauPhi,
Float_t  TauSignalPFChargedCands,
Float_t  TauHasGsf,
Float_t  ElecEta,
Float_t  ElecPhi,
Float_t  ElecPt,
Float_t  ElecEe,
Float_t  ElecEgamma,
Float_t  ElecPin,
Float_t  ElecPout,
Float_t  ElecFbrem,
Float_t  ElecChi2GSF,
Float_t  ElecGSFNumHits,
Float_t  ElecGSFTrackResol,
Float_t  ElecGSFTracklnPt,
Float_t  ElecGSFTrackEta 
)

Definition at line 609 of file AntiElectronIDMVA3.cc.

References deltaR(), Elec_Chi2GSF_, Elec_EgammaOverPdif_, Elec_EtotOverPin_, Elec_Fbrem_, Elec_GSFNumHits_, Elec_GSFTrackEta_, Elec_GSFTracklnPt_, Elec_GSFTrackResol_, fTMVAReader_, isInitialized_, k_NoEleMatchwGwGSF_BL, k_NoEleMatchwGwGSF_EC, k_NoEleMatchwGwoGSF_BL, k_NoEleMatchwGwoGSF_EC, k_NoEleMatchwoGwGSF_BL, k_NoEleMatchwoGwGSF_EC, k_NoEleMatchwoGwoGSF_BL, k_NoEleMatchwoGwoGSF_EC, k_wGwGSF_BL, k_wGwGSF_EC, k_wGwoGSF_BL, k_wGwoGSF_EC, k_woGwGSF_BL, k_woGwGSF_EC, k_woGwoGSF_BL, k_woGwoGSF_EC, siStripFEDMonitor_P5_cff::Max, methodName_, Tau_dCrackEta_, Tau_dCrackPhi_, Tau_EmFraction_, Tau_EtaAtEcalEntrance_, Tau_GammaEnFrac_, Tau_GammaEtaMom_, Tau_GammaPhiMom_, Tau_GSFChi2_, Tau_GSFTrackEta_, Tau_GSFTracklnPt_, Tau_GSFTrackResol_, Tau_HadrEoP_, Tau_HadrHoP_, Tau_HadrMva_, Tau_NumGammaCands_, Tau_NumHitsVariable_, Tau_Pt_, and Tau_VisMass_.

644 {
645 
646  if ( !isInitialized_ ) {
647  //std::cout << "Error: AntiElectronMVA not properly initialized.\n";
648  return -99.;
649  }
650 
651  Tau_EtaAtEcalEntrance_ = TauEtaAtEcalEntrance;
652  Tau_Pt_ = TauPt;
653  Tau_dCrackEta_ = TaudCrackEta;
654  Tau_dCrackPhi_ = TaudCrackPhi;
655  Tau_EmFraction_ = TMath::Max(TauEmFraction,float(0.0));
656  Tau_NumGammaCands_ = TauSignalPFGammaCands;
657  Tau_HadrHoP_ = TauLeadPFChargedHadrHoP;
658  Tau_HadrEoP_ = TauLeadPFChargedHadrEoP;
659  Tau_VisMass_ = TauVisMass;
660  Tau_HadrMva_ = TauHadrMva;
661  Tau_GammaEtaMom_ = TauGammaEtaMom;
662  Tau_GammaPhiMom_ = TauGammaPhiMom;
663  Tau_GammaEnFrac_ = TauGammaEnFrac;
664  Tau_GSFChi2_ = TauGSFChi2;
665  Tau_NumHitsVariable_ = (TauGSFNumHits-TauKFNumHits)/(TauGSFNumHits+TauKFNumHits);
666  Tau_GSFTrackResol_ = TauGSFTrackResol;
667  Tau_GSFTracklnPt_ = TauGSFTracklnPt;
668  Tau_GSFTrackEta_ = TauGSFTrackEta;
669 
670  Elec_EtotOverPin_ = (ElecEe + ElecEgamma)/ElecPin;
671  Elec_EgammaOverPdif_ = ElecEgamma/(ElecPin - ElecPout);
672  Elec_Fbrem_ = ElecFbrem;
673  Elec_Chi2GSF_ = ElecChi2GSF;
674  Elec_GSFNumHits_ = ElecGSFNumHits;
675  Elec_GSFTrackResol_ = ElecGSFTrackResol;
676  Elec_GSFTracklnPt_ = ElecGSFTracklnPt;
677  Elec_GSFTrackEta_ = ElecGSFTrackEta;
678 
679  double mva = -99.;
680  if ( TauSignalPFChargedCands == 3 ) mva = 1.0;
681  else if ( deltaR(TauEtaAtEcalEntrance, TauPhi, ElecEta, ElecPhi) > 0.3 && TauSignalPFGammaCands == 0 && TauHasGsf < 0.5) {
682  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_NoEleMatchwoGwoGSF_BL]->EvaluateMVA(methodName_);
683  else mva = fTMVAReader_[k_NoEleMatchwoGwoGSF_EC]->EvaluateMVA(methodName_);
684  }
685  else if ( deltaR(TauEtaAtEcalEntrance, TauPhi, ElecEta, ElecPhi) > 0.3 && TauSignalPFGammaCands == 0 && TauHasGsf > 0.5) {
686  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_NoEleMatchwoGwGSF_BL]->EvaluateMVA(methodName_);
687  else mva = fTMVAReader_[k_NoEleMatchwoGwGSF_EC]->EvaluateMVA(methodName_);
688  }
689  else if ( deltaR(TauEtaAtEcalEntrance, TauPhi, ElecEta, ElecPhi) > 0.3 && TauSignalPFGammaCands > 0 && TauHasGsf < 0.5) {
690  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_NoEleMatchwGwoGSF_BL]->EvaluateMVA(methodName_);
691  else mva = fTMVAReader_[k_NoEleMatchwGwoGSF_EC]->EvaluateMVA(methodName_);
692  }
693  else if ( deltaR(TauEtaAtEcalEntrance, TauPhi, ElecEta, ElecPhi) > 0.3 && TauSignalPFGammaCands > 0 && TauHasGsf > 0.5) {
694  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_NoEleMatchwGwGSF_BL]->EvaluateMVA(methodName_);
695  else mva = fTMVAReader_[k_NoEleMatchwGwGSF_EC]->EvaluateMVA(methodName_);
696  }
697  else if ( TauSignalPFGammaCands == 0 && TauHasGsf < 0.5) {
698  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_woGwoGSF_BL]->EvaluateMVA(methodName_);
699  else mva = fTMVAReader_[k_woGwoGSF_EC]->EvaluateMVA(methodName_);
700  }
701  else if ( TauSignalPFGammaCands == 0 && TauHasGsf > 0.5) {
702  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_woGwGSF_BL]->EvaluateMVA(methodName_);
703  else mva = fTMVAReader_[k_woGwGSF_EC]->EvaluateMVA(methodName_);
704  }
705  else if ( TauSignalPFGammaCands > 0 && TauHasGsf < 0.5) {
706  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_wGwoGSF_BL]->EvaluateMVA(methodName_);
707  else mva = fTMVAReader_[k_wGwoGSF_EC]->EvaluateMVA(methodName_);
708  }
709  else if ( TauSignalPFGammaCands > 0 && TauHasGsf > 0.5) {
710  if ( TMath::Abs(TauEtaAtEcalEntrance) < 1.5 ) mva = fTMVAReader_[k_wGwGSF_BL]->EvaluateMVA(methodName_);
711  else mva = fTMVAReader_[k_wGwGSF_EC]->EvaluateMVA(methodName_);
712  }
713  return mva;
714 }
siStripFEDMonitor_P5_cff.Max
tuple Max
Definition: siStripFEDMonitor_P5_cff.py:56
deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
AntiElectronIDMVA3::fTMVAReader_
TMVA::Reader * fTMVAReader_[16]
Definition: AntiElectronIDMVA3.h:156
double AntiElectronIDMVA3::MVAValue ( const reco::PFTau &  thePFTau,
const reco::GsfElectron &  theGsfEle 
)

Definition at line 717 of file AntiElectronIDMVA3.cc.

References edm::RefVector< C, T, F >::at(), asciidump::at, dCrackEta(), dCrackPhi(), reco::PFTau::electronPreIDOutput(), reco::PFTau::emFraction(), reco::LeafCandidate::eta(), reco::GsfElectron::fbrem(), reco::GsfElectron::gsfTrack(), i, edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), j, reco::PFTau::leadPFChargedHadrCand(), create_public_lumi_plots::log, reco::LeafCandidate::mass(), siStripFEDMonitor_P5_cff::Max, MVAValue(), reco::GsfElectron::pflowSuperCluster(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::PFTau::signalPFCands(), reco::PFTau::signalPFChargedHadrCands(), reco::PFTau::signalPFGammaCands(), edm::RefVector< C, T, F >::size(), findQualityFiles::size, reco::GsfElectron::trackMomentumAtVtx(), and reco::GsfElectron::trackMomentumOut().

720 {
721  Float_t TauEtaAtEcalEntrance = -99.;
722  float sumEtaTimesEnergy = 0;
723  float sumEnergy = 0;
724  for(unsigned int j = 0 ; j < (thePFTau.signalPFCands()).size() ; j++){
725  reco::PFCandidateRef pfcandidate = (thePFTau.signalPFCands()).at(j);
726  sumEtaTimesEnergy += pfcandidate->positionAtECALEntrance().eta()*pfcandidate->energy();
727  sumEnergy += pfcandidate->energy();
728  }
729  if(sumEnergy>0)TauEtaAtEcalEntrance = sumEtaTimesEnergy/sumEnergy;
730 
731  Float_t TauPt = thePFTau.pt();
732  Float_t TauEmFraction = TMath::Max(thePFTau.emFraction(), (Float_t)0.);
733  Float_t TauSignalPFGammaCands = thePFTau.signalPFGammaCands().size();
734  Float_t TauLeadPFChargedHadrHoP = 0.;
735  Float_t TauLeadPFChargedHadrEoP = 0.;
736  if ( thePFTau.leadPFChargedHadrCand()->p() > 0. ) {
737  TauLeadPFChargedHadrHoP = thePFTau.leadPFChargedHadrCand()->hcalEnergy()/thePFTau.leadPFChargedHadrCand()->p();
738  TauLeadPFChargedHadrEoP = thePFTau.leadPFChargedHadrCand()->ecalEnergy()/thePFTau.leadPFChargedHadrCand()->p();
739  }
740  Float_t TauVisMass = thePFTau.mass();
741  Float_t TauHadrMva = TMath::Max(thePFTau.electronPreIDOutput(),float(-1.0));
742  std::vector<Float_t> GammasdEta;
743  std::vector<Float_t> GammasdPhi;
744  std::vector<Float_t> GammasPt;
745  for ( unsigned i = 0 ; i < thePFTau.signalPFGammaCands().size(); ++i ) {
746  reco::PFCandidateRef gamma = thePFTau.signalPFGammaCands().at(i);
747  if ( thePFTau.leadPFChargedHadrCand().isNonnull() ) {
748  GammasdEta.push_back(gamma->eta() - thePFTau.leadPFChargedHadrCand()->eta());
749  GammasdPhi.push_back(gamma->phi() - thePFTau.leadPFChargedHadrCand()->phi());
750  } else {
751  GammasdEta.push_back(gamma->eta() - thePFTau.eta());
752  GammasdPhi.push_back(gamma->phi() - thePFTau.phi());
753  }
754  GammasPt.push_back(gamma->pt());
755  }
756  Float_t TauKFNumHits = -99.;
757  if((thePFTau.leadPFChargedHadrCand()->trackRef()).isNonnull()){
758  TauKFNumHits = thePFTau.leadPFChargedHadrCand()->trackRef()->numberOfValidHits();
759  }
760  Float_t TauGSFNumHits = -99.;
761  Float_t TauGSFChi2 = -99.;
762  Float_t TauGSFTrackResol = -99.;
763  Float_t TauGSFTracklnPt = -99.;
764  Float_t TauGSFTrackEta = -99.;
765  if((thePFTau.leadPFChargedHadrCand()->gsfTrackRef()).isNonnull()){
766  TauGSFChi2 = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->normalizedChi2();
767  TauGSFNumHits = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->numberOfValidHits();
768  if ( thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->pt() > 0. ) {
769  TauGSFTrackResol = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->ptError()/thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->pt();
770  TauGSFTracklnPt = log(thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->pt())*TMath::Ln10();
771  }
772  TauGSFTrackEta = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->eta();
773  }
774  Float_t TauPhi = thePFTau.phi();
775  float sumPhiTimesEnergy = 0;
776  float sumEnergyPhi = 0;
777  for(unsigned int j = 0 ; j < (thePFTau.signalPFCands()).size() ; j++){
778  reco::PFCandidateRef pfcandidate = (thePFTau.signalPFCands()).at(j);
779  sumPhiTimesEnergy += pfcandidate->positionAtECALEntrance().phi()*pfcandidate->energy();
780  sumEnergyPhi += pfcandidate->energy();
781  }
782  if(sumEnergy>0)TauPhi = sumPhiTimesEnergy/sumEnergyPhi;
783  Float_t TaudCrackPhi = dCrackPhi(TauPhi,TauEtaAtEcalEntrance) ;
784  Float_t TaudCrackEta = dCrackEta(TauEtaAtEcalEntrance) ;
785  Float_t TauSignalPFChargedCands = thePFTau.signalPFChargedHadrCands().size();
786  Float_t TauHasGsf = thePFTau.leadPFChargedHadrCand()->gsfTrackRef().isNonnull();
787 
788  Float_t ElecEta = theGsfEle.eta();
789  Float_t ElecPhi = theGsfEle.phi();
790  Float_t ElecPt = theGsfEle.pt();
791  //Variables related to the electron Cluster
792  Float_t ElecEe = 0.;
793  Float_t ElecEgamma = 0.;
794  reco::SuperClusterRef pfSuperCluster = theGsfEle.pflowSuperCluster();
795  if ( pfSuperCluster.isNonnull() && pfSuperCluster.isAvailable() ) {
796  for ( reco::CaloCluster_iterator pfCluster = pfSuperCluster->clustersBegin();
797  pfCluster != pfSuperCluster->clustersEnd(); ++pfCluster ) {
798  double pfClusterEn = (*pfCluster)->energy();
799  if ( pfCluster == pfSuperCluster->clustersBegin() ) ElecEe += pfClusterEn;
800  else ElecEgamma += pfClusterEn;
801  }
802  }
803  Float_t ElecPin = TMath::Sqrt(theGsfEle.trackMomentumAtVtx().Mag2());
804  Float_t ElecPout = TMath::Sqrt(theGsfEle.trackMomentumOut().Mag2());
805  Float_t ElecFbrem = theGsfEle.fbrem();
806  //Variables related to the GsfTrack
807  Float_t ElecChi2GSF = -99.;
808  Float_t ElecGSFNumHits = -99.;
809  Float_t ElecGSFTrackResol = -99.;
810  Float_t ElecGSFTracklnPt = -99.;
811  Float_t ElecGSFTrackEta = -99.;
812  if ( theGsfEle.gsfTrack().isNonnull() ) {
813  ElecChi2GSF = (theGsfEle).gsfTrack()->normalizedChi2();
814  ElecGSFNumHits = (theGsfEle).gsfTrack()->numberOfValidHits();
815  if ( theGsfEle.gsfTrack()->pt() > 0. ) {
816  ElecGSFTrackResol = theGsfEle.gsfTrack()->ptError()/theGsfEle.gsfTrack()->pt();
817  ElecGSFTracklnPt = log(theGsfEle.gsfTrack()->pt())*TMath::Ln10();
818  }
819  ElecGSFTrackEta = theGsfEle.gsfTrack()->eta();
820  }
821 
822  return MVAValue(TauEtaAtEcalEntrance,
823  TauPt,
824  TaudCrackEta,
825  TaudCrackPhi,
826  TauEmFraction,
827  TauSignalPFGammaCands,
828  TauLeadPFChargedHadrHoP,
829  TauLeadPFChargedHadrEoP,
830  TauVisMass,
831  TauHadrMva,
832  GammasdEta,
833  GammasdPhi,
834  GammasPt,
835  TauKFNumHits,
836  TauGSFNumHits,
837  TauGSFChi2,
838  TauGSFTrackResol,
839  TauGSFTracklnPt,
840  TauGSFTrackEta,
841  TauPhi,
842  TauSignalPFChargedCands,
843  TauHasGsf,
844  ElecEta,
845  ElecPhi,
846  ElecPt,
847  ElecEe,
848  ElecEgamma,
849  ElecPin,
850  ElecPout,
851  ElecFbrem,
852  ElecChi2GSF,
853  ElecGSFNumHits,
854  ElecGSFTrackResol,
855  ElecGSFTracklnPt,
856  ElecGSFTrackEta);
857 }
i
int i
Definition: DBlmapReader.cc:9
reco::GsfElectron::trackMomentumAtVtx
math::XYZVectorF trackMomentumAtVtx() const
Definition: GsfElectron.h:272
create_public_lumi_plots.log
log
Definition: create_public_lumi_plots.py:1108
reco::PFTau::emFraction
float emFraction() const
Definition: PFTau.cc:158
reco::PFTau::signalPFChargedHadrCands
const PFCandidateRefVector & signalPFChargedHadrCands() const
Charged hadrons in signal region.
Definition: PFTau.cc:75
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:639
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:276
siStripFEDMonitor_P5_cff.Max
tuple Max
Definition: siStripFEDMonitor_P5_cff.py:56
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:143
reco::PFTau::signalPFCands
const PFCandidateRefVector & signalPFCands() const
PFCandidates in signal region.
Definition: PFTau.cc:73
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:250
reco::LeafCandidate::mass
virtual double mass() const
mass
Definition: LeafCandidate.h:123
reco::GsfElectron::trackMomentumOut
math::XYZVectorF trackMomentumOut() const
Definition: GsfElectron.h:274
AntiElectronIDMVA3::dCrackEta
double dCrackEta(double eta)
Definition: AntiElectronIDMVA3.cc:1035
j
int j
Definition: DBlmapReader.cc:9
AntiElectronIDMVA3::dCrackPhi
double dCrackPhi(double phi, double eta)
Definition: AntiElectronIDMVA3.cc:976
reco::PFTau::leadPFChargedHadrCand
const PFCandidateRef & leadPFChargedHadrCand() const
Definition: PFTau.cc:62
reco::PFTau::signalPFGammaCands
const PFCandidateRefVector & signalPFGammaCands() const
Gamma candidates in signal region.
Definition: PFTau.cc:79
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.cc:127
edm::RefVector::at
value_type const at(size_type idx) const
Retrieve an element of the RefVector.
Definition: RefVector.h:76
edm::RefVector::size
size_type size() const
Size of the RefVector.
Definition: RefVector.h:89
AntiElectronIDMVA3::MVAValue
double MVAValue(Float_t TauEtaAtEcalEntrance, Float_t TauPt, Float_t TaudCrackEta, Float_t TaudCrackPhi, Float_t TauEmFraction, Float_t TauSignalPFGammaCands, Float_t TauLeadPFChargedHadrHoP, Float_t TauLeadPFChargedHadrEoP, Float_t TauVisMass, Float_t TauHadrMva, const std::vector< Float_t > &GammasdEta, const std::vector< Float_t > &GammasdPhi, const std::vector< Float_t > &GammasPt, Float_t TauKFNumHits, Float_t TauGSFNumHits, Float_t TauGSFChi2, Float_t TauGSFTrackResol, Float_t TauGSFTracklnPt, Float_t TauGSFTrackEta, Float_t TauPhi, Float_t TauSignalPFChargedCands, Float_t TauHasGsf, Float_t ElecEta, Float_t ElecPhi, Float_t ElecPt, Float_t ElecEe, Float_t ElecEgamma, Float_t ElecPin, Float_t ElecPout, Float_t ElecFbrem, Float_t ElecChi2GSF, Float_t ElecGSFNumHits, Float_t ElecGSFTrackResol, Float_t ElecGSFTracklnPt, Float_t ElecGSFTrackEta)
Definition: AntiElectronIDMVA3.cc:504
reco::PFTau::electronPreIDOutput
float electronPreIDOutput() const
Definition: PFTau.cc:165
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:139
reco::GsfElectron::pflowSuperCluster
SuperClusterRef pflowSuperCluster() const
Definition: GsfElectron.h:174
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
asciidump.at
list at
Definition: asciidump.py:428
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:169
double AntiElectronIDMVA3::MVAValue ( const reco::PFTau &  thePFTau)

Definition at line 859 of file AntiElectronIDMVA3.cc.

References edm::RefVector< C, T, F >::at(), asciidump::at, dCrackEta(), dCrackPhi(), reco::PFTau::electronPreIDOutput(), reco::PFTau::emFraction(), reco::LeafCandidate::eta(), i, edm::Ref< C, T, F >::isNonnull(), j, reco::PFTau::leadPFChargedHadrCand(), create_public_lumi_plots::log, reco::LeafCandidate::mass(), siStripFEDMonitor_P5_cff::Max, MVAValue(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::PFTau::signalPFCands(), reco::PFTau::signalPFChargedHadrCands(), reco::PFTau::signalPFGammaCands(), edm::RefVector< C, T, F >::size(), and findQualityFiles::size.

860 {
861  Float_t TauEtaAtEcalEntrance = -99.;
862  float sumEtaTimesEnergy = 0;
863  float sumEnergy = 0;
864  for(unsigned int j = 0 ; j < (thePFTau.signalPFCands()).size() ; j++){
865  reco::PFCandidateRef pfcandidate = (thePFTau.signalPFCands()).at(j);
866  sumEtaTimesEnergy += pfcandidate->positionAtECALEntrance().eta()*pfcandidate->energy();
867  sumEnergy += pfcandidate->energy();
868  }
869  if(sumEnergy>0)TauEtaAtEcalEntrance = sumEtaTimesEnergy/sumEnergy;
870 
871  Float_t TauPt = thePFTau.pt();
872  Float_t TauEmFraction = TMath::Max(thePFTau.emFraction(), (Float_t)0.);
873  Float_t TauSignalPFGammaCands = thePFTau.signalPFGammaCands().size();
874  Float_t TauLeadPFChargedHadrHoP = 0.;
875  Float_t TauLeadPFChargedHadrEoP = 0.;
876  if ( thePFTau.leadPFChargedHadrCand()->p() > 0. ) {
877  TauLeadPFChargedHadrHoP = thePFTau.leadPFChargedHadrCand()->hcalEnergy()/thePFTau.leadPFChargedHadrCand()->p();
878  TauLeadPFChargedHadrEoP = thePFTau.leadPFChargedHadrCand()->ecalEnergy()/thePFTau.leadPFChargedHadrCand()->p();
879  }
880  Float_t TauVisMass = thePFTau.mass();
881  Float_t TauHadrMva = TMath::Max(thePFTau.electronPreIDOutput(),float(-1.0));
882  std::vector<Float_t> GammasdEta;
883  std::vector<Float_t> GammasdPhi;
884  std::vector<Float_t> GammasPt;
885  for ( unsigned i = 0 ; i < thePFTau.signalPFGammaCands().size(); ++i ) {
886  reco::PFCandidateRef gamma = thePFTau.signalPFGammaCands().at(i);
887  if ( thePFTau.leadPFChargedHadrCand().isNonnull() ) {
888  GammasdEta.push_back(gamma->eta() - thePFTau.leadPFChargedHadrCand()->eta());
889  GammasdPhi.push_back(gamma->phi() - thePFTau.leadPFChargedHadrCand()->phi());
890  } else {
891  GammasdEta.push_back(gamma->eta() - thePFTau.eta());
892  GammasdPhi.push_back(gamma->phi() - thePFTau.phi());
893  }
894  GammasPt.push_back(gamma->pt());
895  }
896  Float_t TauKFNumHits = -99.;
897  if((thePFTau.leadPFChargedHadrCand()->trackRef()).isNonnull()){
898  TauKFNumHits = thePFTau.leadPFChargedHadrCand()->trackRef()->numberOfValidHits();
899  }
900  Float_t TauGSFNumHits = -99.;
901  Float_t TauGSFChi2 = -99.;
902  Float_t TauGSFTrackResol = -99.;
903  Float_t TauGSFTracklnPt = -99.;
904  Float_t TauGSFTrackEta = -99.;
905  if((thePFTau.leadPFChargedHadrCand()->gsfTrackRef()).isNonnull()){
906  TauGSFChi2 = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->normalizedChi2();
907  TauGSFNumHits = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->numberOfValidHits();
908  if ( thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->pt() > 0. ) {
909  TauGSFTrackResol = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->ptError()/thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->pt();
910  TauGSFTracklnPt = log(thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->pt())*TMath::Ln10();
911  }
912  TauGSFTrackEta = thePFTau.leadPFChargedHadrCand()->gsfTrackRef()->eta();
913  }
914  Float_t TauPhi = thePFTau.phi();
915  float sumPhiTimesEnergy = 0;
916  float sumEnergyPhi = 0;
917  for(unsigned int j = 0 ; j < (thePFTau.signalPFCands()).size() ; j++){
918  reco::PFCandidateRef pfcandidate = (thePFTau.signalPFCands()).at(j);
919  sumPhiTimesEnergy += pfcandidate->positionAtECALEntrance().phi()*pfcandidate->energy();
920  sumEnergyPhi += pfcandidate->energy();
921  }
922  if(sumEnergy>0)TauPhi = sumPhiTimesEnergy/sumEnergyPhi;
923  Float_t TaudCrackPhi = dCrackPhi(TauPhi,TauEtaAtEcalEntrance) ;
924  Float_t TaudCrackEta = dCrackEta(TauEtaAtEcalEntrance) ;
925  Float_t TauSignalPFChargedCands = thePFTau.signalPFChargedHadrCands().size();
926  Float_t TauHasGsf = thePFTau.leadPFChargedHadrCand()->gsfTrackRef().isNonnull();
927 
928  Float_t dummyElecEta = 9.9;
929 
930  return MVAValue(TauEtaAtEcalEntrance,
931  TauPt,
932  TaudCrackEta,
933  TaudCrackPhi,
934  TauEmFraction,
935  TauSignalPFGammaCands,
936  TauLeadPFChargedHadrHoP,
937  TauLeadPFChargedHadrEoP,
938  TauVisMass,
939  TauHadrMva,
940  GammasdEta,
941  GammasdPhi,
942  GammasPt,
943  TauKFNumHits,
944  TauGSFNumHits,
945  TauGSFChi2,
946  TauGSFTrackResol,
947  TauGSFTracklnPt,
948  TauGSFTrackEta,
949  TauPhi,
950  TauSignalPFChargedCands,
951  TauHasGsf,
952  dummyElecEta,
953  0.,
954  0.,
955  0.,
956  0.,
957  0.,
958  0.,
959  0.,
960  0.,
961  0.,
962  0.,
963  0.,
964  0.);
965 }
i
int i
Definition: DBlmapReader.cc:9
create_public_lumi_plots.log
log
Definition: create_public_lumi_plots.py:1108
reco::PFTau::emFraction
float emFraction() const
Definition: PFTau.cc:158
reco::PFTau::signalPFChargedHadrCands
const PFCandidateRefVector & signalPFChargedHadrCands() const
Charged hadrons in signal region.
Definition: PFTau.cc:75
siStripFEDMonitor_P5_cff.Max
tuple Max
Definition: siStripFEDMonitor_P5_cff.py:56
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:143
reco::PFTau::signalPFCands
const PFCandidateRefVector & signalPFCands() const
PFCandidates in signal region.
Definition: PFTau.cc:73
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:250
reco::LeafCandidate::mass
virtual double mass() const
mass
Definition: LeafCandidate.h:123
AntiElectronIDMVA3::dCrackEta
double dCrackEta(double eta)
Definition: AntiElectronIDMVA3.cc:1035
j
int j
Definition: DBlmapReader.cc:9
AntiElectronIDMVA3::dCrackPhi
double dCrackPhi(double phi, double eta)
Definition: AntiElectronIDMVA3.cc:976
reco::PFTau::leadPFChargedHadrCand
const PFCandidateRef & leadPFChargedHadrCand() const
Definition: PFTau.cc:62
reco::PFTau::signalPFGammaCands
const PFCandidateRefVector & signalPFGammaCands() const
Gamma candidates in signal region.
Definition: PFTau.cc:79
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.cc:127
edm::RefVector::at
value_type const at(size_type idx) const
Retrieve an element of the RefVector.
Definition: RefVector.h:76
edm::RefVector::size
size_type size() const
Size of the RefVector.
Definition: RefVector.h:89
AntiElectronIDMVA3::MVAValue
double MVAValue(Float_t TauEtaAtEcalEntrance, Float_t TauPt, Float_t TaudCrackEta, Float_t TaudCrackPhi, Float_t TauEmFraction, Float_t TauSignalPFGammaCands, Float_t TauLeadPFChargedHadrHoP, Float_t TauLeadPFChargedHadrEoP, Float_t TauVisMass, Float_t TauHadrMva, const std::vector< Float_t > &GammasdEta, const std::vector< Float_t > &GammasdPhi, const std::vector< Float_t > &GammasPt, Float_t TauKFNumHits, Float_t TauGSFNumHits, Float_t TauGSFChi2, Float_t TauGSFTrackResol, Float_t TauGSFTracklnPt, Float_t TauGSFTrackEta, Float_t TauPhi, Float_t TauSignalPFChargedCands, Float_t TauHasGsf, Float_t ElecEta, Float_t ElecPhi, Float_t ElecPt, Float_t ElecEe, Float_t ElecEgamma, Float_t ElecPin, Float_t ElecPout, Float_t ElecFbrem, Float_t ElecChi2GSF, Float_t ElecGSFNumHits, Float_t ElecGSFTrackResol, Float_t ElecGSFTracklnPt, Float_t ElecGSFTrackEta)
Definition: AntiElectronIDMVA3.cc:504
reco::PFTau::electronPreIDOutput
float electronPreIDOutput() const
Definition: PFTau.cc:165
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:139
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
asciidump.at
list at
Definition: asciidump.py:428

Member Data Documentation

Float_t AntiElectronIDMVA3::Elec_Chi2GSF_
private

Definition at line 184 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Elec_EgammaOverPdif_
private

Definition at line 182 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Elec_EtotOverPin_
private

Definition at line 181 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Elec_Fbrem_
private

Definition at line 183 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Elec_GSFNumHits_
private

Definition at line 185 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Elec_GSFTrackEta_
private

Definition at line 188 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Elec_GSFTracklnPt_
private

Definition at line 187 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Elec_GSFTrackResol_
private

Definition at line 186 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

TMVA::Reader* AntiElectronIDMVA3::fTMVAReader_[16]
private

Definition at line 156 of file AntiElectronIDMVA3.h.

Referenced by AntiElectronIDMVA3(), bookMVAs(), Initialize_from_file(), Initialize_from_string(), MVAValue(), and ~AntiElectronIDMVA3().

Float_t AntiElectronIDMVA3::GammadEta_
private

Definition at line 158 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs().

Float_t AntiElectronIDMVA3::GammadPhi_
private

Definition at line 159 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs().

Float_t AntiElectronIDMVA3::GammadPt_
private

Definition at line 160 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs().

Bool_t AntiElectronIDMVA3::isInitialized_
private

Definition at line 154 of file AntiElectronIDMVA3.h.

Referenced by Initialize_from_file(), Initialize_from_string(), and MVAValue().

std::string AntiElectronIDMVA3::methodName_
private

Definition at line 155 of file AntiElectronIDMVA3.h.

Referenced by Initialize_from_file(), Initialize_from_string(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_dCrackEta_
private

Definition at line 164 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_dCrackPhi_
private

Definition at line 165 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_EmFraction_
private

Definition at line 166 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_EtaAtEcalEntrance_
private

Definition at line 162 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_GammaEnFrac_
private

Definition at line 174 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_GammaEtaMom_
private

Definition at line 172 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_GammaPhiMom_
private

Definition at line 173 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_GSFChi2_
private

Definition at line 175 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_GSFTrackEta_
private

Definition at line 179 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_GSFTracklnPt_
private

Definition at line 178 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_GSFTrackResol_
private

Definition at line 177 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_HadrEoP_
private

Definition at line 169 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_HadrHoP_
private

Definition at line 168 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_HadrMva_
private

Definition at line 171 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_NumGammaCands_
private

Definition at line 167 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_NumHitsVariable_
private

Definition at line 176 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_Pt_
private

Definition at line 163 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

Float_t AntiElectronIDMVA3::Tau_VisMass_
private

Definition at line 170 of file AntiElectronIDMVA3.h.

Referenced by bookMVAs(), and MVAValue().

int AntiElectronIDMVA3::verbosity_
private

Definition at line 190 of file AntiElectronIDMVA3.h.

Referenced by AntiElectronIDMVA3(), and bookMVAs().

Generated for CMSSW Reference Manual by doxygen 1.8.5