MuonIDMVA Class Reference

#include <MuonIDMVA.h>

Public Types
enum	EMuonEffectiveAreaType {   kMuChargedIso03, kMuNeutralIso03, kMuChargedIso04, kMuNeutralIso04,   kMuHadEnergy, kMuHoEnergy, kMuEmEnergy, kMuHadS9Energy,   kMuHoS9Energy, kMuEmS9Energy, kMuTrkIso03, kMuEMIso03,   kMuHadIso03, kMuTrkIso05, kMuEMIso05, kMuHadIso05 }

Public Member Functions
void	Initialize (TString methodName, unsigned int version, TString Subdet0Pt10To14p5Weights, TString Subdet1Pt10To14p5Weights, TString Subdet0Pt14p5To20Weights, TString Subdet1Pt14p5To20Weights, TString Subdet0Pt20ToInfWeights, TString Subdet1Pt20ToInfWeights)
Bool_t	IsInitialized () const
	MuonIDMVA ()
Double_t	MVAValue (HWW &, const unsigned int mu, const unsigned int vertex)
Double_t	MVAValue (Double_t MuPt, Double_t MuEta, Double_t MuTkNchi2, Double_t MuGlobalNchi2, Double_t MuNValidHits, Double_t MuNTrackerHits, Double_t MuNPixelHits, Double_t MuNMatches, Double_t MuD0, Double_t MuIP3d, Double_t MuIP3dSig, Double_t MuTrkKink, Double_t MuSegmentCompatibility, Double_t MuCaloCompatibility, Double_t MuHadEnergyOverPt, Double_t MuHoEnergyOverPt, Double_t MuEmEnergyOverPt, Double_t MuHadS9EnergyOverPt, Double_t MuHoS9EnergyOverPt, Double_t MuEmS9EnergyOverPt, Double_t MuTrkIso03OverPt, Double_t MuEMIso03OverPt, Double_t MuHadIso03OverPt, Double_t MuTrkIso05OverPt, Double_t MuEMIso05OverPt, Double_t MuHadIso05OverPt, Bool_t printDebug=kFALSE)
	~MuonIDMVA ()

Static Public Member Functions
static const Double_t	MuonEffectiveArea (EMuonEffectiveAreaType type, Double_t Eta)

Protected Attributes
Bool_t	fIsInitialized
TString	fMethodname
Float_t	fMVAVar_MuCaloCompatibility
Float_t	fMVAVar_MuChargedIso03OverPt
Float_t	fMVAVar_MuChargedIso04OverPt
Float_t	fMVAVar_MuD0
Float_t	fMVAVar_MuEmEnergyOverPt
Float_t	fMVAVar_MuEMIso03OverPt
Float_t	fMVAVar_MuEMIso05OverPt
Float_t	fMVAVar_MuEmS9EnergyOverPt
Float_t	fMVAVar_MuGlobalNchi2
Float_t	fMVAVar_MuHadEnergyOverPt
Float_t	fMVAVar_MuHadIso03OverPt
Float_t	fMVAVar_MuHadIso05OverPt
Float_t	fMVAVar_MuHadS9EnergyOverPt
Float_t	fMVAVar_MuHoEnergyOverPt
Float_t	fMVAVar_MuHoS9EnergyOverPt
Float_t	fMVAVar_MuIP3d
Float_t	fMVAVar_MuIP3dSig
Float_t	fMVAVar_MuNeutralIso03OverPt
Float_t	fMVAVar_MuNeutralIso04OverPt
Float_t	fMVAVar_MuNMatches
Float_t	fMVAVar_MuNPixelHits
Float_t	fMVAVar_MuNTrackerHits
Float_t	fMVAVar_MuNValidHits
Float_t	fMVAVar_MuSegmentCompatibility
Float_t	fMVAVar_MuTkNchi2
Float_t	fMVAVar_MuTrkIso03OverPt
Float_t	fMVAVar_MuTrkIso05OverPt
Float_t	fMVAVar_MuTrkKink
TMVA::Reader *	fTMVAReader [6]

Detailed Description

Definition at line 24 of file MuonIDMVA.h.

Member Enumeration Documentation

enum MuonIDMVA::EMuonEffectiveAreaType

Enumerator
kMuChargedIso03
kMuNeutralIso03
kMuChargedIso04
kMuNeutralIso04
kMuHadEnergy
kMuHoEnergy
kMuEmEnergy
kMuHadS9Energy
kMuHoS9Energy
kMuEmS9Energy
kMuTrkIso03
kMuEMIso03
kMuHadIso03
kMuTrkIso05
kMuEMIso05
kMuHadIso05

Definition at line 68 of file MuonIDMVA.h.

                                {
      kMuChargedIso03, 
      kMuNeutralIso03, 
      kMuChargedIso04, 
      kMuNeutralIso04, 
      kMuHadEnergy, 
      kMuHoEnergy, 
      kMuEmEnergy, 
      kMuHadS9Energy, 
      kMuHoS9Energy, 
      kMuEmS9Energy,
      kMuTrkIso03, 
      kMuEMIso03, 
      kMuHadIso03, 
      kMuTrkIso05, 
      kMuEMIso05, 
      kMuHadIso05 
    };

Constructor & Destructor Documentation

MuonIDMVA::MuonIDMVA

(

)

Definition at line 15 of file MuonIDMVA.cc.

References fTMVAReader, and i.

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

MuonIDMVA::~MuonIDMVA

(

)

Definition at line 27 of file MuonIDMVA.cc.

References fTMVAReader, and i.

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

Member Function Documentation

void MuonIDMVA::Initialize	(	TString	methodName,
		unsigned int	version,
		TString	Subdet0Pt10To14p5Weights,
		TString	Subdet1Pt10To14p5Weights,
		TString	Subdet0Pt14p5To20Weights,
		TString	Subdet1Pt14p5To20Weights,
		TString	Subdet0Pt20ToInfWeights,
		TString	Subdet1Pt20ToInfWeights
	)

Definition at line 35 of file MuonIDMVA.cc.

References fIsInitialized, fMethodname, fMVAVar_MuCaloCompatibility, fMVAVar_MuD0, fMVAVar_MuEmEnergyOverPt, fMVAVar_MuEMIso03OverPt, fMVAVar_MuEMIso05OverPt, fMVAVar_MuEmS9EnergyOverPt, fMVAVar_MuGlobalNchi2, fMVAVar_MuHadEnergyOverPt, fMVAVar_MuHadIso03OverPt, fMVAVar_MuHadIso05OverPt, fMVAVar_MuHadS9EnergyOverPt, fMVAVar_MuIP3d, fMVAVar_MuIP3dSig, fMVAVar_MuNMatches, fMVAVar_MuNPixelHits, fMVAVar_MuNTrackerHits, fMVAVar_MuNValidHits, fMVAVar_MuSegmentCompatibility, fMVAVar_MuTkNchi2, fMVAVar_MuTrkIso03OverPt, fMVAVar_MuTrkIso05OverPt, fMVAVar_MuTrkKink, fTMVAReader, and i.

                                                             {

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

    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]->AddVariable( "TkNchi2",              &fMVAVar_MuTkNchi2               );
      fTMVAReader[i]->AddVariable( "GlobalNchi2",          &fMVAVar_MuGlobalNchi2           );
      fTMVAReader[i]->AddVariable( "NValidHits",           &fMVAVar_MuNValidHits            );
      fTMVAReader[i]->AddVariable( "NTrackerHits",         &fMVAVar_MuNTrackerHits          );
      fTMVAReader[i]->AddVariable( "NPixelHits",           &fMVAVar_MuNPixelHits            );
      fTMVAReader[i]->AddVariable( "NMatches",             &fMVAVar_MuNMatches              );
      fTMVAReader[i]->AddVariable( "D0",                   &fMVAVar_MuD0                    );      
      fTMVAReader[i]->AddVariable( "IP3d",                 &fMVAVar_MuIP3d                  );      
      fTMVAReader[i]->AddVariable( "IP3dSig",              &fMVAVar_MuIP3dSig               );      
      fTMVAReader[i]->AddVariable( "TrkKink",              &fMVAVar_MuTrkKink               );      
      fTMVAReader[i]->AddVariable( "SegmentCompatibility", &fMVAVar_MuSegmentCompatibility  );      
      fTMVAReader[i]->AddVariable( "CaloCompatibility",    &fMVAVar_MuCaloCompatibility     );      
      fTMVAReader[i]->AddVariable( "HadEnergyOverPt",      &fMVAVar_MuHadEnergyOverPt       );      
      fTMVAReader[i]->AddVariable( "EmEnergyOverPt",       &fMVAVar_MuEmEnergyOverPt        );      
      fTMVAReader[i]->AddVariable( "HadS9EnergyOverPt",    &fMVAVar_MuHadS9EnergyOverPt     );      
      fTMVAReader[i]->AddVariable( "EmS9EnergyOverPt",     &fMVAVar_MuEmS9EnergyOverPt      );      
      fTMVAReader[i]->AddVariable( "TrkIso03OverPt",       &fMVAVar_MuTrkIso03OverPt        );
      fTMVAReader[i]->AddVariable( "EMIso03OverPt",        &fMVAVar_MuEMIso03OverPt         );
      fTMVAReader[i]->AddVariable( "HadIso03OverPt",       &fMVAVar_MuHadIso03OverPt        );
      fTMVAReader[i]->AddVariable( "TrkIso05OverPt",       &fMVAVar_MuTrkIso05OverPt        );
      fTMVAReader[i]->AddVariable( "EMIso05OverPt",        &fMVAVar_MuEMIso05OverPt         );
      fTMVAReader[i]->AddVariable( "HadIso05OverPt",       &fMVAVar_MuHadIso05OverPt        ); 
    }
    
    if (i==0) fTMVAReader[i]->BookMVA(fMethodname , Subdet0Pt10To14p5Weights );
    if (i==1) fTMVAReader[i]->BookMVA(fMethodname , Subdet1Pt10To14p5Weights );
    if (i==2) fTMVAReader[i]->BookMVA(fMethodname , Subdet0Pt14p5To20Weights );
    if (i==3) fTMVAReader[i]->BookMVA(fMethodname , Subdet1Pt14p5To20Weights );
    if (i==4) fTMVAReader[i]->BookMVA(fMethodname , Subdet0Pt20ToInfWeights  );
    if (i==5) fTMVAReader[i]->BookMVA(fMethodname , Subdet1Pt20ToInfWeights  );
    
    }
}

Bool_t MuonIDMVA::IsInitialized ( ) const

inline

Definition at line 37 of file MuonIDMVA.h.

References fIsInitialized.

{ return fIsInitialized; }

static const Double_t MuonIDMVA::MuonEffectiveArea ( EMuonEffectiveAreaType  type, Double_t  Eta  )

inlinestatic

Definition at line 87 of file MuonIDMVA.h.

References kMuChargedIso03, kMuChargedIso04, kMuEmEnergy, kMuEMIso03, kMuEMIso05, kMuEmS9Energy, kMuHadEnergy, kMuHadIso03, kMuHadIso05, kMuHadS9Energy, kMuHoEnergy, kMuHoS9Energy, kMuNeutralIso03, kMuNeutralIso04, kMuTrkIso03, and kMuTrkIso05.

Referenced by MVAValue().

                                                                                     {

    Double_t EffectiveArea = 0;
    if (fabs(Eta) < 1.0) {
      if (type == kMuChargedIso03) EffectiveArea = 0.000;
      if (type == kMuNeutralIso03) EffectiveArea = 0.080;
      if (type == kMuChargedIso04) EffectiveArea = 0.000;
      if (type == kMuNeutralIso04) EffectiveArea = 0.163;
      if (type == kMuHadEnergy)    EffectiveArea = 0.000;
      if (type == kMuHoEnergy)     EffectiveArea = 0.000;
      if (type == kMuEmEnergy)     EffectiveArea = 0.000;
      if (type == kMuHadS9Energy)  EffectiveArea = 0.016;
      if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
      if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
      if (type == kMuTrkIso03)     EffectiveArea = 0.000;
      if (type == kMuEMIso03)      EffectiveArea = 0.080;
      if (type == kMuHadIso03)     EffectiveArea = 0.025;
      if (type == kMuTrkIso05)     EffectiveArea = 0.000;
      if (type == kMuEMIso05)      EffectiveArea = 0.290;
      if (type == kMuHadIso05)     EffectiveArea = 0.091;
    } else if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) {
      if (type == kMuChargedIso03) EffectiveArea = 0.000;
      if (type == kMuNeutralIso03) EffectiveArea = 0.083;
      if (type == kMuChargedIso04) EffectiveArea = 0.000;
      if (type == kMuNeutralIso04) EffectiveArea = 0.168;
      if (type == kMuHadEnergy)    EffectiveArea = 0.005;
      if (type == kMuHoEnergy)     EffectiveArea = 0.000;
      if (type == kMuEmEnergy)     EffectiveArea = 0.000;
      if (type == kMuHadS9Energy)  EffectiveArea = 0.041;
      if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
      if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
      if (type == kMuTrkIso03)     EffectiveArea = 0.000;
      if (type == kMuEMIso03)      EffectiveArea = 0.043;
      if (type == kMuHadIso03)     EffectiveArea = 0.028;
      if (type == kMuTrkIso05)     EffectiveArea = 0.000;
      if (type == kMuEMIso05)      EffectiveArea = 0.184;
      if (type == kMuHadIso05)     EffectiveArea = 0.106;
    } else if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) {
      if (type == kMuChargedIso03) EffectiveArea = 0.000;
      if (type == kMuNeutralIso03) EffectiveArea = 0.060;
      if (type == kMuChargedIso04) EffectiveArea = 0.000;
      if (type == kMuNeutralIso04) EffectiveArea = 0.131;
      if (type == kMuHadEnergy)    EffectiveArea = 0.020;
      if (type == kMuHoEnergy)     EffectiveArea = 0.000;
      if (type == kMuEmEnergy)     EffectiveArea = 0.000;
      if (type == kMuHadS9Energy)  EffectiveArea = 0.072;
      if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
      if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
      if (type == kMuTrkIso03)     EffectiveArea = 0.000;
      if (type == kMuEMIso03)      EffectiveArea = 0.025;
      if (type == kMuHadIso03)     EffectiveArea = 0.036;
      if (type == kMuTrkIso05)     EffectiveArea = 0.000;
      if (type == kMuEMIso05)      EffectiveArea = 0.124;
      if (type == kMuHadIso05)     EffectiveArea = 0.140;
    } else if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.25 ) {
      if (type == kMuChargedIso03) EffectiveArea = 0.000;
      if (type == kMuNeutralIso03) EffectiveArea = 0.066;
      if (type == kMuChargedIso04) EffectiveArea = 0.000;
      if (type == kMuNeutralIso04) EffectiveArea = 0.149;
      if (type == kMuHadEnergy)    EffectiveArea = 0.056;
      if (type == kMuHoEnergy)     EffectiveArea = 0.000;
      if (type == kMuEmEnergy)     EffectiveArea = 0.000;
      if (type == kMuHadS9Energy)  EffectiveArea = 0.148;
      if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
      if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
      if (type == kMuTrkIso03)     EffectiveArea = 0.000;
      if (type == kMuEMIso03)      EffectiveArea = 0.025;
      if (type == kMuHadIso03)     EffectiveArea = 0.050;
      if (type == kMuTrkIso05)     EffectiveArea = 0.000;
      if (type == kMuEMIso05)      EffectiveArea = 0.120;
      if (type == kMuHadIso05)     EffectiveArea = 0.186;
    } else if (fabs(Eta) >= 2.25 && fabs(Eta) < 2.4 ) {
      if (type == kMuChargedIso03) EffectiveArea = 0.000;
      if (type == kMuNeutralIso03) EffectiveArea = 0.098;
      if (type == kMuChargedIso04) EffectiveArea = 0.000;
      if (type == kMuNeutralIso04) EffectiveArea = 0.200;
      if (type == kMuHadEnergy)    EffectiveArea = 0.093;
      if (type == kMuHoEnergy)     EffectiveArea = 0.000;
      if (type == kMuEmEnergy)     EffectiveArea = 0.000;
      if (type == kMuHadS9Energy)  EffectiveArea = 0.260;
      if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
      if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
      if (type == kMuTrkIso03)     EffectiveArea = 0.000;
      if (type == kMuEMIso03)      EffectiveArea = 0.027;
      if (type == kMuHadIso03)     EffectiveArea = 0.060;
      if (type == kMuTrkIso05)     EffectiveArea = 0.000;
      if (type == kMuEMIso05)      EffectiveArea = 0.139;
      if (type == kMuHadIso05)     EffectiveArea = 0.228;
    }
    return EffectiveArea;
  }

Double_t MuonIDMVA::MVAValue	(	HWW &	hww,
		const unsigned int	mu,
		const unsigned int	vertex
	)

Definition at line 211 of file MuonIDMVA.cc.

References assert(), HWW::evt_event(), HWW::evt_lumiBlock(), HWW::evt_run(), HWW::evt_ww_rho_vor(), plotBeamSpotDB::first, fIsInitialized, fMethodname, fMVAVar_MuCaloCompatibility, fMVAVar_MuD0, fMVAVar_MuEmEnergyOverPt, fMVAVar_MuEMIso03OverPt, fMVAVar_MuEMIso05OverPt, fMVAVar_MuEmS9EnergyOverPt, fMVAVar_MuGlobalNchi2, fMVAVar_MuHadEnergyOverPt, fMVAVar_MuHadIso03OverPt, fMVAVar_MuHadIso05OverPt, fMVAVar_MuHadS9EnergyOverPt, fMVAVar_MuHoEnergyOverPt, fMVAVar_MuHoS9EnergyOverPt, fMVAVar_MuIP3d, fMVAVar_MuIP3dSig, fMVAVar_MuNMatches, fMVAVar_MuNPixelHits, fMVAVar_MuNTrackerHits, fMVAVar_MuNValidHits, fMVAVar_MuSegmentCompatibility, fMVAVar_MuTkNchi2, fMVAVar_MuTrkIso03OverPt, fMVAVar_MuTrkIso05OverPt, fMVAVar_MuTrkKink, fTMVAReader, kMuEmEnergy, kMuEMIso03, kMuEMIso05, kMuEmS9Energy, kMuHadEnergy, kMuHadIso03, kMuHadIso05, kMuHadS9Energy, kMuHoEnergy, kMuHoS9Energy, kMuTrkIso03, kMuTrkIso05, LogDebug, HWWFunctions::mud0PV_smurfV3(), MuonEffectiveArea(), HWW::mus_caloCompatibility(), HWW::mus_chi2(), HWW::mus_e_em(), HWW::mus_e_emS9(), HWW::mus_e_had(), HWW::mus_e_hadS9(), HWW::mus_e_ho(), HWW::mus_e_hoS9(), HWW::mus_gfit_chi2(), HWW::mus_gfit_ndof(), HWW::mus_gfit_validSTAHits(), HWW::mus_ip3d(), HWW::mus_ip3derr(), HWW::mus_iso03_emEt(), HWW::mus_iso03_hadEt(), HWW::mus_iso03_sumPt(), HWW::mus_iso05_emEt(), HWW::mus_iso05_hadEt(), HWW::mus_iso05_sumPt(), HWW::mus_ndof(), HWW::mus_nmatches(), HWW::mus_p4(), HWW::mus_segmCompatibility(), HWW::mus_trk_p4(), HWW::mus_trkidx(), HWW::mus_trkKink(), HWW::mus_validHits(), matplotRender::reader, HWWFunctions::trks_d0_pv(), and HWW::trks_valid_pixelhits().

Referenced by HWWFunctions::goodMuonTMVA().

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

  Double_t MuPt  = hww.mus_trk_p4().at(mu).pt();
  Double_t MuEta = hww.mus_trk_p4().at(mu).eta();
  Double_t Rho = hww.evt_ww_rho_vor();

  Int_t subdet = 0;
  if (fabs(MuEta) < 1.479) subdet = 0;
  else subdet = 1;
  Int_t ptBin = 0;
  if (MuPt > 14.5) ptBin = 1;
  if (MuPt > 20.0) ptBin = 2;

  //set all input variables
  fMVAVar_MuTkNchi2              = hww.mus_chi2().at(mu) / hww.mus_ndof().at(mu);
  fMVAVar_MuGlobalNchi2          = hww.mus_gfit_chi2().at(mu) / hww.mus_gfit_ndof().at(mu);
  fMVAVar_MuNValidHits           = hww.mus_gfit_validSTAHits().at(mu);
  fMVAVar_MuNTrackerHits         = hww.mus_validHits().at(mu);
  fMVAVar_MuNPixelHits           = hww.trks_valid_pixelhits().at(hww.mus_trkidx().at(mu));
  fMVAVar_MuNMatches             = hww.mus_nmatches().at(mu);
  fMVAVar_MuD0                   = mud0PV_smurfV3(hww, mu);

  const double mud0sign   = ( (trks_d0_pv(hww, hww.mus_trkidx().at(mu),0).first)   >=0 ) ? 1. : -1.;
  fMVAVar_MuIP3d =                   hww.mus_ip3d().at(mu)*mud0sign; 
  if (hww.mus_ip3derr().at(mu) == 0.0) fMVAVar_MuIP3dSig = 0.0;
  else fMVAVar_MuIP3dSig =           hww.mus_ip3d().at(mu)*mud0sign / hww.mus_ip3derr().at(mu); 
  fMVAVar_MuTrkKink              = hww.mus_trkKink().at(mu);
  fMVAVar_MuSegmentCompatibility = hww.mus_segmCompatibility().at(mu);
  fMVAVar_MuCaloCompatibility    = hww.mus_caloCompatibility().at(mu);
  fMVAVar_MuHadEnergyOverPt      = (hww.mus_e_had().at(mu)   - Rho*MuonEffectiveArea(MuonIDMVA::kMuHadEnergy,MuEta))/MuPt;
  fMVAVar_MuHoEnergyOverPt       = (hww.mus_e_ho().at(mu)    - Rho*MuonEffectiveArea(MuonIDMVA::kMuHoEnergy,MuEta))/MuPt;
  fMVAVar_MuEmEnergyOverPt       = (hww.mus_e_em().at(mu)    - Rho*MuonEffectiveArea(MuonIDMVA::kMuEmEnergy,MuEta))/MuPt;
  fMVAVar_MuHadS9EnergyOverPt    = (hww.mus_e_hadS9().at(mu) - Rho*MuonEffectiveArea(MuonIDMVA::kMuHadS9Energy,MuEta))/MuPt;
  fMVAVar_MuHoS9EnergyOverPt     = (hww.mus_e_hoS9().at(mu)  - Rho*MuonEffectiveArea(MuonIDMVA::kMuHoS9Energy,MuEta))/MuPt;
  fMVAVar_MuEmS9EnergyOverPt     = (hww.mus_e_emS9().at(mu)  - Rho*MuonEffectiveArea(MuonIDMVA::kMuEmS9Energy,MuEta))/MuPt;
  fMVAVar_MuTrkIso03OverPt       = (hww.mus_iso03_sumPt().at(mu) - Rho*MuonEffectiveArea(MuonIDMVA::kMuTrkIso03,MuEta))/MuPt;
  fMVAVar_MuEMIso03OverPt        = (hww.mus_iso03_emEt().at(mu)  - Rho*MuonEffectiveArea(MuonIDMVA::kMuEMIso03,MuEta))/MuPt;
  fMVAVar_MuHadIso03OverPt       = (hww.mus_iso03_hadEt().at(mu) - Rho*MuonEffectiveArea(MuonIDMVA::kMuHadIso03,MuEta))/MuPt;
  fMVAVar_MuTrkIso05OverPt       = (hww.mus_iso05_sumPt().at(mu) - Rho*MuonEffectiveArea(MuonIDMVA::kMuTrkIso05,MuEta))/MuPt;
  fMVAVar_MuEMIso05OverPt        = (hww.mus_iso05_emEt().at(mu)  - Rho*MuonEffectiveArea(MuonIDMVA::kMuEMIso05,MuEta))/MuPt;
  fMVAVar_MuHadIso05OverPt       = (hww.mus_iso05_hadEt().at(mu) - Rho*MuonEffectiveArea(MuonIDMVA::kMuHadIso05,MuEta))/MuPt;

  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 == 0 && ptBin == 1) MVABin = 2;
  if (subdet == 1 && ptBin == 1) MVABin = 3;
  if (subdet == 0 && ptBin == 2) MVABin = 4;
  if (subdet == 1 && ptBin == 2) MVABin = 5;
  assert(MVABin >= 0 && MVABin <= 5);
  reader = fTMVAReader[MVABin];

  mva = reader->EvaluateMVA( fMethodname );

  if (0) {
    LogDebug("MuonIDMVA") << "Debug Muon MVA: "
                          << hww.evt_run() << " " << hww.evt_lumiBlock() << " " << hww.evt_event() << " " << Rho << " "
                          << hww.mus_p4().at(mu).pt() << " " << hww.mus_p4().at(mu).eta() << " " << hww.mus_p4().at(mu).phi() << " : "
                          << MuPt << " " << MuEta << " --> MVABin " << MVABin << " : "     
                          << fMVAVar_MuTkNchi2              << " " 
                          << fMVAVar_MuGlobalNchi2          << " " 
                          << fMVAVar_MuNValidHits           << " " 
                          << fMVAVar_MuNTrackerHits         << " " 
                          << fMVAVar_MuNPixelHits           << " "  
                          << fMVAVar_MuNMatches             << " " 
                          << fMVAVar_MuD0                   << " " 
                          << fMVAVar_MuIP3d                 << " " 
                          << fMVAVar_MuIP3dSig              << " " 
                          << fMVAVar_MuTrkKink              << " " 
                          << fMVAVar_MuSegmentCompatibility << " " 
                          << fMVAVar_MuCaloCompatibility    << " " 
                          << fMVAVar_MuHadEnergyOverPt      << " " 
                          << fMVAVar_MuHoEnergyOverPt       << " " 
                          << fMVAVar_MuEmEnergyOverPt       << " " 
                          << fMVAVar_MuHadS9EnergyOverPt    << " " 
                          << fMVAVar_MuHoS9EnergyOverPt     << " " 
                          << fMVAVar_MuEmS9EnergyOverPt     << " " 
                          << fMVAVar_MuTrkIso03OverPt   << " " 
                          << fMVAVar_MuEMIso03OverPt   << " " 
                          << fMVAVar_MuHadIso03OverPt   << " " 
                          << fMVAVar_MuTrkIso05OverPt   << " " 
                          << fMVAVar_MuEMIso05OverPt   << " " 
                          << fMVAVar_MuHadIso05OverPt   << " " 
                          << " === : === "
                          << mva;
  }

  return mva;
}

Double_t MuonIDMVA::MVAValue	(	Double_t	MuPt,
		Double_t	MuEta,
		Double_t	MuTkNchi2,
		Double_t	MuGlobalNchi2,
		Double_t	MuNValidHits,
		Double_t	MuNTrackerHits,
		Double_t	MuNPixelHits,
		Double_t	MuNMatches,
		Double_t	MuD0,
		Double_t	MuIP3d,
		Double_t	MuIP3dSig,
		Double_t	MuTrkKink,
		Double_t	MuSegmentCompatibility,
		Double_t	MuCaloCompatibility,
		Double_t	MuHadEnergyOverPt,
		Double_t	MuHoEnergyOverPt,
		Double_t	MuEmEnergyOverPt,
		Double_t	MuHadS9EnergyOverPt,
		Double_t	MuHoS9EnergyOverPt,
		Double_t	MuEmS9EnergyOverPt,
		Double_t	MuTrkIso03OverPt,
		Double_t	MuEMIso03OverPt,
		Double_t	MuHadIso03OverPt,
		Double_t	MuTrkIso05OverPt,
		Double_t	MuEMIso05OverPt,
		Double_t	MuHadIso05OverPt,
		Bool_t	printDebug = kFALSE
	)

Definition at line 94 of file MuonIDMVA.cc.

References assert(), fIsInitialized, fMethodname, fMVAVar_MuCaloCompatibility, fMVAVar_MuD0, fMVAVar_MuEmEnergyOverPt, fMVAVar_MuEMIso03OverPt, fMVAVar_MuEMIso05OverPt, fMVAVar_MuEmS9EnergyOverPt, fMVAVar_MuGlobalNchi2, fMVAVar_MuHadEnergyOverPt, fMVAVar_MuHadIso03OverPt, fMVAVar_MuHadIso05OverPt, fMVAVar_MuHadS9EnergyOverPt, fMVAVar_MuHoEnergyOverPt, fMVAVar_MuHoS9EnergyOverPt, fMVAVar_MuIP3d, fMVAVar_MuIP3dSig, fMVAVar_MuNMatches, fMVAVar_MuNPixelHits, fMVAVar_MuNTrackerHits, fMVAVar_MuNValidHits, fMVAVar_MuSegmentCompatibility, fMVAVar_MuTkNchi2, fMVAVar_MuTrkIso03OverPt, fMVAVar_MuTrkIso05OverPt, fMVAVar_MuTrkKink, fTMVAReader, LogDebug, and matplotRender::reader.

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

  Int_t subdet = 0;
  if (fabs(MuEta) < 1.479) subdet = 0;
  else subdet = 1;
  Int_t ptBin = 0;
  if (MuPt > 14.5) ptBin = 1;
  if (MuPt > 20.0) ptBin = 2;

  
  //set all input variables
  fMVAVar_MuTkNchi2              = MuTkNchi2; 
  fMVAVar_MuGlobalNchi2          = MuGlobalNchi2; 
  fMVAVar_MuNValidHits           = MuNValidHits; 
  fMVAVar_MuNTrackerHits         = MuNTrackerHits; 
  fMVAVar_MuNPixelHits           = MuNPixelHits;  
  fMVAVar_MuNMatches             = MuNMatches; 
  fMVAVar_MuD0                   = MuD0; 
  fMVAVar_MuIP3d                 = MuIP3d; 
  fMVAVar_MuIP3dSig              = MuIP3dSig; 
  fMVAVar_MuTrkKink              = MuTrkKink; 
  fMVAVar_MuSegmentCompatibility = MuSegmentCompatibility; 
  fMVAVar_MuCaloCompatibility    = MuCaloCompatibility; 
  fMVAVar_MuHadEnergyOverPt      = MuHadEnergyOverPt; 
  fMVAVar_MuHoEnergyOverPt       = MuHoEnergyOverPt; 
  fMVAVar_MuEmEnergyOverPt       = MuEmEnergyOverPt; 
  fMVAVar_MuHadS9EnergyOverPt    = MuHadS9EnergyOverPt; 
  fMVAVar_MuHoS9EnergyOverPt     = MuHoS9EnergyOverPt; 
  fMVAVar_MuEmS9EnergyOverPt     = MuEmS9EnergyOverPt; 
  fMVAVar_MuTrkIso03OverPt       = MuTrkIso03OverPt; 
  fMVAVar_MuEMIso03OverPt        = MuEMIso03OverPt; 
  fMVAVar_MuHadIso03OverPt       = MuHadIso03OverPt; 
  fMVAVar_MuTrkIso05OverPt       = MuTrkIso05OverPt; 
  fMVAVar_MuEMIso05OverPt        = MuEMIso05OverPt; 
  fMVAVar_MuHadIso05OverPt       = MuHadIso05OverPt; 

  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 == 0 && ptBin == 1) MVABin = 2;
  if (subdet == 1 && ptBin == 1) MVABin = 3;
  if (subdet == 0 && ptBin == 2) MVABin = 4;
  if (subdet == 1 && ptBin == 2) MVABin = 5;
  assert(MVABin >= 0 && MVABin <= 5);
  reader = fTMVAReader[MVABin];
                                                
  mva = reader->EvaluateMVA( fMethodname );

  if (printDebug) {
    LogDebug("MuonIDMVA") << "Debug Muon MVA: "
                            << MuPt << " " << MuEta << " --> MVABin " << MVABin << " : "     
                            << fMVAVar_MuTkNchi2              << " " 
                            << fMVAVar_MuGlobalNchi2          << " " 
                            << fMVAVar_MuNValidHits           << " " 
                            << fMVAVar_MuNTrackerHits         << " " 
                            << fMVAVar_MuNPixelHits           << " "  
                            << fMVAVar_MuNMatches             << " " 
                            << fMVAVar_MuD0                   << " " 
                            << fMVAVar_MuIP3d                 << " " 
                            << fMVAVar_MuIP3dSig              << " " 
                            << fMVAVar_MuTrkKink              << " " 
                            << fMVAVar_MuSegmentCompatibility << " " 
                            << fMVAVar_MuCaloCompatibility    << " " 
                            << fMVAVar_MuHadEnergyOverPt      << " " 
                            << fMVAVar_MuHoEnergyOverPt       << " " 
                            << fMVAVar_MuEmEnergyOverPt       << " " 
                            << fMVAVar_MuHadS9EnergyOverPt    << " " 
                            << fMVAVar_MuHoS9EnergyOverPt     << " " 
                            << fMVAVar_MuEmS9EnergyOverPt     << " " 
                            << fMVAVar_MuTrkIso03OverPt   << " " 
                            << fMVAVar_MuEMIso03OverPt   << " " 
                            << fMVAVar_MuHadIso03OverPt   << " " 
                            << fMVAVar_MuTrkIso05OverPt   << " " 
                            << fMVAVar_MuEMIso05OverPt   << " " 
                            << fMVAVar_MuHadIso05OverPt   << " " 
                            << " === : === "
                            << mva;
  }

  return mva;
}

Member Data Documentation

Bool_t MuonIDMVA::fIsInitialized

protected

Definition at line 183 of file MuonIDMVA.h.

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

TString MuonIDMVA::fMethodname

protected

Definition at line 182 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuCaloCompatibility

protected

Definition at line 196 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuChargedIso03OverPt

protected

Definition at line 203 of file MuonIDMVA.h.

Float_t MuonIDMVA::fMVAVar_MuChargedIso04OverPt

protected

Definition at line 205 of file MuonIDMVA.h.

Float_t MuonIDMVA::fMVAVar_MuD0

protected

Definition at line 191 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuEmEnergyOverPt

protected

Definition at line 199 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuEMIso03OverPt

protected

Definition at line 208 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuEMIso05OverPt

protected

Definition at line 211 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuEmS9EnergyOverPt

protected

Definition at line 202 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuGlobalNchi2

protected

Definition at line 186 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuHadEnergyOverPt

protected

Definition at line 197 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuHadIso03OverPt

protected

Definition at line 209 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuHadIso05OverPt

protected

Definition at line 212 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuHadS9EnergyOverPt

protected

Definition at line 200 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuHoEnergyOverPt

protected

Definition at line 198 of file MuonIDMVA.h.

Referenced by MVAValue().

Float_t MuonIDMVA::fMVAVar_MuHoS9EnergyOverPt

protected

Definition at line 201 of file MuonIDMVA.h.

Referenced by MVAValue().

Float_t MuonIDMVA::fMVAVar_MuIP3d

protected

Definition at line 192 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuIP3dSig

protected

Definition at line 193 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuNeutralIso03OverPt

protected

Definition at line 204 of file MuonIDMVA.h.

Float_t MuonIDMVA::fMVAVar_MuNeutralIso04OverPt

protected

Definition at line 206 of file MuonIDMVA.h.

Float_t MuonIDMVA::fMVAVar_MuNMatches

protected

Definition at line 190 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuNPixelHits

protected

Definition at line 189 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuNTrackerHits

protected

Definition at line 188 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuNValidHits

protected

Definition at line 187 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuSegmentCompatibility

protected

Definition at line 195 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuTkNchi2

protected

Definition at line 185 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuTrkIso03OverPt

protected

Definition at line 207 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuTrkIso05OverPt

protected

Definition at line 210 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

Float_t MuonIDMVA::fMVAVar_MuTrkKink

protected

Definition at line 194 of file MuonIDMVA.h.

Referenced by Initialize(), and MVAValue().

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

protected

Definition at line 181 of file MuonIDMVA.h.

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