ElectronIDMVA Class Reference

#include <ElectronIDMVA.h>

Public Types
enum	EElectronEffectiveAreaType {   kEleChargedIso03, kEleNeutralHadronIso03, kEleGammaIso03, kEleGammaIsoVetoEtaStrip03,   kEleChargedIso04, kEleNeutralHadronIso04, kEleGammaIso04, kEleGammaIsoVetoEtaStrip04,   kEleNeutralHadronIso007, kEleHoverE, kEleHcalDepth1OverEcal, kEleHcalDepth2OverEcal }

Public Member Functions
Double_t	ElectronEffectiveArea (EElectronEffectiveAreaType type, Double_t Eta)
	ElectronIDMVA ()
void	Initialize (TString methodName, unsigned int version, TString Subdet0Pt10To20Weights, TString Subdet1Pt10To20Weights, TString Subdet2Pt10To20Weights, TString Subdet0Pt20ToInfWeights, TString Subdet1Pt20ToInfWeights, TString Subdet2Pt20ToInfWeights)
Bool_t	IsInitialized () const
Double_t	MVAValue (HWW &, const unsigned int ele, const unsigned int vertex)
Double_t	MVAValue (Double_t ElePt, Double_t EleSCEta, Double_t EleSigmaIEtaIEta, Double_t EleDEtaIn, Double_t EleDPhiIn, Double_t EleD0, Double_t EleDZ, Double_t EleFBrem, Double_t EleEOverP, Double_t EleESeedClusterOverPout, Double_t EleSigmaIPhiIPhi, Double_t EleNBrem, Double_t EleOneOverEMinusOneOverP, Double_t EleESeedClusterOverPIn, Double_t EleIP3d, Double_t EleIP3dSig)
Double_t	MVAValue (Double_t ElePt, Double_t EleEta, Double_t PileupEnergyDensity, Double_t EleSigmaIEtaIEta, Double_t EleDEtaIn, Double_t EleDPhiIn, Double_t EleHoverE, Double_t EleD0, Double_t EleDZ, Double_t EleFBrem, Double_t EleEOverP, Double_t EleESeedClusterOverPout, Double_t EleSigmaIPhiIPhi, Double_t EleNBrem, Double_t EleOneOverEMinusOneOverP, Double_t EleESeedClusterOverPIn, Double_t EleIP3d, Double_t EleIP3dSig, Double_t EleGsfTrackChi2OverNdof, Double_t EledEtaCalo, Double_t EledPhiCalo, Double_t EleR9, Double_t EleSCEtaWidth, Double_t EleSCPhiWidth, Double_t EleCovIEtaIPhi, Double_t ElePreShowerOverRaw, Double_t EleChargedIso03, Double_t EleNeutralHadronIso03, Double_t EleGammaIso03, Double_t EleChargedIso04, Double_t EleNeutralHadronIso04, Double_t EleGammaIso04, Bool_t printDebug=kFALSE)
	~ElectronIDMVA ()

Protected Attributes
Bool_t	fIsInitialized
TString	fMethodname
Float_t	fMVAVar_EleChargedIso03OverPt
Float_t	fMVAVar_EleChargedIso04OverPt
Float_t	fMVAVar_EleCovIEtaIPhi
Float_t	fMVAVar_EleD0
Float_t	fMVAVar_EledEtaCalo
Float_t	fMVAVar_EleDEtaIn
Float_t	fMVAVar_EledPhiCalo
Float_t	fMVAVar_EleDPhiIn
Float_t	fMVAVar_EleDZ
Float_t	fMVAVar_EleEOverP
Float_t	fMVAVar_EleESeedClusterOverPIn
Float_t	fMVAVar_EleESeedClusterOverPout
Float_t	fMVAVar_EleFBrem
Float_t	fMVAVar_EleGammaIso03OverPt
Float_t	fMVAVar_EleGammaIso04OverPt
Float_t	fMVAVar_EleGsfTrackChi2OverNdof
Float_t	fMVAVar_EleHoverE
Float_t	fMVAVar_EleIP3d
Float_t	fMVAVar_EleIP3dSig
Float_t	fMVAVar_EleNBrem
Float_t	fMVAVar_EleNeutralHadronIso03OverPt
Float_t	fMVAVar_EleNeutralHadronIso04OverPt
Float_t	fMVAVar_EleOneOverEMinusOneOverP
Float_t	fMVAVar_ElePreShowerOverRaw
Float_t	fMVAVar_EleR9
Float_t	fMVAVar_EleSCEtaWidth
Float_t	fMVAVar_EleSCPhiWidth
Float_t	fMVAVar_EleSigmaIEtaIEta
Float_t	fMVAVar_EleSigmaIPhiIPhi
TMVA::Reader *	fTMVAReader [6]
unsigned int	version_

Detailed Description

Definition at line 24 of file ElectronIDMVA.h.

Member Enumeration Documentation

enum ElectronIDMVA::EElectronEffectiveAreaType

Enumerator
kEleChargedIso03
kEleNeutralHadronIso03
kEleGammaIso03
kEleGammaIsoVetoEtaStrip03
kEleChargedIso04
kEleNeutralHadronIso04
kEleGammaIso04
kEleGammaIsoVetoEtaStrip04
kEleNeutralHadronIso007
kEleHoverE
kEleHcalDepth1OverEcal
kEleHcalDepth2OverEcal

Definition at line 87 of file ElectronIDMVA.h.

                                    {
      kEleChargedIso03, 
      kEleNeutralHadronIso03, 
      kEleGammaIso03, 
      kEleGammaIsoVetoEtaStrip03, 
      kEleChargedIso04, 
      kEleNeutralHadronIso04, 
      kEleGammaIso04, 
      kEleGammaIsoVetoEtaStrip04, 
      kEleNeutralHadronIso007, 
      kEleHoverE, 
      kEleHcalDepth1OverEcal, 
      kEleHcalDepth2OverEcal    
    };

Constructor & Destructor Documentation

ElectronIDMVA::ElectronIDMVA

(

)

Definition at line 15 of file ElectronIDMVA.cc.

References fTMVAReader, and i.

                             :
    fMethodname("BDTG method"),
    fIsInitialized(kFALSE)
{
    // Constructor.
    for(UInt_t i=0; i<6; ++i) {
        fTMVAReader[i] = 0;
    }
}

ElectronIDMVA::~ElectronIDMVA

(

)

Definition at line 27 of file ElectronIDMVA.cc.

References fTMVAReader, and i.

{
    for(UInt_t i=0; i<6; ++i) {
        if (fTMVAReader[i]) delete fTMVAReader[i];
    }
}

Member Function Documentation

Double_t ElectronIDMVA::ElectronEffectiveArea	(	EElectronEffectiveAreaType	type,
		Double_t	Eta
	)

Definition at line 469 of file ElectronIDMVA.cc.

References kEleChargedIso03, kEleChargedIso04, kEleGammaIso03, kEleGammaIso04, kEleGammaIsoVetoEtaStrip03, kEleGammaIsoVetoEtaStrip04, kEleHcalDepth1OverEcal, kEleHcalDepth2OverEcal, kEleHoverE, kEleNeutralHadronIso007, kEleNeutralHadronIso03, and kEleNeutralHadronIso04.

Referenced by MVAValue().

                                                                                           {

  Double_t EffectiveArea = 0;

  if (fabs(Eta) < 1.0) {
    if (type == ElectronIDMVA::kEleChargedIso03) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso03) EffectiveArea = 0.017;
    if (type == ElectronIDMVA::kEleGammaIso03) EffectiveArea = 0.045;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.014;
    if (type == ElectronIDMVA::kEleChargedIso04) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso04) EffectiveArea = 0.034;
    if (type == ElectronIDMVA::kEleGammaIso04) EffectiveArea = 0.079;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.014;
    if (type == ElectronIDMVA::kEleNeutralHadronIso007) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleHoverE) EffectiveArea = 0.00016;
    if (type == ElectronIDMVA::kEleHcalDepth1OverEcal) EffectiveArea = 0.00016;
    if (type == ElectronIDMVA::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
  } else if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) {
    if (type == ElectronIDMVA::kEleChargedIso03) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso03) EffectiveArea = 0.025;
    if (type == ElectronIDMVA::kEleGammaIso03) EffectiveArea = 0.052;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.030;
    if (type == ElectronIDMVA::kEleChargedIso04) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso04) EffectiveArea = 0.050;
    if (type == ElectronIDMVA::kEleGammaIso04) EffectiveArea = 0.073;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.030;
    if (type == ElectronIDMVA::kEleNeutralHadronIso007) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleHoverE) EffectiveArea = 0.00022;
    if (type == ElectronIDMVA::kEleHcalDepth1OverEcal) EffectiveArea = 0.00022;
    if (type == ElectronIDMVA::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
  } else if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) {
    if (type == ElectronIDMVA::kEleChargedIso03) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso03) EffectiveArea = 0.030;
    if (type == ElectronIDMVA::kEleGammaIso03) EffectiveArea = 0.170;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.134;
    if (type == ElectronIDMVA::kEleChargedIso04) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso04) EffectiveArea = 0.060;
    if (type == ElectronIDMVA::kEleGammaIso04) EffectiveArea = 0.187;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.134;
    if (type == ElectronIDMVA::kEleNeutralHadronIso007) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleHoverE) EffectiveArea = 0.00030;
    if (type == ElectronIDMVA::kEleHcalDepth1OverEcal) EffectiveArea = 0.00026;
    if (type == ElectronIDMVA::kEleHcalDepth2OverEcal) EffectiveArea = 0.00002;        
  } else if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.25 ) {
    if (type == ElectronIDMVA::kEleChargedIso03) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso03) EffectiveArea = 0.022;
    if (type == ElectronIDMVA::kEleGammaIso03) EffectiveArea = 0.623;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.516;
    if (type == ElectronIDMVA::kEleChargedIso04) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso04) EffectiveArea = 0.055;
    if (type == ElectronIDMVA::kEleGammaIso04) EffectiveArea = 0.659;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.517;
    if (type == ElectronIDMVA::kEleNeutralHadronIso007) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleHoverE) EffectiveArea = 0.00054;
    if (type == ElectronIDMVA::kEleHcalDepth1OverEcal) EffectiveArea = 0.00045;
    if (type == ElectronIDMVA::kEleHcalDepth2OverEcal) EffectiveArea = 0.00003;
  } else if (fabs(Eta) >= 2.25 && fabs(Eta) < 2.5 ) {
    if (type == ElectronIDMVA::kEleChargedIso03) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso03) EffectiveArea = 0.018;
    if (type == ElectronIDMVA::kEleGammaIso03) EffectiveArea = 1.198;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 1.049;
    if (type == ElectronIDMVA::kEleChargedIso04) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleNeutralHadronIso04) EffectiveArea = 0.073;
    if (type == ElectronIDMVA::kEleGammaIso04) EffectiveArea = 1.258;
    if (type == ElectronIDMVA::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 1.051;
    if (type == ElectronIDMVA::kEleNeutralHadronIso007) EffectiveArea = 0.000;
    if (type == ElectronIDMVA::kEleHoverE) EffectiveArea = 0.00082;
    if (type == ElectronIDMVA::kEleHcalDepth1OverEcal) EffectiveArea = 0.00066;
    if (type == ElectronIDMVA::kEleHcalDepth2OverEcal) EffectiveArea = 0.00004;
  }
    
  return EffectiveArea;  
}

void ElectronIDMVA::Initialize	(	TString	methodName,
		unsigned int	version,
		TString	Subdet0Pt10To20Weights,
		TString	Subdet1Pt10To20Weights,
		TString	Subdet2Pt10To20Weights,
		TString	Subdet0Pt20ToInfWeights,
		TString	Subdet1Pt20ToInfWeights,
		TString	Subdet2Pt20ToInfWeights
	)

Definition at line 35 of file ElectronIDMVA.cc.

References fIsInitialized, fMethodname, fMVAVar_EleChargedIso03OverPt, fMVAVar_EleChargedIso04OverPt, fMVAVar_EleCovIEtaIPhi, fMVAVar_EleD0, fMVAVar_EledEtaCalo, fMVAVar_EleDEtaIn, fMVAVar_EledPhiCalo, fMVAVar_EleDPhiIn, fMVAVar_EleEOverP, fMVAVar_EleESeedClusterOverPIn, fMVAVar_EleESeedClusterOverPout, fMVAVar_EleFBrem, fMVAVar_EleGammaIso03OverPt, fMVAVar_EleGammaIso04OverPt, fMVAVar_EleGsfTrackChi2OverNdof, fMVAVar_EleIP3d, fMVAVar_EleIP3dSig, fMVAVar_EleNBrem, fMVAVar_EleNeutralHadronIso03OverPt, fMVAVar_EleNeutralHadronIso04OverPt, fMVAVar_EleOneOverEMinusOneOverP, fMVAVar_ElePreShowerOverRaw, fMVAVar_EleR9, fMVAVar_EleSCEtaWidth, fMVAVar_EleSCPhiWidth, fMVAVar_EleSigmaIEtaIEta, fMVAVar_EleSigmaIPhiIPhi, fTMVAReader, i, BeamSplash_cfg::version, and version_.

                                         {

    if (version != 1 && version != 2 && version != 3) {
        edm::LogError("InvalidInput") << "[ElectronIDMVA::Initialize] Version must be 1 or 2 or 3.  Aborting.";
        return;
    }

    version_ = version;

    fIsInitialized = kTRUE;
    fMethodname = methodName;

    for(UInt_t i=0; i<6; ++i) {
        if (fTMVAReader[i]) delete fTMVAReader[i];

    fTMVAReader[i] = new TMVA::Reader( "!Color:!Silent:Error" );  
    fTMVAReader[i]->SetVerbose(kTRUE);
        // order matters!

        if (version == 1) {
            fTMVAReader[i] = new TMVA::Reader( "!Color:!Silent:Error" );  
            fTMVAReader[i]->SetVerbose(kTRUE);
            fTMVAReader[i]->AddVariable( "SigmaIEtaIEta",         &fMVAVar_EleSigmaIEtaIEta         );
            fTMVAReader[i]->AddVariable( "DEtaIn",                &fMVAVar_EleDEtaIn                );
            fTMVAReader[i]->AddVariable( "DPhiIn",                &fMVAVar_EleDPhiIn                );
            fTMVAReader[i]->AddVariable( "FBrem",                 &fMVAVar_EleFBrem                 );
            fTMVAReader[i]->AddVariable( "EOverP",                &fMVAVar_EleEOverP                );
            fTMVAReader[i]->AddVariable( "ESeedClusterOverPout",  &fMVAVar_EleESeedClusterOverPout  );
            fTMVAReader[i]->AddVariable( "SigmaIPhiIPhi",         &fMVAVar_EleSigmaIPhiIPhi         );
            fTMVAReader[i]->AddVariable( "NBrem",                 &fMVAVar_EleNBrem                 );
            fTMVAReader[i]->AddVariable( "OneOverEMinusOneOverP", &fMVAVar_EleOneOverEMinusOneOverP );
            fTMVAReader[i]->AddVariable( "ESeedClusterOverPIn",   &fMVAVar_EleESeedClusterOverPIn   );
        } else if (version == 2) {
            fTMVAReader[i] = new TMVA::Reader( "!Color:!Silent:Error" );
            fTMVAReader[i]->AddVariable( "SigmaIEtaIEta",         &fMVAVar_EleSigmaIEtaIEta         );
            fTMVAReader[i]->AddVariable( "DEtaIn",                &fMVAVar_EleDEtaIn                );
            fTMVAReader[i]->AddVariable( "DPhiIn",                &fMVAVar_EleDPhiIn                );
            fTMVAReader[i]->AddVariable( "D0",                    &fMVAVar_EleD0                    );
            fTMVAReader[i]->AddVariable( "FBrem",                 &fMVAVar_EleFBrem                 );
            fTMVAReader[i]->AddVariable( "EOverP",                &fMVAVar_EleEOverP                );
            fTMVAReader[i]->AddVariable( "ESeedClusterOverPout",  &fMVAVar_EleESeedClusterOverPout  );
            fTMVAReader[i]->AddVariable( "SigmaIPhiIPhi",         &fMVAVar_EleSigmaIPhiIPhi         );
            fTMVAReader[i]->AddVariable( "NBrem",                 &fMVAVar_EleNBrem                 );
            fTMVAReader[i]->AddVariable( "OneOverEMinusOneOverP", &fMVAVar_EleOneOverEMinusOneOverP );
            fTMVAReader[i]->AddVariable( "ESeedClusterOverPIn",   &fMVAVar_EleESeedClusterOverPIn   );
            fTMVAReader[i]->AddVariable( "IP3d",                  &fMVAVar_EleIP3d                  );
            fTMVAReader[i]->AddVariable( "IP3dSig",               &fMVAVar_EleIP3dSig               );
        } else if (version == 3) {
        fTMVAReader[i]->AddVariable( "SigmaIEtaIEta",         &fMVAVar_EleSigmaIEtaIEta            );
        fTMVAReader[i]->AddVariable( "DEtaIn",                &fMVAVar_EleDEtaIn                   );
        fTMVAReader[i]->AddVariable( "DPhiIn",                &fMVAVar_EleDPhiIn                   );
        //fTMVAReader[i]->AddVariable( "HoverE",                &fMVAVar_EleHoverE                   );
        fTMVAReader[i]->AddVariable( "D0",                    &fMVAVar_EleD0                       );
        fTMVAReader[i]->AddVariable( "FBrem",                 &fMVAVar_EleFBrem                    );
        fTMVAReader[i]->AddVariable( "EOverP",                &fMVAVar_EleEOverP                   );
        fTMVAReader[i]->AddVariable( "ESeedClusterOverPout",  &fMVAVar_EleESeedClusterOverPout     );
        fTMVAReader[i]->AddVariable( "SigmaIPhiIPhi",         &fMVAVar_EleSigmaIPhiIPhi            );
        fTMVAReader[i]->AddVariable( "OneOverEMinusOneOverP", &fMVAVar_EleOneOverEMinusOneOverP    );      
        fTMVAReader[i]->AddVariable( "ESeedClusterOverPIn",   &fMVAVar_EleESeedClusterOverPIn      );
        fTMVAReader[i]->AddVariable( "IP3d",                  &fMVAVar_EleIP3d                     );
        fTMVAReader[i]->AddVariable( "IP3dSig",               &fMVAVar_EleIP3dSig                  );
        
        fTMVAReader[i]->AddVariable( "GsfTrackChi2OverNdof",  &fMVAVar_EleGsfTrackChi2OverNdof     );
        fTMVAReader[i]->AddVariable( "dEtaCalo",              &fMVAVar_EledEtaCalo                 );
        fTMVAReader[i]->AddVariable( "dPhiCalo",              &fMVAVar_EledPhiCalo                 );
        fTMVAReader[i]->AddVariable( "R9",                    &fMVAVar_EleR9                       );
        fTMVAReader[i]->AddVariable( "SCEtaWidth",            &fMVAVar_EleSCEtaWidth               );
        fTMVAReader[i]->AddVariable( "SCPhiWidth",            &fMVAVar_EleSCPhiWidth               );
        fTMVAReader[i]->AddVariable( "CovIEtaIPhi",           &fMVAVar_EleCovIEtaIPhi              );
        if (i == 2 || i == 5) {
          fTMVAReader[i]->AddVariable( "PreShowerOverRaw",      &fMVAVar_ElePreShowerOverRaw       );
        }
        fTMVAReader[i]->AddVariable( "ChargedIso03",          &fMVAVar_EleChargedIso03OverPt       );
        fTMVAReader[i]->AddVariable( "NeutralHadronIso03",    &fMVAVar_EleNeutralHadronIso03OverPt );
        fTMVAReader[i]->AddVariable( "GammaIso03",            &fMVAVar_EleGammaIso03OverPt         );
        fTMVAReader[i]->AddVariable( "ChargedIso04",          &fMVAVar_EleChargedIso04OverPt       );
        fTMVAReader[i]->AddVariable( "NeutralHadronIso04",    &fMVAVar_EleNeutralHadronIso04OverPt );
        fTMVAReader[i]->AddVariable( "GammaIso04",            &fMVAVar_EleGammaIso04OverPt         );
    }

        if (i==0) fTMVAReader[i]->BookMVA(fMethodname , Subdet0Pt10To20Weights );
        if (i==1) fTMVAReader[i]->BookMVA(fMethodname , Subdet1Pt10To20Weights );
        if (i==2) fTMVAReader[i]->BookMVA(fMethodname , Subdet2Pt10To20Weights );
        if (i==3) fTMVAReader[i]->BookMVA(fMethodname , Subdet0Pt20ToInfWeights );
        if (i==4) fTMVAReader[i]->BookMVA(fMethodname , Subdet1Pt20ToInfWeights );
        if (i==5) fTMVAReader[i]->BookMVA(fMethodname , Subdet2Pt20ToInfWeights );

    }

}

Bool_t ElectronIDMVA::IsInitialized ( ) const

inline

Definition at line 37 of file ElectronIDMVA.h.

References fIsInitialized.

{ return fIsInitialized; }

Double_t ElectronIDMVA::MVAValue	(	HWW &	hww,
		const unsigned int	ele,
		const unsigned int	vertex
	)

Definition at line 133 of file ElectronIDMVA.cc.

References HWWFunctions::electron_d0PV_wwV1(), HWWFunctions::electron_dzPV_wwV1(), ElectronEffectiveArea(), HWW::els_chi2(), HWW::els_dEtaIn(), HWW::els_dEtaOut(), HWW::els_dPhiIn(), HWW::els_dPhiOut(), HWW::els_e3x3(), HWW::els_eOverPIn(), HWW::els_eSC(), HWW::els_eSCPresh(), HWW::els_eSCRaw(), HWW::els_eSeedOverPIn(), HWW::els_eSeedOverPOut(), HWW::els_etaSC(), HWW::els_etaSCwidth(), HWW::els_fbrem(), HWW::els_gsftrkidx(), HWW::els_hOverE(), HWW::els_ip3d(), HWW::els_ip3derr(), HWW::els_iso03_pf_ch(), HWW::els_iso03_pf_gamma05(), HWW::els_iso03_pf_nhad05(), HWW::els_iso04_pf_ch(), HWW::els_iso04_pf_gamma05(), HWW::els_iso04_pf_nhad05(), HWW::els_ndof(), HWW::els_nSeed(), HWW::els_p4(), HWW::els_phiSCwidth(), HWW::els_scindex(), HWW::els_sigmaIEtaIEta(), HWW::els_sigmaIPhiIPhi(), HWW::els_trk_p4(), HWW::evt_event(), HWW::evt_lumiBlock(), HWW::evt_run(), HWW::evt_ww_rho_vor(), first, fIsInitialized, fMethodname, fMVAVar_EleChargedIso03OverPt, fMVAVar_EleChargedIso04OverPt, fMVAVar_EleCovIEtaIPhi, fMVAVar_EleD0, fMVAVar_EledEtaCalo, fMVAVar_EleDEtaIn, fMVAVar_EledPhiCalo, fMVAVar_EleDPhiIn, fMVAVar_EleDZ, fMVAVar_EleEOverP, fMVAVar_EleESeedClusterOverPIn, fMVAVar_EleESeedClusterOverPout, fMVAVar_EleFBrem, fMVAVar_EleGammaIso03OverPt, fMVAVar_EleGammaIso04OverPt, fMVAVar_EleGsfTrackChi2OverNdof, fMVAVar_EleHoverE, fMVAVar_EleIP3d, fMVAVar_EleIP3dSig, fMVAVar_EleNBrem, fMVAVar_EleNeutralHadronIso03OverPt, fMVAVar_EleNeutralHadronIso04OverPt, fMVAVar_EleOneOverEMinusOneOverP, fMVAVar_ElePreShowerOverRaw, fMVAVar_EleR9, fMVAVar_EleSCEtaWidth, fMVAVar_EleSCPhiWidth, fMVAVar_EleSigmaIEtaIEta, fMVAVar_EleSigmaIPhiIPhi, fTMVAReader, HWWFunctions::gsftrks_d0_pv(), kEleChargedIso03, kEleChargedIso04, kEleGammaIso03, kEleGammaIso04, kEleGammaIsoVetoEtaStrip03, kEleGammaIsoVetoEtaStrip04, kEleNeutralHadronIso007, kEleNeutralHadronIso03, kEleNeutralHadronIso04, LogDebug, funct::pow(), matplotRender::reader, HWW::scs_sigmaIEtaIPhi(), and version_.

                                                                                            {

    if (!fIsInitialized) { 
        edm::LogError("NotInitialized") << "Error: ElectronIDMVA not properly initialized."; 
        return -9999;
    }

    Int_t subdet = 0;
    if (fabs(hww.els_etaSC().at(ele)) < 1.0) subdet = 0;
    else if (fabs(hww.els_etaSC().at(ele)) < 1.479) subdet = 1;
    else subdet = 2;
    Int_t ptBin = 0;
    if (hww.els_p4().at(ele).Pt() > 20.0) ptBin = 1;

    //set all input variables
    fMVAVar_EleSigmaIEtaIEta =          hww.els_sigmaIEtaIEta().at(ele); 
    fMVAVar_EleDEtaIn =                 hww.els_dEtaIn().at(ele); 
    fMVAVar_EleDPhiIn =                 hww.els_dPhiIn().at(ele); 
    if (version_ == 3) fMVAVar_EleHoverE =                 hww.els_hOverE().at(ele); //this is new
    fMVAVar_EleD0 =                     electron_d0PV_wwV1(hww, ele);
    fMVAVar_EleDZ =                     electron_dzPV_wwV1(hww, ele);
    fMVAVar_EleFBrem =                  hww.els_fbrem().at(ele); 
    fMVAVar_EleEOverP =                 hww.els_eOverPIn().at(ele);
    fMVAVar_EleESeedClusterOverPout =   hww.els_eSeedOverPOut().at(ele);
    fMVAVar_EleSigmaIPhiIPhi =          hww.els_sigmaIPhiIPhi().at(ele);
    fMVAVar_EleNBrem =                  hww.els_nSeed().at(ele);
    TVector3 pIn(hww.els_trk_p4().at(ele).px(), hww.els_trk_p4().at(ele).py(), hww.els_trk_p4().at(ele).pz());
    fMVAVar_EleOneOverEMinusOneOverP =  1./hww.els_eSC().at(ele) - 1./pIn.Mag();
    fMVAVar_EleESeedClusterOverPIn =    hww.els_eSeedOverPIn().at(ele);
    const double gsfsign   = ( (gsftrks_d0_pv(hww, hww.els_gsftrkidx().at(ele),0).first)   >=0 ) ? 1. : -1.;
    fMVAVar_EleIP3d =                   hww.els_ip3d().at(ele)*gsfsign; 
    if (hww.els_ip3derr().at(ele) == 0.0) fMVAVar_EleIP3dSig = 0.0;
    else fMVAVar_EleIP3dSig =           hww.els_ip3d().at(ele)*gsfsign / hww.els_ip3derr().at(ele); 

    if (version_ == 3) {    //these are new
      Double_t ElePt = hww.els_p4().at(ele).pt();
      Double_t EleEta = hww.els_etaSC().at(ele);
      Double_t Rho = hww.evt_ww_rho_vor();
      fMVAVar_EleGsfTrackChi2OverNdof = hww.els_chi2().at(ele) / hww.els_ndof().at(ele);
      fMVAVar_EledEtaCalo = hww.els_dEtaOut().at(ele); 
      fMVAVar_EledPhiCalo = hww.els_dPhiOut().at(ele);
      fMVAVar_EleR9 = hww.els_e3x3().at(ele) / hww.els_eSCRaw().at(ele);
      fMVAVar_EleSCEtaWidth = hww.els_etaSCwidth().at(ele);
      fMVAVar_EleSCPhiWidth = hww.els_phiSCwidth().at(ele);
      fMVAVar_EleCovIEtaIPhi = hww.scs_sigmaIEtaIPhi().at(hww.els_scindex().at(ele))>0 ? pow(hww.scs_sigmaIEtaIPhi().at(hww.els_scindex().at(ele)),2) : -1.* pow(hww.scs_sigmaIEtaIPhi().at(hww.els_scindex().at(ele)),2);
      fMVAVar_ElePreShowerOverRaw = hww.els_eSCPresh().at(ele) / hww.els_eSCRaw().at(ele);
      fMVAVar_EleChargedIso03OverPt 
    = (hww.els_iso03_pf_ch().at(ele) //electronIsoValuePF(ele, vertex, 0.3, 99999., 0.1, 0.07, 0.025, -999., 0)*ElePt
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleChargedIso03, EleEta)) / ElePt;
      fMVAVar_EleNeutralHadronIso03OverPt 
    = (hww.els_iso03_pf_nhad05().at(ele) //electronIsoValuePF(ele, vertex, 0.3, 0.5, 0.1, 0.07, 0.025, -999., 130)*ElePt
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso03, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso007,EleEta)) / ElePt;
      fMVAVar_EleGammaIso03OverPt 
    = (hww.els_iso03_pf_gamma05().at(ele) //electronIsoValuePF(ele, vertex, 0.3, 0.5, 0.1, 0.07, 0.025, -999., 22)*ElePt
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIso03, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIsoVetoEtaStrip03,EleEta))/ElePt;      
      fMVAVar_EleChargedIso04OverPt 
    = (hww.els_iso04_pf_ch().at(ele) //electronIsoValuePF(ele, vertex, 0.4, 99999., 0.1, 0.07, 0.025, -999., 0)*ElePt
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleChargedIso04, EleEta)) / ElePt;
      fMVAVar_EleNeutralHadronIso04OverPt 
    = (hww.els_iso04_pf_nhad05().at(ele) //electronIsoValuePF(ele, vertex, 0.4, 0.5, 0.1, 0.07, 0.025, -999., 130)*ElePt
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso04, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso007,EleEta)) / ElePt;
      fMVAVar_EleGammaIso04OverPt 
    = (hww.els_iso04_pf_gamma05().at(ele) //electronIsoValuePF(ele, vertex, 0.4, 0.5, 0.1, 0.07, 0.025, -999., 22)*ElePt
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIso04, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIsoVetoEtaStrip04,EleEta))/ElePt;
    }

    Double_t mva = -9999;  
    TMVA::Reader *reader = 0;
    Int_t MVABin = -1;
    if (subdet == 0 && ptBin == 0) MVABin = 0;
    if (subdet == 1 && ptBin == 0) MVABin = 1;
    if (subdet == 2 && ptBin == 0) MVABin = 2;
    if (subdet == 0 && ptBin == 1) MVABin = 3;
    if (subdet == 1 && ptBin == 1) MVABin = 4;
    if (subdet == 2 && ptBin == 1) MVABin = 5;
    assert(MVABin >= 0 && MVABin <= 5);
    reader = fTMVAReader[MVABin];

    mva = reader->EvaluateMVA( fMethodname );

    
    //DEBUG
    if (0) {
      LogDebug("ElectronIDMVA") << hww.evt_run() << " " << hww.evt_lumiBlock() << " " << hww.evt_event() << " " << hww.evt_ww_rho_vor()
                                    << hww.els_p4().at(ele).pt() << " " << hww.els_etaSC().at(ele) << " --> MVABin " << MVABin << " : "     
                                    << fMVAVar_EleSigmaIEtaIEta << " " 
                                    << fMVAVar_EleDEtaIn << " " 
                                    << fMVAVar_EleDPhiIn << " " 
                                    << fMVAVar_EleHoverE << " " 
                                    << fMVAVar_EleD0 << " " 
                                    << fMVAVar_EleDZ << " " 
                                    << fMVAVar_EleFBrem << " " 
                                    << fMVAVar_EleEOverP << " " 
                                    << fMVAVar_EleESeedClusterOverPout << " " 
                                    << fMVAVar_EleSigmaIPhiIPhi << " " 
                                    << fMVAVar_EleNBrem << " " 
                                    << fMVAVar_EleOneOverEMinusOneOverP << " " 
                                    << fMVAVar_EleESeedClusterOverPIn << " " 
                                    << fMVAVar_EleIP3d << " " 
                                    << fMVAVar_EleIP3dSig << " " 
                                    << fMVAVar_EleGsfTrackChi2OverNdof << " "
                                    << fMVAVar_EledEtaCalo << " "
                                    << fMVAVar_EledPhiCalo << " "
                                    << fMVAVar_EleR9 << " "
                                    << fMVAVar_EleSCEtaWidth << " "
                                    << fMVAVar_EleSCPhiWidth << " "
                                    << fMVAVar_EleCovIEtaIPhi << " "
                                    << fMVAVar_ElePreShowerOverRaw << " "
                                    << fMVAVar_EleChargedIso03OverPt  << " "
                                    << fMVAVar_EleNeutralHadronIso03OverPt  << " "
                                    << fMVAVar_EleGammaIso03OverPt  << " "
                                    << fMVAVar_EleChargedIso04OverPt  << " "
                                    << fMVAVar_EleNeutralHadronIso04OverPt  << " "
                                    << fMVAVar_EleGammaIso04OverPt  << " "
                                    << " === : === "
                                    << mva;
    }
    return mva;

}

Double_t ElectronIDMVA::MVAValue	(	Double_t	ElePt,
		Double_t	EleSCEta,
		Double_t	EleSigmaIEtaIEta,
		Double_t	EleDEtaIn,
		Double_t	EleDPhiIn,
		Double_t	EleD0,
		Double_t	EleDZ,
		Double_t	EleFBrem,
		Double_t	EleEOverP,
		Double_t	EleESeedClusterOverPout,
		Double_t	EleSigmaIPhiIPhi,
		Double_t	EleNBrem,
		Double_t	EleOneOverEMinusOneOverP,
		Double_t	EleESeedClusterOverPIn,
		Double_t	EleIP3d,
		Double_t	EleIP3dSig
	)

Definition at line 259 of file ElectronIDMVA.cc.

References fIsInitialized, fMethodname, fMVAVar_EleD0, fMVAVar_EleDEtaIn, fMVAVar_EleDPhiIn, fMVAVar_EleDZ, fMVAVar_EleEOverP, fMVAVar_EleESeedClusterOverPIn, fMVAVar_EleESeedClusterOverPout, fMVAVar_EleFBrem, fMVAVar_EleIP3d, fMVAVar_EleIP3dSig, fMVAVar_EleNBrem, fMVAVar_EleOneOverEMinusOneOverP, fMVAVar_EleSigmaIEtaIEta, fMVAVar_EleSigmaIPhiIPhi, fTMVAReader, and matplotRender::reader.

          {

    if (!fIsInitialized) { 
        edm::LogError("NotInitialized") << "Error: ElectronIDMVA not properly initialized."; 
        return -9999;
    }

    Int_t subdet = 0;
    if (fabs(EleSCEta) < 1.0) subdet = 0;
    else if (fabs(EleSCEta) < 1.479) subdet = 1;
    else subdet = 2;
    Int_t ptBin = 0;
    if (ElePt > 20.0) ptBin = 1;

    //set all input variables
    fMVAVar_EleSigmaIEtaIEta = EleSigmaIEtaIEta;
    fMVAVar_EleDEtaIn = EleDEtaIn;
    fMVAVar_EleDPhiIn = EleDPhiIn;
    fMVAVar_EleD0 = EleD0;
    fMVAVar_EleDZ = EleDZ;
    fMVAVar_EleFBrem = EleFBrem;
    fMVAVar_EleEOverP = EleEOverP;
    fMVAVar_EleESeedClusterOverPout = EleESeedClusterOverPout;
    fMVAVar_EleSigmaIPhiIPhi = EleSigmaIPhiIPhi;
    fMVAVar_EleNBrem = EleNBrem;
    fMVAVar_EleOneOverEMinusOneOverP = EleOneOverEMinusOneOverP;
    fMVAVar_EleESeedClusterOverPIn = EleESeedClusterOverPIn;
    fMVAVar_EleIP3d = EleIP3d;
    fMVAVar_EleIP3dSig = EleIP3dSig;

    Double_t mva = -9999;  
    TMVA::Reader *reader = 0;
    Int_t MVABin = -1;
    if (subdet == 0 && ptBin == 0) MVABin = 0;
    if (subdet == 1 && ptBin == 0) MVABin = 1;
    if (subdet == 2 && ptBin == 0) MVABin = 2;
    if (subdet == 0 && ptBin == 1) MVABin = 3;
    if (subdet == 1 && ptBin == 1) MVABin = 4;
    if (subdet == 2 && ptBin == 1) MVABin = 5;
    assert(MVABin >= 0 && MVABin <= 5);
    reader = fTMVAReader[MVABin];

    mva = reader->EvaluateMVA( fMethodname );

    return mva;
}

Double_t ElectronIDMVA::MVAValue	(	Double_t	ElePt,
		Double_t	EleEta,
		Double_t	PileupEnergyDensity,
		Double_t	EleSigmaIEtaIEta,
		Double_t	EleDEtaIn,
		Double_t	EleDPhiIn,
		Double_t	EleHoverE,
		Double_t	EleD0,
		Double_t	EleDZ,
		Double_t	EleFBrem,
		Double_t	EleEOverP,
		Double_t	EleESeedClusterOverPout,
		Double_t	EleSigmaIPhiIPhi,
		Double_t	EleNBrem,
		Double_t	EleOneOverEMinusOneOverP,
		Double_t	EleESeedClusterOverPIn,
		Double_t	EleIP3d,
		Double_t	EleIP3dSig,
		Double_t	EleGsfTrackChi2OverNdof,
		Double_t	EledEtaCalo,
		Double_t	EledPhiCalo,
		Double_t	EleR9,
		Double_t	EleSCEtaWidth,
		Double_t	EleSCPhiWidth,
		Double_t	EleCovIEtaIPhi,
		Double_t	ElePreShowerOverRaw,
		Double_t	EleChargedIso03,
		Double_t	EleNeutralHadronIso03,
		Double_t	EleGammaIso03,
		Double_t	EleChargedIso04,
		Double_t	EleNeutralHadronIso04,
		Double_t	EleGammaIso04,
		Bool_t	printDebug = kFALSE
	)

Definition at line 321 of file ElectronIDMVA.cc.

References ElectronEffectiveArea(), fIsInitialized, fMethodname, fMVAVar_EleChargedIso03OverPt, fMVAVar_EleChargedIso04OverPt, fMVAVar_EleCovIEtaIPhi, fMVAVar_EleD0, fMVAVar_EledEtaCalo, fMVAVar_EleDEtaIn, fMVAVar_EledPhiCalo, fMVAVar_EleDPhiIn, fMVAVar_EleDZ, fMVAVar_EleEOverP, fMVAVar_EleESeedClusterOverPIn, fMVAVar_EleESeedClusterOverPout, fMVAVar_EleFBrem, fMVAVar_EleGammaIso03OverPt, fMVAVar_EleGammaIso04OverPt, fMVAVar_EleGsfTrackChi2OverNdof, fMVAVar_EleHoverE, fMVAVar_EleIP3d, fMVAVar_EleIP3dSig, fMVAVar_EleNBrem, fMVAVar_EleNeutralHadronIso03OverPt, fMVAVar_EleNeutralHadronIso04OverPt, fMVAVar_EleOneOverEMinusOneOverP, fMVAVar_ElePreShowerOverRaw, fMVAVar_EleR9, fMVAVar_EleSCEtaWidth, fMVAVar_EleSCPhiWidth, fMVAVar_EleSigmaIEtaIEta, fMVAVar_EleSigmaIPhiIPhi, fTMVAReader, kEleChargedIso03, kEleChargedIso04, kEleGammaIso03, kEleGammaIso04, kEleGammaIsoVetoEtaStrip03, kEleGammaIsoVetoEtaStrip04, kEleNeutralHadronIso007, kEleNeutralHadronIso03, kEleNeutralHadronIso04, LogDebug, and matplotRender::reader.

    {
  
  if (!fIsInitialized) { 
    edm::LogError("NotInitialized") << "Error: ElectronIDMVA not properly initialized."; 
    return -9999;
  }

  Double_t Rho = 0;
  //jgran if (!(TMath::IsNaN(PileupEnergyDensity) || isinf(PileupEnergyDensity))) Rho = PileupEnergyDensity;
  Rho = PileupEnergyDensity; //jgran

  Int_t subdet = 0;
  if (fabs(EleEta) < 1.0) subdet = 0;
  else if (fabs(EleEta) < 1.479) subdet = 1;
  else subdet = 2;
  Int_t ptBin = 0;
  if (ElePt > 20.0) ptBin = 1;
  
  //set all input variables
  fMVAVar_EleSigmaIEtaIEta = EleSigmaIEtaIEta;
  fMVAVar_EleDEtaIn = EleDEtaIn;
  fMVAVar_EleDPhiIn = EleDPhiIn;
  fMVAVar_EleHoverE = EleHoverE;
  fMVAVar_EleD0 = EleD0;
  fMVAVar_EleDZ = EleDZ;
  fMVAVar_EleFBrem = EleFBrem;
  fMVAVar_EleEOverP = EleEOverP;
  fMVAVar_EleESeedClusterOverPout = EleESeedClusterOverPout;
  fMVAVar_EleSigmaIPhiIPhi = EleSigmaIPhiIPhi;
  fMVAVar_EleNBrem = EleNBrem;
  fMVAVar_EleOneOverEMinusOneOverP = EleOneOverEMinusOneOverP;
  fMVAVar_EleESeedClusterOverPIn = EleESeedClusterOverPIn;
  fMVAVar_EleIP3d = EleIP3d;
  fMVAVar_EleIP3dSig = EleIP3dSig;
  fMVAVar_EleGsfTrackChi2OverNdof = EleGsfTrackChi2OverNdof;
  fMVAVar_EledEtaCalo = EledEtaCalo;
  fMVAVar_EledPhiCalo = EledPhiCalo;
  fMVAVar_EleR9 = EleR9;
  fMVAVar_EleSCEtaWidth = EleSCEtaWidth;
  fMVAVar_EleSCPhiWidth = EleSCPhiWidth;
  fMVAVar_EleCovIEtaIPhi = EleCovIEtaIPhi;
  fMVAVar_ElePreShowerOverRaw = ElePreShowerOverRaw;
  fMVAVar_EleChargedIso03OverPt 
    = (EleChargedIso03 
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleChargedIso03, EleEta)) / ElePt;
  fMVAVar_EleNeutralHadronIso03OverPt 
    = (EleNeutralHadronIso03
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso03, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso007,EleEta)) / ElePt;
  fMVAVar_EleGammaIso03OverPt 
    = (EleGammaIso03 
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIso03, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIsoVetoEtaStrip03,EleEta))/ElePt;
  fMVAVar_EleChargedIso04OverPt 
    = (EleChargedIso04 
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleChargedIso04, EleEta))/ElePt;
  fMVAVar_EleNeutralHadronIso04OverPt
    = (EleNeutralHadronIso04 
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso04, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleNeutralHadronIso007,EleEta))/ElePt;
  fMVAVar_EleGammaIso04OverPt 
    = (EleGammaIso04 
       - Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIso04, EleEta) 
       + Rho * ElectronEffectiveArea(ElectronIDMVA::kEleGammaIsoVetoEtaStrip04,EleEta))/ElePt;
  
  Double_t mva = -9999;  
  TMVA::Reader *reader = 0;
  Int_t MVABin = -1;
  if (subdet == 0 && ptBin == 0) MVABin = 0;
  if (subdet == 1 && ptBin == 0) MVABin = 1;
  if (subdet == 2 && ptBin == 0) MVABin = 2;
  if (subdet == 0 && ptBin == 1) MVABin = 3;
  if (subdet == 1 && ptBin == 1) MVABin = 4;
  if (subdet == 2 && ptBin == 1) MVABin = 5;
  assert(MVABin >= 0 && MVABin <= 5);
  reader = fTMVAReader[MVABin];
                                                
  mva = reader->EvaluateMVA( fMethodname );

  if (printDebug == kTRUE) {
      LogDebug("ElectronIDMVA") << ElePt << " " << EleEta << " " << " --> MVABin " << MVABin << " : "     
                                << fMVAVar_EleSigmaIEtaIEta << " " 
                                << fMVAVar_EleDEtaIn << " " 
                                << fMVAVar_EleDPhiIn << " " 
                                << fMVAVar_EleHoverE << " " 
                                << fMVAVar_EleD0 << " " 
                                << fMVAVar_EleDZ << " " 
                                << fMVAVar_EleFBrem << " " 
                                << fMVAVar_EleEOverP << " " 
                                << fMVAVar_EleESeedClusterOverPout << " " 
                                << fMVAVar_EleSigmaIPhiIPhi << " " 
                                << fMVAVar_EleNBrem << " " 
                                << fMVAVar_EleOneOverEMinusOneOverP << " " 
                                << fMVAVar_EleESeedClusterOverPIn << " " 
                                << fMVAVar_EleIP3d << " " 
                                << fMVAVar_EleIP3dSig << " " 
                                << fMVAVar_EleGsfTrackChi2OverNdof << " "
                                << fMVAVar_EledEtaCalo << " "
                                << fMVAVar_EledPhiCalo << " "
                                << fMVAVar_EleR9 << " "
                                << fMVAVar_EleSCEtaWidth << " "
                                << fMVAVar_EleSCPhiWidth << " "
                                << fMVAVar_EleCovIEtaIPhi << " "
                                << fMVAVar_ElePreShowerOverRaw << " "
                                << fMVAVar_EleChargedIso03OverPt  << " "
                                << fMVAVar_EleNeutralHadronIso03OverPt  << " "
                                << fMVAVar_EleGammaIso03OverPt  << " "
                                << fMVAVar_EleChargedIso04OverPt  << " "
                                << fMVAVar_EleNeutralHadronIso04OverPt  << " "
                                << fMVAVar_EleGammaIso04OverPt  << " "
                                << " === : === "
                                << mva;
  }

  return mva;
}

Member Data Documentation

Bool_t ElectronIDMVA::fIsInitialized

protected

Definition at line 107 of file ElectronIDMVA.h.

Referenced by Initialize(), IsInitialized(), and MVAValue().

TString ElectronIDMVA::fMethodname

protected

Definition at line 106 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleChargedIso03OverPt

protected

Definition at line 132 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleChargedIso04OverPt

protected

Definition at line 135 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleCovIEtaIPhi

protected

Definition at line 130 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleD0

protected

Definition at line 113 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EledEtaCalo

protected

Definition at line 125 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleDEtaIn

protected

Definition at line 110 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EledPhiCalo

protected

Definition at line 126 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleDPhiIn

protected

Definition at line 111 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleDZ

protected

Definition at line 114 of file ElectronIDMVA.h.

Referenced by MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleEOverP

protected

Definition at line 116 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleESeedClusterOverPIn

protected

Definition at line 121 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleESeedClusterOverPout

protected

Definition at line 117 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleFBrem

protected

Definition at line 115 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleGammaIso03OverPt

protected

Definition at line 134 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleGammaIso04OverPt

protected

Definition at line 137 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleGsfTrackChi2OverNdof

protected

Definition at line 124 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleHoverE

protected

Definition at line 112 of file ElectronIDMVA.h.

Referenced by MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleIP3d

protected

Definition at line 122 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleIP3dSig

protected

Definition at line 123 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleNBrem

protected

Definition at line 119 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleNeutralHadronIso03OverPt

protected

Definition at line 133 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleNeutralHadronIso04OverPt

protected

Definition at line 136 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleOneOverEMinusOneOverP

protected

Definition at line 120 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_ElePreShowerOverRaw

protected

Definition at line 131 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleR9

protected

Definition at line 127 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleSCEtaWidth

protected

Definition at line 128 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleSCPhiWidth

protected

Definition at line 129 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleSigmaIEtaIEta

protected

Definition at line 109 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t ElectronIDMVA::fMVAVar_EleSigmaIPhiIPhi

protected

Definition at line 118 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().

TMVA::Reader* ElectronIDMVA::fTMVAReader[6]

protected

Definition at line 105 of file ElectronIDMVA.h.

Referenced by ElectronIDMVA(), Initialize(), MVAValue(), and ~ElectronIDMVA().

unsigned int ElectronIDMVA::version_

protected

Definition at line 139 of file ElectronIDMVA.h.

Referenced by Initialize(), and MVAValue().