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

#include <EGammaMvaEleEstimator.h>

Public Types

enum  MVAType {
  kTrig = 0, kTrigNoIP = 1, kNonTrig = 2, kIsoRings,
  kTrigIDIsoCombined, kTrigIDIsoCombinedPUCorrected
}
 

Public Member Functions

void bindVariables ()
 
 EGammaMvaEleEstimator ()
 
UInt_t GetMVABin (double eta, double pt) const
 
Double_t IDIsoCombinedMvaValue (const reco::GsfElectron &ele, const reco::Vertex &vertex, const TransientTrackBuilder &transientTrackBuilder, EcalClusterLazyTools myEcalCluster, const reco::PFCandidateCollection &PFCandidates, double Rho, ElectronEffectiveArea::ElectronEffectiveAreaTarget EATarget, bool printDebug=kFALSE)
 
Double_t IDIsoCombinedMvaValue (Double_t fbrem, Double_t kfchi2, Int_t kfhits, Double_t gsfchi2, Double_t deta, Double_t dphi, Double_t detacalo, Double_t see, Double_t spp, Double_t etawidth, Double_t phiwidth, Double_t OneMinusE1x5E5x5, Double_t R9, Double_t HoE, Double_t EoP, Double_t IoEmIoP, Double_t eleEoPout, Double_t PreShowerOverRaw, Double_t d0, Double_t ip3d, Double_t ChargedIso_DR0p0To0p1, Double_t ChargedIso_DR0p1To0p2, Double_t ChargedIso_DR0p2To0p3, Double_t ChargedIso_DR0p3To0p4, Double_t ChargedIso_DR0p4To0p5, Double_t GammaIso_DR0p0To0p1, Double_t GammaIso_DR0p1To0p2, Double_t GammaIso_DR0p2To0p3, Double_t GammaIso_DR0p3To0p4, Double_t GammaIso_DR0p4To0p5, Double_t NeutralHadronIso_DR0p0To0p1, Double_t NeutralHadronIso_DR0p1To0p2, Double_t NeutralHadronIso_DR0p2To0p3, Double_t NeutralHadronIso_DR0p3To0p4, Double_t NeutralHadronIso_DR0p4To0p5, Double_t Rho, Double_t eta, Double_t pt, Bool_t printDebug=kFALSE)
 
void initialize (std::string methodName, std::string weightsfile, EGammaMvaEleEstimator::MVAType type)
 
void initialize (std::string methodName, EGammaMvaEleEstimator::MVAType type, Bool_t useBinnedVersion, std::vector< std::string > weightsfiles)
 
Bool_t isInitialized () const
 
Double_t isoMvaValue (const reco::GsfElectron &ele, const reco::Vertex &vertex, const reco::PFCandidateCollection &PFCandidates, double Rho, ElectronEffectiveArea::ElectronEffectiveAreaTarget EATarget, const reco::GsfElectronCollection &IdentifiedElectrons, const reco::MuonCollection &IdentifiedMuons, bool printDebug=kFALSE)
 
Double_t isoMvaValue (Double_t Pt, Double_t Eta, Double_t Rho, ElectronEffectiveArea::ElectronEffectiveAreaTarget EATarget, Double_t ChargedIso_DR0p0To0p1, Double_t ChargedIso_DR0p1To0p2, Double_t ChargedIso_DR0p2To0p3, Double_t ChargedIso_DR0p3To0p4, Double_t ChargedIso_DR0p4To0p5, Double_t GammaIso_DR0p0To0p1, Double_t GammaIso_DR0p1To0p2, Double_t GammaIso_DR0p2To0p3, Double_t GammaIso_DR0p3To0p4, Double_t GammaIso_DR0p4To0p5, Double_t NeutralHadronIso_DR0p0To0p1, Double_t NeutralHadronIso_DR0p1To0p2, Double_t NeutralHadronIso_DR0p2To0p3, Double_t NeutralHadronIso_DR0p3To0p4, Double_t NeutralHadronIso_DR0p4To0p5, Bool_t printDebug=kFALSE)
 
Double_t mvaValue (const reco::GsfElectron &ele, const reco::Vertex &vertex, const TransientTrackBuilder &transientTrackBuilder, EcalClusterLazyTools myEcalCluster, bool printDebug=kFALSE)
 
Double_t mvaValue (const reco::GsfElectron &ele, const reco::Vertex &vertex, double rho, EcalClusterLazyTools myEcalCluster, bool printDebug=kFALSE)
 
Double_t mvaValue (const pat::Electron &ele, double rho, bool printDebug=kFALSE)
 
Double_t mvaValue (const pat::Electron &ele, const reco::Vertex &vertex, double rho, bool useFull5x5=kFALSE, bool printDebug=kFALSE)
 
Double_t mvaValue (Double_t fbrem, Double_t kfchi2, Int_t kfhits, Double_t gsfchi2, Double_t deta, Double_t dphi, Double_t detacalo, Double_t see, Double_t spp, Double_t etawidth, Double_t phiwidth, Double_t e1x5e5x5, Double_t R9, Double_t HoE, Double_t EoP, Double_t IoEmIoP, Double_t eleEoPout, Double_t PreShowerOverRaw, Double_t d0, Double_t ip3d, Double_t eta, Double_t pt, Bool_t printDebug=kFALSE)
 
Double_t mvaValue (Double_t fbrem, Double_t kfchi2, Int_t kfhits, Double_t gsfchi2, Double_t deta, Double_t dphi, Double_t detacalo, Double_t see, Double_t spp, Double_t etawidth, Double_t phiwidth, Double_t e1x5e5x5, Double_t R9, Double_t HoE, Double_t EoP, Double_t IoEmIoP, Double_t eleEoPout, Double_t rho, Double_t PreShowerOverRaw, Double_t eta, Double_t pt, Bool_t printDebug=kFALSE)
 
Double_t mvaValue (Double_t fbrem, Double_t kfchi2, Int_t kfhits, Double_t gsfchi2, Double_t deta, Double_t dphi, Double_t detacalo, Double_t see, Double_t spp, Double_t etawidth, Double_t phiwidth, Double_t e1x5e5x5, Double_t R9, Double_t HoE, Double_t EoP, Double_t IoEmIoP, Double_t eleEoPout, Double_t PreShowerOverRaw, Double_t eta, Double_t pt, Bool_t printDebug=kFALSE)
 
 ~EGammaMvaEleEstimator ()
 

Private Attributes

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

Detailed Description

–> NOTE if you want to use this class as standalone without the CMSSW part you need to uncomment the below line and compile normally with scramv1 b Then you need just to load it in your root macro the lib with the correct path, eg: gSystem->Load("/data/benedet/CMSSW_5_2_2/lib/slc5_amd64_gcc462/pluginEGammaEGammaAnalysisTools.so");

Definition at line 39 of file EGammaMvaEleEstimator.h.

Member Enumeration Documentation

enum EGammaMvaEleEstimator::MVAType
Enumerator
kTrig 
kTrigNoIP 
kNonTrig 
kIsoRings 
kTrigIDIsoCombined 
kTrigIDIsoCombinedPUCorrected 

Definition at line 44 of file EGammaMvaEleEstimator.h.

44  {
45  kTrig = 0, // MVA for triggering electrons
46  kTrigNoIP = 1, // MVA for triggering electrons without IP info
47  kNonTrig = 2, // MVA for non-triggering electrons
48  kIsoRings, // Isolation MVA for non-trigger electrons
49  kTrigIDIsoCombined, // ID+Iso Combined MVA for triggering electrons
50  kTrigIDIsoCombinedPUCorrected // ID+Iso Combined MVA for triggering electrons
51  };

Constructor & Destructor Documentation

EGammaMvaEleEstimator::EGammaMvaEleEstimator ( )

Definition at line 29 of file EGammaMvaEleEstimator.cc.

29  :
30 fMethodname("BDTG method"),
31 fisInitialized(kFALSE),
32 fMVAType(kTrig),
33 fUseBinnedVersion(kTRUE),
34 fNMVABins(0)
35 {
36  // Constructor.
37 }
EGammaMvaEleEstimator::~EGammaMvaEleEstimator ( )

Definition at line 40 of file EGammaMvaEleEstimator.cc.

References fTMVAMethod, fTMVAReader, and i.

41 {
42  for (unsigned int i=0;i<fTMVAReader.size(); ++i) {
43  if (fTMVAMethod[i]) delete fTMVAMethod[i];
44  if (fTMVAReader[i]) delete fTMVAReader[i];
45  }
46 }
i
int i
Definition: DBlmapReader.cc:9
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250

Member Function Documentation

void EGammaMvaEleEstimator::bindVariables ( )

Definition at line 1947 of file EGammaMvaEleEstimator.cc.

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

Referenced by isoMvaValue(), and mvaValue().

1947  {
1948 
1949  // this binding is needed for variables that sometime diverge.
1950 
1951 
1952  if(fMVAVar_fbrem < -1.)
1953  fMVAVar_fbrem = -1.;
1954 
1955  fMVAVar_deta = fabs(fMVAVar_deta);
1956  if(fMVAVar_deta > 0.06)
1957  fMVAVar_deta = 0.06;
1958 
1959 
1960  fMVAVar_dphi = fabs(fMVAVar_dphi);
1961  if(fMVAVar_dphi > 0.6)
1962  fMVAVar_dphi = 0.6;
1963 
1964 
1965  if(fMVAVar_EoP > 20.)
1966  fMVAVar_EoP = 20.;
1967 
1968  if(fMVAVar_eleEoPout > 20.)
1969  fMVAVar_eleEoPout = 20.;
1970 
1971 
1972  fMVAVar_detacalo = fabs(fMVAVar_detacalo);
1973  if(fMVAVar_detacalo > 0.2)
1974  fMVAVar_detacalo = 0.2;
1975 
1976  if(fMVAVar_OneMinusE1x5E5x5 < -1.)
1977  fMVAVar_OneMinusE1x5E5x5 = -1;
1978 
1979  if(fMVAVar_OneMinusE1x5E5x5 > 2.)
1980  fMVAVar_OneMinusE1x5E5x5 = 2.;
1981 
1982 
1983 
1984  if(fMVAVar_R9 > 5)
1985  fMVAVar_R9 = 5;
1986 
1987  if(fMVAVar_gsfchi2 > 200.)
1988  fMVAVar_gsfchi2 = 200;
1989 
1990 
1991  if(fMVAVar_kfchi2 > 10.)
1992  fMVAVar_kfchi2 = 10.;
1993 
1994 
1995  // Needed for a bug in CMSSW_420, fixed in more recent CMSSW versions
1996  if(std::isnan(fMVAVar_spp))
1997  fMVAVar_spp = 0.;
1998 
1999 
2000  return;
2001 }
edm::detail::isnan
bool isnan(float x)
Definition: math.h:13
UInt_t EGammaMvaEleEstimator::GetMVABin ( double  eta,
double  pt 
) const

Definition at line 411 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, fMVAType, kIsoRings, kNonTrig, kTrig, kTrigIDIsoCombined, kTrigIDIsoCombinedPUCorrected, and kTrigNoIP.

Referenced by IDIsoCombinedMvaValue(), isoMvaValue(), and mvaValue().

411  {
412 
413  //Default is to return the first bin
414  unsigned int bin = 0;
415 
416  if (fMVAType == EGammaMvaEleEstimator::kIsoRings) {
417  if (pt < 10 && fabs(eta) < 1.479) bin = 0;
418  if (pt < 10 && fabs(eta) >= 1.479) bin = 1;
419  if (pt >= 10 && fabs(eta) < 1.479) bin = 2;
420  if (pt >= 10 && fabs(eta) >= 1.479) bin = 3;
421  }
422 
423  if (fMVAType == EGammaMvaEleEstimator::kNonTrig ) {
424  bin = 0;
425  if (pt < 10 && fabs(eta) < 0.8) bin = 0;
426  if (pt < 10 && fabs(eta) >= 0.8 && fabs(eta) < 1.479 ) bin = 1;
427  if (pt < 10 && fabs(eta) >= 1.479) bin = 2;
428  if (pt >= 10 && fabs(eta) < 0.8) bin = 3;
429  if (pt >= 10 && fabs(eta) >= 0.8 && fabs(eta) < 1.479 ) bin = 4;
430  if (pt >= 10 && fabs(eta) >= 1.479) bin = 5;
431  }
432 
433 
434  if (fMVAType == EGammaMvaEleEstimator::kTrig ||
435  fMVAType == EGammaMvaEleEstimator::kTrigIDIsoCombined ||
436  fMVAType == EGammaMvaEleEstimator::kTrigIDIsoCombinedPUCorrected ||
437  fMVAType == EGammaMvaEleEstimator::kTrigNoIP
438  ) {
439  bin = 0;
440  if (pt < 20 && fabs(eta) < 0.8) bin = 0;
441  if (pt < 20 && fabs(eta) >= 0.8 && fabs(eta) < 1.479 ) bin = 1;
442  if (pt < 20 && fabs(eta) >= 1.479) bin = 2;
443  if (pt >= 20 && fabs(eta) < 0.8) bin = 3;
444  if (pt >= 20 && fabs(eta) >= 0.8 && fabs(eta) < 1.479 ) bin = 4;
445  if (pt >= 20 && fabs(eta) >= 1.479) bin = 5;
446  }
447 
448 
449 
450  return bin;
451 }
Double_t EGammaMvaEleEstimator::IDIsoCombinedMvaValue ( const reco::GsfElectron &  ele,
const reco::Vertex &  vertex,
const TransientTrackBuilder &  transientTrackBuilder,
EcalClusterLazyTools  myEcalCluster,
const reco::PFCandidateCollection &  PFCandidates,
double  Rho,
ElectronEffectiveArea::ElectronEffectiveAreaTarget  EATarget,
bool  printDebug = kFALSE 
)

Definition at line 1684 of file EGammaMvaEleEstimator.cc.

References IPTools::absoluteImpactParameter3D(), newFWLiteAna::bin, TransientTrackBuilder::build(), reco::GsfElectron::closestCtfTrackRef(), funct::cos(), gather_cfg::cout, reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), reco::PFCandidate::e, reco::GsfElectron::e1x5(), reco::GsfElectron::e5x5(), reco::GsfElectron::eEleClusterOverPout(), reco::GsfElectron::eSuperClusterOverP(), reco::LeafCandidate::eta(), reco::GsfElectron::fbrem(), fisInitialized, fMethodname, fMVAType, fMVAVar_ChargedIso_DR0p0To0p1, fMVAVar_ChargedIso_DR0p1To0p2, fMVAVar_ChargedIso_DR0p2To0p3, fMVAVar_ChargedIso_DR0p3To0p4, fMVAVar_ChargedIso_DR0p4To0p5, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_GammaIso_DR0p0To0p1, fMVAVar_GammaIso_DR0p1To0p2, fMVAVar_GammaIso_DR0p2To0p3, fMVAVar_GammaIso_DR0p3To0p4, fMVAVar_GammaIso_DR0p4To0p5, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_ip3d, fMVAVar_ip3dSig, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_kfhitsall, fMVAVar_NeutralHadronIso_DR0p0To0p1, fMVAVar_NeutralHadronIso_DR0p1To0p2, fMVAVar_NeutralHadronIso_DR0p2To0p3, fMVAVar_NeutralHadronIso_DR0p3To0p4, fMVAVar_NeutralHadronIso_DR0p4To0p5, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_rho, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, reco::PFCandidate::gamma, ElectronEffectiveArea::GetElectronEffectiveArea(), GetMVABin(), reco::GsfElectron::gsfTrack(), reco::GsfElectron::hadronicOverEm(), edm::Ref< C, T, F >::isAvailable(), edm::detail::isnan(), edm::Ref< C, T, F >::isNonnull(), ElectronEffectiveArea::kEleGammaIsoDR0p0To0p1, ElectronEffectiveArea::kEleGammaIsoDR0p1To0p2, ElectronEffectiveArea::kEleGammaIsoDR0p2To0p3, ElectronEffectiveArea::kEleGammaIsoDR0p3To0p4, ElectronEffectiveArea::kEleGammaIsoDR0p4To0p5, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p0To0p1, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p1To0p2, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p2To0p3, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p3To0p4, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p4To0p5, kTrig, kTrigIDIsoCombined, kTrigIDIsoCombinedPUCorrected, Max(), bookConverter::max, Min(), min(), reco::PFCandidate::mu, reco::LeafCandidate::phi(), reco::Vertex::position(), funct::pow(), reco::LeafCandidate::pt(), reco::GsfElectron::sigmaIetaIeta(), mathSSE::sqrt(), reco::GsfElectron::superCluster(), reco::GsfElectron::trackMomentumAtVtx(), and groupFilesInBlocks::tt.

1691  {
1692 
1693  if (!fisInitialized) {
1694  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
1695  return -9999;
1696  }
1697 
1698  bool validKF= false;
1699  reco::TrackRef myTrackRef = ele.closestCtfTrackRef();
1700  validKF = (myTrackRef.isAvailable());
1701  validKF = (myTrackRef.isNonnull());
1702 
1703  // Pure tracking variables
1704  fMVAVar_fbrem = (ele.fbrem() < -1. ) ? -1. : ele.fbrem();
1705  fMVAVar_kfchi2 = (validKF) ? myTrackRef->normalizedChi2() : 0 ;
1706  if (fMVAVar_kfchi2 > 10) fMVAVar_kfchi2 = 10;
1707  fMVAVar_kfhits = (validKF) ? myTrackRef->hitPattern().trackerLayersWithMeasurement() : -1. ;
1708  fMVAVar_kfhitsall = (validKF) ? myTrackRef->numberOfValidHits() : -1. ; // save also this in your ntuple as possible alternative
1709  fMVAVar_gsfchi2 = ele.gsfTrack()->normalizedChi2();
1710  if (fMVAVar_gsfchi2 > 200) fMVAVar_gsfchi2 = 200;
1711 
1712 
1713  // Geometrical matchings
1714  fMVAVar_deta = ( fabs(ele.deltaEtaSuperClusterTrackAtVtx()) > 0.06 ) ? 0.06 : fabs(ele.deltaEtaSuperClusterTrackAtVtx());
1715  fMVAVar_dphi = ele.deltaPhiSuperClusterTrackAtVtx();
1716  fMVAVar_detacalo = ele.deltaEtaSeedClusterTrackAtCalo();
1717 
1718 
1719  // Pure ECAL -> shower shapes
1720  fMVAVar_see = ele.sigmaIetaIeta(); //EleSigmaIEtaIEta
1721  std::vector<float> vCov = myEcalCluster.localCovariances(*(ele.superCluster()->seed())) ;
1722  if (!isnan(vCov[2])) fMVAVar_spp = sqrt (vCov[2]); //EleSigmaIPhiIPhi
1723  else fMVAVar_spp = 0.;
1724 
1725  fMVAVar_etawidth = ele.superCluster()->etaWidth();
1726  fMVAVar_phiwidth = ele.superCluster()->phiWidth();
1727  fMVAVar_OneMinusE1x5E5x5 = (ele.e5x5()) !=0. ? 1.-(ele.e1x5()/ele.e5x5()) : -1. ;
1728  fMVAVar_OneMinusE1x5E5x5 = max(min(double(fMVAVar_OneMinusE1x5E5x5),2.0),-1.0);
1729  fMVAVar_R9 = myEcalCluster.e3x3(*(ele.superCluster()->seed())) / ele.superCluster()->rawEnergy();
1730  if (fMVAVar_R9 > 5) fMVAVar_R9 = 5;
1731 
1732  // Energy matching
1733  fMVAVar_HoE = ele.hadronicOverEm();
1734  fMVAVar_EoP = ( ele.eSuperClusterOverP() > 20 ) ? 20 : ele.eSuperClusterOverP();
1735  fMVAVar_IoEmIoP = (1.0/ele.superCluster()->energy()) - (1.0 / ele.trackMomentumAtVtx().R()); //this is the proper variable
1736  fMVAVar_eleEoPout = ( ele.eEleClusterOverPout() > 20 ) ? 20 : ele.eEleClusterOverPout();
1737  fMVAVar_PreShowerOverRaw= ele.superCluster()->preshowerEnergy() / ele.superCluster()->rawEnergy();
1738 
1739  // Spectators
1740  fMVAVar_eta = ele.superCluster()->eta();
1741  fMVAVar_pt = ele.pt();
1742 
1743 
1744 
1745  // for triggering electrons get the impact parameteres
1746  if(fMVAType == kTrig) {
1747  //d0
1748  if (ele.gsfTrack().isNonnull()) {
1749  fMVAVar_d0 = (-1.0)*ele.gsfTrack()->dxy(vertex.position());
1750  } else if (ele.closestCtfTrackRef().isNonnull()) {
1751  fMVAVar_d0 = (-1.0)*ele.closestCtfTrackRef()->dxy(vertex.position());
1752  } else {
1753  fMVAVar_d0 = -9999.0;
1754  }
1755 
1756  //default values for IP3D
1757  fMVAVar_ip3d = -999.0;
1758  fMVAVar_ip3dSig = 0.0;
1759  if (ele.gsfTrack().isNonnull()) {
1760  const double gsfsign = ( (-ele.gsfTrack()->dxy(vertex.position())) >=0 ) ? 1. : -1.;
1761 
1762  const reco::TransientTrack &tt = transientTrackBuilder.build(ele.gsfTrack());
1763  const std::pair<bool,Measurement1D> &ip3dpv = IPTools::absoluteImpactParameter3D(tt,vertex);
1764  if (ip3dpv.first) {
1765  double ip3d = gsfsign*ip3dpv.second.value();
1766  double ip3derr = ip3dpv.second.error();
1767  fMVAVar_ip3d = ip3d;
1768  fMVAVar_ip3dSig = ip3d/ip3derr;
1769  }
1770  }
1771  }
1772 
1773  //**********************************************************
1774  //Isolation variables
1775  //**********************************************************
1776  Double_t tmpChargedIso_DR0p0To0p1 = 0;
1777  Double_t tmpChargedIso_DR0p1To0p2 = 0;
1778  Double_t tmpChargedIso_DR0p2To0p3 = 0;
1779  Double_t tmpChargedIso_DR0p3To0p4 = 0;
1780  Double_t tmpChargedIso_DR0p4To0p5 = 0;
1781  Double_t tmpGammaIso_DR0p0To0p1 = 0;
1782  Double_t tmpGammaIso_DR0p1To0p2 = 0;
1783  Double_t tmpGammaIso_DR0p2To0p3 = 0;
1784  Double_t tmpGammaIso_DR0p3To0p4 = 0;
1785  Double_t tmpGammaIso_DR0p4To0p5 = 0;
1786  Double_t tmpNeutralHadronIso_DR0p0To0p1 = 0;
1787  Double_t tmpNeutralHadronIso_DR0p1To0p2 = 0;
1788  Double_t tmpNeutralHadronIso_DR0p2To0p3 = 0;
1789  Double_t tmpNeutralHadronIso_DR0p3To0p4 = 0;
1790  Double_t tmpNeutralHadronIso_DR0p4To0p5 = 0;
1791 
1792  for (reco::PFCandidateCollection::const_iterator iP = PFCandidates.begin();
1793  iP != PFCandidates.end(); ++iP) {
1794 
1795  double dr = sqrt(pow(iP->eta() - ele.eta(),2) + pow(acos(cos(iP->phi() - ele.phi())),2));
1796 
1797  Bool_t passVeto = kTRUE;
1798  //Charged
1799  if(iP->trackRef().isNonnull()) {
1800 
1801  //make sure charged pf candidates pass the PFNoPU condition (assumed)
1802 
1803  //************************************************************
1804  // Veto any PFmuon, or PFEle
1805  if (iP->particleId() == reco::PFCandidate::e || iP->particleId() == reco::PFCandidate::mu) passVeto = kFALSE;
1806  //************************************************************
1807  //************************************************************
1808  // Footprint Veto
1809  if (fabs(fMVAVar_eta) > 1.479 && dr < 0.015) passVeto = kFALSE;
1810  //************************************************************
1811  if (passVeto) {
1812  if (dr < 0.1) tmpChargedIso_DR0p0To0p1 += iP->pt();
1813  if (dr >= 0.1 && dr < 0.2) tmpChargedIso_DR0p1To0p2 += iP->pt();
1814  if (dr >= 0.2 && dr < 0.3) tmpChargedIso_DR0p2To0p3 += iP->pt();
1815  if (dr >= 0.3 && dr < 0.4) tmpChargedIso_DR0p3To0p4 += iP->pt();
1816  if (dr >= 0.4 && dr < 0.5) tmpChargedIso_DR0p4To0p5 += iP->pt();
1817  } //pass veto
1818  }
1819  //Gamma
1820  else if (iP->particleId() == reco::PFCandidate::gamma) {
1821  //************************************************************
1822  // Footprint Veto
1823  if (fabs(fMVAVar_eta) > 1.479 && dr < 0.08) passVeto = kFALSE;
1824  //************************************************************
1825  if (passVeto) {
1826  if (dr < 0.1) tmpGammaIso_DR0p0To0p1 += iP->pt();
1827  if (dr >= 0.1 && dr < 0.2) tmpGammaIso_DR0p1To0p2 += iP->pt();
1828  if (dr >= 0.2 && dr < 0.3) tmpGammaIso_DR0p2To0p3 += iP->pt();
1829  if (dr >= 0.3 && dr < 0.4) tmpGammaIso_DR0p3To0p4 += iP->pt();
1830  if (dr >= 0.4 && dr < 0.5) tmpGammaIso_DR0p4To0p5 += iP->pt();
1831  }
1832  }
1833  //NeutralHadron
1834  else {
1835  if (dr < 0.1) tmpNeutralHadronIso_DR0p0To0p1 += iP->pt();
1836  if (dr >= 0.1 && dr < 0.2) tmpNeutralHadronIso_DR0p1To0p2 += iP->pt();
1837  if (dr >= 0.2 && dr < 0.3) tmpNeutralHadronIso_DR0p2To0p3 += iP->pt();
1838  if (dr >= 0.3 && dr < 0.4) tmpNeutralHadronIso_DR0p3To0p4 += iP->pt();
1839  if (dr >= 0.4 && dr < 0.5) tmpNeutralHadronIso_DR0p4To0p5 += iP->pt();
1840  }
1841  } //loop over PF candidates
1842 
1843  if (fMVAType == kTrigIDIsoCombinedPUCorrected) {
1844  fMVAVar_ChargedIso_DR0p0To0p1 = TMath::Min((tmpChargedIso_DR0p0To0p1)/ele.pt(), 2.5);
1845  fMVAVar_ChargedIso_DR0p1To0p2 = TMath::Min((tmpChargedIso_DR0p1To0p2)/ele.pt(), 2.5);
1846  fMVAVar_ChargedIso_DR0p2To0p3 = TMath::Min((tmpChargedIso_DR0p2To0p3)/ele.pt(), 2.5);
1847  fMVAVar_ChargedIso_DR0p3To0p4 = TMath::Min((tmpChargedIso_DR0p3To0p4)/ele.pt(), 2.5);
1848  fMVAVar_ChargedIso_DR0p4To0p5 = TMath::Min((tmpChargedIso_DR0p4To0p5)/ele.pt(), 2.5);
1849  fMVAVar_GammaIso_DR0p0To0p1 = TMath::Max(TMath::Min((tmpGammaIso_DR0p0To0p1 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p0To0p1, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1850  fMVAVar_GammaIso_DR0p1To0p2 = TMath::Max(TMath::Min((tmpGammaIso_DR0p1To0p2 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p1To0p2, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1851  fMVAVar_GammaIso_DR0p2To0p3 = TMath::Max(TMath::Min((tmpGammaIso_DR0p2To0p3 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p2To0p3, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1852  fMVAVar_GammaIso_DR0p3To0p4 = TMath::Max(TMath::Min((tmpGammaIso_DR0p3To0p4 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p3To0p4, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1853  fMVAVar_GammaIso_DR0p4To0p5 = TMath::Max(TMath::Min((tmpGammaIso_DR0p4To0p5 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p4To0p5, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1854  fMVAVar_NeutralHadronIso_DR0p0To0p1 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p0To0p1 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p0To0p1, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1855  fMVAVar_NeutralHadronIso_DR0p1To0p2 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p1To0p2 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p1To0p2, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1856  fMVAVar_NeutralHadronIso_DR0p2To0p3 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p2To0p3 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p2To0p3, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1857  fMVAVar_NeutralHadronIso_DR0p3To0p4 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p3To0p4 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p3To0p4, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1858  fMVAVar_NeutralHadronIso_DR0p4To0p5 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p4To0p5 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p4To0p5, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1859  } else if (fMVAType == kTrigIDIsoCombined) {
1860  fMVAVar_ChargedIso_DR0p0To0p1 = TMath::Min((tmpChargedIso_DR0p0To0p1)/ele.pt(), 2.5);
1861  fMVAVar_ChargedIso_DR0p1To0p2 = TMath::Min((tmpChargedIso_DR0p1To0p2)/ele.pt(), 2.5);
1862  fMVAVar_ChargedIso_DR0p2To0p3 = TMath::Min((tmpChargedIso_DR0p2To0p3)/ele.pt(), 2.5);
1863  fMVAVar_ChargedIso_DR0p3To0p4 = TMath::Min((tmpChargedIso_DR0p3To0p4)/ele.pt(), 2.5);
1864  fMVAVar_ChargedIso_DR0p4To0p5 = TMath::Min((tmpChargedIso_DR0p4To0p5)/ele.pt(), 2.5);
1865  fMVAVar_GammaIso_DR0p0To0p1 = TMath::Max(TMath::Min((tmpGammaIso_DR0p0To0p1)/ele.pt(), 2.5), 0.0);
1866  fMVAVar_GammaIso_DR0p1To0p2 = TMath::Max(TMath::Min((tmpGammaIso_DR0p1To0p2)/ele.pt(), 2.5), 0.0);
1867  fMVAVar_GammaIso_DR0p2To0p3 = TMath::Max(TMath::Min((tmpGammaIso_DR0p2To0p3)/ele.pt(), 2.5), 0.0);
1868  fMVAVar_GammaIso_DR0p3To0p4 = TMath::Max(TMath::Min((tmpGammaIso_DR0p3To0p4)/ele.pt(), 2.5), 0.0);
1869  fMVAVar_GammaIso_DR0p4To0p5 = TMath::Max(TMath::Min((tmpGammaIso_DR0p4To0p5)/ele.pt(), 2.5), 0.0);
1870  fMVAVar_NeutralHadronIso_DR0p0To0p1 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p0To0p1)/ele.pt(), 2.5), 0.0);
1871  fMVAVar_NeutralHadronIso_DR0p1To0p2 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p1To0p2)/ele.pt(), 2.5), 0.0);
1872  fMVAVar_NeutralHadronIso_DR0p2To0p3 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p2To0p3)/ele.pt(), 2.5), 0.0);
1873  fMVAVar_NeutralHadronIso_DR0p3To0p4 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p3To0p4)/ele.pt(), 2.5), 0.0);
1874  fMVAVar_NeutralHadronIso_DR0p4To0p5 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p4To0p5)/ele.pt(), 2.5), 0.0);
1875  fMVAVar_rho = Rho;
1876  } else {
1877  cout << "Warning: Type " << fMVAType << " is not supported.\n";
1878  }
1879 
1880  // evaluate
1881  Double_t mva = -9999;
1882  if (fUseBinnedVersion) {
1883  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
1884  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
1885  } else {
1886  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
1887  }
1888 
1889 
1890 
1891  if(printDebug) {
1892  cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType << endl;
1893  cout << " fbrem " << fMVAVar_fbrem
1894  << " kfchi2 " << fMVAVar_kfchi2
1895  << " mykfhits " << fMVAVar_kfhits
1896  << " gsfchi2 " << fMVAVar_gsfchi2
1897  << " deta " << fMVAVar_deta
1898  << " dphi " << fMVAVar_dphi
1899  << " detacalo " << fMVAVar_detacalo
1900  << " see " << fMVAVar_see
1901  << " spp " << fMVAVar_spp
1902  << " etawidth " << fMVAVar_etawidth
1903  << " phiwidth " << fMVAVar_phiwidth
1904  << " OneMinusE1x5E5x5 " << fMVAVar_OneMinusE1x5E5x5
1905  << " R9 " << fMVAVar_R9
1906  << " HoE " << fMVAVar_HoE
1907  << " EoP " << fMVAVar_EoP
1908  << " IoEmIoP " << fMVAVar_IoEmIoP
1909  << " eleEoPout " << fMVAVar_eleEoPout
1910  << " d0 " << fMVAVar_d0
1911  << " ip3d " << fMVAVar_ip3d
1912  << " eta " << fMVAVar_eta
1913  << " pt " << fMVAVar_pt << endl;
1914  cout << "ChargedIso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
1915  << fMVAVar_ChargedIso_DR0p0To0p1 << " "
1916  << fMVAVar_ChargedIso_DR0p1To0p2 << " "
1917  << fMVAVar_ChargedIso_DR0p2To0p3 << " "
1918  << fMVAVar_ChargedIso_DR0p3To0p4 << " "
1919  << fMVAVar_ChargedIso_DR0p4To0p5 << endl;
1920  cout << "PF Gamma Iso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
1921  << fMVAVar_GammaIso_DR0p0To0p1 << " "
1922  << fMVAVar_GammaIso_DR0p1To0p2 << " "
1923  << fMVAVar_GammaIso_DR0p2To0p3 << " "
1924  << fMVAVar_GammaIso_DR0p3To0p4 << " "
1925  << fMVAVar_GammaIso_DR0p4To0p5 << endl;
1926  cout << "PF Neutral Hadron Iso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
1927  << fMVAVar_NeutralHadronIso_DR0p0To0p1 << " "
1928  << fMVAVar_NeutralHadronIso_DR0p1To0p2 << " "
1929  << fMVAVar_NeutralHadronIso_DR0p2To0p3 << " "
1930  << fMVAVar_NeutralHadronIso_DR0p3To0p4 << " "
1931  << fMVAVar_NeutralHadronIso_DR0p4To0p5 << " "
1932  << endl;
1933  cout << "Rho : " << Rho << endl;
1934  cout << " ### MVA " << mva << endl;
1935  }
1936 
1937 
1938 
1939  return mva;
1940 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:576
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:243
reco::GsfElectron::trackMomentumAtVtx
math::XYZVectorF trackMomentumAtVtx() const
Definition: GsfElectron.h:289
ElectronEffectiveArea::GetElectronEffectiveArea
static Double_t GetElectronEffectiveArea(ElectronEffectiveAreaType type, Double_t SCEta, ElectronEffectiveAreaTarget EffectiveAreaTarget=kEleEAData2011)
Definition: ElectronEffectiveArea.h:64
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:10
IPTools::absoluteImpactParameter3D
std::pair< bool, Measurement1D > absoluteImpactParameter3D(const reco::TransientTrack &transientTrack, const reco::Vertex &vertex)
Definition: IPTools.cc:37
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:698
Min
T Min(T a, T b)
Definition: MathUtil.h:39
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:106
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:137
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
reco::GsfElectron::closestCtfTrackRef
TrackRef closestCtfTrackRef() const
Definition: GsfElectron.h:199
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:247
reco::GsfElectron::sigmaIetaIeta
float sigmaIetaIeta() const
Definition: GsfElectron.h:416
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:457
edm::detail::isnan
bool isnan(float x)
Definition: math.h:13
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::GsfElectron::superCluster
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
Definition: GsfElectron.h:182
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:250
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:246
min
T min(T a, T b)
Definition: MathUtil.h:58
Max
T Max(T a, T b)
Definition: MathUtil.h:44
reco::GsfElectron::e1x5
float e1x5() const
Definition: GsfElectron.h:418
reco::GsfElectron::e5x5
float e5x5() const
Definition: GsfElectron.h:420
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:248
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:133
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:183
benchmark_cfg.PFCandidates
tuple PFCandidates
Definition: benchmark_cfg.py:32
Double_t EGammaMvaEleEstimator::IDIsoCombinedMvaValue ( Double_t  fbrem,
Double_t  kfchi2,
Int_t  kfhits,
Double_t  gsfchi2,
Double_t  deta,
Double_t  dphi,
Double_t  detacalo,
Double_t  see,
Double_t  spp,
Double_t  etawidth,
Double_t  phiwidth,
Double_t  OneMinusE1x5E5x5,
Double_t  R9,
Double_t  HoE,
Double_t  EoP,
Double_t  IoEmIoP,
Double_t  eleEoPout,
Double_t  PreShowerOverRaw,
Double_t  d0,
Double_t  ip3d,
Double_t  ChargedIso_DR0p0To0p1,
Double_t  ChargedIso_DR0p1To0p2,
Double_t  ChargedIso_DR0p2To0p3,
Double_t  ChargedIso_DR0p3To0p4,
Double_t  ChargedIso_DR0p4To0p5,
Double_t  GammaIso_DR0p0To0p1,
Double_t  GammaIso_DR0p1To0p2,
Double_t  GammaIso_DR0p2To0p3,
Double_t  GammaIso_DR0p3To0p4,
Double_t  GammaIso_DR0p4To0p5,
Double_t  NeutralHadronIso_DR0p0To0p1,
Double_t  NeutralHadronIso_DR0p1To0p2,
Double_t  NeutralHadronIso_DR0p2To0p3,
Double_t  NeutralHadronIso_DR0p3To0p4,
Double_t  NeutralHadronIso_DR0p4To0p5,
Double_t  Rho,
Double_t  eta,
Double_t  pt,
Bool_t  printDebug = kFALSE 
)

Definition at line 773 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, gather_cfg::cout, eta, fisInitialized, fMethodname, fMVAVar_ChargedIso_DR0p0To0p1, fMVAVar_ChargedIso_DR0p1To0p2, fMVAVar_ChargedIso_DR0p2To0p3, fMVAVar_ChargedIso_DR0p3To0p4, fMVAVar_ChargedIso_DR0p4To0p5, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_GammaIso_DR0p0To0p1, fMVAVar_GammaIso_DR0p1To0p2, fMVAVar_GammaIso_DR0p2To0p3, fMVAVar_GammaIso_DR0p3To0p4, fMVAVar_GammaIso_DR0p4To0p5, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_ip3d, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_NeutralHadronIso_DR0p0To0p1, fMVAVar_NeutralHadronIso_DR0p1To0p2, fMVAVar_NeutralHadronIso_DR0p2To0p3, fMVAVar_NeutralHadronIso_DR0p3To0p4, fMVAVar_NeutralHadronIso_DR0p4To0p5, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_rho, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, GetMVABin(), bookConverter::max, min(), and EnergyCorrector::pt.

811  {
812 
813  if (!fisInitialized) {
814  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
815  return -9999;
816  }
817 
818  fMVAVar_fbrem = ( fbrem < -1.0 ) ? -1.0 : fbrem;
819  fMVAVar_kfchi2 = ( kfchi2 > 10 ) ? 10 : kfchi2;
820  fMVAVar_kfhits = float(kfhits); // BTD does not support int variables
821  fMVAVar_gsfchi2 = ( gsfchi2 > 200 ) ? 200 : gsfchi2;
822  fMVAVar_deta = ( fabs(deta) > 0.06 ) ? 0.06 : fabs(deta);
823  fMVAVar_dphi = dphi;
824  fMVAVar_detacalo = detacalo;
825 
826  fMVAVar_see = see;
827  fMVAVar_spp = spp;
828  fMVAVar_etawidth = etawidth;
829  fMVAVar_phiwidth = phiwidth;
830  fMVAVar_OneMinusE1x5E5x5= max(min(double(OneMinusE1x5E5x5),2.0),-1.0);
831  fMVAVar_R9 = (R9 > 5) ? 5: R9;
832 
833  fMVAVar_HoE = HoE;
834  fMVAVar_EoP = (EoP > 20) ? 20 : EoP;
835  fMVAVar_IoEmIoP = IoEmIoP;
836  fMVAVar_eleEoPout = (eleEoPout > 20) ? 20 : eleEoPout;
837  fMVAVar_PreShowerOverRaw= PreShowerOverRaw;
838 
839  fMVAVar_d0 = d0;
840  fMVAVar_ip3d = ip3d;
841 
842  fMVAVar_ChargedIso_DR0p0To0p1 = ChargedIso_DR0p0To0p1;
843  fMVAVar_ChargedIso_DR0p1To0p2 = ChargedIso_DR0p1To0p2;
844  fMVAVar_ChargedIso_DR0p2To0p3 = ChargedIso_DR0p2To0p3;
845  fMVAVar_ChargedIso_DR0p3To0p4 = ChargedIso_DR0p3To0p4;
846  fMVAVar_ChargedIso_DR0p4To0p5 = ChargedIso_DR0p4To0p5;
847  fMVAVar_GammaIso_DR0p0To0p1 = GammaIso_DR0p0To0p1;
848  fMVAVar_GammaIso_DR0p1To0p2 = GammaIso_DR0p1To0p2;
849  fMVAVar_GammaIso_DR0p2To0p3 = GammaIso_DR0p2To0p3;
850  fMVAVar_GammaIso_DR0p3To0p4 = GammaIso_DR0p3To0p4;
851  fMVAVar_GammaIso_DR0p4To0p5 = GammaIso_DR0p4To0p5;
852  fMVAVar_NeutralHadronIso_DR0p0To0p1 = NeutralHadronIso_DR0p0To0p1;
853  fMVAVar_NeutralHadronIso_DR0p1To0p2 = NeutralHadronIso_DR0p1To0p2;
854  fMVAVar_NeutralHadronIso_DR0p2To0p3 = NeutralHadronIso_DR0p2To0p3;
855  fMVAVar_NeutralHadronIso_DR0p3To0p4 = NeutralHadronIso_DR0p3To0p4;
856  fMVAVar_NeutralHadronIso_DR0p4To0p5 = NeutralHadronIso_DR0p4To0p5;
857 
858  fMVAVar_rho = Rho;
859  fMVAVar_eta = eta;
860  fMVAVar_pt = pt;
861 
862  Double_t mva = -9999;
863  if (fUseBinnedVersion) {
864  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
865  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
866  } else {
867  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
868  }
869 
870  if(printDebug) {
871  cout << " *** Inside the class fMethodname " << fMethodname << endl;
872  cout << " fbrem " << fMVAVar_fbrem
873  << " kfchi2 " << fMVAVar_kfchi2
874  << " mykfhits " << fMVAVar_kfhits
875  << " gsfchi2 " << fMVAVar_gsfchi2
876  << " deta " << fMVAVar_deta
877  << " dphi " << fMVAVar_dphi
878  << " detacalo " << fMVAVar_detacalo
879  << " see " << fMVAVar_see
880  << " spp " << fMVAVar_spp
881  << " etawidth " << fMVAVar_etawidth
882  << " phiwidth " << fMVAVar_phiwidth
883  << " OneMinusE1x5E5x5 " << fMVAVar_OneMinusE1x5E5x5
884  << " R9 " << fMVAVar_R9
885  << " HoE " << fMVAVar_HoE
886  << " EoP " << fMVAVar_EoP
887  << " IoEmIoP " << fMVAVar_IoEmIoP
888  << " eleEoPout " << fMVAVar_eleEoPout
889  << " PreShowerOverRaw " << fMVAVar_PreShowerOverRaw
890  << " d0 " << fMVAVar_d0
891  << " ip3d " << fMVAVar_ip3d
892  << " ChargedIso_DR0p0To0p1 " << ChargedIso_DR0p0To0p1
893  << " ChargedIso_DR0p1To0p2 " << ChargedIso_DR0p1To0p2
894  << " ChargedIso_DR0p2To0p3 " << ChargedIso_DR0p2To0p3
895  << " ChargedIso_DR0p3To0p4 " << ChargedIso_DR0p3To0p4
896  << " ChargedIso_DR0p4To0p5 " << ChargedIso_DR0p4To0p5
897  << " GammaIso_DR0p0To0p1 " << GammaIso_DR0p0To0p1
898  << " GammaIso_DR0p1To0p2 " << GammaIso_DR0p1To0p2
899  << " GammaIso_DR0p2To0p3 " << GammaIso_DR0p2To0p3
900  << " GammaIso_DR0p3To0p4 " << GammaIso_DR0p3To0p4
901  << " GammaIso_DR0p4To0p5 " << GammaIso_DR0p4To0p5
902  << " NeutralHadronIso_DR0p0To0p1 " << NeutralHadronIso_DR0p0To0p1
903  << " NeutralHadronIso_DR0p1To0p2 " << NeutralHadronIso_DR0p1To0p2
904  << " NeutralHadronIso_DR0p2To0p3 " << NeutralHadronIso_DR0p2To0p3
905  << " NeutralHadronIso_DR0p3To0p4 " << NeutralHadronIso_DR0p3To0p4
906  << " NeutralHadronIso_DR0p4To0p5 " << NeutralHadronIso_DR0p4To0p5
907  << " Rho " << Rho
908  << " eta " << fMVAVar_eta
909  << " pt " << fMVAVar_pt << endl;
910  cout << " ### MVA " << mva << endl;
911  }
912 
913  return mva;
914 }
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
min
T min(T a, T b)
Definition: MathUtil.h:58
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
void EGammaMvaEleEstimator::initialize ( std::string  methodName,
std::string  weightsfile,
EGammaMvaEleEstimator::MVAType  type 
)

Definition at line 49 of file EGammaMvaEleEstimator.cc.

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

52 {
53 
54  std::vector<std::string> tempWeightFileVector;
55  tempWeightFileVector.push_back(weightsfile);
56  initialize(methodName,type,kFALSE,tempWeightFileVector);
57 }
type
type
Definition: HCALResponse.h:21
EGammaMvaEleEstimator::initialize
void initialize(std::string methodName, std::string weightsfile, EGammaMvaEleEstimator::MVAType type)
Definition: EGammaMvaEleEstimator.cc:49
void EGammaMvaEleEstimator::initialize ( std::string  methodName,
EGammaMvaEleEstimator::MVAType  type,
Bool_t  useBinnedVersion,
std::vector< std::string >  weightsfiles 
)

Definition at line 61 of file EGammaMvaEleEstimator.cc.

References assert(), gather_cfg::cout, data, Exception, mergeVDriftHistosByStation::file, fisInitialized, fMethodname, fMVAType, fMVAVar_ChargedIso_DR0p0To0p1, fMVAVar_ChargedIso_DR0p1To0p2, fMVAVar_ChargedIso_DR0p2To0p3, fMVAVar_ChargedIso_DR0p3To0p4, fMVAVar_ChargedIso_DR0p4To0p5, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_GammaIso_DR0p0To0p1, fMVAVar_GammaIso_DR0p1To0p2, fMVAVar_GammaIso_DR0p2To0p3, fMVAVar_GammaIso_DR0p3To0p4, fMVAVar_GammaIso_DR0p4To0p5, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_ip3d, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_NeutralHadronIso_DR0p0To0p1, fMVAVar_NeutralHadronIso_DR0p1To0p2, fMVAVar_NeutralHadronIso_DR0p2To0p3, fMVAVar_NeutralHadronIso_DR0p3To0p4, fMVAVar_NeutralHadronIso_DR0p4To0p5, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_rho, fMVAVar_see, fMVAVar_spp, fNMVABins, fTMVAMethod, fTMVAReader, fUseBinnedVersion, i, kIsoRings, kNonTrig, kTrig, kTrigIDIsoCombined, kTrigIDIsoCombinedPUCorrected, kTrigNoIP, and SiPixelLorentzAngle_cfi::read.

65  {
66 
67  //clean up first
68  for (unsigned int i=0;i<fTMVAReader.size(); ++i) {
69  if (fTMVAReader[i]) delete fTMVAReader[i];
70  if (fTMVAMethod[i]) delete fTMVAMethod[i];
71  }
72  fTMVAReader.clear();
73  fTMVAMethod.clear();
74 
75  //initialize
76  fisInitialized = kTRUE;
77  fMVAType = type;
78  fMethodname = methodName;
79  fUseBinnedVersion = useBinnedVersion;
80 
81  //Define expected number of bins
82  UInt_t ExpectedNBins = 0;
83  if (!fUseBinnedVersion) {
84  ExpectedNBins = 1;
85  } else if (type == kTrig) {
86  ExpectedNBins = 6;
87  }else if (type == kTrigNoIP) {
88  ExpectedNBins = 6;
89  } else if (type == kNonTrig) {
90  ExpectedNBins = 6;
91  } else if (type == kIsoRings) {
92  ExpectedNBins = 4;
93  } else if (type == kTrigIDIsoCombined) {
94  ExpectedNBins = 6;
95  } else if (type == kTrigIDIsoCombinedPUCorrected) {
96  ExpectedNBins = 6;
97  }
98 
99  fNMVABins = ExpectedNBins;
100 
101  //Check number of weight files given
102  if (fNMVABins != weightsfiles.size() ) {
103  std::cout << "Error: Expected Number of bins = " << fNMVABins << " does not equal to weightsfiles.size() = "
104  << weightsfiles.size() << std::endl;
105 
106  #ifndef STANDALONE
107  assert(fNMVABins == weightsfiles.size());
108  #endif
109  }
110 
111  //Loop over all bins
112  for (unsigned int i=0;i<fNMVABins; ++i) {
113 
114  TMVA::Reader *tmpTMVAReader = new TMVA::Reader( "!Color:!Silent:Error" );
115  tmpTMVAReader->SetVerbose(kTRUE);
116 
117  if (type == kTrig) {
118  // Pure tracking variables
119  tmpTMVAReader->AddVariable("fbrem", &fMVAVar_fbrem);
120  tmpTMVAReader->AddVariable("kfchi2", &fMVAVar_kfchi2);
121  tmpTMVAReader->AddVariable("kfhits", &fMVAVar_kfhits);
122  tmpTMVAReader->AddVariable("gsfchi2", &fMVAVar_gsfchi2);
123 
124  // Geometrical matchings
125  tmpTMVAReader->AddVariable("deta", &fMVAVar_deta);
126  tmpTMVAReader->AddVariable("dphi", &fMVAVar_dphi);
127  tmpTMVAReader->AddVariable("detacalo", &fMVAVar_detacalo);
128 
129  // Pure ECAL -> shower shapes
130  tmpTMVAReader->AddVariable("see", &fMVAVar_see);
131  tmpTMVAReader->AddVariable("spp", &fMVAVar_spp);
132  tmpTMVAReader->AddVariable("etawidth", &fMVAVar_etawidth);
133  tmpTMVAReader->AddVariable("phiwidth", &fMVAVar_phiwidth);
134  tmpTMVAReader->AddVariable("e1x5e5x5", &fMVAVar_OneMinusE1x5E5x5);
135  tmpTMVAReader->AddVariable("R9", &fMVAVar_R9);
136 
137  // Energy matching
138  tmpTMVAReader->AddVariable("HoE", &fMVAVar_HoE);
139  tmpTMVAReader->AddVariable("EoP", &fMVAVar_EoP);
140  tmpTMVAReader->AddVariable("IoEmIoP", &fMVAVar_IoEmIoP);
141  tmpTMVAReader->AddVariable("eleEoPout", &fMVAVar_eleEoPout);
142  if(i == 2 || i == 5)
143  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
144 
145  if(!fUseBinnedVersion)
146  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
147 
148  // IP
149  tmpTMVAReader->AddVariable("d0", &fMVAVar_d0);
150  tmpTMVAReader->AddVariable("ip3d", &fMVAVar_ip3d);
151 
152  tmpTMVAReader->AddSpectator("eta", &fMVAVar_eta);
153  tmpTMVAReader->AddSpectator("pt", &fMVAVar_pt);
154  }
155 
156  if (type == kTrigNoIP) {
157 
158 
159  // Pure tracking variables
160  tmpTMVAReader->AddVariable("fbrem", &fMVAVar_fbrem);
161  tmpTMVAReader->AddVariable("kfchi2", &fMVAVar_kfchi2);
162  tmpTMVAReader->AddVariable("kfhits", &fMVAVar_kfhits);
163  tmpTMVAReader->AddVariable("gsfchi2", &fMVAVar_gsfchi2);
164 
165  // Geometrical matchings
166  tmpTMVAReader->AddVariable("deta", &fMVAVar_deta);
167  tmpTMVAReader->AddVariable("dphi", &fMVAVar_dphi);
168  tmpTMVAReader->AddVariable("detacalo", &fMVAVar_detacalo);
169 
170  // Pure ECAL -> shower shapes
171  tmpTMVAReader->AddVariable("see", &fMVAVar_see);
172  tmpTMVAReader->AddVariable("spp", &fMVAVar_spp);
173  tmpTMVAReader->AddVariable("etawidth", &fMVAVar_etawidth);
174  tmpTMVAReader->AddVariable("phiwidth", &fMVAVar_phiwidth);
175  tmpTMVAReader->AddVariable("e1x5e5x5", &fMVAVar_OneMinusE1x5E5x5);
176  tmpTMVAReader->AddVariable("R9", &fMVAVar_R9);
177 
178  // Energy matching
179  tmpTMVAReader->AddVariable("HoE", &fMVAVar_HoE);
180  tmpTMVAReader->AddVariable("EoP", &fMVAVar_EoP);
181  tmpTMVAReader->AddVariable("IoEmIoP", &fMVAVar_IoEmIoP);
182  tmpTMVAReader->AddVariable("eleEoPout", &fMVAVar_eleEoPout);
183  tmpTMVAReader->AddVariable("rho", &fMVAVar_rho);
184 
185  if(i == 2 || i == 5)
186  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
187 
188  if(!fUseBinnedVersion)
189  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
190 
191  tmpTMVAReader->AddSpectator("eta", &fMVAVar_eta);
192  tmpTMVAReader->AddSpectator("pt", &fMVAVar_pt);
193 
194  }
195 
196  if (type == kNonTrig) {
197 
198  // Pure tracking variables
199  tmpTMVAReader->AddVariable("fbrem", &fMVAVar_fbrem);
200  tmpTMVAReader->AddVariable("kfchi2", &fMVAVar_kfchi2);
201  tmpTMVAReader->AddVariable("kfhits", &fMVAVar_kfhits);
202  tmpTMVAReader->AddVariable("gsfchi2", &fMVAVar_gsfchi2);
203 
204  // Geometrical matchings
205  tmpTMVAReader->AddVariable("deta", &fMVAVar_deta);
206  tmpTMVAReader->AddVariable("dphi", &fMVAVar_dphi);
207  tmpTMVAReader->AddVariable("detacalo", &fMVAVar_detacalo);
208 
209  // Pure ECAL -> shower shapes
210  tmpTMVAReader->AddVariable("see", &fMVAVar_see);
211  tmpTMVAReader->AddVariable("spp", &fMVAVar_spp);
212  tmpTMVAReader->AddVariable("etawidth", &fMVAVar_etawidth);
213  tmpTMVAReader->AddVariable("phiwidth", &fMVAVar_phiwidth);
214  tmpTMVAReader->AddVariable("e1x5e5x5", &fMVAVar_OneMinusE1x5E5x5);
215  tmpTMVAReader->AddVariable("R9", &fMVAVar_R9);
216 
217  // Energy matching
218  tmpTMVAReader->AddVariable("HoE", &fMVAVar_HoE);
219  tmpTMVAReader->AddVariable("EoP", &fMVAVar_EoP);
220  tmpTMVAReader->AddVariable("IoEmIoP", &fMVAVar_IoEmIoP);
221  tmpTMVAReader->AddVariable("eleEoPout", &fMVAVar_eleEoPout);
222  if(i == 2 || i == 5)
223  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
224 
225  if(!fUseBinnedVersion)
226  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
227 
228  tmpTMVAReader->AddSpectator("eta", &fMVAVar_eta);
229  tmpTMVAReader->AddSpectator("pt", &fMVAVar_pt);
230 
231  }
232 
233  if (type == kIsoRings) {
234  tmpTMVAReader->AddVariable( "ChargedIso_DR0p0To0p1", &fMVAVar_ChargedIso_DR0p0To0p1 );
235  tmpTMVAReader->AddVariable( "ChargedIso_DR0p1To0p2", &fMVAVar_ChargedIso_DR0p1To0p2 );
236  tmpTMVAReader->AddVariable( "ChargedIso_DR0p2To0p3", &fMVAVar_ChargedIso_DR0p2To0p3 );
237  tmpTMVAReader->AddVariable( "ChargedIso_DR0p3To0p4", &fMVAVar_ChargedIso_DR0p3To0p4 );
238  tmpTMVAReader->AddVariable( "ChargedIso_DR0p4To0p5", &fMVAVar_ChargedIso_DR0p4To0p5 );
239  tmpTMVAReader->AddVariable( "GammaIso_DR0p0To0p1", &fMVAVar_GammaIso_DR0p0To0p1 );
240  tmpTMVAReader->AddVariable( "GammaIso_DR0p1To0p2", &fMVAVar_GammaIso_DR0p1To0p2 );
241  tmpTMVAReader->AddVariable( "GammaIso_DR0p2To0p3", &fMVAVar_GammaIso_DR0p2To0p3 );
242  tmpTMVAReader->AddVariable( "GammaIso_DR0p3To0p4", &fMVAVar_GammaIso_DR0p3To0p4 );
243  tmpTMVAReader->AddVariable( "GammaIso_DR0p4To0p5", &fMVAVar_GammaIso_DR0p4To0p5 );
244  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p0To0p1", &fMVAVar_NeutralHadronIso_DR0p0To0p1 );
245  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p1To0p2", &fMVAVar_NeutralHadronIso_DR0p1To0p2 );
246  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p2To0p3", &fMVAVar_NeutralHadronIso_DR0p2To0p3 );
247  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p3To0p4", &fMVAVar_NeutralHadronIso_DR0p3To0p4 );
248  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p4To0p5", &fMVAVar_NeutralHadronIso_DR0p4To0p5 );
249  tmpTMVAReader->AddSpectator("eta", &fMVAVar_eta);
250  tmpTMVAReader->AddSpectator("pt", &fMVAVar_pt);
251  }
252 
253 
254  if (type == kTrigIDIsoCombinedPUCorrected) {
255 
256  // Pure tracking variables
257  tmpTMVAReader->AddVariable("fbrem", &fMVAVar_fbrem);
258  tmpTMVAReader->AddVariable("kfchi2", &fMVAVar_kfchi2);
259  tmpTMVAReader->AddVariable("kflayers", &fMVAVar_kfhits);
260  tmpTMVAReader->AddVariable("gsfchi2", &fMVAVar_gsfchi2);
261 
262  // Geometrical matchings
263  tmpTMVAReader->AddVariable("deta", &fMVAVar_deta);
264  tmpTMVAReader->AddVariable("dphi", &fMVAVar_dphi);
265  tmpTMVAReader->AddVariable("detacalo", &fMVAVar_detacalo);
266 
267  // Pure ECAL -> shower shapes
268  tmpTMVAReader->AddVariable("see", &fMVAVar_see);
269  tmpTMVAReader->AddVariable("spp", &fMVAVar_spp);
270  tmpTMVAReader->AddVariable("etawidth", &fMVAVar_etawidth);
271  tmpTMVAReader->AddVariable("phiwidth", &fMVAVar_phiwidth);
272  tmpTMVAReader->AddVariable("OneMinusSeedE1x5OverE5x5", &fMVAVar_OneMinusE1x5E5x5);
273  tmpTMVAReader->AddVariable("R9", &fMVAVar_R9);
274 
275  // Energy matching
276  tmpTMVAReader->AddVariable("HoE", &fMVAVar_HoE);
277  tmpTMVAReader->AddVariable("EoP", &fMVAVar_EoP);
278  tmpTMVAReader->AddVariable("IoEmIoP", &fMVAVar_IoEmIoP);
279  tmpTMVAReader->AddVariable("EEleoPout", &fMVAVar_eleEoPout);
280  if(i == 2 || i == 5) {
281  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
282  }
283  if(!fUseBinnedVersion) {
284  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
285  }
286 
287  // IP
288  tmpTMVAReader->AddVariable("d0", &fMVAVar_d0);
289  tmpTMVAReader->AddVariable("ip3d", &fMVAVar_ip3d);
290 
291  //isolation variables
292  tmpTMVAReader->AddVariable( "ChargedIso_DR0p0To0p1", &fMVAVar_ChargedIso_DR0p0To0p1 );
293  tmpTMVAReader->AddVariable( "ChargedIso_DR0p1To0p2", &fMVAVar_ChargedIso_DR0p1To0p2 );
294  tmpTMVAReader->AddVariable( "ChargedIso_DR0p2To0p3", &fMVAVar_ChargedIso_DR0p2To0p3 );
295  tmpTMVAReader->AddVariable( "ChargedIso_DR0p3To0p4", &fMVAVar_ChargedIso_DR0p3To0p4 );
296  tmpTMVAReader->AddVariable( "ChargedIso_DR0p4To0p5", &fMVAVar_ChargedIso_DR0p4To0p5 );
297  tmpTMVAReader->AddVariable( "GammaIso_DR0p0To0p1", &fMVAVar_GammaIso_DR0p0To0p1 );
298  tmpTMVAReader->AddVariable( "GammaIso_DR0p1To0p2", &fMVAVar_GammaIso_DR0p1To0p2 );
299  tmpTMVAReader->AddVariable( "GammaIso_DR0p2To0p3", &fMVAVar_GammaIso_DR0p2To0p3 );
300  tmpTMVAReader->AddVariable( "GammaIso_DR0p3To0p4", &fMVAVar_GammaIso_DR0p3To0p4 );
301  tmpTMVAReader->AddVariable( "GammaIso_DR0p4To0p5", &fMVAVar_GammaIso_DR0p4To0p5 );
302  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p0To0p1", &fMVAVar_NeutralHadronIso_DR0p0To0p1 );
303  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p1To0p2", &fMVAVar_NeutralHadronIso_DR0p1To0p2 );
304  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p2To0p3", &fMVAVar_NeutralHadronIso_DR0p2To0p3 );
305  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p3To0p4", &fMVAVar_NeutralHadronIso_DR0p3To0p4 );
306  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p4To0p5", &fMVAVar_NeutralHadronIso_DR0p4To0p5 );
307 
308  //spectators
309  tmpTMVAReader->AddSpectator("eta", &fMVAVar_eta);
310  tmpTMVAReader->AddSpectator("pt", &fMVAVar_pt);
311 
312  }
313 
314  if (type == kTrigIDIsoCombined) {
315 
316  // Pure tracking variables
317  tmpTMVAReader->AddVariable("fbrem", &fMVAVar_fbrem);
318  tmpTMVAReader->AddVariable("kfchi2", &fMVAVar_kfchi2);
319  tmpTMVAReader->AddVariable("kflayers", &fMVAVar_kfhits);
320  tmpTMVAReader->AddVariable("gsfchi2", &fMVAVar_gsfchi2);
321 
322  // Geometrical matchings
323  tmpTMVAReader->AddVariable("deta", &fMVAVar_deta);
324  tmpTMVAReader->AddVariable("dphi", &fMVAVar_dphi);
325  tmpTMVAReader->AddVariable("detacalo", &fMVAVar_detacalo);
326 
327  // Pure ECAL -> shower shapes
328  tmpTMVAReader->AddVariable("see", &fMVAVar_see);
329  tmpTMVAReader->AddVariable("spp", &fMVAVar_spp);
330  tmpTMVAReader->AddVariable("etawidth", &fMVAVar_etawidth);
331  tmpTMVAReader->AddVariable("phiwidth", &fMVAVar_phiwidth);
332  tmpTMVAReader->AddVariable("OneMinusSeedE1x5OverE5x5", &fMVAVar_OneMinusE1x5E5x5);
333  tmpTMVAReader->AddVariable("R9", &fMVAVar_R9);
334 
335  // Energy matching
336  tmpTMVAReader->AddVariable("HoE", &fMVAVar_HoE);
337  tmpTMVAReader->AddVariable("EoP", &fMVAVar_EoP);
338  tmpTMVAReader->AddVariable("IoEmIoP", &fMVAVar_IoEmIoP);
339  tmpTMVAReader->AddVariable("EEleoPout", &fMVAVar_eleEoPout);
340  if(i == 2 || i == 5) {
341  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
342  }
343  if(!fUseBinnedVersion) {
344  tmpTMVAReader->AddVariable("PreShowerOverRaw",&fMVAVar_PreShowerOverRaw);
345  }
346 
347  // IP
348  tmpTMVAReader->AddVariable("d0", &fMVAVar_d0);
349  tmpTMVAReader->AddVariable("ip3d", &fMVAVar_ip3d);
350 
351  //isolation variables
352  tmpTMVAReader->AddVariable( "ChargedIso_DR0p0To0p1", &fMVAVar_ChargedIso_DR0p0To0p1 );
353  tmpTMVAReader->AddVariable( "ChargedIso_DR0p1To0p2", &fMVAVar_ChargedIso_DR0p1To0p2 );
354  tmpTMVAReader->AddVariable( "ChargedIso_DR0p2To0p3", &fMVAVar_ChargedIso_DR0p2To0p3 );
355  tmpTMVAReader->AddVariable( "ChargedIso_DR0p3To0p4", &fMVAVar_ChargedIso_DR0p3To0p4 );
356  tmpTMVAReader->AddVariable( "ChargedIso_DR0p4To0p5", &fMVAVar_ChargedIso_DR0p4To0p5 );
357  tmpTMVAReader->AddVariable( "GammaIso_DR0p0To0p1", &fMVAVar_GammaIso_DR0p0To0p1 );
358  tmpTMVAReader->AddVariable( "GammaIso_DR0p1To0p2", &fMVAVar_GammaIso_DR0p1To0p2 );
359  tmpTMVAReader->AddVariable( "GammaIso_DR0p2To0p3", &fMVAVar_GammaIso_DR0p2To0p3 );
360  tmpTMVAReader->AddVariable( "GammaIso_DR0p3To0p4", &fMVAVar_GammaIso_DR0p3To0p4 );
361  tmpTMVAReader->AddVariable( "GammaIso_DR0p4To0p5", &fMVAVar_GammaIso_DR0p4To0p5 );
362  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p0To0p1", &fMVAVar_NeutralHadronIso_DR0p0To0p1 );
363  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p1To0p2", &fMVAVar_NeutralHadronIso_DR0p1To0p2 );
364  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p2To0p3", &fMVAVar_NeutralHadronIso_DR0p2To0p3 );
365  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p3To0p4", &fMVAVar_NeutralHadronIso_DR0p3To0p4 );
366  tmpTMVAReader->AddVariable( "NeutralHadronIso_DR0p4To0p5", &fMVAVar_NeutralHadronIso_DR0p4To0p5 );
367  tmpTMVAReader->AddVariable( "rho", &fMVAVar_rho );
368 
369  //spectators
370  tmpTMVAReader->AddSpectator("eta", &fMVAVar_eta);
371  tmpTMVAReader->AddSpectator("pt", &fMVAVar_pt);
372 
373  }
374 
375 #ifndef STANDALONE
376  if ((fMethodname == "BDT") && (weightsfiles[i].rfind(".xml.gz") == weightsfiles[i].length()-strlen(".xml.gz"))) {
377  gzFile file = gzopen(weightsfiles[i].c_str(), "rb");
378  if (file == nullptr) { std::cout << "Error opening gzip file associated to " << weightsfiles[i] << std::endl; throw cms::Exception("Configuration","Error reading zipped XML file"); }
379  std::vector<char> data;
380  data.reserve(1024*1024*10);
381  unsigned int bufflen = 32*1024;
382  char *buff = reinterpret_cast<char *>(malloc(bufflen));
383  if (buff == nullptr) { std::cout << "Error creating buffer for " << weightsfiles[i] << std::endl; gzclose(file); throw cms::Exception("Configuration","Error reading zipped XML file"); }
384  int read;
385  while ((read = gzread(file, buff, bufflen)) != 0) {
386  if (read == -1) { std::cout << "Error reading gzip file associated to " << weightsfiles[i] << ": " << gzerror(file,&read) << std::endl; gzclose(file); free(buff); throw cms::Exception("Configuration","Error reading zipped XML file"); }
387  data.insert(data.end(), buff, buff+read);
388  }
389  if (gzclose(file) != Z_OK) { std::cout << "Error closing gzip file associated to " << weightsfiles[i] << std::endl; }
390  free(buff);
391  data.push_back('\0'); // IMPORTANT
392  fTMVAMethod.push_back(dynamic_cast<TMVA::MethodBase*>(tmpTMVAReader->BookMVA(TMVA::Types::kBDT, &data[0])));
393  } else {
394  fTMVAMethod.push_back(dynamic_cast<TMVA::MethodBase*>(tmpTMVAReader->BookMVA(fMethodname , weightsfiles[i])));
395  }
396 #else
397  fTMVAMethod.push_back(dynamic_cast<TMVA::MethodBase*>(tmpTMVAReader->BookMVA(fMethodname , weightsfiles[i])));
398 #endif
399  std::cout << "MVABin " << i << " : MethodName = " << fMethodname
400  << " , type == " << type << " , "
401  << "Load weights file : " << weightsfiles[i]
402  << std::endl;
403  fTMVAReader.push_back(tmpTMVAReader);
404  }
405  std::cout << "Electron ID MVA Completed\n";
406 
407 }
type
type
Definition: HCALResponse.h:21
i
int i
Definition: DBlmapReader.cc:9
assert
assert(m_qm.get())
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
SiPixelLorentzAngle_cfi.read
tuple read
Definition: SiPixelLorentzAngle_cfi.py:7
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
Bool_t EGammaMvaEleEstimator::isInitialized ( ) const
inline

Definition at line 61 of file EGammaMvaEleEstimator.h.

References fisInitialized.

61 { return fisInitialized; }
Double_t EGammaMvaEleEstimator::isoMvaValue ( const reco::GsfElectron &  ele,
const reco::Vertex &  vertex,
const reco::PFCandidateCollection &  PFCandidates,
double  Rho,
ElectronEffectiveArea::ElectronEffectiveAreaTarget  EATarget,
const reco::GsfElectronCollection &  IdentifiedElectrons,
const reco::MuonCollection &  IdentifiedMuons,
bool  printDebug = kFALSE 
)

Definition at line 1483 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, bindVariables(), reco::GsfElectron::closestCtfTrackRef(), funct::cos(), gather_cfg::cout, reco::PFCandidate::e, reco::LeafCandidate::eta(), fisInitialized, fMethodname, fMVAType, fMVAVar_ChargedIso_DR0p0To0p1, fMVAVar_ChargedIso_DR0p1To0p2, fMVAVar_ChargedIso_DR0p2To0p3, fMVAVar_ChargedIso_DR0p3To0p4, fMVAVar_ChargedIso_DR0p4To0p5, fMVAVar_eta, fMVAVar_GammaIso_DR0p0To0p1, fMVAVar_GammaIso_DR0p1To0p2, fMVAVar_GammaIso_DR0p2To0p3, fMVAVar_GammaIso_DR0p3To0p4, fMVAVar_GammaIso_DR0p4To0p5, fMVAVar_NeutralHadronIso_DR0p0To0p1, fMVAVar_NeutralHadronIso_DR0p1To0p2, fMVAVar_NeutralHadronIso_DR0p2To0p3, fMVAVar_NeutralHadronIso_DR0p3To0p4, fMVAVar_NeutralHadronIso_DR0p4To0p5, fMVAVar_pt, fTMVAMethod, fTMVAReader, fUseBinnedVersion, reco::PFCandidate::gamma, ElectronEffectiveArea::GetElectronEffectiveArea(), GetMVABin(), reco::GsfElectron::gsfTrack(), edm::Ref< C, T, F >::isNonnull(), ElectronEffectiveArea::kEleGammaIsoDR0p0To0p1, ElectronEffectiveArea::kEleGammaIsoDR0p1To0p2, ElectronEffectiveArea::kEleGammaIsoDR0p2To0p3, ElectronEffectiveArea::kEleGammaIsoDR0p3To0p4, ElectronEffectiveArea::kEleGammaIsoDR0p4To0p5, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p0To0p1, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p1To0p2, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p2To0p3, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p3To0p4, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p4To0p5, Max(), Min(), reco::PFCandidate::mu, reco::LeafCandidate::phi(), reco::Vertex::position(), funct::pow(), reco::LeafCandidate::pt(), edm::refToPtr(), mathSSE::sqrt(), and reco::GsfElectron::superCluster().

1490  {
1491 
1492  if (!fisInitialized) {
1493  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
1494  return -9999;
1495  }
1496 
1497  // Spectators
1498  fMVAVar_eta = ele.superCluster()->eta();
1499  fMVAVar_pt = ele.pt();
1500 
1501  //**********************************************************
1502  //Isolation variables
1503  //**********************************************************
1504  Double_t tmpChargedIso_DR0p0To0p1 = 0;
1505  Double_t tmpChargedIso_DR0p1To0p2 = 0;
1506  Double_t tmpChargedIso_DR0p2To0p3 = 0;
1507  Double_t tmpChargedIso_DR0p3To0p4 = 0;
1508  Double_t tmpChargedIso_DR0p4To0p5 = 0;
1509  Double_t tmpGammaIso_DR0p0To0p1 = 0;
1510  Double_t tmpGammaIso_DR0p1To0p2 = 0;
1511  Double_t tmpGammaIso_DR0p2To0p3 = 0;
1512  Double_t tmpGammaIso_DR0p3To0p4 = 0;
1513  Double_t tmpGammaIso_DR0p4To0p5 = 0;
1514  Double_t tmpNeutralHadronIso_DR0p0To0p1 = 0;
1515  Double_t tmpNeutralHadronIso_DR0p1To0p2 = 0;
1516  Double_t tmpNeutralHadronIso_DR0p2To0p3 = 0;
1517  Double_t tmpNeutralHadronIso_DR0p3To0p4 = 0;
1518  Double_t tmpNeutralHadronIso_DR0p4To0p5 = 0;
1519 
1520  double electronTrackZ = 0;
1521  if (ele.gsfTrack().isNonnull()) {
1522  electronTrackZ = ele.gsfTrack()->dz(vertex.position());
1523  } else if (ele.closestCtfTrackRef().isNonnull()) {
1524  electronTrackZ = ele.closestCtfTrackRef()->dz(vertex.position());
1525  }
1526 
1527  for (reco::PFCandidateCollection::const_iterator iP = PFCandidates.begin();
1528  iP != PFCandidates.end(); ++iP) {
1529 
1530  //exclude the electron itself
1531  if(iP->gsfTrackRef().isNonnull() && ele.gsfTrack().isNonnull() &&
1532  refToPtr(iP->gsfTrackRef()) == refToPtr(ele.gsfTrack())) continue;
1533  if(iP->trackRef().isNonnull() && ele.closestCtfTrackRef().isNonnull() &&
1534  refToPtr(iP->trackRef()) == refToPtr(ele.closestCtfTrackRef())) continue;
1535 
1536  //************************************************************
1537  // New Isolation Calculations
1538  //************************************************************
1539  double dr = sqrt(pow(iP->eta() - ele.eta(),2) + pow(acos(cos(iP->phi() - ele.phi())),2));
1540  //Double_t deta = (iP->eta() - ele.eta());
1541 
1542  if (dr < 1.0) {
1543  Bool_t IsLeptonFootprint = kFALSE;
1544  //************************************************************
1545  // Lepton Footprint Removal
1546  //************************************************************
1547  for (reco::GsfElectronCollection::const_iterator iE = IdentifiedElectrons.begin();
1548  iE != IdentifiedElectrons.end(); ++iE) {
1549  //if pf candidate matches an electron passing ID cuts, then veto it
1550  if(iP->gsfTrackRef().isNonnull() && iE->gsfTrack().isNonnull() &&
1551  refToPtr(iP->gsfTrackRef()) == refToPtr(iE->gsfTrack())) IsLeptonFootprint = kTRUE;
1552  if(iP->trackRef().isNonnull() && iE->closestCtfTrackRef().isNonnull() &&
1553  refToPtr(iP->trackRef()) == refToPtr(iE->closestCtfTrackRef())) IsLeptonFootprint = kTRUE;
1554 
1555  //if pf candidate lies in veto regions of electron passing ID cuts, then veto it
1556  double tmpDR = sqrt(pow(iP->eta() - iE->eta(),2) + pow(acos(cos(iP->phi() - iE->phi())),2));
1557  if(iP->trackRef().isNonnull() && fabs(iE->superCluster()->eta()) >= 1.479
1558  && tmpDR < 0.015) IsLeptonFootprint = kTRUE;
1559  if(iP->particleId() == reco::PFCandidate::gamma && fabs(iE->superCluster()->eta()) >= 1.479
1560  && tmpDR < 0.08) IsLeptonFootprint = kTRUE;
1561  }
1562  for (reco::MuonCollection::const_iterator iM = IdentifiedMuons.begin();
1563  iM != IdentifiedMuons.end(); ++iM) {
1564  //if pf candidate matches an muon passing ID cuts, then veto it
1565  if(iP->trackRef().isNonnull() && iM->innerTrack().isNonnull() &&
1566  refToPtr(iP->trackRef()) == refToPtr(iM->innerTrack())) IsLeptonFootprint = kTRUE;
1567 
1568  //if pf candidate lies in veto regions of muon passing ID cuts, then veto it
1569  double tmpDR = sqrt(pow(iP->eta() - iM->eta(),2) + pow(acos(cos(iP->phi() - iM->phi())),2));
1570  if(iP->trackRef().isNonnull() && tmpDR < 0.01) IsLeptonFootprint = kTRUE;
1571  }
1572 
1573  if (!IsLeptonFootprint) {
1574  Bool_t passVeto = kTRUE;
1575  //Charged
1576  if(iP->trackRef().isNonnull()) {
1577  if (!(fabs(iP->trackRef()->dz(vertex.position()) - electronTrackZ) < 0.2)) passVeto = kFALSE;
1578  //************************************************************
1579  // Veto any PFmuon, or PFEle
1580  if (iP->particleId() == reco::PFCandidate::e || iP->particleId() == reco::PFCandidate::mu) passVeto = kFALSE;
1581  //************************************************************
1582  //************************************************************
1583  // Footprint Veto
1584  if (fabs(fMVAVar_eta) > 1.479 && dr < 0.015) passVeto = kFALSE;
1585  //************************************************************
1586  if (passVeto) {
1587  if (dr < 0.1) tmpChargedIso_DR0p0To0p1 += iP->pt();
1588  if (dr >= 0.1 && dr < 0.2) tmpChargedIso_DR0p1To0p2 += iP->pt();
1589  if (dr >= 0.2 && dr < 0.3) tmpChargedIso_DR0p2To0p3 += iP->pt();
1590  if (dr >= 0.3 && dr < 0.4) tmpChargedIso_DR0p3To0p4 += iP->pt();
1591  if (dr >= 0.4 && dr < 0.5) tmpChargedIso_DR0p4To0p5 += iP->pt();
1592  } //pass veto
1593  }
1594  //Gamma
1595  else if (iP->particleId() == reco::PFCandidate::gamma) {
1596  //************************************************************
1597  // Footprint Veto
1598  if (fabs(fMVAVar_eta) > 1.479 && dr < 0.08) passVeto = kFALSE;
1599  //************************************************************
1600  if (passVeto) {
1601  if (dr < 0.1) tmpGammaIso_DR0p0To0p1 += iP->pt();
1602  if (dr >= 0.1 && dr < 0.2) tmpGammaIso_DR0p1To0p2 += iP->pt();
1603  if (dr >= 0.2 && dr < 0.3) tmpGammaIso_DR0p2To0p3 += iP->pt();
1604  if (dr >= 0.3 && dr < 0.4) tmpGammaIso_DR0p3To0p4 += iP->pt();
1605  if (dr >= 0.4 && dr < 0.5) tmpGammaIso_DR0p4To0p5 += iP->pt();
1606  }
1607  }
1608  //NeutralHadron
1609  else {
1610  if (dr < 0.1) tmpNeutralHadronIso_DR0p0To0p1 += iP->pt();
1611  if (dr >= 0.1 && dr < 0.2) tmpNeutralHadronIso_DR0p1To0p2 += iP->pt();
1612  if (dr >= 0.2 && dr < 0.3) tmpNeutralHadronIso_DR0p2To0p3 += iP->pt();
1613  if (dr >= 0.3 && dr < 0.4) tmpNeutralHadronIso_DR0p3To0p4 += iP->pt();
1614  if (dr >= 0.4 && dr < 0.5) tmpNeutralHadronIso_DR0p4To0p5 += iP->pt();
1615  }
1616  } //not lepton footprint
1617  } //in 1.0 dr cone
1618  } //loop over PF candidates
1619 
1620  fMVAVar_ChargedIso_DR0p0To0p1 = TMath::Min((tmpChargedIso_DR0p0To0p1)/ele.pt(), 2.5);
1621  fMVAVar_ChargedIso_DR0p1To0p2 = TMath::Min((tmpChargedIso_DR0p1To0p2)/ele.pt(), 2.5);
1622  fMVAVar_ChargedIso_DR0p2To0p3 = TMath::Min((tmpChargedIso_DR0p2To0p3)/ele.pt(), 2.5);
1623  fMVAVar_ChargedIso_DR0p3To0p4 = TMath::Min((tmpChargedIso_DR0p3To0p4)/ele.pt(), 2.5);
1624  fMVAVar_ChargedIso_DR0p4To0p5 = TMath::Min((tmpChargedIso_DR0p4To0p5)/ele.pt(), 2.5);
1625  fMVAVar_GammaIso_DR0p0To0p1 = TMath::Max(TMath::Min((tmpGammaIso_DR0p0To0p1 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p0To0p1, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1626  fMVAVar_GammaIso_DR0p1To0p2 = TMath::Max(TMath::Min((tmpGammaIso_DR0p1To0p2 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p1To0p2, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1627  fMVAVar_GammaIso_DR0p2To0p3 = TMath::Max(TMath::Min((tmpGammaIso_DR0p2To0p3 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p2To0p3, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1628  fMVAVar_GammaIso_DR0p3To0p4 = TMath::Max(TMath::Min((tmpGammaIso_DR0p3To0p4 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p3To0p4, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1629  fMVAVar_GammaIso_DR0p4To0p5 = TMath::Max(TMath::Min((tmpGammaIso_DR0p4To0p5 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p4To0p5, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1630  fMVAVar_NeutralHadronIso_DR0p0To0p1 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p0To0p1 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p0To0p1, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1631  fMVAVar_NeutralHadronIso_DR0p1To0p2 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p1To0p2 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p1To0p2, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1632  fMVAVar_NeutralHadronIso_DR0p2To0p3 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p2To0p3 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p2To0p3, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1633  fMVAVar_NeutralHadronIso_DR0p3To0p4 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p3To0p4 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p3To0p4, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1634  fMVAVar_NeutralHadronIso_DR0p4To0p5 = TMath::Max(TMath::Min((tmpNeutralHadronIso_DR0p4To0p5 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p4To0p5, fMVAVar_eta, EATarget))/ele.pt(), 2.5), 0.0);
1635 
1636  if (printDebug) {
1637  cout << "UseBinnedVersion=" << fUseBinnedVersion << " -> BIN: " << fMVAVar_eta << " " << fMVAVar_pt << " : " << GetMVABin(fMVAVar_eta,fMVAVar_pt) << endl;
1638  }
1639 
1640  // evaluate
1641  bindVariables();
1642  Double_t mva = -9999;
1643 
1644 // mva = fTMVAReader[0]->EvaluateMVA(fMethodname);
1645  if (fUseBinnedVersion) {
1646  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
1647  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
1648  } else {
1649  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
1650  }
1651 
1652 
1653  if(printDebug) {
1654  cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType << endl;
1655  cout << "ChargedIso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
1656  << fMVAVar_ChargedIso_DR0p0To0p1 << " "
1657  << fMVAVar_ChargedIso_DR0p1To0p2 << " "
1658  << fMVAVar_ChargedIso_DR0p2To0p3 << " "
1659  << fMVAVar_ChargedIso_DR0p3To0p4 << " "
1660  << fMVAVar_ChargedIso_DR0p4To0p5 << endl;
1661  cout << "PF Gamma Iso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
1662  << fMVAVar_GammaIso_DR0p0To0p1 << " "
1663  << fMVAVar_GammaIso_DR0p1To0p2 << " "
1664  << fMVAVar_GammaIso_DR0p2To0p3 << " "
1665  << fMVAVar_GammaIso_DR0p3To0p4 << " "
1666  << fMVAVar_GammaIso_DR0p4To0p5 << endl;
1667  cout << "PF Neutral Hadron Iso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
1668  << fMVAVar_NeutralHadronIso_DR0p0To0p1 << " "
1669  << fMVAVar_NeutralHadronIso_DR0p1To0p2 << " "
1670  << fMVAVar_NeutralHadronIso_DR0p2To0p3 << " "
1671  << fMVAVar_NeutralHadronIso_DR0p3To0p4 << " "
1672  << fMVAVar_NeutralHadronIso_DR0p4To0p5 << " "
1673  << endl;
1674  cout << " ### MVA " << mva << endl;
1675  }
1676 
1677 
1678  return mva;
1679 }
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
edm::refToPtr
Ptr< typename C::value_type > refToPtr(Ref< C, typename C::value_type, refhelper::FindUsingAdvance< C, typename C::value_type > > const &ref)
Definition: RefToPtr.h:18
ElectronEffectiveArea::GetElectronEffectiveArea
static Double_t GetElectronEffectiveArea(ElectronEffectiveAreaType type, Double_t SCEta, ElectronEffectiveAreaTarget EffectiveAreaTarget=kEleEAData2011)
Definition: ElectronEffectiveArea.h:64
Min
T Min(T a, T b)
Definition: MathUtil.h:39
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:106
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:137
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
reco::GsfElectron::closestCtfTrackRef
TrackRef closestCtfTrackRef() const
Definition: GsfElectron.h:199
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::GsfElectron::superCluster
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
Definition: GsfElectron.h:182
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
Max
T Max(T a, T b)
Definition: MathUtil.h:44
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:133
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:183
benchmark_cfg.PFCandidates
tuple PFCandidates
Definition: benchmark_cfg.py:32
Double_t EGammaMvaEleEstimator::isoMvaValue ( Double_t  Pt,
Double_t  Eta,
Double_t  Rho,
ElectronEffectiveArea::ElectronEffectiveAreaTarget  EATarget,
Double_t  ChargedIso_DR0p0To0p1,
Double_t  ChargedIso_DR0p1To0p2,
Double_t  ChargedIso_DR0p2To0p3,
Double_t  ChargedIso_DR0p3To0p4,
Double_t  ChargedIso_DR0p4To0p5,
Double_t  GammaIso_DR0p0To0p1,
Double_t  GammaIso_DR0p1To0p2,
Double_t  GammaIso_DR0p2To0p3,
Double_t  GammaIso_DR0p3To0p4,
Double_t  GammaIso_DR0p4To0p5,
Double_t  NeutralHadronIso_DR0p0To0p1,
Double_t  NeutralHadronIso_DR0p1To0p2,
Double_t  NeutralHadronIso_DR0p2To0p3,
Double_t  NeutralHadronIso_DR0p3To0p4,
Double_t  NeutralHadronIso_DR0p4To0p5,
Bool_t  printDebug = kFALSE 
)

Definition at line 921 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, gather_cfg::cout, fisInitialized, fMethodname, fMVAType, fMVAVar_ChargedIso_DR0p0To0p1, fMVAVar_ChargedIso_DR0p1To0p2, fMVAVar_ChargedIso_DR0p2To0p3, fMVAVar_ChargedIso_DR0p3To0p4, fMVAVar_ChargedIso_DR0p4To0p5, fMVAVar_GammaIso_DR0p0To0p1, fMVAVar_GammaIso_DR0p1To0p2, fMVAVar_GammaIso_DR0p2To0p3, fMVAVar_GammaIso_DR0p3To0p4, fMVAVar_GammaIso_DR0p4To0p5, fMVAVar_NeutralHadronIso_DR0p0To0p1, fMVAVar_NeutralHadronIso_DR0p1To0p2, fMVAVar_NeutralHadronIso_DR0p2To0p3, fMVAVar_NeutralHadronIso_DR0p3To0p4, fMVAVar_NeutralHadronIso_DR0p4To0p5, fTMVAMethod, fTMVAReader, ElectronEffectiveArea::GetElectronEffectiveArea(), GetMVABin(), ElectronEffectiveArea::kEleGammaIsoDR0p0To0p1, ElectronEffectiveArea::kEleGammaIsoDR0p1To0p2, ElectronEffectiveArea::kEleGammaIsoDR0p2To0p3, ElectronEffectiveArea::kEleGammaIsoDR0p3To0p4, ElectronEffectiveArea::kEleGammaIsoDR0p4To0p5, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p0To0p1, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p1To0p2, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p2To0p3, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p3To0p4, ElectronEffectiveArea::kEleNeutralHadronIsoDR0p4To0p5, Max(), and Min().

940  {
941 
942  if (!fisInitialized) {
943  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
944  return -9999;
945  }
946 
947  fMVAVar_ChargedIso_DR0p0To0p1 = TMath::Min((ChargedIso_DR0p0To0p1)/Pt, 2.5);
948  fMVAVar_ChargedIso_DR0p1To0p2 = TMath::Min((ChargedIso_DR0p1To0p2)/Pt, 2.5);
949  fMVAVar_ChargedIso_DR0p2To0p3 = TMath::Min((ChargedIso_DR0p2To0p3)/Pt, 2.5);
950  fMVAVar_ChargedIso_DR0p3To0p4 = TMath::Min((ChargedIso_DR0p3To0p4)/Pt, 2.5);
951  fMVAVar_ChargedIso_DR0p4To0p5 = TMath::Min((ChargedIso_DR0p4To0p5)/Pt, 2.5);
952  fMVAVar_GammaIso_DR0p0To0p1 = TMath::Max(TMath::Min((GammaIso_DR0p0To0p1 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p0To0p1, Eta, EATarget))/Pt, 2.5), 0.0);
953  fMVAVar_GammaIso_DR0p1To0p2 = TMath::Max(TMath::Min((GammaIso_DR0p1To0p2 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p1To0p2, Eta, EATarget))/Pt, 2.5), 0.0);
954  fMVAVar_GammaIso_DR0p2To0p3 = TMath::Max(TMath::Min((GammaIso_DR0p2To0p3 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p2To0p3, Eta, EATarget))/Pt, 2.5), 0.0);
955  fMVAVar_GammaIso_DR0p3To0p4 = TMath::Max(TMath::Min((GammaIso_DR0p3To0p4 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p3To0p4, Eta, EATarget))/Pt, 2.5), 0.0);
956  fMVAVar_GammaIso_DR0p4To0p5 = TMath::Max(TMath::Min((GammaIso_DR0p4To0p5 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaIsoDR0p4To0p5, Eta, EATarget))/Pt, 2.5), 0.0);
957  fMVAVar_NeutralHadronIso_DR0p0To0p1 = TMath::Max(TMath::Min((NeutralHadronIso_DR0p0To0p1 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p0To0p1, Eta, EATarget))/Pt, 2.5), 0.0);
958  fMVAVar_NeutralHadronIso_DR0p1To0p2 = TMath::Max(TMath::Min((NeutralHadronIso_DR0p1To0p2 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p1To0p2, Eta, EATarget))/Pt, 2.5), 0.0);
959  fMVAVar_NeutralHadronIso_DR0p2To0p3 = TMath::Max(TMath::Min((NeutralHadronIso_DR0p2To0p3 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p2To0p3, Eta, EATarget))/Pt, 2.5), 0.0);
960  fMVAVar_NeutralHadronIso_DR0p3To0p4 = TMath::Max(TMath::Min((NeutralHadronIso_DR0p3To0p4 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p3To0p4, Eta, EATarget))/Pt, 2.5), 0.0);
961  fMVAVar_NeutralHadronIso_DR0p4To0p5 = TMath::Max(TMath::Min((NeutralHadronIso_DR0p4To0p5 - Rho*ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleNeutralHadronIsoDR0p4To0p5, Eta, EATarget))/Pt, 2.5), 0.0);
962 
963  // evaluate
964  int bin = GetMVABin(Eta,Pt);
965  Double_t mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
966 
967  if(printDebug) {
968  cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType << endl;
969  cout << "ChargedIso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
970  << fMVAVar_ChargedIso_DR0p0To0p1 << " "
971  << fMVAVar_ChargedIso_DR0p1To0p2 << " "
972  << fMVAVar_ChargedIso_DR0p2To0p3 << " "
973  << fMVAVar_ChargedIso_DR0p3To0p4 << " "
974  << fMVAVar_ChargedIso_DR0p4To0p5 << endl;
975  cout << "PF Gamma Iso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
976  << fMVAVar_GammaIso_DR0p0To0p1 << " "
977  << fMVAVar_GammaIso_DR0p1To0p2 << " "
978  << fMVAVar_GammaIso_DR0p2To0p3 << " "
979  << fMVAVar_GammaIso_DR0p3To0p4 << " "
980  << fMVAVar_GammaIso_DR0p4To0p5 << endl;
981  cout << "PF Neutral Hadron Iso ( 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 ): "
982  << fMVAVar_NeutralHadronIso_DR0p0To0p1 << " "
983  << fMVAVar_NeutralHadronIso_DR0p1To0p2 << " "
984  << fMVAVar_NeutralHadronIso_DR0p2To0p3 << " "
985  << fMVAVar_NeutralHadronIso_DR0p3To0p4 << " "
986  << fMVAVar_NeutralHadronIso_DR0p4To0p5 << " "
987  << endl;
988  cout << " ### MVA " << mva << endl;
989  }
990 
991  return mva;
992 
993 }
ElectronEffectiveArea::GetElectronEffectiveArea
static Double_t GetElectronEffectiveArea(ElectronEffectiveAreaType type, Double_t SCEta, ElectronEffectiveAreaTarget EffectiveAreaTarget=kEleEAData2011)
Definition: ElectronEffectiveArea.h:64
Min
T Min(T a, T b)
Definition: MathUtil.h:39
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
Max
T Max(T a, T b)
Definition: MathUtil.h:44
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
Double_t EGammaMvaEleEstimator::mvaValue ( const reco::GsfElectron &  ele,
const reco::Vertex &  vertex,
const TransientTrackBuilder &  transientTrackBuilder,
EcalClusterLazyTools  myEcalCluster,
bool  printDebug = kFALSE 
)

Definition at line 999 of file EGammaMvaEleEstimator.cc.

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

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

1003  {
1004 
1005  if (!fisInitialized) {
1006  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
1007  return -9999;
1008  }
1009 
1010  if ( (fMVAType != EGammaMvaEleEstimator::kTrig) && (fMVAType != EGammaMvaEleEstimator::kNonTrig )) {
1011  std::cout << "Error: This method should be called for kTrig or kNonTrig MVA only" << endl;
1012  return -9999;
1013  }
1014 
1015  bool validKF= false;
1016  reco::TrackRef myTrackRef = ele.closestCtfTrackRef();
1017  validKF = (myTrackRef.isAvailable());
1018  validKF = (myTrackRef.isNonnull());
1019 
1020  // Pure tracking variables
1021  fMVAVar_fbrem = ele.fbrem();
1022  fMVAVar_kfchi2 = (validKF) ? myTrackRef->normalizedChi2() : 0 ;
1023  fMVAVar_kfhits = (validKF) ? myTrackRef->hitPattern().trackerLayersWithMeasurement() : -1. ;
1024  fMVAVar_kfhitsall = (validKF) ? myTrackRef->numberOfValidHits() : -1. ; // save also this in your ntuple as possible alternative
1025  fMVAVar_gsfchi2 = ele.gsfTrack()->normalizedChi2();
1026 
1027 
1028  // Geometrical matchings
1029  fMVAVar_deta = ele.deltaEtaSuperClusterTrackAtVtx();
1030  fMVAVar_dphi = ele.deltaPhiSuperClusterTrackAtVtx();
1031  fMVAVar_detacalo = ele.deltaEtaSeedClusterTrackAtCalo();
1032 
1033 
1034  // Pure ECAL -> shower shapes
1035  fMVAVar_see = ele.sigmaIetaIeta(); //EleSigmaIEtaIEta
1036  std::vector<float> vCov = myEcalCluster.localCovariances(*(ele.superCluster()->seed())) ;
1037  if (!isnan(vCov[2])) fMVAVar_spp = sqrt (vCov[2]); //EleSigmaIPhiIPhi
1038  else fMVAVar_spp = 0.;
1039 
1040  fMVAVar_etawidth = ele.superCluster()->etaWidth();
1041  fMVAVar_phiwidth = ele.superCluster()->phiWidth();
1042  fMVAVar_OneMinusE1x5E5x5 = (ele.e5x5()) !=0. ? 1.-(ele.e1x5()/ele.e5x5()) : -1. ;
1043  fMVAVar_R9 = myEcalCluster.e3x3(*(ele.superCluster()->seed())) / ele.superCluster()->rawEnergy();
1044 
1045  // Energy matching
1046  fMVAVar_HoE = ele.hadronicOverEm();
1047  fMVAVar_EoP = ele.eSuperClusterOverP();
1048  fMVAVar_IoEmIoP = (1.0/ele.ecalEnergy()) - (1.0 / ele.p()); // in the future to be changed with ele.gsfTrack()->p()
1049  fMVAVar_eleEoPout = ele.eEleClusterOverPout();
1050  fMVAVar_PreShowerOverRaw= ele.superCluster()->preshowerEnergy() / ele.superCluster()->rawEnergy();
1051 
1052 
1053  // Spectators
1054  fMVAVar_eta = ele.superCluster()->eta();
1055  fMVAVar_pt = ele.pt();
1056 
1057 
1058 
1059  // for triggering electrons get the impact parameteres
1060  if(fMVAType == kTrig) {
1061  //d0
1062  if (ele.gsfTrack().isNonnull()) {
1063  fMVAVar_d0 = (-1.0)*ele.gsfTrack()->dxy(vertex.position());
1064  } else if (ele.closestCtfTrackRef().isNonnull()) {
1065  fMVAVar_d0 = (-1.0)*ele.closestCtfTrackRef()->dxy(vertex.position());
1066  } else {
1067  fMVAVar_d0 = -9999.0;
1068  }
1069 
1070  //default values for IP3D
1071  fMVAVar_ip3d = -999.0;
1072  fMVAVar_ip3dSig = 0.0;
1073  if (ele.gsfTrack().isNonnull()) {
1074  const double gsfsign = ( (-ele.gsfTrack()->dxy(vertex.position())) >=0 ) ? 1. : -1.;
1075 
1076  const reco::TransientTrack &tt = transientTrackBuilder.build(ele.gsfTrack());
1077  const std::pair<bool,Measurement1D> &ip3dpv = IPTools::absoluteImpactParameter3D(tt,vertex);
1078  if (ip3dpv.first) {
1079  double ip3d = gsfsign*ip3dpv.second.value();
1080  double ip3derr = ip3dpv.second.error();
1081  fMVAVar_ip3d = ip3d;
1082  fMVAVar_ip3dSig = ip3d/ip3derr;
1083  }
1084  }
1085  }
1086 
1087 
1088  // evaluate
1089  bindVariables();
1090  Double_t mva = -9999;
1091  if (fUseBinnedVersion) {
1092  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
1093  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
1094  } else {
1095  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
1096  }
1097 
1098 
1099 
1100  if(printDebug) {
1101  cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType << endl;
1102  cout << " fbrem " << fMVAVar_fbrem
1103  << " kfchi2 " << fMVAVar_kfchi2
1104  << " mykfhits " << fMVAVar_kfhits
1105  << " gsfchi2 " << fMVAVar_gsfchi2
1106  << " deta " << fMVAVar_deta
1107  << " dphi " << fMVAVar_dphi
1108  << " detacalo " << fMVAVar_detacalo
1109  << " see " << fMVAVar_see
1110  << " spp " << fMVAVar_spp
1111  << " etawidth " << fMVAVar_etawidth
1112  << " phiwidth " << fMVAVar_phiwidth
1113  << " OneMinusE1x5E5x5 " << fMVAVar_OneMinusE1x5E5x5
1114  << " R9 " << fMVAVar_R9
1115  << " HoE " << fMVAVar_HoE
1116  << " EoP " << fMVAVar_EoP
1117  << " IoEmIoP " << fMVAVar_IoEmIoP
1118  << " eleEoPout " << fMVAVar_eleEoPout
1119  << " d0 " << fMVAVar_d0
1120  << " ip3d " << fMVAVar_ip3d
1121  << " eta " << fMVAVar_eta
1122  << " pt " << fMVAVar_pt << endl;
1123  cout << " ### MVA " << mva << endl;
1124  }
1125 
1126 
1127 
1128  return mva;
1129 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:576
reco::LeafCandidate::p
virtual double p() const
magnitude of momentum vector
Definition: LeafCandidate.h:108
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:243
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:10
IPTools::absoluteImpactParameter3D
std::pair< bool, Measurement1D > absoluteImpactParameter3D(const reco::TransientTrack &transientTrack, const reco::Vertex &vertex)
Definition: IPTools.cc:37
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:698
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:106
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
reco::GsfElectron::closestCtfTrackRef
TrackRef closestCtfTrackRef() const
Definition: GsfElectron.h:199
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:247
reco::GsfElectron::sigmaIetaIeta
float sigmaIetaIeta() const
Definition: GsfElectron.h:416
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:457
edm::detail::isnan
bool isnan(float x)
Definition: math.h:13
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::GsfElectron::superCluster
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
Definition: GsfElectron.h:182
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:250
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:246
reco::GsfElectron::e1x5
float e1x5() const
Definition: GsfElectron.h:418
reco::GsfElectron::ecalEnergy
float ecalEnergy() const
Definition: GsfElectron.h:785
reco::GsfElectron::e5x5
float e5x5() const
Definition: GsfElectron.h:420
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:248
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:183
Double_t EGammaMvaEleEstimator::mvaValue ( const reco::GsfElectron &  ele,
const reco::Vertex &  vertex,
double  rho,
EcalClusterLazyTools  myEcalCluster,
bool  printDebug = kFALSE 
)

Definition at line 1133 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, bindVariables(), reco::GsfElectron::closestCtfTrackRef(), gather_cfg::cout, reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), reco::GsfElectron::e1x5(), reco::GsfElectron::e5x5(), reco::GsfElectron::eEleClusterOverPout(), reco::GsfElectron::eSuperClusterOverP(), reco::GsfElectron::fbrem(), fisInitialized, fMethodname, fMVAType, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_rho, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, GetMVABin(), reco::GsfElectron::gsfTrack(), reco::GsfElectron::hadronicOverEm(), edm::Ref< C, T, F >::isAvailable(), edm::detail::isnan(), edm::Ref< C, T, F >::isNonnull(), kTrigNoIP, reco::LeafCandidate::pt(), rho, reco::GsfElectron::sigmaIetaIeta(), mathSSE::sqrt(), and reco::GsfElectron::superCluster().

1138  {
1139 
1140  if (!fisInitialized) {
1141  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
1142  return -9999;
1143  }
1144 
1145  if (fMVAType != EGammaMvaEleEstimator::kTrigNoIP) {
1146  std::cout << "Error: This method should be called for kTrigNoIP MVA only" << endl;
1147  return -9999;
1148  }
1149 
1150  bool validKF= false;
1151  reco::TrackRef myTrackRef = ele.closestCtfTrackRef();
1152  validKF = (myTrackRef.isAvailable());
1153  validKF = (myTrackRef.isNonnull());
1154 
1155  // Pure tracking variables
1156  fMVAVar_fbrem = ele.fbrem();
1157  fMVAVar_kfchi2 = (validKF) ? myTrackRef->normalizedChi2() : 0 ;
1158  fMVAVar_kfhits = (validKF) ? myTrackRef->hitPattern().trackerLayersWithMeasurement() : -1. ;
1159  fMVAVar_gsfchi2 = ele.gsfTrack()->normalizedChi2();
1160 
1161 
1162  // Geometrical matchings
1163  fMVAVar_deta = ele.deltaEtaSuperClusterTrackAtVtx();
1164  fMVAVar_dphi = ele.deltaPhiSuperClusterTrackAtVtx();
1165  fMVAVar_detacalo = ele.deltaEtaSeedClusterTrackAtCalo();
1166 
1167 
1168  // Pure ECAL -> shower shapes
1169  fMVAVar_see = ele.sigmaIetaIeta(); //EleSigmaIEtaIEta
1170  std::vector<float> vCov = myEcalCluster.localCovariances(*(ele.superCluster()->seed())) ;
1171  if (!isnan(vCov[2])) fMVAVar_spp = sqrt (vCov[2]); //EleSigmaIPhiIPhi
1172  else fMVAVar_spp = 0.;
1173 
1174 
1175  fMVAVar_etawidth = ele.superCluster()->etaWidth();
1176  fMVAVar_phiwidth = ele.superCluster()->phiWidth();
1177  fMVAVar_OneMinusE1x5E5x5 = (ele.e5x5()) !=0. ? 1.-(ele.e1x5()/ele.e5x5()) : -1. ;
1178  fMVAVar_R9 = myEcalCluster.e3x3(*(ele.superCluster()->seed())) / ele.superCluster()->rawEnergy();
1179 
1180 
1181  // Energy matching
1182  fMVAVar_HoE = ele.hadronicOverEm();
1183  fMVAVar_EoP = ele.eSuperClusterOverP();
1184  fMVAVar_IoEmIoP = (1.0/ele.superCluster()->energy()) - (1.0 / ele.gsfTrack()->p()); // in the future to be changed with ele.gsfTrack()->p()
1185  fMVAVar_eleEoPout = ele.eEleClusterOverPout();
1186  fMVAVar_rho = rho;
1187  fMVAVar_PreShowerOverRaw= ele.superCluster()->preshowerEnergy() / ele.superCluster()->rawEnergy();
1188 
1189 
1190  // Spectators
1191  fMVAVar_eta = ele.superCluster()->eta();
1192  fMVAVar_pt = ele.pt();
1193 
1194 
1195 
1196 
1197 
1198  // evaluate
1199  bindVariables();
1200  Double_t mva = -9999;
1201  if (fUseBinnedVersion) {
1202  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
1203  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
1204  } else {
1205  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
1206  }
1207 
1208 
1209 
1210  if(printDebug) {
1211  cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType << endl;
1212  cout << " fbrem " << fMVAVar_fbrem
1213  << " kfchi2 " << fMVAVar_kfchi2
1214  << " mykfhits " << fMVAVar_kfhits
1215  << " gsfchi2 " << fMVAVar_gsfchi2
1216  << " deta " << fMVAVar_deta
1217  << " dphi " << fMVAVar_dphi
1218  << " detacalo " << fMVAVar_detacalo
1219  // << " dphicalo " << fMVAVar_dphicalo
1220  << " see " << fMVAVar_see
1221  << " spp " << fMVAVar_spp
1222  << " etawidth " << fMVAVar_etawidth
1223  << " phiwidth " << fMVAVar_phiwidth
1224  << " e1x5e5x5 " << fMVAVar_OneMinusE1x5E5x5
1225  << " R9 " << fMVAVar_R9
1226  // << " mynbrems " << fMVAVar_nbrems
1227  << " HoE " << fMVAVar_HoE
1228  << " EoP " << fMVAVar_EoP
1229  << " IoEmIoP " << fMVAVar_IoEmIoP
1230  << " eleEoPout " << fMVAVar_eleEoPout
1231  << " rho " << fMVAVar_rho
1232  // << " EoPout " << fMVAVar_EoPout
1233  << " eta " << fMVAVar_eta
1234  << " pt " << fMVAVar_pt << endl;
1235  cout << " ### MVA " << mva << endl;
1236  }
1237 
1238 
1239 
1240  return mva;
1241 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:576
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:243
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:698
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
reco::GsfElectron::closestCtfTrackRef
TrackRef closestCtfTrackRef() const
Definition: GsfElectron.h:199
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:247
reco::GsfElectron::sigmaIetaIeta
float sigmaIetaIeta() const
Definition: GsfElectron.h:416
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:457
edm::detail::isnan
bool isnan(float x)
Definition: math.h:13
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::GsfElectron::superCluster
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
Definition: GsfElectron.h:182
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:250
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:246
reco::GsfElectron::e1x5
float e1x5() const
Definition: GsfElectron.h:418
reco::GsfElectron::e5x5
float e5x5() const
Definition: GsfElectron.h:420
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:248
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:183
Double_t EGammaMvaEleEstimator::mvaValue ( const pat::Electron &  ele,
double  rho,
bool  printDebug = kFALSE 
)

Definition at line 1370 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, bindVariables(), pat::Electron::closestCtfTrackRef(), gather_cfg::cout, reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), reco::GsfElectron::e1x5(), reco::GsfElectron::e5x5(), reco::GsfElectron::eEleClusterOverPout(), reco::GsfElectron::eSuperClusterOverP(), reco::GsfElectron::fbrem(), fisInitialized, fMethodname, fMVAType, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_rho, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, GetMVABin(), pat::Electron::gsfTrack(), reco::GsfElectron::hadronicOverEm(), edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), kTrigNoIP, reco::LeafCandidate::pt(), reco::GsfElectron::r9(), rho, reco::GsfElectron::sigmaIetaIeta(), reco::GsfElectron::sigmaIphiIphi(), and pat::Electron::superCluster().

1372  {
1373 
1374  if (!fisInitialized) {
1375  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
1376  return -9999;
1377  }
1378 
1379  if ( (fMVAType != EGammaMvaEleEstimator::kTrigNoIP) ) {
1380  std::cout << "Error: This method should be called for kTrigNoIP mva only" << endl;
1381  return -9999;
1382  }
1383 
1384 
1385  bool validKF= false;
1386  reco::TrackRef myTrackRef = ele.closestCtfTrackRef();
1387  validKF = (myTrackRef.isAvailable());
1388  validKF = (myTrackRef.isNonnull());
1389 
1390  // Pure tracking variables
1391  fMVAVar_fbrem = ele.fbrem();
1392  fMVAVar_kfchi2 = (validKF) ? myTrackRef->normalizedChi2() : 0 ;
1393  fMVAVar_kfhits = (validKF) ? myTrackRef->hitPattern().trackerLayersWithMeasurement() : -1. ;
1394  fMVAVar_gsfchi2 = ele.gsfTrack()->normalizedChi2();
1395 
1396 
1397  // Geometrical matchings
1398  fMVAVar_deta = ele.deltaEtaSuperClusterTrackAtVtx();
1399  fMVAVar_dphi = ele.deltaPhiSuperClusterTrackAtVtx();
1400  fMVAVar_detacalo = ele.deltaEtaSeedClusterTrackAtCalo();
1401 
1402 
1403  // Pure ECAL -> shower shapes
1404  fMVAVar_see = ele.sigmaIetaIeta(); //EleSigmaIEtaIEta
1405 
1406  fMVAVar_spp = ele.sigmaIphiIphi();
1407 
1408  fMVAVar_etawidth = ele.superCluster()->etaWidth();
1409  fMVAVar_phiwidth = ele.superCluster()->phiWidth();
1410  fMVAVar_OneMinusE1x5E5x5 = (ele.e5x5()) !=0. ? 1.-(ele.e1x5()/ele.e5x5()) : -1. ;
1411  fMVAVar_R9 = ele.r9();
1412 
1413  // Energy matching
1414  fMVAVar_HoE = ele.hadronicOverEm();
1415  fMVAVar_EoP = ele.eSuperClusterOverP();
1416  fMVAVar_IoEmIoP = (1.0/ele.superCluster()->energy()) - (1.0 / ele.gsfTrack()->p()); // in the future to be changed with ele.gsfTrack()->p()
1417  fMVAVar_eleEoPout = ele.eEleClusterOverPout();
1418  fMVAVar_rho = rho;
1419  fMVAVar_PreShowerOverRaw= ele.superCluster()->preshowerEnergy() / ele.superCluster()->rawEnergy();
1420 
1421 
1422  // Spectators
1423  fMVAVar_eta = ele.superCluster()->eta();
1424  fMVAVar_pt = ele.pt();
1425 
1426 
1427 
1428 
1429 
1430  // evaluate
1431  bindVariables();
1432  Double_t mva = -9999;
1433  if (fUseBinnedVersion) {
1434  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
1435  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
1436  } else {
1437  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
1438  }
1439 
1440 
1441 
1442  if(printDebug) {
1443  cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType << endl;
1444  cout << " fbrem " << fMVAVar_fbrem
1445  << " kfchi2 " << fMVAVar_kfchi2
1446  << " mykfhits " << fMVAVar_kfhits
1447  << " gsfchi2 " << fMVAVar_gsfchi2
1448  << " deta " << fMVAVar_deta
1449  << " dphi " << fMVAVar_dphi
1450  << " detacalo " << fMVAVar_detacalo
1451  // << " dphicalo " << fMVAVar_dphicalo
1452  << " see " << fMVAVar_see
1453  << " spp " << fMVAVar_spp
1454  << " etawidth " << fMVAVar_etawidth
1455  << " phiwidth " << fMVAVar_phiwidth
1456  << " e1x5e5x5 " << fMVAVar_OneMinusE1x5E5x5
1457  << " R9 " << fMVAVar_R9
1458  // << " mynbrems " << fMVAVar_nbrems
1459  << " HoE " << fMVAVar_HoE
1460  << " EoP " << fMVAVar_EoP
1461  << " IoEmIoP " << fMVAVar_IoEmIoP
1462  << " eleEoPout " << fMVAVar_eleEoPout
1463  << " rho " << fMVAVar_rho
1464  // << " EoPout " << fMVAVar_EoPout
1465  << " eta " << fMVAVar_eta
1466  << " pt " << fMVAVar_pt << endl;
1467  cout << " ### MVA " << mva << endl;
1468  }
1469 
1470 
1471 
1472  return mva;
1473 }
reco::GsfElectron::sigmaIphiIphi
float sigmaIphiIphi() const
Definition: GsfElectron.h:417
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:576
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:243
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:698
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
pat::Electron::superCluster
reco::SuperClusterRef superCluster() const
override the reco::GsfElectron::superCluster method, to access the internal storage of the superclust...
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:247
reco::GsfElectron::sigmaIetaIeta
float sigmaIetaIeta() const
Definition: GsfElectron.h:416
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:457
pat::Electron::gsfTrack
reco::GsfTrackRef gsfTrack() const
override the reco::GsfElectron::gsfTrack method, to access the internal storage of the supercluster ...
pat::Electron::closestCtfTrackRef
reco::TrackRef closestCtfTrackRef() const
override the reco::GsfElectron::closestCtfTrackRef method, to access the internal storage of the trac...
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:250
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:246
reco::GsfElectron::e1x5
float e1x5() const
Definition: GsfElectron.h:418
reco::GsfElectron::e5x5
float e5x5() const
Definition: GsfElectron.h:420
reco::GsfElectron::r9
float r9() const
Definition: GsfElectron.h:421
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:248
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
Double_t EGammaMvaEleEstimator::mvaValue ( const pat::Electron &  ele,
const reco::Vertex &  vertex,
double  rho,
bool  useFull5x5 = kFALSE,
bool  printDebug = kFALSE 
)

Definition at line 1244 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, bindVariables(), pat::Electron::closestCtfTrackRef(), gather_cfg::cout, pat::Electron::dB(), reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo(), reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx(), reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx(), reco::GsfElectron::e1x5(), reco::GsfElectron::e5x5(), reco::GsfElectron::ecalEnergy(), pat::Electron::edB(), reco::GsfElectron::eEleClusterOverPout(), reco::GsfElectron::eSuperClusterOverP(), reco::GsfElectron::fbrem(), fisInitialized, fMethodname, fMVAType, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_ip3d, fMVAVar_ip3dSig, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_rho, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, reco::GsfElectron::full5x5_e1x5(), reco::GsfElectron::full5x5_e5x5(), reco::GsfElectron::full5x5_r9(), reco::GsfElectron::full5x5_sigmaIetaIeta(), reco::GsfElectron::full5x5_sigmaIphiIphi(), fUseBinnedVersion, GetMVABin(), pat::Electron::gsfTrack(), reco::GsfElectron::hadronicOverEm(), edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), kNonTrig, kTrig, reco::LeafCandidate::p(), reco::Vertex::position(), reco::LeafCandidate::pt(), reco::GsfElectron::r9(), rho, reco::GsfElectron::sigmaIetaIeta(), reco::GsfElectron::sigmaIphiIphi(), and pat::Electron::superCluster().

1248  {
1249 
1250  if (!fisInitialized) {
1251  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
1252  return -9999;
1253  }
1254 
1255  bool validKF= false;
1256  reco::TrackRef myTrackRef = ele.closestCtfTrackRef();
1257  validKF = (myTrackRef.isAvailable());
1258  validKF = (myTrackRef.isNonnull());
1259 
1260  // Pure tracking variables
1261  fMVAVar_fbrem = ele.fbrem();
1262  fMVAVar_kfchi2 = (validKF) ? myTrackRef->normalizedChi2() : 0 ;
1263  fMVAVar_kfhits = (validKF) ? myTrackRef->hitPattern().trackerLayersWithMeasurement() : -1. ;
1264  fMVAVar_gsfchi2 = ele.gsfTrack()->normalizedChi2();
1265 
1266 
1267  // Geometrical matchings
1268  fMVAVar_deta = ele.deltaEtaSuperClusterTrackAtVtx();
1269  fMVAVar_dphi = ele.deltaPhiSuperClusterTrackAtVtx();
1270  fMVAVar_detacalo = ele.deltaEtaSeedClusterTrackAtCalo();
1271 
1272  // Pure ECAL -> shower shapes
1273  fMVAVar_see = useFull5x5 ? ele.full5x5_sigmaIetaIeta() : ele.sigmaIetaIeta(); //EleSigmaIEtaIEta
1274  fMVAVar_spp = useFull5x5 ? ele.full5x5_sigmaIphiIphi() : ele.sigmaIphiIphi(); //EleSigmaIEtaIEta
1275 
1276  fMVAVar_etawidth = ele.superCluster()->etaWidth();
1277  fMVAVar_phiwidth = ele.superCluster()->phiWidth();
1278  fMVAVar_OneMinusE1x5E5x5 = useFull5x5 ? ( (ele.full5x5_e5x5()) !=0. ? 1.-(ele.full5x5_e1x5()/ele.full5x5_e5x5()) : -1. )
1279  : ( (ele.e5x5() ) !=0. ? 1.-(ele.e1x5() /ele.e5x5() ) : -1. ) ;
1280  fMVAVar_R9 = useFull5x5 ? ele.full5x5_r9() : ele.r9();
1281 
1282  // Energy matching
1283  fMVAVar_HoE = ele.hadronicOverEm(); // this is identical for both
1284  fMVAVar_EoP = ele.eSuperClusterOverP();
1285 
1286  // unify this in the future?
1287  if( fMVAType == kTrig || fMVAType == kNonTrig){
1288  fMVAVar_IoEmIoP = (1.0/ele.ecalEnergy()) - (1.0 / ele.p()); // in the future to be changed with ele.gsfTrack()->p()
1289  }else{
1290  fMVAVar_IoEmIoP = (1.0/ele.superCluster()->energy()) - (1.0 / ele.gsfTrack()->p()); // in the future to be changed with ele.gsfTrack()->p()
1291  }
1292  fMVAVar_eleEoPout = ele.eEleClusterOverPout();
1293  fMVAVar_rho = rho;
1294  fMVAVar_PreShowerOverRaw= ele.superCluster()->preshowerEnergy() / ele.superCluster()->rawEnergy();
1295 
1296  // for triggering electrons get the impact parameteres
1297  if(fMVAType == kTrig) {
1298  //d0
1299  if (ele.gsfTrack().isNonnull()) {
1300  fMVAVar_d0 = (-1.0)*ele.gsfTrack()->dxy(vertex.position());
1301  } else if (ele.closestCtfTrackRef().isNonnull()) {
1302  fMVAVar_d0 = (-1.0)*ele.closestCtfTrackRef()->dxy(vertex.position());
1303  } else {
1304  fMVAVar_d0 = -9999.0;
1305  }
1306 
1307  //default values for IP3D
1308  fMVAVar_ip3d = -999.0;
1309  fMVAVar_ip3dSig = 0.0;
1310  if (ele.gsfTrack().isNonnull()) {
1311  const double gsfsign = ( (-ele.gsfTrack()->dxy(vertex.position())) >=0 ) ? 1. : -1.;
1312 
1313  //std::cout << "Warning : if usePV = false when pat-electrons were produced, dB() was calculated with respect to beamspot while original mva uses primary vertex" << std::endl;
1314  double ip3d = gsfsign*ele.dB();
1315  double ip3derr = ele.edB();
1316  fMVAVar_ip3d = ip3d;
1317  fMVAVar_ip3dSig = ip3d/ip3derr;
1318  }
1319  }
1320 
1321  // Spectators
1322  fMVAVar_eta = ele.superCluster()->eta();
1323  fMVAVar_pt = ele.pt();
1324 
1325  // evaluate
1326  bindVariables();
1327  Double_t mva = -9999;
1328  if (fUseBinnedVersion) {
1329  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
1330  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
1331  } else {
1332  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
1333  }
1334 
1335  if(printDebug) {
1336  cout << " *** Inside the class fMethodname " << fMethodname << " fMVAType " << fMVAType << endl;
1337  cout << " fbrem " << fMVAVar_fbrem
1338  << " kfchi2 " << fMVAVar_kfchi2
1339  << " mykfhits " << fMVAVar_kfhits
1340  << " gsfchi2 " << fMVAVar_gsfchi2
1341  << " deta " << fMVAVar_deta
1342  << " dphi " << fMVAVar_dphi
1343  << " detacalo " << fMVAVar_detacalo
1344  // << " dphicalo " << fMVAVar_dphicalo
1345  << " see " << fMVAVar_see
1346  << " spp " << fMVAVar_spp
1347  << " etawidth " << fMVAVar_etawidth
1348  << " phiwidth " << fMVAVar_phiwidth
1349  << " e1x5e5x5 " << fMVAVar_OneMinusE1x5E5x5
1350  << " R9 " << fMVAVar_R9
1351  // << " mynbrems " << fMVAVar_nbrems
1352  << " HoE " << fMVAVar_HoE
1353  << " EoP " << fMVAVar_EoP
1354  << " IoEmIoP " << fMVAVar_IoEmIoP
1355  << " eleEoPout " << fMVAVar_eleEoPout
1356  << " rho " << fMVAVar_rho
1357  // << " EoPout " << fMVAVar_EoPout
1358  << " d0 " << fMVAVar_d0
1359  << " ip3d " << fMVAVar_ip3d
1360  << " eta " << fMVAVar_eta
1361  << " pt " << fMVAVar_pt << endl;
1362  cout << " ### MVA " << mva << endl;
1363  }
1364 
1365  return mva;
1366 }
reco::GsfElectron::sigmaIphiIphi
float sigmaIphiIphi() const
Definition: GsfElectron.h:417
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:576
reco::LeafCandidate::p
virtual double p() const
magnitude of momentum vector
Definition: LeafCandidate.h:108
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:243
reco::GsfElectron::full5x5_e5x5
float full5x5_e5x5() const
Definition: GsfElectron.h:437
reco::GsfElectron::full5x5_e1x5
float full5x5_e1x5() const
Definition: GsfElectron.h:435
reco::GsfElectron::full5x5_sigmaIphiIphi
float full5x5_sigmaIphiIphi() const
Definition: GsfElectron.h:434
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:698
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:106
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
reco::GsfElectron::full5x5_sigmaIetaIeta
float full5x5_sigmaIetaIeta() const
Definition: GsfElectron.h:433
pat::Electron::superCluster
reco::SuperClusterRef superCluster() const
override the reco::GsfElectron::superCluster method, to access the internal storage of the superclust...
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:247
reco::GsfElectron::sigmaIetaIeta
float sigmaIetaIeta() const
Definition: GsfElectron.h:416
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:457
pat::Electron::gsfTrack
reco::GsfTrackRef gsfTrack() const
override the reco::GsfElectron::gsfTrack method, to access the internal storage of the supercluster ...
pat::Electron::closestCtfTrackRef
reco::TrackRef closestCtfTrackRef() const
override the reco::GsfElectron::closestCtfTrackRef method, to access the internal storage of the trac...
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:250
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:246
reco::GsfElectron::ecalEnergy
float ecalEnergy() const
Definition: GsfElectron.h:785
reco::GsfElectron::full5x5_r9
float full5x5_r9() const
Definition: GsfElectron.h:438
pat::Electron::dB
double dB(IPTYPE type) const
Impact parameter wrt primary vertex or beamspot.
reco::GsfElectron::r9
float r9() const
Definition: GsfElectron.h:421
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:248
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
pat::Electron::edB
double edB(IPTYPE type) const
Uncertainty on the corresponding impact parameter.
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
Double_t EGammaMvaEleEstimator::mvaValue ( Double_t  fbrem,
Double_t  kfchi2,
Int_t  kfhits,
Double_t  gsfchi2,
Double_t  deta,
Double_t  dphi,
Double_t  detacalo,
Double_t  see,
Double_t  spp,
Double_t  etawidth,
Double_t  phiwidth,
Double_t  e1x5e5x5,
Double_t  R9,
Double_t  HoE,
Double_t  EoP,
Double_t  IoEmIoP,
Double_t  eleEoPout,
Double_t  PreShowerOverRaw,
Double_t  d0,
Double_t  ip3d,
Double_t  eta,
Double_t  pt,
Bool_t  printDebug = kFALSE 
)

Definition at line 456 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, bindVariables(), gather_cfg::cout, eta, fisInitialized, fMethodname, fMVAType, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_ip3d, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, GetMVABin(), kTrig, and EnergyCorrector::pt.

481  {
482 
483  if (!fisInitialized) {
484  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
485  return -9999;
486  }
487 
488  if (fMVAType != EGammaMvaEleEstimator::kTrig) {
489  std::cout << "Error: This method should be called for kTrig MVA only" << endl;
490  return -9999;
491  }
492 
493  fMVAVar_fbrem = fbrem;
494  fMVAVar_kfchi2 = kfchi2;
495  fMVAVar_kfhits = float(kfhits); // BTD does not support int variables
496  fMVAVar_gsfchi2 = gsfchi2;
497 
498  fMVAVar_deta = deta;
499  fMVAVar_dphi = dphi;
500  fMVAVar_detacalo = detacalo;
501 
502 
503  fMVAVar_see = see;
504  fMVAVar_spp = spp;
505  fMVAVar_etawidth = etawidth;
506  fMVAVar_phiwidth = phiwidth;
507  fMVAVar_OneMinusE1x5E5x5 = OneMinusE1x5E5x5;
508  fMVAVar_R9 = R9;
509 
510 
511  fMVAVar_HoE = HoE;
512  fMVAVar_EoP = EoP;
513  fMVAVar_IoEmIoP = IoEmIoP;
514  fMVAVar_eleEoPout = eleEoPout;
515  fMVAVar_PreShowerOverRaw= PreShowerOverRaw;
516 
517  fMVAVar_d0 = d0;
518  fMVAVar_ip3d = ip3d;
519  fMVAVar_eta = eta;
520  fMVAVar_pt = pt;
521 
522 
523  bindVariables();
524  Double_t mva = -9999;
525  if (fUseBinnedVersion) {
526  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
527  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
528  } else {
529  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
530  }
531 
532  if(printDebug) {
533  cout << " *** Inside the class fMethodname " << fMethodname << endl;
534  cout << " fbrem " << fMVAVar_fbrem
535  << " kfchi2 " << fMVAVar_kfchi2
536  << " mykfhits " << fMVAVar_kfhits
537  << " gsfchi2 " << fMVAVar_gsfchi2
538  << " deta " << fMVAVar_deta
539  << " dphi " << fMVAVar_dphi
540  << " detacalo " << fMVAVar_detacalo
541  << " see " << fMVAVar_see
542  << " spp " << fMVAVar_spp
543  << " etawidth " << fMVAVar_etawidth
544  << " phiwidth " << fMVAVar_phiwidth
545  << " OneMinusE1x5E5x5 " << fMVAVar_OneMinusE1x5E5x5
546  << " R9 " << fMVAVar_R9
547  << " HoE " << fMVAVar_HoE
548  << " EoP " << fMVAVar_EoP
549  << " IoEmIoP " << fMVAVar_IoEmIoP
550  << " eleEoPout " << fMVAVar_eleEoPout
551  << " PreShowerOverRaw " << fMVAVar_PreShowerOverRaw
552  << " d0 " << fMVAVar_d0
553  << " ip3d " << fMVAVar_ip3d
554  << " eta " << fMVAVar_eta
555  << " pt " << fMVAVar_pt << endl;
556  cout << " ### MVA " << mva << endl;
557  }
558 
559 
560  return mva;
561 }
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
Double_t EGammaMvaEleEstimator::mvaValue ( Double_t  fbrem,
Double_t  kfchi2,
Int_t  kfhits,
Double_t  gsfchi2,
Double_t  deta,
Double_t  dphi,
Double_t  detacalo,
Double_t  see,
Double_t  spp,
Double_t  etawidth,
Double_t  phiwidth,
Double_t  e1x5e5x5,
Double_t  R9,
Double_t  HoE,
Double_t  EoP,
Double_t  IoEmIoP,
Double_t  eleEoPout,
Double_t  rho,
Double_t  PreShowerOverRaw,
Double_t  eta,
Double_t  pt,
Bool_t  printDebug = kFALSE 
)

Definition at line 565 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, bindVariables(), gather_cfg::cout, eta, fisInitialized, fMethodname, fMVAType, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_rho, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, GetMVABin(), kTrigNoIP, EnergyCorrector::pt, and rho.

586  {
587 
588  if (!fisInitialized) {
589  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
590  return -9999;
591  }
592 
593  if (fMVAType != EGammaMvaEleEstimator::kTrigNoIP) {
594  std::cout << "Error: This method should be called for kTrigNoIP MVA only" << endl;
595  return -9999;
596  }
597 
598  fMVAVar_fbrem = fbrem;
599  fMVAVar_kfchi2 = kfchi2;
600  fMVAVar_kfhits = float(kfhits); // BTD does not support int variables
601  fMVAVar_gsfchi2 = gsfchi2;
602 
603  fMVAVar_deta = deta;
604  fMVAVar_dphi = dphi;
605  fMVAVar_detacalo = detacalo;
606 
607  fMVAVar_see = see;
608  fMVAVar_spp = spp;
609  fMVAVar_etawidth = etawidth;
610  fMVAVar_phiwidth = phiwidth;
611  fMVAVar_OneMinusE1x5E5x5 = e1x5e5x5;
612  fMVAVar_R9 = R9;
613 
614  fMVAVar_HoE = HoE;
615  fMVAVar_EoP = EoP;
616  fMVAVar_IoEmIoP = IoEmIoP;
617  fMVAVar_eleEoPout = eleEoPout;
618  fMVAVar_rho = rho;
619  fMVAVar_PreShowerOverRaw= PreShowerOverRaw;
620 
621 
622  fMVAVar_eta = eta;
623  fMVAVar_pt = pt;
624 
625 
626  bindVariables();
627  Double_t mva = -9999;
628  if (fUseBinnedVersion) {
629  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
630  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
631  } else {
632  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
633  }
634 
635  if(printDebug) {
636  cout << " *** Inside the class fMethodname " << fMethodname << endl;
637  cout << " fbrem " << fMVAVar_fbrem
638  << " kfchi2 " << fMVAVar_kfchi2
639  << " mykfhits " << fMVAVar_kfhits
640  << " gsfchi2 " << fMVAVar_gsfchi2
641  << " deta " << fMVAVar_deta
642  << " dphi " << fMVAVar_dphi
643  << " detacalo " << fMVAVar_detacalo
644  << " see " << fMVAVar_see
645  << " spp " << fMVAVar_spp
646  << " etawidth " << fMVAVar_etawidth
647  << " phiwidth " << fMVAVar_phiwidth
648  << " e1x5e5x5 " << fMVAVar_OneMinusE1x5E5x5
649  << " R9 " << fMVAVar_R9
650  << " HoE " << fMVAVar_HoE
651  << " EoP " << fMVAVar_EoP
652  << " IoEmIoP " << fMVAVar_IoEmIoP
653  << " eleEoPout " << fMVAVar_eleEoPout
654  << " rho " << fMVAVar_rho
655  << " PreShowerOverRaw " << fMVAVar_PreShowerOverRaw
656  << " eta " << fMVAVar_eta
657  << " pt " << fMVAVar_pt << endl;
658  cout << " ### MVA " << mva << endl;
659  }
660 
661 
662  return mva;
663 }
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37
Double_t EGammaMvaEleEstimator::mvaValue ( Double_t  fbrem,
Double_t  kfchi2,
Int_t  kfhits,
Double_t  gsfchi2,
Double_t  deta,
Double_t  dphi,
Double_t  detacalo,
Double_t  see,
Double_t  spp,
Double_t  etawidth,
Double_t  phiwidth,
Double_t  e1x5e5x5,
Double_t  R9,
Double_t  HoE,
Double_t  EoP,
Double_t  IoEmIoP,
Double_t  eleEoPout,
Double_t  PreShowerOverRaw,
Double_t  eta,
Double_t  pt,
Bool_t  printDebug = kFALSE 
)

Definition at line 667 of file EGammaMvaEleEstimator.cc.

References newFWLiteAna::bin, bindVariables(), gather_cfg::cout, eta, fisInitialized, fMethodname, fMVAType, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_eleEoPout, fMVAVar_EoP, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_OneMinusE1x5E5x5, fMVAVar_phiwidth, fMVAVar_PreShowerOverRaw, fMVAVar_pt, fMVAVar_R9, fMVAVar_see, fMVAVar_spp, fTMVAMethod, fTMVAReader, fUseBinnedVersion, GetMVABin(), kNonTrig, and EnergyCorrector::pt.

690  {
691 
692  if (!fisInitialized) {
693  std::cout << "Error: EGammaMvaEleEstimator not properly initialized.\n";
694  return -9999;
695  }
696 
697  if (fMVAType != EGammaMvaEleEstimator::kNonTrig) {
698  std::cout << "Error: This method should be called for kNonTrig MVA only" << endl;
699  return -9999;
700  }
701 
702  fMVAVar_fbrem = fbrem;
703  fMVAVar_kfchi2 = kfchi2;
704  fMVAVar_kfhits = float(kfhits); // BTD does not support int variables
705  fMVAVar_gsfchi2 = gsfchi2;
706 
707  fMVAVar_deta = deta;
708  fMVAVar_dphi = dphi;
709  fMVAVar_detacalo = detacalo;
710 
711 
712  fMVAVar_see = see;
713  fMVAVar_spp = spp;
714  fMVAVar_etawidth = etawidth;
715  fMVAVar_phiwidth = phiwidth;
716  fMVAVar_OneMinusE1x5E5x5 = OneMinusE1x5E5x5;
717  fMVAVar_R9 = R9;
718 
719 
720  fMVAVar_HoE = HoE;
721  fMVAVar_EoP = EoP;
722  fMVAVar_IoEmIoP = IoEmIoP;
723  fMVAVar_eleEoPout = eleEoPout;
724  fMVAVar_PreShowerOverRaw= PreShowerOverRaw;
725 
726  fMVAVar_eta = eta;
727  fMVAVar_pt = pt;
728 
729 
730  bindVariables();
731  Double_t mva = -9999;
732  if (fUseBinnedVersion) {
733  int bin = GetMVABin(fMVAVar_eta,fMVAVar_pt);
734  mva = fTMVAReader[bin]->EvaluateMVA(fTMVAMethod[bin]);
735  } else {
736  mva = fTMVAReader[0]->EvaluateMVA(fTMVAMethod[0]);
737  }
738 
739 
740 
741  if(printDebug) {
742  cout << " *** Inside the class fMethodname " << fMethodname << endl;
743  cout << " fbrem " << fMVAVar_fbrem
744  << " kfchi2 " << fMVAVar_kfchi2
745  << " mykfhits " << fMVAVar_kfhits
746  << " gsfchi2 " << fMVAVar_gsfchi2
747  << " deta " << fMVAVar_deta
748  << " dphi " << fMVAVar_dphi
749  << " detacalo " << fMVAVar_detacalo
750  << " see " << fMVAVar_see
751  << " spp " << fMVAVar_spp
752  << " etawidth " << fMVAVar_etawidth
753  << " phiwidth " << fMVAVar_phiwidth
754  << " OneMinusE1x5E5x5 " << fMVAVar_OneMinusE1x5E5x5
755  << " R9 " << fMVAVar_R9
756  << " HoE " << fMVAVar_HoE
757  << " EoP " << fMVAVar_EoP
758  << " IoEmIoP " << fMVAVar_IoEmIoP
759  << " eleEoPout " << fMVAVar_eleEoPout
760  << " PreShowerOverRaw " << fMVAVar_PreShowerOverRaw
761  << " eta " << fMVAVar_eta
762  << " pt " << fMVAVar_pt << endl;
763  cout << " ### MVA " << mva << endl;
764  }
765 
766 
767  return mva;
768 }
EGammaMvaEleEstimator::GetMVABin
UInt_t GetMVABin(double eta, double pt) const
Definition: EGammaMvaEleEstimator.cc:411
EGammaMvaEleEstimator::fTMVAReader
std::vector< TMVA::Reader * > fTMVAReader
Definition: EGammaMvaEleEstimator.h:249
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EGammaMvaEleEstimator::fTMVAMethod
std::vector< TMVA::MethodBase * > fTMVAMethod
Definition: EGammaMvaEleEstimator.h:250
SQLiteCheck_cfg.printDebug
tuple printDebug
Definition: SQLiteCheck_cfg.py:37

Member Data Documentation

Bool_t EGammaMvaEleEstimator::fisInitialized
private

Definition at line 252 of file EGammaMvaEleEstimator.h.

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

std::string EGammaMvaEleEstimator::fMethodname
private

Definition at line 251 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), isoMvaValue(), and mvaValue().

MVAType EGammaMvaEleEstimator::fMVAType
private

Definition at line 253 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_ChargedIso_DR0p0To0p1
private

Definition at line 289 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_ChargedIso_DR0p1To0p2
private

Definition at line 290 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_ChargedIso_DR0p2To0p3
private

Definition at line 291 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_ChargedIso_DR0p3To0p4
private

Definition at line 292 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_ChargedIso_DR0p4To0p5
private

Definition at line 293 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_d0
private

Definition at line 281 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_deta
private

Definition at line 263 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_detacalo
private

Definition at line 265 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_dphi
private

Definition at line 264 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_eleEoPout
private

Definition at line 277 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_EoP
private

Definition at line 275 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_EoPout
private

Definition at line 278 of file EGammaMvaEleEstimator.h.

Float_t EGammaMvaEleEstimator::fMVAVar_eta
private

Definition at line 285 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), isoMvaValue(), and mvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_etawidth
private

Definition at line 269 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_fbrem
private

Definition at line 257 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_GammaIso_DR0p0To0p1
private

Definition at line 294 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_GammaIso_DR0p1To0p2
private

Definition at line 295 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_GammaIso_DR0p2To0p3
private

Definition at line 296 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_GammaIso_DR0p3To0p4
private

Definition at line 297 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_GammaIso_DR0p4To0p5
private

Definition at line 298 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_gsfchi2
private

Definition at line 261 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_HoE
private

Definition at line 274 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_IoEmIoP
private

Definition at line 276 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_ip3d
private

Definition at line 282 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_ip3dSig
private

Definition at line 283 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), and mvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_kfchi2
private

Definition at line 258 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_kfhits
private

Definition at line 259 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_kfhitsall
private

Definition at line 260 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), and mvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_NeutralHadronIso_DR0p0To0p1
private

Definition at line 299 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_NeutralHadronIso_DR0p1To0p2
private

Definition at line 300 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_NeutralHadronIso_DR0p2To0p3
private

Definition at line 301 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_NeutralHadronIso_DR0p3To0p4
private

Definition at line 302 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_NeutralHadronIso_DR0p4To0p5
private

Definition at line 303 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), and isoMvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_OneMinusE1x5E5x5
private

Definition at line 271 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_phiwidth
private

Definition at line 270 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_PreShowerOverRaw
private

Definition at line 279 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_pt
private

Definition at line 286 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), isoMvaValue(), and mvaValue().

Float_t EGammaMvaEleEstimator::fMVAVar_R9
private

Definition at line 272 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_rho
private

Definition at line 287 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_see
private

Definition at line 267 of file EGammaMvaEleEstimator.h.

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

Float_t EGammaMvaEleEstimator::fMVAVar_spp
private

Definition at line 268 of file EGammaMvaEleEstimator.h.

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

UInt_t EGammaMvaEleEstimator::fNMVABins
private

Definition at line 255 of file EGammaMvaEleEstimator.h.

Referenced by initialize().

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

Definition at line 250 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), isoMvaValue(), mvaValue(), and ~EGammaMvaEleEstimator().

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

Definition at line 249 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), isoMvaValue(), mvaValue(), and ~EGammaMvaEleEstimator().

Bool_t EGammaMvaEleEstimator::fUseBinnedVersion
private

Definition at line 254 of file EGammaMvaEleEstimator.h.

Referenced by IDIsoCombinedMvaValue(), initialize(), isoMvaValue(), and mvaValue().

Generated for CMSSW Reference Manual by doxygen 1.8.5