#include <AntiElectronIDMVA6.h>

Public Types
typedef std::vector< ElectronType >	ElectronCollection

typedef edm::Ref< ElectronCollection >	ElectronRef

Public Member Functions
	AntiElectronIDMVA6 (const edm::ParameterSet &, edm::ConsumesCollector &&)

void	beginEvent (const edm::Event &, const edm::EventSetup &)

antiElecIDMVA6_blocks::ElecVars	getElecVars (const ElectronRef &theEleRef)

void	getElecVarsHGCalTypeSpecific (const reco::GsfElectronRef &theEleRef, antiElecIDMVA6_blocks::ElecVars &elecVars)

void	getElecVarsHGCalTypeSpecific (const pat::ElectronRef &theEleRef, antiElecIDMVA6_blocks::ElecVars &elecVars)

antiElecIDMVA6_blocks::TauGammaVecs	getTauGammaVecs (const TauType &theTau)

antiElecIDMVA6_blocks::TauVars	getTauVars (const TauType &theTau)

antiElecIDMVA6_blocks::TauVars	getTauVarsTypeSpecific (const reco::PFTau &theTau)

antiElecIDMVA6_blocks::TauVars	getTauVarsTypeSpecific (const pat::Tau &theTau)

double	mvaValue (const antiElecIDMVA6_blocks::TauVars &tauVars, const antiElecIDMVA6_blocks::TauGammaVecs &tauGammaVecs, const antiElecIDMVA6_blocks::ElecVars &elecVars)

double	mvaValue (const antiElecIDMVA6_blocks::TauVars &tauVars, const antiElecIDMVA6_blocks::TauGammaMoms &tauGammaMoms, const antiElecIDMVA6_blocks::ElecVars &elecVars)

double	mvaValue (const TauType &theTau, const ElectronRef &theEleRef)

double	mvaValue (const TauType &theTau)

double	mvaValuePhase2 (const antiElecIDMVA6_blocks::TauVars &tauVars, const antiElecIDMVA6_blocks::TauGammaMoms &tauGammaMoms, const antiElecIDMVA6_blocks::ElecVars &elecVars)

	~AntiElectronIDMVA6 ()

Private Member Functions
double	dCrackEta (double eta)

double	dCrackPhi (double phi, double eta)

bool	energyWeightedEtaAndPhiAtECal (const pat::Tau &theTau, float &eta, float &phi)

double	minimum (double a, double b)

Private Attributes
std::map< std::string, edm::Handle< edm::ValueMap< float > > >	electronIds_

std::map< std::string, edm::EDGetTokenT< edm::ValueMap< float > > >	electronIds_tokens_

edm::FileInPath	inputFileName_

std::vector< TFile * >	inputFilesToDelete_

bool	isInitialized_

const bool	isPhase2_

bool	loadMVAfromDB_

const GBRForest *	mva_NoEleMatch_wGwoGSF_BL_

const GBRForest *	mva_NoEleMatch_wGwoGSF_EC_

const GBRForest *	mva_NoEleMatch_wGwoGSF_VFEC_

const GBRForest *	mva_NoEleMatch_woGwoGSF_BL_

const GBRForest *	mva_NoEleMatch_woGwoGSF_EC_

const GBRForest *	mva_NoEleMatch_woGwoGSF_VFEC_

const GBRForest *	mva_wGwGSF_BL_

const GBRForest *	mva_wGwGSF_EC_

const GBRForest *	mva_wGwGSF_VFEC_

const GBRForest *	mva_woGwGSF_BL_

const GBRForest *	mva_woGwGSF_EC_

const GBRForest *	mva_woGwGSF_VFEC_

std::string	mvaName_NoEleMatch_wGwoGSF_BL_

std::string	mvaName_NoEleMatch_wGwoGSF_EC_

std::string	mvaName_NoEleMatch_wGwoGSF_VFEC_

std::string	mvaName_NoEleMatch_woGwoGSF_BL_

std::string	mvaName_NoEleMatch_woGwoGSF_EC_

std::string	mvaName_NoEleMatch_woGwoGSF_VFEC_

std::string	mvaName_wGwGSF_BL_

std::string	mvaName_wGwGSF_EC_

std::string	mvaName_wGwGSF_VFEC_

std::string	mvaName_woGwGSF_BL_

std::string	mvaName_woGwGSF_EC_

std::string	mvaName_woGwGSF_VFEC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_NoEleMatch_wGwoGSF_BL_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_NoEleMatch_wGwoGSF_EC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_NoEleMatch_wGwoGSF_VFEC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_NoEleMatch_woGwoGSF_BL_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_NoEleMatch_woGwoGSF_EC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_NoEleMatch_woGwoGSF_VFEC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_wGwGSF_BL_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_wGwGSF_EC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_wGwGSF_VFEC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_woGwGSF_BL_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_woGwGSF_EC_

edm::ESGetToken< GBRForest, GBRWrapperRcd >	mvaToken_woGwGSF_VFEC_

PositionAtECalEntranceComputer	positionAtECalEntrance_

bool	usePhiAtEcalEntranceExtrapolation_

std::vector< float >	var_NoEleMatch_wGwoGSF_Barrel_

std::vector< float >	var_NoEleMatch_wGwoGSF_Endcap_

std::vector< float >	var_NoEleMatch_wGwoGSF_VFEndcap_

std::vector< float >	var_NoEleMatch_woGwoGSF_Barrel_

std::vector< float >	var_NoEleMatch_woGwoGSF_Endcap_

std::vector< float >	var_NoEleMatch_woGwoGSF_VFEndcap_

std::vector< float >	var_wGwGSF_Barrel_

std::vector< float >	var_wGwGSF_Endcap_

std::vector< float >	var_wGwGSF_VFEndcap_

std::vector< float >	var_woGwGSF_Barrel_

std::vector< float >	var_woGwGSF_Endcap_

std::vector< float >	var_woGwGSF_VFEndcap_

const int	verbosity_

Static Private Attributes
static constexpr float	ecalBarrelEndcapEtaBorder_ = 1.479

static constexpr float	ecalEndcapVFEndcapEtaBorder_ = 2.4

Detailed Description

template<class TauType, class ElectronType>
class AntiElectronIDMVA6< TauType, ElectronType >

Definition at line 116 of file AntiElectronIDMVA6.h.

Member Typedef Documentation

◆ ElectronCollection

template<class TauType , class ElectronType >

typedef std::vector<ElectronType> AntiElectronIDMVA6< TauType, ElectronType >::ElectronCollection

Definition at line 118 of file AntiElectronIDMVA6.h.

◆ ElectronRef

template<class TauType , class ElectronType >

typedef edm::Ref<ElectronCollection> AntiElectronIDMVA6< TauType, ElectronType >::ElectronRef

Definition at line 119 of file AntiElectronIDMVA6.h.

Constructor & Destructor Documentation

◆ AntiElectronIDMVA6()

template<class TauType , class ElectronType >

AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	cc
	)

Definition at line 36 of file AntiElectronIDMVA6.cc.

     : isInitialized_(false),
       mva_NoEleMatch_woGwoGSF_BL_(nullptr),
       mva_NoEleMatch_wGwoGSF_BL_(nullptr),
       mva_woGwGSF_BL_(nullptr),
       mva_wGwGSF_BL_(nullptr),
       mva_NoEleMatch_woGwoGSF_EC_(nullptr),
       mva_NoEleMatch_wGwoGSF_EC_(nullptr),
       mva_woGwGSF_EC_(nullptr),
       mva_wGwGSF_EC_(nullptr),
       mva_NoEleMatch_woGwoGSF_VFEC_(nullptr),
       mva_NoEleMatch_wGwoGSF_VFEC_(nullptr),
       mva_woGwGSF_VFEC_(nullptr),
       mva_wGwGSF_VFEC_(nullptr),
       isPhase2_(cfg.getParameter<bool>("isPhase2")),
       positionAtECalEntrance_(PositionAtECalEntranceComputer(cc, isPhase2_)),
       verbosity_(cfg.getParameter<int>("verbosity")) {
   loadMVAfromDB_ = cfg.exists("loadMVAfromDB") ? cfg.getParameter<bool>("loadMVAfromDB") : false;
   if (!loadMVAfromDB_) {
     if (cfg.exists("inputFileName")) {
       inputFileName_ = cfg.getParameter<edm::FileInPath>("inputFileName");
     } else
       throw cms::Exception("MVA input not defined")
           << "Requested to load tau MVA input from ROOT file but no file provided in cfg file";
   }
 
   mvaName_NoEleMatch_woGwoGSF_BL_ = cfg.getParameter<std::string>("mvaName_NoEleMatch_woGwoGSF_BL");
   mvaName_NoEleMatch_wGwoGSF_BL_ = cfg.getParameter<std::string>("mvaName_NoEleMatch_wGwoGSF_BL");
   mvaName_woGwGSF_BL_ = cfg.getParameter<std::string>("mvaName_woGwGSF_BL");
   mvaName_wGwGSF_BL_ = cfg.getParameter<std::string>("mvaName_wGwGSF_BL");
   mvaName_NoEleMatch_woGwoGSF_EC_ = cfg.getParameter<std::string>("mvaName_NoEleMatch_woGwoGSF_EC");
   mvaName_NoEleMatch_wGwoGSF_EC_ = cfg.getParameter<std::string>("mvaName_NoEleMatch_wGwoGSF_EC");
   mvaName_woGwGSF_EC_ = cfg.getParameter<std::string>("mvaName_woGwGSF_EC");
   mvaName_wGwGSF_EC_ = cfg.getParameter<std::string>("mvaName_wGwGSF_EC");
   if (isPhase2_) {
     mvaName_NoEleMatch_woGwoGSF_VFEC_ = cfg.getParameter<std::string>("mvaName_NoEleMatch_woGwoGSF_VFEC");
     mvaName_NoEleMatch_wGwoGSF_VFEC_ = cfg.getParameter<std::string>("mvaName_NoEleMatch_wGwoGSF_VFEC");
     mvaName_woGwGSF_VFEC_ = cfg.getParameter<std::string>("mvaName_woGwGSF_VFEC");
     mvaName_wGwGSF_VFEC_ = cfg.getParameter<std::string>("mvaName_wGwGSF_VFEC");
   }
   if (loadMVAfromDB_) {
     mvaToken_NoEleMatch_woGwoGSF_BL_ = cc.esConsumes(edm::ESInputTag{"", mvaName_NoEleMatch_woGwoGSF_BL_});
     mvaToken_NoEleMatch_wGwoGSF_BL_ = cc.esConsumes(edm::ESInputTag{"", mvaName_NoEleMatch_wGwoGSF_BL_});
     mvaToken_woGwGSF_BL_ = cc.esConsumes(edm::ESInputTag{"", mvaName_woGwGSF_BL_});
     mvaToken_wGwGSF_BL_ = cc.esConsumes(edm::ESInputTag{"", mvaName_wGwGSF_BL_});
     mvaToken_NoEleMatch_woGwoGSF_EC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_NoEleMatch_woGwoGSF_EC_});
     mvaToken_NoEleMatch_wGwoGSF_EC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_NoEleMatch_wGwoGSF_EC_});
     mvaToken_woGwGSF_EC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_woGwGSF_EC_});
     mvaToken_wGwGSF_EC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_wGwGSF_EC_});
     if (isPhase2_) {
       mvaToken_NoEleMatch_woGwoGSF_VFEC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_NoEleMatch_woGwoGSF_VFEC_});
       mvaToken_NoEleMatch_wGwoGSF_VFEC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_NoEleMatch_wGwoGSF_VFEC_});
       mvaToken_woGwGSF_VFEC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_woGwGSF_VFEC_});
       mvaToken_wGwGSF_VFEC_ = cc.esConsumes(edm::ESInputTag{"", mvaName_wGwGSF_VFEC_});
     }
   }
   usePhiAtEcalEntranceExtrapolation_ = cfg.getParameter<bool>("usePhiAtEcalEntranceExtrapolation");
 
   if (!isPhase2_) {
     var_NoEleMatch_woGwoGSF_Barrel_ = std::vector<float>(10);
     var_NoEleMatch_wGwoGSF_Barrel_ = std::vector<float>(18);
     var_woGwGSF_Barrel_ = std::vector<float>(24);
     var_wGwGSF_Barrel_ = std::vector<float>(32);
     var_NoEleMatch_woGwoGSF_Endcap_ = std::vector<float>(9);
     var_NoEleMatch_wGwoGSF_Endcap_ = std::vector<float>(17);
     var_woGwGSF_Endcap_ = std::vector<float>(23);
     var_wGwGSF_Endcap_ = std::vector<float>(31);
   } else {
     var_NoEleMatch_woGwoGSF_Barrel_ = std::vector<float>(9);
     var_NoEleMatch_wGwoGSF_Barrel_ = std::vector<float>(17);
     var_woGwGSF_Barrel_ = std::vector<float>(27);
     var_wGwGSF_Barrel_ = std::vector<float>(36);
     var_NoEleMatch_woGwoGSF_Endcap_ = std::vector<float>(6);
     var_NoEleMatch_wGwoGSF_Endcap_ = std::vector<float>(14);
     var_woGwGSF_Endcap_ = std::vector<float>(31);
     var_wGwGSF_Endcap_ = std::vector<float>(38);
     var_NoEleMatch_woGwoGSF_VFEndcap_ = std::vector<float>(6);
     var_NoEleMatch_wGwoGSF_VFEndcap_ = std::vector<float>(14);
     var_woGwGSF_VFEndcap_ = std::vector<float>(32);
     var_wGwGSF_VFEndcap_ = std::vector<float>(40);
 
     //MB: Tokens for additional inputs (HGCal EleID variables) only for Phase2 and reco::GsfElectrons
     if (std::is_same<ElectronType, reco::GsfElectron>::value) {
       for (const auto& inputTag : cfg.getParameter<std::vector<edm::InputTag>>("hgcalElectronIDs")) {
         std::string elecIdLabel = "hgcElectronID:" + inputTag.instance();
         electronIds_tokens_[elecIdLabel] = cc.consumes<edm::ValueMap<float>>(
             inputTag);  //MB: It assumes that instances are not empty and meaningful (as for userData in patElectrons)
         electronIds_[elecIdLabel] = edm::Handle<edm::ValueMap<float>>();
       }
     }
   }
 }

◆ ~AntiElectronIDMVA6()

template<class TauType , class ElectronType >

AntiElectronIDMVA6< TauType, ElectronType >::~AntiElectronIDMVA6 ( )

Definition at line 130 of file AntiElectronIDMVA6.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

                                                                {
   if (!loadMVAfromDB_) {
     delete mva_NoEleMatch_woGwoGSF_BL_;
     delete mva_NoEleMatch_wGwoGSF_BL_;
     delete mva_woGwGSF_BL_;
     delete mva_wGwGSF_BL_;
     delete mva_NoEleMatch_woGwoGSF_EC_;
     delete mva_NoEleMatch_wGwoGSF_EC_;
     delete mva_woGwGSF_EC_;
     delete mva_wGwGSF_EC_;
     if (isPhase2_) {
       delete mva_NoEleMatch_woGwoGSF_VFEC_;
       delete mva_NoEleMatch_wGwoGSF_VFEC_;
       delete mva_woGwGSF_VFEC_;
       delete mva_wGwGSF_VFEC_;
     }
   }
 
   for (std::vector<TFile*>::iterator it = inputFilesToDelete_.begin(); it != inputFilesToDelete_.end(); ++it) {
     delete (*it);
   }
 }

Member Function Documentation

◆ beginEvent()

template<class TauType , class ElectronType >

void AntiElectronIDMVA6< TauType, ElectronType >::beginEvent	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

Definition at line 166 of file AntiElectronIDMVA6.cc.

References edm::EventSetup::getData(), edm::Event::getHandle(), makeListRunsInFiles::inputFile, edm::FileInPath::Unknown, and relativeConstraints::value.

                                                                                                        {
   if (!isInitialized_) {
     if (loadMVAfromDB_) {
       mva_NoEleMatch_woGwoGSF_BL_ = &es.getData(mvaToken_NoEleMatch_woGwoGSF_BL_);
       mva_NoEleMatch_wGwoGSF_BL_ = &es.getData(mvaToken_NoEleMatch_wGwoGSF_BL_);
       mva_woGwGSF_BL_ = &es.getData(mvaToken_woGwGSF_BL_);
       mva_wGwGSF_BL_ = &es.getData(mvaToken_wGwGSF_BL_);
       mva_NoEleMatch_woGwoGSF_EC_ = &es.getData(mvaToken_NoEleMatch_woGwoGSF_EC_);
       mva_NoEleMatch_wGwoGSF_EC_ = &es.getData(mvaToken_NoEleMatch_wGwoGSF_EC_);
       mva_woGwGSF_EC_ = &es.getData(mvaToken_woGwGSF_EC_);
       mva_wGwGSF_EC_ = &es.getData(mvaToken_wGwGSF_EC_);
       if (isPhase2_) {
         mva_NoEleMatch_woGwoGSF_VFEC_ = &es.getData(mvaToken_NoEleMatch_woGwoGSF_VFEC_);
         mva_NoEleMatch_wGwoGSF_VFEC_ = &es.getData(mvaToken_NoEleMatch_wGwoGSF_VFEC_);
         mva_woGwGSF_VFEC_ = &es.getData(mvaToken_woGwGSF_VFEC_);
         mva_wGwGSF_VFEC_ = &es.getData(mvaToken_wGwGSF_VFEC_);
       }
     } else {
       if (inputFileName_.location() == edm::FileInPath::Unknown)
         throw cms::Exception("PFRecoTauDiscriminationAgainstElectronMVA6::loadMVA")
             << " Failed to find File = " << inputFileName_ << " !!\n";
       TFile* inputFile = new TFile(inputFileName_.fullPath().data());
 
       mva_NoEleMatch_woGwoGSF_BL_ = loadMVAfromFile(inputFile, mvaName_NoEleMatch_woGwoGSF_BL_);
       mva_NoEleMatch_wGwoGSF_BL_ = loadMVAfromFile(inputFile, mvaName_NoEleMatch_wGwoGSF_BL_);
       mva_woGwGSF_BL_ = loadMVAfromFile(inputFile, mvaName_woGwGSF_BL_);
       mva_wGwGSF_BL_ = loadMVAfromFile(inputFile, mvaName_wGwGSF_BL_);
       mva_NoEleMatch_woGwoGSF_EC_ = loadMVAfromFile(inputFile, mvaName_NoEleMatch_woGwoGSF_EC_);
       mva_NoEleMatch_wGwoGSF_EC_ = loadMVAfromFile(inputFile, mvaName_NoEleMatch_wGwoGSF_EC_);
       mva_woGwGSF_EC_ = loadMVAfromFile(inputFile, mvaName_woGwGSF_EC_);
       mva_wGwGSF_EC_ = loadMVAfromFile(inputFile, mvaName_wGwGSF_EC_);
       if (isPhase2_) {
         mva_NoEleMatch_woGwoGSF_VFEC_ = loadMVAfromFile(inputFile, mvaName_NoEleMatch_woGwoGSF_VFEC_);
         mva_NoEleMatch_wGwoGSF_VFEC_ = loadMVAfromFile(inputFile, mvaName_NoEleMatch_wGwoGSF_VFEC_);
         mva_woGwGSF_VFEC_ = loadMVAfromFile(inputFile, mvaName_woGwGSF_VFEC_);
         mva_wGwGSF_VFEC_ = loadMVAfromFile(inputFile, mvaName_wGwGSF_VFEC_);
       }
       inputFilesToDelete_.push_back(inputFile);
     }
     isInitialized_ = true;
   }
   positionAtECalEntrance_.beginEvent(es);
   //MB: Handle additional inputs (HGCal EleID variables) only for Phase2 and reco::GsfElectrons
   if (isPhase2_ && std::is_same<ElectronType, reco::GsfElectron>::value) {
     for (const auto& eId_token : electronIds_tokens_) {
       electronIds_[eId_token.first] = evt.getHandle(eId_token.second);
     }
   }
 }

◆ dCrackEta()

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::dCrackEta ( double eta )

private

Definition at line 1082 of file AntiElectronIDMVA6.cc.

References funct::abs(), ztail::d, and PVValHelper::eta.

                                                                       {
   //--- compute the (unsigned) distance to the closest eta-crack in the ECAL barrel
 
   // IN: define locations of the eta-cracks
   double cracks[5] = {0., 4.44747e-01, 7.92824e-01, 1.14090e+00, 1.47464e+00};
 
   double retVal = 99.;
 
   for (int iCrack = 0; iCrack < 5; ++iCrack) {
     double d = minimum(eta - cracks[iCrack], eta + cracks[iCrack]);
     if (std::abs(d) < std::abs(retVal)) {
       retVal = d;
     }
   }
 
   return std::abs(retVal);
 }

◆ dCrackPhi()

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::dCrackPhi	(	double	phi,
		double	eta
	)

private

Definition at line 1031 of file AntiElectronIDMVA6.cc.

References funct::abs(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), PVValHelper::eta, M_PI, PVValHelper::phi, and pi.

                                                                                   {
   //--- compute the (unsigned) distance to the closest phi-crack in the ECAL barrel
 
   constexpr double pi = M_PI;  // 3.14159265358979323846;
 
   // IN: shift of this location if eta < 0
   constexpr double delta_cPhi = 0.00638;
 
   double retVal = 99.;
 
   if (eta >= -1.47464 && eta <= 1.47464) {
     // the location is shifted
     if (eta < 0.)
       phi += delta_cPhi;
 
     // CV: need to bring-back phi into interval [-pi,+pi]
     if (phi > pi)
       phi -= 2. * pi;
     if (phi < -pi)
       phi += 2. * pi;
 
     if (phi >= -pi && phi <= pi) {
       // the problem of the extrema:
       if (phi < cPhi[17] || phi >= cPhi[0]) {
         if (phi < 0.)
           phi += 2. * pi;
         retVal = minimum(phi - cPhi[0], phi - cPhi[17] - 2. * pi);
       } else {
         // between these extrema...
         bool OK = false;
         unsigned iCrack = 16;
         while (!OK) {
           if (phi < cPhi[iCrack]) {
             retVal = minimum(phi - cPhi[iCrack + 1], phi - cPhi[iCrack]);
             OK = true;
           } else {
             iCrack -= 1;
           }
         }
       }
     } else {
       retVal = 0.;  // IN: if there is a problem, we assume that we are in a crack
     }
   } else {
     return -99.;
   }
 
   return std::abs(retVal);
 }

◆ energyWeightedEtaAndPhiAtECal()

template<class TauType , class ElectronType >

bool AntiElectronIDMVA6< TauType, ElectronType >::energyWeightedEtaAndPhiAtECal	(	const pat::Tau &	theTau,
		float &	eta,
		float &	phi
	)

private

Definition at line 1329 of file AntiElectronIDMVA6.cc.

References PVValHelper::eta, PVValHelper::phi, pat::Tau::signalCands(), and summarizeEdmComparisonLogfiles::success.

                                                                                           {
   eta = 0.;
   phi = 0.;
   float sumEnergy = 0.;
   for (const auto& candidate : theTau.signalCands()) {
     float etaAtECalEntrance = candidate->eta();
     float phiAtECalEntrance = candidate->phi();
     bool success = false;
     reco::Candidate::Point posAtECal = positionAtECalEntrance_(candidate.get(), success);
     if (success) {
       etaAtECalEntrance = posAtECal.eta();
       phiAtECalEntrance = posAtECal.phi();
     }
     eta += etaAtECalEntrance * candidate->energy();
     phi += phiAtECalEntrance * candidate->energy();
     sumEnergy += candidate->energy();
   }
   if (sumEnergy > 0.) {
     eta = eta / sumEnergy;
     phi = phi / sumEnergy;
   } else {
     eta = -99.;
     phi = -99.;
   }
   return (sumEnergy > 0.);
 }

◆ getElecVars()

template<class TauType , class ElectronType >

ElecVars AntiElectronIDMVA6< TauType, ElectronType >::getElecVars ( const ElectronRef & theEleRef )

Definition at line 924 of file AntiElectronIDMVA6.cc.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::mvaValue().

                                                                                             {
   ElecVars elecVars;
 
   elecVars.eta = theEleRef->eta();
   elecVars.phi = theEleRef->phi();
 
   // Variables related to the electron Cluster
   float elecEe = 0.;
   float elecEgamma = 0.;
   reco::SuperClusterRef pfSuperCluster = theEleRef->superCluster();
   if (pfSuperCluster.isNonnull() && pfSuperCluster.isAvailable()) {
     if (!isPhase2_) {
       for (reco::CaloCluster_iterator pfCluster = pfSuperCluster->clustersBegin();
            pfCluster != pfSuperCluster->clustersEnd();
            ++pfCluster) {
         double pfClusterEn = (*pfCluster)->energy();
         if (pfCluster == pfSuperCluster->clustersBegin())
           elecEe += pfClusterEn;
         else
           elecEgamma += pfClusterEn;
       }
     }
     elecVars.superClusterEtaWidth = pfSuperCluster->etaWidth();
     elecVars.superClusterPhiWidth = pfSuperCluster->phiWidth();
   }
   elecVars.eSeedClusterOverPout = theEleRef->eSeedClusterOverPout();
   elecVars.showerCircularity = 1. - theEleRef->e1x5() / theEleRef->e5x5();
   elecVars.r9 = theEleRef->r9();
   elecVars.sigmaIEtaIEta5x5 = theEleRef->full5x5_sigmaIetaIeta();
   elecVars.sigmaIPhiIPhi5x5 = theEleRef->full5x5_sigmaIphiIphi();
 
   elecVars.pIn = std::sqrt(theEleRef->trackMomentumAtVtx().Mag2());
   elecVars.pOut = std::sqrt(theEleRef->trackMomentumOut().Mag2());
   elecVars.eTotOverPin = (elecVars.pIn > 0.0) ? ((elecEe + elecEgamma) / elecVars.pIn) : -0.1;
   elecVars.eEcal = theEleRef->ecalEnergy();
   if (!isPhase2_) {
     elecVars.deltaEta = theEleRef->deltaEtaSeedClusterTrackAtCalo();
     elecVars.deltaPhi = theEleRef->deltaPhiSeedClusterTrackAtCalo();
   } else {
     elecVars.deltaEta = theEleRef->deltaEtaEleClusterTrackAtCalo();
     elecVars.deltaPhi = theEleRef->deltaPhiEleClusterTrackAtCalo();
   }
   elecVars.mvaInSigmaEtaEta = theEleRef->mvaInput().sigmaEtaEta;
   elecVars.mvaInHadEnergy = theEleRef->mvaInput().hadEnergy;
   elecVars.mvaInDeltaEta = theEleRef->mvaInput().deltaEta;
 
   // Variables related to the GsfTrack
   elecVars.chi2NormGSF = -99.;
   elecVars.gsfNumHits = -99.;
   elecVars.gsfTrackResol = -99.;
   elecVars.gsfTracklnPt = -99.;
   if (theEleRef->gsfTrack().isNonnull()) {
     elecVars.chi2NormGSF = theEleRef->gsfTrack()->normalizedChi2();
     elecVars.gsfNumHits = theEleRef->gsfTrack()->numberOfValidHits();
     if (theEleRef->gsfTrack()->pt() > 0.) {
       elecVars.gsfTrackResol = theEleRef->gsfTrack()->ptError() / theEleRef->gsfTrack()->pt();
       elecVars.gsfTracklnPt = log(theEleRef->gsfTrack()->pt()) * M_LN10;
     }
   }
 
   // Variables related to the CtfTrack
   elecVars.chi2NormKF = -99.;
   elecVars.kfNumHits = -99.;
   if (theEleRef->closestCtfTrackRef().isNonnull()) {
     elecVars.chi2NormKF = theEleRef->closestCtfTrackRef()->normalizedChi2();
     elecVars.kfNumHits = theEleRef->closestCtfTrackRef()->numberOfValidHits();
   }
 
   // Variables related to HGCal
   if (isPhase2_ && !theEleRef->isEB()) {
     if (std::is_same<ElectronType, reco::GsfElectron>::value || std::is_same<ElectronType, pat::Electron>::value)
       getElecVarsHGCalTypeSpecific(theEleRef, elecVars);
     else
       throw cms::Exception("AntiElectronIDMVA6")
           << "Unsupported ElectronType used. You must use either reco::GsfElectron or pat::Electron.";
   }
 
   return elecVars;
 }

◆ getElecVarsHGCalTypeSpecific() [1/2]

template<class TauType , class ElectronType >

void AntiElectronIDMVA6< TauType, ElectronType >::getElecVarsHGCalTypeSpecific	(	const reco::GsfElectronRef &	theEleRef,
		antiElecIDMVA6_blocks::ElecVars &	elecVars
	)

Definition at line 1279 of file AntiElectronIDMVA6.cc.

                                                                                   {
   //MB: Assumed that presence of one of the HGCal EleID variables guarantee presence of all
   if (!(electronIds_.find("hgcElectronID:sigmaUU") != electronIds_.end() &&
         electronIds_.at("hgcElectronID:sigmaUU").isValid()))
     return;
 
   elecVars.hgcalSigmaUU = (*electronIds_.at("hgcElectronID:sigmaUU"))[theEleRef];
   elecVars.hgcalSigmaVV = (*electronIds_.at("hgcElectronID:sigmaVV"))[theEleRef];
   elecVars.hgcalSigmaEE = (*electronIds_.at("hgcElectronID:sigmaEE"))[theEleRef];
   elecVars.hgcalSigmaPP = (*electronIds_.at("hgcElectronID:sigmaPP"))[theEleRef];
   elecVars.hgcalNLayers = (*electronIds_.at("hgcElectronID:nLayers"))[theEleRef];
   elecVars.hgcalFirstLayer = (*electronIds_.at("hgcElectronID:firstLayer"))[theEleRef];
   elecVars.hgcalLastLayer = (*electronIds_.at("hgcElectronID:lastLayer"))[theEleRef];
   elecVars.hgcalLayerEfrac10 = (*electronIds_.at("hgcElectronID:layerEfrac10"))[theEleRef];
   elecVars.hgcalLayerEfrac90 = (*electronIds_.at("hgcElectronID:layerEfrac90"))[theEleRef];
   elecVars.hgcalEcEnergyEE = (*electronIds_.at("hgcElectronID:ecEnergyEE"))[theEleRef];
   elecVars.hgcalEcEnergyFH = (*electronIds_.at("hgcElectronID:ecEnergyFH"))[theEleRef];
   elecVars.hgcalMeasuredDepth = (*electronIds_.at("hgcElectronID:measuredDepth"))[theEleRef];
   elecVars.hgcalExpectedDepth = (*electronIds_.at("hgcElectronID:expectedDepth"))[theEleRef];
   elecVars.hgcalExpectedSigma = (*electronIds_.at("hgcElectronID:expectedSigma"))[theEleRef];
   elecVars.hgcalDepthCompatibility = (*electronIds_.at("hgcElectronID:depthCompatibility"))[theEleRef];
 }

◆ getElecVarsHGCalTypeSpecific() [2/2]

template<class TauType , class ElectronType >

void AntiElectronIDMVA6< TauType, ElectronType >::getElecVarsHGCalTypeSpecific	(	const pat::ElectronRef &	theEleRef,
		antiElecIDMVA6_blocks::ElecVars &	elecVars
	)

Definition at line 1305 of file AntiElectronIDMVA6.cc.

                                                                               {
   //MB: Assumed that presence of one of the HGCal EleID variables guarantee presence of all
   if (!theEleRef->hasUserFloat("hgcElectronID:sigmaUU"))
     return;
 
   elecVars.hgcalSigmaUU = theEleRef->userFloat("hgcElectronID:sigmaUU");
   elecVars.hgcalSigmaVV = theEleRef->userFloat("hgcElectronID:sigmaVV");
   elecVars.hgcalSigmaEE = theEleRef->userFloat("hgcElectronID:sigmaEE");
   elecVars.hgcalSigmaPP = theEleRef->userFloat("hgcElectronID:sigmaPP");
   elecVars.hgcalNLayers = theEleRef->userFloat("hgcElectronID:nLayers");
   elecVars.hgcalFirstLayer = theEleRef->userFloat("hgcElectronID:firstLayer");
   elecVars.hgcalLastLayer = theEleRef->userFloat("hgcElectronID:lastLayer");
   elecVars.hgcalLayerEfrac10 = theEleRef->userFloat("hgcElectronID:layerEfrac10");
   elecVars.hgcalLayerEfrac90 = theEleRef->userFloat("hgcElectronID:layerEfrac90");
   elecVars.hgcalEcEnergyEE = theEleRef->userFloat("hgcElectronID:ecEnergyEE");
   elecVars.hgcalEcEnergyFH = theEleRef->userFloat("hgcElectronID:ecEnergyFH");
   elecVars.hgcalMeasuredDepth = theEleRef->userFloat("hgcElectronID:measuredDepth");
   elecVars.hgcalExpectedDepth = theEleRef->userFloat("hgcElectronID:expectedDepth");
   elecVars.hgcalExpectedSigma = theEleRef->userFloat("hgcElectronID:expectedSigma");
   elecVars.hgcalDepthCompatibility = theEleRef->userFloat("hgcElectronID:depthCompatibility");
 }

◆ getTauGammaVecs()

template<class TauType , class ElectronType >

TauGammaVecs AntiElectronIDMVA6< TauType, ElectronType >::getTauGammaVecs ( const TauType & theTau )

Definition at line 901 of file AntiElectronIDMVA6.cc.

References PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, CustomPhysics_cfi::gamma, antiElecIDMVA6_blocks::TauGammaVecs::gammasdEtaInSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasdEtaOutSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasdPhiInSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasdPhiOutSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasPtInSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasPtOutSigCone, and SiStripPI::max.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::mvaValue().

                                                                                              {
   TauGammaVecs tauGammaVecs;
 
   float signalrad = std::clamp(3.0 / std::max(1.0, theTau.pt()), 0.05, 0.10);
   for (const auto& gamma : theTau.signalGammaCands()) {
     float dR = deltaR(gamma->p4(), theTau.leadChargedHadrCand()->p4());
     // pfGammas inside the tau signal cone
     if (dR < signalrad) {
       tauGammaVecs.gammasdEtaInSigCone.push_back(gamma->eta() - theTau.leadChargedHadrCand()->eta());
       tauGammaVecs.gammasdPhiInSigCone.push_back(gamma->phi() - theTau.leadChargedHadrCand()->phi());
       tauGammaVecs.gammasPtInSigCone.push_back(gamma->pt());
     }
     // pfGammas outside the tau signal cone
     else {
       tauGammaVecs.gammasdEtaOutSigCone.push_back(gamma->eta() - theTau.leadChargedHadrCand()->eta());
       tauGammaVecs.gammasdPhiOutSigCone.push_back(gamma->phi() - theTau.leadChargedHadrCand()->phi());
       tauGammaVecs.gammasPtOutSigCone.push_back(gamma->pt());
     }
   }
   return tauGammaVecs;
 }

◆ getTauVars()

template<class TauType , class ElectronType >

TauVars AntiElectronIDMVA6< TauType, ElectronType >::getTauVars ( const TauType & theTau )

Definition at line 851 of file AntiElectronIDMVA6.cc.

References funct::abs(), PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, Exception, CustomPhysics_cfi::gamma, reco::PFCandidate::gsfTrackRef(), edm::Ref< C, T, F >::isNonnull(), EgHLTOffHistBins_cfi::mass, SiStripPI::max, pat::PackedCandidate::pdgId(), leptonTimeLifeInfo_common_cff::tauVars, and relativeConstraints::value.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::mvaValue().

                                                                                    {
   TauVars tauVars;
   if (std::is_same<TauType, reco::PFTau>::value || std::is_same<TauType, pat::Tau>::value)
     tauVars = getTauVarsTypeSpecific(theTau);
   else
     throw cms::Exception("AntiElectronIDMVA6")
         << "Unsupported TauType used. You must use either reco::PFTau or pat::Tau.";
 
   tauVars.pt = theTau.pt();
 
   reco::Candidate::LorentzVector pfGammaSum(0, 0, 0, 0);
   reco::Candidate::LorentzVector pfChargedSum(0, 0, 0, 0);
   float signalrad = std::clamp(3.0 / std::max(1.0, theTau.pt()), 0.05, 0.10);
   for (const auto& gamma : theTau.signalGammaCands()) {
     float dR = deltaR(gamma->p4(), theTau.leadChargedHadrCand()->p4());
     // pfGammas inside the tau signal cone
     if (dR < signalrad) {
       pfGammaSum += gamma->p4();
     }
   }
   for (const auto& charged : theTau.signalChargedHadrCands()) {
     float dR = deltaR(charged->p4(), theTau.leadChargedHadrCand()->p4());
     // charged particles inside the tau signal cone
     if (dR < signalrad) {
       pfChargedSum += charged->p4();
     }
   }
   tauVars.visMassIn = (pfGammaSum + pfChargedSum).mass();
 
   tauVars.hasGsf = 0;
   if (theTau.leadChargedHadrCand().isNonnull()) {
     const pat::PackedCandidate* packedLeadChCand =
         dynamic_cast<const pat::PackedCandidate*>(theTau.leadChargedHadrCand().get());
     if (packedLeadChCand != nullptr) {
       if (std::abs(packedLeadChCand->pdgId()) == 11)
         tauVars.hasGsf = 1;
     } else {
       const reco::PFCandidate* pfLeadChCand =
           dynamic_cast<const reco::PFCandidate*>(theTau.leadChargedHadrCand().get());
       if (pfLeadChCand != nullptr && pfLeadChCand->gsfTrackRef().isNonnull())
         tauVars.hasGsf = 1;
     }
   }
   tauVars.dCrackPhi = dCrackPhi(tauVars.phi, tauVars.etaAtEcalEntrance);
   tauVars.dCrackEta = dCrackEta(tauVars.etaAtEcalEntrance);
 
   return tauVars;
 }

◆ getTauVarsTypeSpecific() [1/2]

template<class TauType , class ElectronType >

TauVars AntiElectronIDMVA6< TauType, ElectronType >::getTauVarsTypeSpecific ( const reco::PFTau & theTau )

Definition at line 1152 of file AntiElectronIDMVA6.cc.

References funct::abs(), pat::PackedCandidate::caloFraction(), reco::PFCandidate::ecalEnergy(), reco::PFTau::emFraction(), pat::PackedCandidate::energy(), reco::LeafCandidate::eta(), edm::Ptr< T >::get(), reco::PFCandidate::hcalEnergy(), pat::PackedCandidate::hcalFraction(), edm::Ptr< T >::isNonnull(), reco::PFTau::leadCand(), reco::PFTau::leadChargedHadrCand(), SiStripPI::max, reco::LeafCandidate::p(), pat::PackedCandidate::p(), reco::LeafCandidate::phi(), reco::PFCandidate::positionAtECALEntrance(), reco::PFTau::signalCands(), summarizeEdmComparisonLogfiles::success, leptonTimeLifeInfo_common_cff::tauVars, and HLT_2024v14_cff::track.

                                                                                                  {
   TauVars tauVars;
   tauVars.etaAtEcalEntrance = -99.;
   tauVars.leadChargedPFCandEtaAtEcalEntrance = -99.;
   tauVars.leadChargedPFCandPt = -99.;
   float sumEtaTimesEnergy = 0.;
   float sumPhiTimesEnergy = 0.;
   float sumEnergy = 0.;
   tauVars.phi = theTau.phi();
   // Check type of candidates building tau to avoid dynamic casts further
   bool isFromPFCands =
       (theTau.leadCand().isNonnull() && dynamic_cast<const reco::PFCandidate*>(theTau.leadCand().get()) != nullptr);
   if (!isPhase2_) {
     for (const auto& candidate : theTau.signalCands()) {
       float etaAtECalEntrance = candidate->eta();
       float phiAtECalEntrance = candidate->phi();
       const reco::Track* track = nullptr;
       if (isFromPFCands) {
         const reco::PFCandidate* pfCandidate = static_cast<const reco::PFCandidate*>(candidate.get());
         etaAtECalEntrance = pfCandidate->positionAtECALEntrance().eta();
         if (!usePhiAtEcalEntranceExtrapolation_) {
           phiAtECalEntrance = pfCandidate->positionAtECALEntrance().phi();
         } else {
           bool success = false;
           reco::Candidate::Point posAtECal = positionAtECalEntrance_(candidate.get(), success);
           if (success) {
             phiAtECalEntrance = posAtECal.phi();
           }
         }
         track = getTrackFromPFCand(pfCandidate);
       } else {
         bool success = false;
         reco::Candidate::Point posAtECal = positionAtECalEntrance_(candidate.get(), success);
         if (success) {
           etaAtECalEntrance = posAtECal.eta();
           phiAtECalEntrance = posAtECal.phi();
         }
         track = candidate->bestTrack();
       }
       if (track != nullptr) {
         if (track->pt() > tauVars.leadChargedPFCandPt) {
           tauVars.leadChargedPFCandEtaAtEcalEntrance = etaAtECalEntrance;
           tauVars.leadChargedPFCandPt = track->pt();
         }
       }
       sumEtaTimesEnergy += etaAtECalEntrance * candidate->energy();
       sumPhiTimesEnergy += phiAtECalEntrance * candidate->energy();
       sumEnergy += candidate->energy();
     }
     if (sumEnergy > 0.) {
       tauVars.etaAtEcalEntrance = sumEtaTimesEnergy / sumEnergy;
       tauVars.phi = sumPhiTimesEnergy / sumEnergy;
     }
     tauVars.emFraction = std::max(theTau.emFraction(), 0.f);
   } else {  // Phase2
     for (const auto& candidate : theTau.signalCands()) {
       float etaAtECalEntrance = candidate->eta();
       const reco::Track* track = nullptr;
       if (isFromPFCands) {
         const reco::PFCandidate* pfCandidate = static_cast<const reco::PFCandidate*>(candidate.get());
         etaAtECalEntrance = pfCandidate->positionAtECALEntrance().eta();
         if (std::abs(theTau.eta()) >= ecalBarrelEndcapEtaBorder_) {  //HGCal
           bool success = false;
           reco::Candidate::Point posAtECal = positionAtECalEntrance_(candidate.get(), success);
           if (success) {
             etaAtECalEntrance = posAtECal.eta();
           }
         }
         track = getTrackFromPFCand(pfCandidate);
       } else {
         bool success = false;
         reco::Candidate::Point posAtECal = positionAtECalEntrance_(candidate.get(), success);
         if (success) {
           etaAtECalEntrance = posAtECal.eta();
         }
         track = candidate->bestTrack();
       }
       if (track != nullptr) {
         if (track->pt() > tauVars.leadChargedPFCandPt) {
           tauVars.leadChargedPFCandEtaAtEcalEntrance = etaAtECalEntrance;
           tauVars.leadChargedPFCandPt = track->pt();
         }
       }
       sumEtaTimesEnergy += etaAtECalEntrance * candidate->energy();
       sumEnergy += candidate->energy();
     }
     if (sumEnergy > 0.) {
       tauVars.etaAtEcalEntrance = sumEtaTimesEnergy / sumEnergy;
     }
     if (isFromPFCands) {
       const reco::PFCandidate* pfLeadCandidte =
           static_cast<const reco::PFCandidate*>(theTau.leadChargedHadrCand().get());
       tauVars.emFraction =
           std::max(pfLeadCandidte->ecalEnergy() / (pfLeadCandidte->ecalEnergy() + pfLeadCandidte->hcalEnergy()), 0.);
     } else {
       const pat::PackedCandidate* patLeadCandiate =
           dynamic_cast<const pat::PackedCandidate*>(theTau.leadChargedHadrCand().get());
       if (patLeadCandiate != nullptr) {
         tauVars.emFraction = std::max(1. - patLeadCandiate->hcalFraction(), 0.);
       }
     }
   }
   tauVars.leadPFChargedHadrHoP = 0.;
   tauVars.leadPFChargedHadrEoP = 0.;
   if (theTau.leadChargedHadrCand()->p() > 0.) {
     if (isFromPFCands) {
       const reco::PFCandidate* pfLeadCandiate =
           static_cast<const reco::PFCandidate*>(theTau.leadChargedHadrCand().get());
       tauVars.leadPFChargedHadrHoP = pfLeadCandiate->hcalEnergy() / pfLeadCandiate->p();
       tauVars.leadPFChargedHadrEoP = pfLeadCandiate->ecalEnergy() / pfLeadCandiate->p();
     } else {
       const pat::PackedCandidate* patLeadCandiate =
           dynamic_cast<const pat::PackedCandidate*>(theTau.leadChargedHadrCand().get());
       if (patLeadCandiate != nullptr) {
         tauVars.leadPFChargedHadrHoP = patLeadCandiate->caloFraction() * patLeadCandiate->energy() *
                                        patLeadCandiate->hcalFraction() / patLeadCandiate->p();
         tauVars.leadPFChargedHadrEoP = patLeadCandiate->caloFraction() * patLeadCandiate->energy() *
                                        (1. - patLeadCandiate->hcalFraction()) / patLeadCandiate->p();
       }
     }
   }
 
   return tauVars;
 }

◆ getTauVarsTypeSpecific() [2/2]

template<class TauType , class ElectronType >

TauVars AntiElectronIDMVA6< TauType, ElectronType >::getTauVarsTypeSpecific ( const pat::Tau & theTau )

Definition at line 1102 of file AntiElectronIDMVA6.cc.

References funct::abs(), pat::Tau::ecalEnergyLeadChargedHadrCand(), pat::Tau::emFraction_MVA(), reco::LeafCandidate::eta(), pat::Tau::etaAtEcalEntrance(), pat::Tau::etaAtEcalEntranceLeadChargedCand(), f, edm::Ptr< T >::get(), pat::Tau::hcalEnergyLeadChargedHadrCand(), pat::Tau::leadChargedHadrCand(), SiStripPI::max, reco::LeafCandidate::phi(), pat::Tau::phiAtEcalEntrance(), pat::Tau::ptLeadChargedCand(), summarizeEdmComparisonLogfiles::success, and leptonTimeLifeInfo_common_cff::tauVars.

                                                                                               {
   TauVars tauVars;
   tauVars.etaAtEcalEntrance = theTau.etaAtEcalEntrance();
   tauVars.leadChargedPFCandEtaAtEcalEntrance = theTau.etaAtEcalEntranceLeadChargedCand();
   tauVars.leadChargedPFCandPt = theTau.ptLeadChargedCand();
   tauVars.phi = theTau.phi();
   if (!isPhase2_) {
     if (!usePhiAtEcalEntranceExtrapolation_) {
       tauVars.phi = theTau.phiAtEcalEntrance();
     } else {
       float etaAtECal = 0.;
       float phiAtECal = 0.;
       bool success = energyWeightedEtaAndPhiAtECal(theTau, etaAtECal, phiAtECal);
       if (success) {
         tauVars.phi = phiAtECal;
       }
     }
     tauVars.emFraction = std::max(theTau.emFraction_MVA(), 0.f);
   } else {
     if (std::abs(theTau.eta()) >= ecalBarrelEndcapEtaBorder_) {  //HGCal
       tauVars.etaAtEcalEntrance = -99.;
       tauVars.leadChargedPFCandEtaAtEcalEntrance = -99.;
       bool success = false;
       reco::Candidate::Point posAtECal = positionAtECalEntrance_(theTau.leadChargedHadrCand().get(), success);
       if (success) {
         tauVars.leadChargedPFCandEtaAtEcalEntrance = posAtECal.eta();
       }
       float phiAtECal = 0.;
       float etaAtECal = 0.;
       success = energyWeightedEtaAndPhiAtECal(theTau, etaAtECal, phiAtECal);
       if (success) {
         tauVars.etaAtEcalEntrance = etaAtECal;
       }
     }
     tauVars.emFraction = std::max(theTau.ecalEnergyLeadChargedHadrCand() /
                                       (theTau.ecalEnergyLeadChargedHadrCand() + theTau.hcalEnergyLeadChargedHadrCand()),
                                   0.f);
   }
   tauVars.leadPFChargedHadrHoP = 0.;
   tauVars.leadPFChargedHadrEoP = 0.;
   if (theTau.leadChargedHadrCand()->p() > 0.) {
     tauVars.leadPFChargedHadrHoP = theTau.hcalEnergyLeadChargedHadrCand() / theTau.leadChargedHadrCand()->p();
     tauVars.leadPFChargedHadrEoP = theTau.ecalEnergyLeadChargedHadrCand() / theTau.leadChargedHadrCand()->p();
   }
 
   return tauVars;
 }

◆ minimum()

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::minimum	(	double	a,
		double	b
	)

private

Definition at line 1005 of file AntiElectronIDMVA6.cc.

References a, funct::abs(), and b.

                                                                             {
   if (std::abs(b) < std::abs(a))
     return b;
   else
     return a;
 }

◆ mvaValue() [1/4]

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::mvaValue	(	const antiElecIDMVA6_blocks::TauVars &	tauVars,
		const antiElecIDMVA6_blocks::TauGammaVecs &	tauGammaVecs,
		const antiElecIDMVA6_blocks::ElecVars &	elecVars
	)

Definition at line 217 of file AntiElectronIDMVA6.cc.

                                                                                      {
   TauGammaMoms tauGammaMoms;
   double sumPt = 0.;
   double dEta2 = 0.;
   double dPhi2 = 0.;
   tauGammaMoms.signalPFGammaCandsIn = tauGammaVecs.gammasPtInSigCone.size();
   for (size_t i = 0; i < tauGammaVecs.gammasPtInSigCone.size(); ++i) {
     double pt_i = tauGammaVecs.gammasPtInSigCone[i];
     double phi_i = tauGammaVecs.gammasdPhiInSigCone[i];
     if (tauGammaVecs.gammasdPhiInSigCone[i] > M_PI)
       phi_i = tauGammaVecs.gammasdPhiInSigCone[i] - 2 * M_PI;
     else if (tauGammaVecs.gammasdPhiInSigCone[i] < -M_PI)
       phi_i = tauGammaVecs.gammasdPhiInSigCone[i] + 2 * M_PI;
     double eta_i = tauGammaVecs.gammasdEtaInSigCone[i];
     sumPt += pt_i;
     dEta2 += (pt_i * eta_i * eta_i);
     dPhi2 += (pt_i * phi_i * phi_i);
   }
 
   tauGammaMoms.gammaEnFracIn = -99.;
   if (tauVars.pt > 0.) {
     tauGammaMoms.gammaEnFracIn = sumPt / tauVars.pt;
   }
   if (sumPt > 0.) {
     dEta2 /= sumPt;
     dPhi2 /= sumPt;
   }
   tauGammaMoms.gammaEtaMomIn = std::sqrt(dEta2 * tauGammaMoms.gammaEnFracIn) * tauVars.pt;
   tauGammaMoms.gammaPhiMomIn = std::sqrt(dPhi2 * tauGammaMoms.gammaEnFracIn) * tauVars.pt;
 
   sumPt = 0.;
   dEta2 = 0.;
   dPhi2 = 0.;
   tauGammaMoms.signalPFGammaCandsOut = tauGammaVecs.gammasPtOutSigCone.size();
   for (size_t i = 0; i < tauGammaVecs.gammasPtOutSigCone.size(); ++i) {
     double pt_i = tauGammaVecs.gammasPtOutSigCone[i];
     double phi_i = tauGammaVecs.gammasdPhiOutSigCone[i];
     if (tauGammaVecs.gammasdPhiOutSigCone[i] > M_PI)
       phi_i = tauGammaVecs.gammasdPhiOutSigCone[i] - 2 * M_PI;
     else if (tauGammaVecs.gammasdPhiOutSigCone[i] < -M_PI)
       phi_i = tauGammaVecs.gammasdPhiOutSigCone[i] + 2 * M_PI;
     double eta_i = tauGammaVecs.gammasdEtaOutSigCone[i];
     sumPt += pt_i;
     dEta2 += (pt_i * eta_i * eta_i);
     dPhi2 += (pt_i * phi_i * phi_i);
   }
 
   tauGammaMoms.gammaEnFracOut = -99.;
   if (tauVars.pt > 0.) {
     tauGammaMoms.gammaEnFracOut = sumPt / tauVars.pt;
   }
   if (sumPt > 0.) {
     dEta2 /= sumPt;
     dPhi2 /= sumPt;
   }
   tauGammaMoms.gammaEtaMomOut = std::sqrt(dEta2 * tauGammaMoms.gammaEnFracOut) * tauVars.pt;
   tauGammaMoms.gammaPhiMomOut = std::sqrt(dPhi2 * tauGammaMoms.gammaEnFracOut) * tauVars.pt;
 
   if (!isPhase2_) {
     return mvaValue(tauVars, tauGammaMoms, elecVars);
   } else {
     return mvaValuePhase2(tauVars, tauGammaMoms, elecVars);
   }
 }

◆ mvaValue() [2/4]

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::mvaValue	(	const antiElecIDMVA6_blocks::TauVars &	tauVars,
		const antiElecIDMVA6_blocks::TauGammaMoms &	tauGammaMoms,
		const antiElecIDMVA6_blocks::ElecVars &	elecVars
	)

Definition at line 285 of file AntiElectronIDMVA6.cc.

                                                                                      {
   if (!isInitialized_) {
     throw cms::Exception("ClassNotInitialized") << " AntiElectronMVA6 not properly initialized !!\n";
   }
 
   double mvaValue = -99.;
 
   float elecDeltaPinPoutOverPin = (elecVars.pIn > 0.0) ? (std::abs(elecVars.pIn - elecVars.pOut) / elecVars.pIn) : 1.0;
   float elecEecalOverPout = (elecVars.pOut > 0.0) ? (elecVars.eEcal / elecVars.pOut) : 20.0;
   float elecNumHitsDiffOverSum =
       ((elecVars.gsfNumHits + elecVars.kfNumHits) > 0.0)
           ? ((elecVars.gsfNumHits - elecVars.kfNumHits) / (elecVars.gsfNumHits + elecVars.kfNumHits))
           : 1.0;
 
   if (deltaR(tauVars.etaAtEcalEntrance, tauVars.phi, elecVars.eta, elecVars.phi) > 0.3 &&
       tauGammaMoms.signalPFGammaCandsIn == 0 && tauVars.hasGsf < 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_NoEleMatch_woGwoGSF_Barrel_[0] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_woGwoGSF_Barrel_[1] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_NoEleMatch_woGwoGSF_Barrel_[2] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Barrel_[3] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Barrel_[4] = tauVars.emFraction;
       var_NoEleMatch_woGwoGSF_Barrel_[5] = tauVars.leadPFChargedHadrHoP;
       var_NoEleMatch_woGwoGSF_Barrel_[6] = tauVars.leadPFChargedHadrEoP;
       var_NoEleMatch_woGwoGSF_Barrel_[7] = tauVars.visMassIn;
       var_NoEleMatch_woGwoGSF_Barrel_[8] = tauVars.dCrackEta;
       var_NoEleMatch_woGwoGSF_Barrel_[9] = tauVars.dCrackPhi;
       mvaValue = mva_NoEleMatch_woGwoGSF_BL_->GetClassifier(var_NoEleMatch_woGwoGSF_Barrel_.data());
     } else {
       var_NoEleMatch_woGwoGSF_Endcap_[0] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_woGwoGSF_Endcap_[1] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_NoEleMatch_woGwoGSF_Endcap_[2] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Endcap_[3] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Endcap_[4] = tauVars.emFraction;
       var_NoEleMatch_woGwoGSF_Endcap_[5] = tauVars.leadPFChargedHadrHoP;
       var_NoEleMatch_woGwoGSF_Endcap_[6] = tauVars.leadPFChargedHadrEoP;
       var_NoEleMatch_woGwoGSF_Endcap_[7] = tauVars.visMassIn;
       var_NoEleMatch_woGwoGSF_Endcap_[8] = tauVars.dCrackEta;
       mvaValue = mva_NoEleMatch_woGwoGSF_EC_->GetClassifier(var_NoEleMatch_woGwoGSF_Endcap_.data());
     }
   } else if (deltaR(tauVars.etaAtEcalEntrance, tauVars.phi, elecVars.eta, elecVars.phi) > 0.3 &&
              tauGammaMoms.signalPFGammaCandsIn > 0 && tauVars.hasGsf < 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_NoEleMatch_wGwoGSF_Barrel_[0] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_wGwoGSF_Barrel_[1] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_NoEleMatch_wGwoGSF_Barrel_[2] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Barrel_[3] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Barrel_[4] = tauVars.emFraction;
       var_NoEleMatch_wGwoGSF_Barrel_[5] = tauGammaMoms.signalPFGammaCandsIn;
       var_NoEleMatch_wGwoGSF_Barrel_[6] = tauGammaMoms.signalPFGammaCandsOut;
       var_NoEleMatch_wGwoGSF_Barrel_[7] = tauVars.leadPFChargedHadrHoP;
       var_NoEleMatch_wGwoGSF_Barrel_[8] = tauVars.leadPFChargedHadrEoP;
       var_NoEleMatch_wGwoGSF_Barrel_[9] = tauVars.visMassIn;
       var_NoEleMatch_wGwoGSF_Barrel_[10] = tauGammaMoms.gammaEtaMomIn;
       var_NoEleMatch_wGwoGSF_Barrel_[11] = tauGammaMoms.gammaEtaMomOut;
       var_NoEleMatch_wGwoGSF_Barrel_[12] = tauGammaMoms.gammaPhiMomIn;
       var_NoEleMatch_wGwoGSF_Barrel_[13] = tauGammaMoms.gammaPhiMomOut;
       var_NoEleMatch_wGwoGSF_Barrel_[14] = tauGammaMoms.gammaEnFracIn;
       var_NoEleMatch_wGwoGSF_Barrel_[15] = tauGammaMoms.gammaEnFracOut;
       var_NoEleMatch_wGwoGSF_Barrel_[16] = tauVars.dCrackEta;
       var_NoEleMatch_wGwoGSF_Barrel_[17] = tauVars.dCrackPhi;
       mvaValue = mva_NoEleMatch_wGwoGSF_BL_->GetClassifier(var_NoEleMatch_wGwoGSF_Barrel_.data());
     } else {
       var_NoEleMatch_wGwoGSF_Endcap_[0] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_wGwoGSF_Endcap_[1] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_NoEleMatch_wGwoGSF_Endcap_[2] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Endcap_[3] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Endcap_[4] = tauVars.emFraction;
       var_NoEleMatch_wGwoGSF_Endcap_[5] = tauGammaMoms.signalPFGammaCandsIn;
       var_NoEleMatch_wGwoGSF_Endcap_[6] = tauGammaMoms.signalPFGammaCandsOut;
       var_NoEleMatch_wGwoGSF_Endcap_[7] = tauVars.leadPFChargedHadrHoP;
       var_NoEleMatch_wGwoGSF_Endcap_[8] = tauVars.leadPFChargedHadrEoP;
       var_NoEleMatch_wGwoGSF_Endcap_[9] = tauVars.visMassIn;
       var_NoEleMatch_wGwoGSF_Endcap_[10] = tauGammaMoms.gammaEtaMomIn;
       var_NoEleMatch_wGwoGSF_Endcap_[11] = tauGammaMoms.gammaEtaMomOut;
       var_NoEleMatch_wGwoGSF_Endcap_[12] = tauGammaMoms.gammaPhiMomIn;
       var_NoEleMatch_wGwoGSF_Endcap_[13] = tauGammaMoms.gammaPhiMomOut;
       var_NoEleMatch_wGwoGSF_Endcap_[14] = tauGammaMoms.gammaEnFracIn;
       var_NoEleMatch_wGwoGSF_Endcap_[15] = tauGammaMoms.gammaEnFracOut;
       var_NoEleMatch_wGwoGSF_Endcap_[16] = tauVars.dCrackEta;
       mvaValue = mva_NoEleMatch_wGwoGSF_EC_->GetClassifier(var_NoEleMatch_wGwoGSF_Endcap_.data());
     }
   } else if (tauGammaMoms.signalPFGammaCandsIn == 0 && tauVars.hasGsf > 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_woGwGSF_Barrel_[0] = std::max(-0.1f, elecVars.eTotOverPin);
       var_woGwGSF_Barrel_[1] = std::log(std::max(0.01f, elecVars.chi2NormGSF));
       var_woGwGSF_Barrel_[2] = elecVars.gsfNumHits;
       var_woGwGSF_Barrel_[3] = std::log(std::max(0.01f, elecVars.gsfTrackResol));
       var_woGwGSF_Barrel_[4] = elecVars.gsfTracklnPt;
       var_woGwGSF_Barrel_[5] = elecNumHitsDiffOverSum;
       var_woGwGSF_Barrel_[6] = std::log(std::max(0.01f, elecVars.chi2NormKF));
       var_woGwGSF_Barrel_[7] = std::min(elecDeltaPinPoutOverPin, 1.f);
       var_woGwGSF_Barrel_[8] = std::min(elecEecalOverPout, 20.f);
       var_woGwGSF_Barrel_[9] = elecVars.deltaEta;
       var_woGwGSF_Barrel_[10] = elecVars.deltaPhi;
       var_woGwGSF_Barrel_[11] = std::min(elecVars.mvaInSigmaEtaEta, 0.01f);
       var_woGwGSF_Barrel_[12] = std::min(elecVars.mvaInHadEnergy, 20.f);
       var_woGwGSF_Barrel_[13] = std::min(elecVars.mvaInDeltaEta, 0.1f);
       var_woGwGSF_Barrel_[14] = tauVars.etaAtEcalEntrance;
       var_woGwGSF_Barrel_[15] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_woGwGSF_Barrel_[16] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_woGwGSF_Barrel_[17] = std::log(std::max(1.f, tauVars.pt));
       var_woGwGSF_Barrel_[18] = tauVars.emFraction;
       var_woGwGSF_Barrel_[19] = tauVars.leadPFChargedHadrHoP;
       var_woGwGSF_Barrel_[20] = tauVars.leadPFChargedHadrEoP;
       var_woGwGSF_Barrel_[21] = tauVars.visMassIn;
       var_woGwGSF_Barrel_[22] = tauVars.dCrackEta;
       var_woGwGSF_Barrel_[23] = tauVars.dCrackPhi;
       mvaValue = mva_woGwGSF_BL_->GetClassifier(var_woGwGSF_Barrel_.data());
     } else {
       var_woGwGSF_Endcap_[0] = std::max(-0.1f, elecVars.eTotOverPin);
       var_woGwGSF_Endcap_[1] = std::log(std::max(0.01f, elecVars.chi2NormGSF));
       var_woGwGSF_Endcap_[2] = elecVars.gsfNumHits;
       var_woGwGSF_Endcap_[3] = std::log(std::max(0.01f, elecVars.gsfTrackResol));
       var_woGwGSF_Endcap_[4] = elecVars.gsfTracklnPt;
       var_woGwGSF_Endcap_[5] = elecNumHitsDiffOverSum;
       var_woGwGSF_Endcap_[6] = std::log(std::max(0.01f, elecVars.chi2NormKF));
       var_woGwGSF_Endcap_[7] = std::min(elecDeltaPinPoutOverPin, 1.f);
       var_woGwGSF_Endcap_[8] = std::min(elecEecalOverPout, 20.f);
       var_woGwGSF_Endcap_[9] = elecVars.deltaEta;
       var_woGwGSF_Endcap_[10] = elecVars.deltaPhi;
       var_woGwGSF_Endcap_[11] = std::min(elecVars.mvaInSigmaEtaEta, 0.01f);
       var_woGwGSF_Endcap_[12] = std::min(elecVars.mvaInHadEnergy, 20.f);
       var_woGwGSF_Endcap_[13] = std::min(elecVars.mvaInDeltaEta, 0.1f);
       var_woGwGSF_Endcap_[14] = tauVars.etaAtEcalEntrance;
       var_woGwGSF_Endcap_[15] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_woGwGSF_Endcap_[16] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_woGwGSF_Endcap_[17] = std::log(std::max(1.f, tauVars.pt));
       var_woGwGSF_Endcap_[18] = tauVars.emFraction;
       var_woGwGSF_Endcap_[19] = tauVars.leadPFChargedHadrHoP;
       var_woGwGSF_Endcap_[20] = tauVars.leadPFChargedHadrEoP;
       var_woGwGSF_Endcap_[21] = tauVars.visMassIn;
       var_woGwGSF_Endcap_[22] = tauVars.dCrackEta;
       mvaValue = mva_woGwGSF_EC_->GetClassifier(var_woGwGSF_Endcap_.data());
     }
   } else if (tauGammaMoms.signalPFGammaCandsIn > 0 && tauVars.hasGsf > 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_wGwGSF_Barrel_[0] = std::max(-0.1f, elecVars.eTotOverPin);
       var_wGwGSF_Barrel_[1] = std::log(std::max(0.01f, elecVars.chi2NormGSF));
       var_wGwGSF_Barrel_[2] = elecVars.gsfNumHits;
       var_wGwGSF_Barrel_[3] = std::log(std::max(0.01f, elecVars.gsfTrackResol));
       var_wGwGSF_Barrel_[4] = elecVars.gsfTracklnPt;
       var_wGwGSF_Barrel_[5] = elecNumHitsDiffOverSum;
       var_wGwGSF_Barrel_[6] = std::log(std::max(0.01f, elecVars.chi2NormKF));
       var_wGwGSF_Barrel_[7] = std::min(elecDeltaPinPoutOverPin, 1.f);
       var_wGwGSF_Barrel_[8] = std::min(elecEecalOverPout, 20.f);
       var_wGwGSF_Barrel_[9] = elecVars.deltaEta;
       var_wGwGSF_Barrel_[10] = elecVars.deltaPhi;
       var_wGwGSF_Barrel_[11] = std::min(elecVars.mvaInSigmaEtaEta, 0.01f);
       var_wGwGSF_Barrel_[12] = std::min(elecVars.mvaInHadEnergy, 20.f);
       var_wGwGSF_Barrel_[13] = std::min(elecVars.mvaInDeltaEta, 0.1f);
       var_wGwGSF_Barrel_[14] = tauVars.etaAtEcalEntrance;
       var_wGwGSF_Barrel_[15] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_wGwGSF_Barrel_[16] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_wGwGSF_Barrel_[17] = std::log(std::max(1.f, tauVars.pt));
       var_wGwGSF_Barrel_[18] = tauVars.emFraction;
       var_wGwGSF_Barrel_[19] = tauGammaMoms.signalPFGammaCandsIn;
       var_wGwGSF_Barrel_[20] = tauGammaMoms.signalPFGammaCandsOut;
       var_wGwGSF_Barrel_[21] = tauVars.leadPFChargedHadrHoP;
       var_wGwGSF_Barrel_[22] = tauVars.leadPFChargedHadrEoP;
       var_wGwGSF_Barrel_[23] = tauVars.visMassIn;
       var_wGwGSF_Barrel_[24] = tauGammaMoms.gammaEtaMomIn;
       var_wGwGSF_Barrel_[25] = tauGammaMoms.gammaEtaMomOut;
       var_wGwGSF_Barrel_[26] = tauGammaMoms.gammaPhiMomIn;
       var_wGwGSF_Barrel_[27] = tauGammaMoms.gammaPhiMomOut;
       var_wGwGSF_Barrel_[28] = tauGammaMoms.gammaEnFracIn;
       var_wGwGSF_Barrel_[29] = tauGammaMoms.gammaEnFracOut;
       var_wGwGSF_Barrel_[30] = tauVars.dCrackEta;
       var_wGwGSF_Barrel_[31] = tauVars.dCrackPhi;
       mvaValue = mva_wGwGSF_BL_->GetClassifier(var_wGwGSF_Barrel_.data());
     } else {
       var_wGwGSF_Endcap_[0] = std::max(-0.1f, elecVars.eTotOverPin);
       var_wGwGSF_Endcap_[1] = std::log(std::max(0.01f, elecVars.chi2NormGSF));
       var_wGwGSF_Endcap_[2] = elecVars.gsfNumHits;
       var_wGwGSF_Endcap_[3] = std::log(std::max(0.01f, elecVars.gsfTrackResol));
       var_wGwGSF_Endcap_[4] = elecVars.gsfTracklnPt;
       var_wGwGSF_Endcap_[5] = elecNumHitsDiffOverSum;
       var_wGwGSF_Endcap_[6] = std::log(std::max(0.01f, elecVars.chi2NormKF));
       var_wGwGSF_Endcap_[7] = std::min(elecDeltaPinPoutOverPin, 1.f);
       var_wGwGSF_Endcap_[8] = std::min(elecEecalOverPout, 20.f);
       var_wGwGSF_Endcap_[9] = elecVars.deltaEta;
       var_wGwGSF_Endcap_[10] = elecVars.deltaPhi;
       var_wGwGSF_Endcap_[11] = std::min(elecVars.mvaInSigmaEtaEta, 0.01f);
       var_wGwGSF_Endcap_[12] = std::min(elecVars.mvaInHadEnergy, 20.f);
       var_wGwGSF_Endcap_[13] = std::min(elecVars.mvaInDeltaEta, 0.1f);
       var_wGwGSF_Endcap_[14] = tauVars.etaAtEcalEntrance;
       var_wGwGSF_Endcap_[15] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_wGwGSF_Endcap_[16] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_wGwGSF_Endcap_[17] = std::log(std::max(1.f, tauVars.pt));
       var_wGwGSF_Endcap_[18] = tauVars.emFraction;
       var_wGwGSF_Endcap_[19] = tauGammaMoms.signalPFGammaCandsIn;
       var_wGwGSF_Endcap_[20] = tauGammaMoms.signalPFGammaCandsOut;
       var_wGwGSF_Endcap_[21] = tauVars.leadPFChargedHadrHoP;
       var_wGwGSF_Endcap_[22] = tauVars.leadPFChargedHadrEoP;
       var_wGwGSF_Endcap_[23] = tauVars.visMassIn;
       var_wGwGSF_Endcap_[24] = tauGammaMoms.gammaEtaMomIn;
       var_wGwGSF_Endcap_[25] = tauGammaMoms.gammaEtaMomOut;
       var_wGwGSF_Endcap_[26] = tauGammaMoms.gammaPhiMomIn;
       var_wGwGSF_Endcap_[27] = tauGammaMoms.gammaPhiMomOut;
       var_wGwGSF_Endcap_[28] = tauGammaMoms.gammaEnFracIn;
       var_wGwGSF_Endcap_[29] = tauGammaMoms.gammaEnFracOut;
       var_wGwGSF_Endcap_[30] = tauVars.dCrackEta;
       mvaValue = mva_wGwGSF_EC_->GetClassifier(var_wGwGSF_Endcap_.data());
     }
   }
   return mvaValue;
 }

◆ mvaValue() [3/4]

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::mvaValue	(	const TauType &	theTau,
		const ElectronRef &	theEleRef
	)

Definition at line 824 of file AntiElectronIDMVA6.cc.

References AntiElectronIDMVA6< TauType, ElectronType >::getElecVars(), AntiElectronIDMVA6< TauType, ElectronType >::getTauGammaVecs(), AntiElectronIDMVA6< TauType, ElectronType >::getTauVars(), and leptonTimeLifeInfo_common_cff::tauVars.

 {
   // === tau variables ===
   TauVars tauVars = AntiElectronIDMVA6<TauType, ElectronType>::getTauVars(theTau);
   TauGammaVecs tauGammaVecs = AntiElectronIDMVA6<TauType, ElectronType>::getTauGammaVecs(theTau);
 
   // === electron variables ===
   ElecVars elecVars = AntiElectronIDMVA6<TauType, ElectronType>::getElecVars(theEleRef);
 
   return mvaValue(tauVars, tauGammaVecs, elecVars);
 }

◆ mvaValue() [4/4]

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::mvaValue ( const TauType & theTau )

Definition at line 838 of file AntiElectronIDMVA6.cc.

References antiElecIDMVA6_blocks::ElecVars::eta, AntiElectronIDMVA6< TauType, ElectronType >::getTauGammaVecs(), AntiElectronIDMVA6< TauType, ElectronType >::getTauVars(), and leptonTimeLifeInfo_common_cff::tauVars.

                                                                                 {
   // === tau variables ===
   TauVars tauVars = AntiElectronIDMVA6<TauType, ElectronType>::getTauVars(theTau);
   TauGammaVecs tauGammaVecs = AntiElectronIDMVA6<TauType, ElectronType>::getTauGammaVecs(theTau);
 
   // === electron variables ===
   ElecVars elecVars;
   elecVars.eta = 9.9;  // Dummy value used in MVA training
 
   return mvaValue(tauVars, tauGammaVecs, elecVars);
 }

◆ mvaValuePhase2()

template<class TauType , class ElectronType >

double AntiElectronIDMVA6< TauType, ElectronType >::mvaValuePhase2	(	const antiElecIDMVA6_blocks::TauVars &	tauVars,
		const antiElecIDMVA6_blocks::TauGammaMoms &	tauGammaMoms,
		const antiElecIDMVA6_blocks::ElecVars &	elecVars
	)

Definition at line 496 of file AntiElectronIDMVA6.cc.

                                                                                            {
   if (!isInitialized_) {
     throw cms::Exception("ClassNotInitialized") << " AntiElectronMVA6 not properly initialized !!\n";
   }
 
   double mvaValue = -99.;
 
   //do not consider tau candidates outside the HGCal border at |eta|=3
   if (std::abs(tauVars.etaAtEcalEntrance) > 3.0) {
     return mvaValue;
   }
 
   float elecDeltaPinPoutOverPin = (elecVars.pIn > 0.0) ? (std::abs(elecVars.pIn - elecVars.pOut) / elecVars.pIn) : 1.0;
   float elecEecalOverPout = (elecVars.pOut > 0.0) ? (elecVars.eEcal / elecVars.pOut) : 20.0;
   float elecNumHitsDiffOverSum =
       ((elecVars.gsfNumHits + elecVars.kfNumHits) > 0.0)
           ? ((elecVars.gsfNumHits - elecVars.kfNumHits) / (elecVars.gsfNumHits + elecVars.kfNumHits))
           : 1.0;
 
   if (tauGammaMoms.signalPFGammaCandsIn == 0 && tauVars.hasGsf < 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_NoEleMatch_woGwoGSF_Barrel_[0] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Barrel_[1] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Barrel_[2] = tauVars.emFraction;
       var_NoEleMatch_woGwoGSF_Barrel_[3] = tauVars.leadPFChargedHadrHoP;
       var_NoEleMatch_woGwoGSF_Barrel_[4] = tauVars.leadPFChargedHadrEoP;
       var_NoEleMatch_woGwoGSF_Barrel_[5] = tauVars.visMassIn;
       var_NoEleMatch_woGwoGSF_Barrel_[6] = tauVars.dCrackEta;
       var_NoEleMatch_woGwoGSF_Barrel_[7] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_woGwoGSF_Barrel_[8] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       mvaValue = mva_NoEleMatch_woGwoGSF_BL_->GetClassifier(var_NoEleMatch_woGwoGSF_Barrel_.data());
     } else if (std::abs(tauVars.etaAtEcalEntrance) < ecalEndcapVFEndcapEtaBorder_) {
       var_NoEleMatch_woGwoGSF_Endcap_[0] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Endcap_[1] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_Endcap_[2] = tauVars.visMassIn;
       var_NoEleMatch_woGwoGSF_Endcap_[3] = tauVars.dCrackEta;
       var_NoEleMatch_woGwoGSF_Endcap_[4] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_woGwoGSF_Endcap_[5] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       mvaValue = mva_NoEleMatch_woGwoGSF_EC_->GetClassifier(var_NoEleMatch_woGwoGSF_Endcap_.data());
     } else {
       var_NoEleMatch_woGwoGSF_VFEndcap_[0] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_VFEndcap_[1] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_woGwoGSF_VFEndcap_[2] = tauVars.visMassIn;
       var_NoEleMatch_woGwoGSF_VFEndcap_[3] = tauVars.dCrackEta;
       var_NoEleMatch_woGwoGSF_VFEndcap_[4] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_woGwoGSF_VFEndcap_[5] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       mvaValue = mva_NoEleMatch_woGwoGSF_VFEC_->GetClassifier(var_NoEleMatch_woGwoGSF_VFEndcap_.data());
     }
   } else if (tauGammaMoms.signalPFGammaCandsIn > 0 && tauVars.hasGsf < 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_NoEleMatch_wGwoGSF_Barrel_[0] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Barrel_[1] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Barrel_[2] = tauVars.emFraction;
       var_NoEleMatch_wGwoGSF_Barrel_[3] = tauGammaMoms.signalPFGammaCandsIn;
       var_NoEleMatch_wGwoGSF_Barrel_[4] = tauGammaMoms.signalPFGammaCandsOut;
       var_NoEleMatch_wGwoGSF_Barrel_[5] = tauVars.leadPFChargedHadrHoP;
       var_NoEleMatch_wGwoGSF_Barrel_[6] = tauVars.leadPFChargedHadrEoP;
       var_NoEleMatch_wGwoGSF_Barrel_[7] = tauVars.visMassIn;
       var_NoEleMatch_wGwoGSF_Barrel_[7] = tauGammaMoms.gammaEtaMomIn;
       var_NoEleMatch_wGwoGSF_Barrel_[9] = tauGammaMoms.gammaEtaMomOut;
       var_NoEleMatch_wGwoGSF_Barrel_[10] = tauGammaMoms.gammaPhiMomIn;
       var_NoEleMatch_wGwoGSF_Barrel_[11] = tauGammaMoms.gammaPhiMomOut;
       var_NoEleMatch_wGwoGSF_Barrel_[12] = tauGammaMoms.gammaEnFracIn;
       var_NoEleMatch_wGwoGSF_Barrel_[13] = tauGammaMoms.gammaEnFracOut;
       var_NoEleMatch_wGwoGSF_Barrel_[14] = tauVars.dCrackEta;
       var_NoEleMatch_wGwoGSF_Barrel_[15] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_wGwoGSF_Barrel_[16] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       mvaValue = mva_NoEleMatch_wGwoGSF_BL_->GetClassifier(var_NoEleMatch_wGwoGSF_Barrel_.data());
     } else if (std::abs(tauVars.etaAtEcalEntrance) < ecalEndcapVFEndcapEtaBorder_) {
       var_NoEleMatch_wGwoGSF_Endcap_[0] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Endcap_[1] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_Endcap_[2] = tauGammaMoms.signalPFGammaCandsIn;
       var_NoEleMatch_wGwoGSF_Endcap_[3] = tauGammaMoms.signalPFGammaCandsOut;
       var_NoEleMatch_wGwoGSF_Endcap_[4] = tauVars.visMassIn;
       var_NoEleMatch_wGwoGSF_Endcap_[5] = tauGammaMoms.gammaEtaMomIn;
       var_NoEleMatch_wGwoGSF_Endcap_[6] = tauGammaMoms.gammaEtaMomOut;
       var_NoEleMatch_wGwoGSF_Endcap_[7] = tauGammaMoms.gammaPhiMomIn;
       var_NoEleMatch_wGwoGSF_Endcap_[8] = tauGammaMoms.gammaPhiMomOut;
       var_NoEleMatch_wGwoGSF_Endcap_[9] = tauGammaMoms.gammaEnFracIn;
       var_NoEleMatch_wGwoGSF_Endcap_[10] = tauGammaMoms.gammaEnFracOut;
       var_NoEleMatch_wGwoGSF_Endcap_[11] = tauVars.dCrackEta;
       var_NoEleMatch_wGwoGSF_Endcap_[12] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_wGwoGSF_Endcap_[13] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       mvaValue = mva_NoEleMatch_wGwoGSF_EC_->GetClassifier(var_NoEleMatch_wGwoGSF_Endcap_.data());
     } else {
       var_NoEleMatch_wGwoGSF_VFEndcap_[0] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_VFEndcap_[1] = std::log(std::max(1.f, tauVars.pt));
       var_NoEleMatch_wGwoGSF_VFEndcap_[2] = tauGammaMoms.signalPFGammaCandsIn;
       var_NoEleMatch_wGwoGSF_VFEndcap_[3] = tauGammaMoms.signalPFGammaCandsOut;
       var_NoEleMatch_wGwoGSF_VFEndcap_[4] = tauVars.visMassIn;
       var_NoEleMatch_wGwoGSF_VFEndcap_[5] = tauGammaMoms.gammaEtaMomIn;
       var_NoEleMatch_wGwoGSF_VFEndcap_[6] = tauGammaMoms.gammaEtaMomOut;
       var_NoEleMatch_wGwoGSF_VFEndcap_[7] = tauGammaMoms.gammaPhiMomIn;
       var_NoEleMatch_wGwoGSF_VFEndcap_[8] = tauGammaMoms.gammaPhiMomOut;
       var_NoEleMatch_wGwoGSF_VFEndcap_[9] = tauGammaMoms.gammaEnFracIn;
       var_NoEleMatch_wGwoGSF_VFEndcap_[10] = tauGammaMoms.gammaEnFracOut;
       var_NoEleMatch_wGwoGSF_VFEndcap_[11] = tauVars.dCrackEta;
       var_NoEleMatch_wGwoGSF_VFEndcap_[12] = tauVars.etaAtEcalEntrance;
       var_NoEleMatch_wGwoGSF_VFEndcap_[13] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       mvaValue = mva_NoEleMatch_wGwoGSF_VFEC_->GetClassifier(var_NoEleMatch_wGwoGSF_VFEndcap_.data());
     }
   } else if (tauGammaMoms.signalPFGammaCandsIn == 0 && tauVars.hasGsf > 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_woGwGSF_Barrel_[0] = std::log(std::max(0.1f, elecVars.chi2NormGSF));
       var_woGwGSF_Barrel_[1] = elecVars.gsfNumHits;
       var_woGwGSF_Barrel_[2] = std::log(std::max(0.1f, elecVars.gsfTrackResol));
       var_woGwGSF_Barrel_[3] = elecVars.gsfTracklnPt;
       var_woGwGSF_Barrel_[4] = elecNumHitsDiffOverSum;
       var_woGwGSF_Barrel_[5] = std::log(std::max(0.1f, elecVars.chi2NormKF));
       var_woGwGSF_Barrel_[6] = std::min(elecDeltaPinPoutOverPin, 1.f);
       var_woGwGSF_Barrel_[7] = std::min(elecEecalOverPout, 20.f);
       var_woGwGSF_Barrel_[8] = std::min(elecVars.mvaInSigmaEtaEta, 0.01f);
       var_woGwGSF_Barrel_[9] = std::min(elecVars.mvaInHadEnergy, 20.f);
       var_woGwGSF_Barrel_[10] = std::min(elecVars.mvaInDeltaEta, 0.1f);
       var_woGwGSF_Barrel_[11] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_woGwGSF_Barrel_[12] = std::log(std::max(1.f, tauVars.pt));
       var_woGwGSF_Barrel_[13] = tauVars.emFraction;
       var_woGwGSF_Barrel_[14] = tauVars.leadPFChargedHadrHoP;
       var_woGwGSF_Barrel_[15] = tauVars.leadPFChargedHadrEoP;
       var_woGwGSF_Barrel_[16] = tauVars.visMassIn;
       var_woGwGSF_Barrel_[17] = tauVars.dCrackEta;
       var_woGwGSF_Barrel_[18] = tauVars.etaAtEcalEntrance;
       var_woGwGSF_Barrel_[19] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_woGwGSF_Barrel_[20] = elecVars.deltaEta;
       var_woGwGSF_Barrel_[21] = elecVars.deltaPhi;
       var_woGwGSF_Barrel_[22] = elecVars.sigmaIEtaIEta5x5;
       var_woGwGSF_Barrel_[23] = elecVars.showerCircularity;
       var_woGwGSF_Barrel_[24] = elecVars.r9;
       var_woGwGSF_Barrel_[25] = elecVars.superClusterEtaWidth;
       var_woGwGSF_Barrel_[26] = elecVars.superClusterPhiWidth;
       mvaValue = mva_woGwGSF_BL_->GetClassifier(var_woGwGSF_Barrel_.data());
     } else if (std::abs(tauVars.etaAtEcalEntrance) < ecalEndcapVFEndcapEtaBorder_) {
       var_woGwGSF_Endcap_[0] = std::log(std::max(0.1f, elecVars.chi2NormGSF));
       var_woGwGSF_Endcap_[1] = elecVars.gsfNumHits;
       var_woGwGSF_Endcap_[2] = std::log(std::max(0.1f, elecVars.gsfTrackResol));
       var_woGwGSF_Endcap_[3] = elecVars.gsfTracklnPt;
       var_woGwGSF_Endcap_[4] = elecNumHitsDiffOverSum;
       var_woGwGSF_Endcap_[5] = std::log(std::max(0.1f, elecVars.chi2NormKF));
       var_woGwGSF_Endcap_[6] = elecVars.eEcal;
       var_woGwGSF_Endcap_[7] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_woGwGSF_Endcap_[8] = std::log(std::max(1.f, tauVars.pt));
       var_woGwGSF_Endcap_[9] = tauVars.visMassIn;
       var_woGwGSF_Endcap_[10] = tauVars.dCrackEta;
       var_woGwGSF_Endcap_[11] = tauVars.etaAtEcalEntrance;
       var_woGwGSF_Endcap_[12] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_woGwGSF_Endcap_[13] = elecVars.hgcalSigmaUU;
       var_woGwGSF_Endcap_[14] = elecVars.hgcalSigmaVV;
       var_woGwGSF_Endcap_[15] = elecVars.hgcalSigmaEE;
       var_woGwGSF_Endcap_[16] = elecVars.hgcalSigmaPP;
       var_woGwGSF_Endcap_[17] = elecVars.hgcalNLayers;
       var_woGwGSF_Endcap_[18] = elecVars.hgcalLastLayer;
       var_woGwGSF_Endcap_[19] = elecVars.hgcalLayerEfrac10;
       var_woGwGSF_Endcap_[20] = elecVars.hgcalLayerEfrac90;
       var_woGwGSF_Endcap_[21] = elecVars.hgcalEcEnergyEE;
       var_woGwGSF_Endcap_[22] = elecVars.hgcalEcEnergyFH;
       var_woGwGSF_Endcap_[23] = elecVars.hgcalMeasuredDepth;
       var_woGwGSF_Endcap_[24] = elecVars.hgcalExpectedDepth;
       var_woGwGSF_Endcap_[25] = elecVars.hgcalDepthCompatibility;
       var_woGwGSF_Endcap_[26] = elecVars.deltaEta;
       var_woGwGSF_Endcap_[27] = elecVars.deltaPhi;
       var_woGwGSF_Endcap_[28] = elecVars.eSeedClusterOverPout;
       var_woGwGSF_Endcap_[29] = elecVars.superClusterEtaWidth;
       var_woGwGSF_Endcap_[30] = elecVars.superClusterPhiWidth;
       mvaValue = mva_woGwGSF_EC_->GetClassifier(var_woGwGSF_Endcap_.data());
     } else {
       var_woGwGSF_VFEndcap_[0] = std::log(std::max(0.1f, elecVars.chi2NormGSF));
       var_woGwGSF_VFEndcap_[1] = elecVars.gsfNumHits;
       var_woGwGSF_VFEndcap_[2] = std::log(std::max(0.1f, elecVars.gsfTrackResol));
       var_woGwGSF_VFEndcap_[3] = elecVars.gsfTracklnPt;
       var_woGwGSF_VFEndcap_[4] = elecNumHitsDiffOverSum;
       var_woGwGSF_VFEndcap_[5] = std::log(std::max(0.1f, elecVars.chi2NormKF));
       var_woGwGSF_VFEndcap_[6] = elecVars.eEcal;
       var_woGwGSF_VFEndcap_[7] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_woGwGSF_VFEndcap_[8] = std::log(std::max(1.f, tauVars.pt));
       var_woGwGSF_VFEndcap_[9] = tauVars.visMassIn;
       var_woGwGSF_VFEndcap_[10] = tauVars.dCrackEta;
       var_woGwGSF_VFEndcap_[11] = tauVars.etaAtEcalEntrance;
       var_woGwGSF_VFEndcap_[12] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_woGwGSF_VFEndcap_[13] = elecVars.hgcalSigmaUU;
       var_woGwGSF_VFEndcap_[14] = elecVars.hgcalSigmaVV;
       var_woGwGSF_VFEndcap_[15] = elecVars.hgcalSigmaEE;
       var_woGwGSF_VFEndcap_[16] = elecVars.hgcalSigmaPP;
       var_woGwGSF_VFEndcap_[17] = elecVars.hgcalNLayers;
       var_woGwGSF_VFEndcap_[18] = elecVars.hgcalLastLayer;
       var_woGwGSF_VFEndcap_[19] = elecVars.hgcalLayerEfrac10;
       var_woGwGSF_VFEndcap_[20] = elecVars.hgcalLayerEfrac90;
       var_woGwGSF_VFEndcap_[21] = elecVars.hgcalEcEnergyEE;
       var_woGwGSF_VFEndcap_[22] = elecVars.hgcalEcEnergyFH;
       var_woGwGSF_VFEndcap_[23] = elecVars.hgcalMeasuredDepth;
       var_woGwGSF_VFEndcap_[24] = elecVars.hgcalExpectedDepth;
       var_woGwGSF_VFEndcap_[25] = elecVars.hgcalExpectedSigma;
       var_woGwGSF_VFEndcap_[26] = elecVars.hgcalDepthCompatibility;
       var_woGwGSF_VFEndcap_[27] = elecVars.deltaEta;
       var_woGwGSF_VFEndcap_[28] = elecVars.deltaPhi;
       var_woGwGSF_VFEndcap_[29] = elecVars.eSeedClusterOverPout;
       var_woGwGSF_VFEndcap_[30] = elecVars.superClusterEtaWidth;
       var_woGwGSF_VFEndcap_[31] = elecVars.superClusterPhiWidth;
       mvaValue = mva_woGwGSF_VFEC_->GetClassifier(var_woGwGSF_VFEndcap_.data());
     }
   } else if (tauGammaMoms.signalPFGammaCandsIn > 0 && tauVars.hasGsf > 0.5) {
     if (std::abs(tauVars.etaAtEcalEntrance) < ecalBarrelEndcapEtaBorder_) {
       var_wGwGSF_Barrel_[0] = std::log(std::max(0.1f, elecVars.chi2NormGSF));
       var_wGwGSF_Barrel_[1] = elecVars.gsfNumHits;
       var_wGwGSF_Barrel_[2] = std::log(std::max(0.1f, elecVars.gsfTrackResol));
       var_wGwGSF_Barrel_[3] = elecVars.gsfTracklnPt;
       var_wGwGSF_Barrel_[4] = elecNumHitsDiffOverSum;
       var_wGwGSF_Barrel_[5] = std::log(std::max(0.1f, elecVars.chi2NormKF));
       var_wGwGSF_Barrel_[6] = std::min(elecDeltaPinPoutOverPin, 1.f);
       var_wGwGSF_Barrel_[7] = std::min(elecEecalOverPout, 20.f);
       var_wGwGSF_Barrel_[8] = std::min(elecVars.mvaInSigmaEtaEta, 0.01f);
       var_wGwGSF_Barrel_[9] = std::min(elecVars.mvaInHadEnergy, 20.f);
       var_wGwGSF_Barrel_[10] = std::min(elecVars.mvaInDeltaEta, 0.1f);
       var_wGwGSF_Barrel_[11] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_wGwGSF_Barrel_[12] = std::log(std::max(1.f, tauVars.pt));
       var_wGwGSF_Barrel_[13] = tauVars.emFraction;
       var_wGwGSF_Barrel_[14] = tauGammaMoms.signalPFGammaCandsIn;
       var_wGwGSF_Barrel_[15] = tauGammaMoms.signalPFGammaCandsOut;
       var_wGwGSF_Barrel_[16] = tauVars.leadPFChargedHadrHoP;
       var_wGwGSF_Barrel_[17] = tauVars.leadPFChargedHadrEoP;
       var_wGwGSF_Barrel_[18] = tauVars.visMassIn;
       var_wGwGSF_Barrel_[19] = tauGammaMoms.gammaEtaMomIn;
       var_wGwGSF_Barrel_[20] = tauGammaMoms.gammaEtaMomOut;
       var_wGwGSF_Barrel_[21] = tauGammaMoms.gammaPhiMomIn;
       var_wGwGSF_Barrel_[22] = tauGammaMoms.gammaPhiMomOut;
       var_wGwGSF_Barrel_[23] = tauGammaMoms.gammaEnFracIn;
       var_wGwGSF_Barrel_[24] = tauGammaMoms.gammaEnFracOut;
       var_wGwGSF_Barrel_[25] = tauVars.dCrackEta;
       var_wGwGSF_Barrel_[26] = tauVars.etaAtEcalEntrance;
       var_wGwGSF_Barrel_[27] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_wGwGSF_Barrel_[28] = elecVars.deltaEta;
       var_wGwGSF_Barrel_[29] = elecVars.deltaPhi;
       var_wGwGSF_Barrel_[30] = elecVars.sigmaIPhiIPhi5x5;
       var_wGwGSF_Barrel_[31] = elecVars.sigmaIEtaIEta5x5;
       var_wGwGSF_Barrel_[32] = elecVars.showerCircularity;
       var_wGwGSF_Barrel_[33] = elecVars.eSeedClusterOverPout;
       var_wGwGSF_Barrel_[34] = elecVars.superClusterEtaWidth;
       var_wGwGSF_Barrel_[35] = elecVars.superClusterPhiWidth;
       mvaValue = mva_wGwGSF_BL_->GetClassifier(var_wGwGSF_Barrel_.data());
     } else if (std::abs(tauVars.etaAtEcalEntrance) < ecalEndcapVFEndcapEtaBorder_) {
       var_wGwGSF_Endcap_[0] = std::log(std::max(0.1f, elecVars.chi2NormGSF));
       var_wGwGSF_Endcap_[1] = elecVars.gsfNumHits;
       var_wGwGSF_Endcap_[2] = std::log(std::max(0.1f, elecVars.gsfTrackResol));
       var_wGwGSF_Endcap_[3] = elecVars.gsfTracklnPt;
       var_wGwGSF_Endcap_[4] = elecNumHitsDiffOverSum;
       var_wGwGSF_Endcap_[5] = std::log(std::max(0.1f, elecVars.chi2NormKF));
       var_wGwGSF_Endcap_[6] = elecVars.eEcal;
       var_wGwGSF_Endcap_[7] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_wGwGSF_Endcap_[8] = std::log(std::max(1.f, tauVars.pt));
       var_wGwGSF_Endcap_[9] = tauGammaMoms.signalPFGammaCandsIn;
       var_wGwGSF_Endcap_[10] = tauGammaMoms.signalPFGammaCandsOut;
       var_wGwGSF_Endcap_[11] = tauVars.visMassIn;
       var_wGwGSF_Endcap_[12] = tauGammaMoms.gammaEtaMomIn;
       var_wGwGSF_Endcap_[13] = tauGammaMoms.gammaEtaMomOut;
       var_wGwGSF_Endcap_[14] = tauGammaMoms.gammaPhiMomIn;
       var_wGwGSF_Endcap_[15] = tauGammaMoms.gammaPhiMomOut;
       var_wGwGSF_Endcap_[16] = tauGammaMoms.gammaEnFracIn;
       var_wGwGSF_Endcap_[17] = tauGammaMoms.gammaEnFracOut;
       var_wGwGSF_Endcap_[18] = tauVars.dCrackEta;
       var_wGwGSF_Endcap_[19] = tauVars.etaAtEcalEntrance;
       var_wGwGSF_Endcap_[20] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_wGwGSF_Endcap_[21] = elecVars.hgcalSigmaVV;
       var_wGwGSF_Endcap_[22] = elecVars.hgcalSigmaEE;
       var_wGwGSF_Endcap_[23] = elecVars.hgcalSigmaPP;
       var_wGwGSF_Endcap_[24] = elecVars.hgcalNLayers;
       var_wGwGSF_Endcap_[25] = elecVars.hgcalFirstLayer;
       var_wGwGSF_Endcap_[26] = elecVars.hgcalLastLayer;
       var_wGwGSF_Endcap_[27] = elecVars.hgcalLayerEfrac10;
       var_wGwGSF_Endcap_[28] = elecVars.hgcalLayerEfrac90;
       var_wGwGSF_Endcap_[29] = elecVars.hgcalEcEnergyEE;
       var_wGwGSF_Endcap_[30] = elecVars.hgcalEcEnergyFH;
       var_wGwGSF_Endcap_[31] = elecVars.hgcalMeasuredDepth;
       var_wGwGSF_Endcap_[32] = elecVars.hgcalExpectedDepth;
       var_wGwGSF_Endcap_[33] = elecVars.deltaEta;
       var_wGwGSF_Endcap_[34] = elecVars.deltaPhi;
       var_wGwGSF_Endcap_[35] = elecVars.eSeedClusterOverPout;
       var_wGwGSF_Endcap_[36] = elecVars.superClusterEtaWidth;
       var_wGwGSF_Endcap_[37] = elecVars.superClusterPhiWidth;
       mvaValue = mva_wGwGSF_EC_->GetClassifier(var_wGwGSF_Endcap_.data());
     } else {
       var_wGwGSF_VFEndcap_[0] = std::log(std::max(0.1f, elecVars.chi2NormGSF));
       var_wGwGSF_VFEndcap_[1] = elecVars.gsfNumHits;
       var_wGwGSF_VFEndcap_[2] = std::log(std::max(0.1f, elecVars.gsfTrackResol));
       var_wGwGSF_VFEndcap_[3] = elecVars.gsfTracklnPt;
       var_wGwGSF_VFEndcap_[4] = elecNumHitsDiffOverSum;
       var_wGwGSF_VFEndcap_[5] = std::log(std::max(0.1f, elecVars.chi2NormKF));
       var_wGwGSF_VFEndcap_[6] = elecVars.eEcal;
       var_wGwGSF_VFEndcap_[7] = std::min(2.f, tauVars.leadChargedPFCandPt / std::max(1.f, tauVars.pt));
       var_wGwGSF_VFEndcap_[8] = std::log(std::max(1.f, tauVars.pt));
       var_wGwGSF_VFEndcap_[9] = tauGammaMoms.signalPFGammaCandsIn;
       var_wGwGSF_VFEndcap_[10] = tauGammaMoms.signalPFGammaCandsOut;
       var_wGwGSF_VFEndcap_[11] = tauVars.visMassIn;
       var_wGwGSF_VFEndcap_[12] = tauGammaMoms.gammaEtaMomIn;
       var_wGwGSF_VFEndcap_[13] = tauGammaMoms.gammaEtaMomOut;
       var_wGwGSF_VFEndcap_[14] = tauGammaMoms.gammaPhiMomIn;
       var_wGwGSF_VFEndcap_[15] = tauGammaMoms.gammaPhiMomOut;
       var_wGwGSF_VFEndcap_[16] = tauGammaMoms.gammaEnFracIn;
       var_wGwGSF_VFEndcap_[17] = tauGammaMoms.gammaEnFracOut;
       var_wGwGSF_VFEndcap_[18] = tauVars.dCrackEta;
       var_wGwGSF_VFEndcap_[19] = tauVars.etaAtEcalEntrance;
       var_wGwGSF_VFEndcap_[20] = tauVars.leadChargedPFCandEtaAtEcalEntrance;
       var_wGwGSF_VFEndcap_[21] = elecVars.hgcalSigmaUU;
       var_wGwGSF_VFEndcap_[22] = elecVars.hgcalSigmaVV;
       var_wGwGSF_VFEndcap_[23] = elecVars.hgcalSigmaEE;
       var_wGwGSF_VFEndcap_[24] = elecVars.hgcalSigmaPP;
       var_wGwGSF_VFEndcap_[25] = elecVars.hgcalNLayers;
       var_wGwGSF_VFEndcap_[26] = elecVars.hgcalLastLayer;
       var_wGwGSF_VFEndcap_[27] = elecVars.hgcalLayerEfrac10;
       var_wGwGSF_VFEndcap_[28] = elecVars.hgcalLayerEfrac90;
       var_wGwGSF_VFEndcap_[29] = elecVars.hgcalEcEnergyEE;
       var_wGwGSF_VFEndcap_[30] = elecVars.hgcalEcEnergyFH;
       var_wGwGSF_VFEndcap_[31] = elecVars.hgcalMeasuredDepth;
       var_wGwGSF_VFEndcap_[32] = elecVars.hgcalExpectedDepth;
       var_wGwGSF_VFEndcap_[33] = elecVars.hgcalExpectedSigma;
       var_wGwGSF_VFEndcap_[34] = elecVars.hgcalDepthCompatibility;
       var_wGwGSF_VFEndcap_[35] = elecVars.deltaEta;
       var_wGwGSF_VFEndcap_[36] = elecVars.deltaPhi;
       var_wGwGSF_VFEndcap_[37] = elecVars.eSeedClusterOverPout;
       var_wGwGSF_VFEndcap_[38] = elecVars.superClusterEtaWidth;
       var_wGwGSF_VFEndcap_[39] = elecVars.superClusterPhiWidth;
       mvaValue = mva_wGwGSF_VFEC_->GetClassifier(var_wGwGSF_VFEndcap_.data());
     }
   }
   return mvaValue;
 }

Member Data Documentation

◆ ecalBarrelEndcapEtaBorder_

template<class TauType , class ElectronType >

constexpr float AntiElectronIDMVA6< TauType, ElectronType >::ecalBarrelEndcapEtaBorder_ = 1.479

staticprivate

Definition at line 162 of file AntiElectronIDMVA6.h.

◆ ecalEndcapVFEndcapEtaBorder_

template<class TauType , class ElectronType >

constexpr float AntiElectronIDMVA6< TauType, ElectronType >::ecalEndcapVFEndcapEtaBorder_ = 2.4

staticprivate

Definition at line 163 of file AntiElectronIDMVA6.h.

◆ electronIds_

template<class TauType , class ElectronType >

std::map<std::string, edm::Handle<edm::ValueMap<float> > > AntiElectronIDMVA6< TauType, ElectronType >::electronIds_

private

Definition at line 230 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ electronIds_tokens_

template<class TauType , class ElectronType >

std::map<std::string, edm::EDGetTokenT<edm::ValueMap<float> > > AntiElectronIDMVA6< TauType, ElectronType >::electronIds_tokens_

private

Definition at line 229 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ inputFileName_

template<class TauType , class ElectronType >

edm::FileInPath AntiElectronIDMVA6< TauType, ElectronType >::inputFileName_

private

Definition at line 167 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ inputFilesToDelete_

template<class TauType , class ElectronType >

std::vector<TFile*> AntiElectronIDMVA6< TauType, ElectronType >::inputFilesToDelete_

private

Definition at line 223 of file AntiElectronIDMVA6.h.

◆ isInitialized_

template<class TauType , class ElectronType >

bool AntiElectronIDMVA6< TauType, ElectronType >::isInitialized_

private

Definition at line 165 of file AntiElectronIDMVA6.h.

◆ isPhase2_

template<class TauType , class ElectronType >

const bool AntiElectronIDMVA6< TauType, ElectronType >::isPhase2_

private

Definition at line 225 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ loadMVAfromDB_

template<class TauType , class ElectronType >

bool AntiElectronIDMVA6< TauType, ElectronType >::loadMVAfromDB_

private

Definition at line 166 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mva_NoEleMatch_wGwoGSF_BL_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_NoEleMatch_wGwoGSF_BL_

private

Definition at line 211 of file AntiElectronIDMVA6.h.

◆ mva_NoEleMatch_wGwoGSF_EC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_NoEleMatch_wGwoGSF_EC_

private

Definition at line 215 of file AntiElectronIDMVA6.h.

◆ mva_NoEleMatch_wGwoGSF_VFEC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_NoEleMatch_wGwoGSF_VFEC_

private

Definition at line 219 of file AntiElectronIDMVA6.h.

◆ mva_NoEleMatch_woGwoGSF_BL_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_NoEleMatch_woGwoGSF_BL_

private

Definition at line 210 of file AntiElectronIDMVA6.h.

◆ mva_NoEleMatch_woGwoGSF_EC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_NoEleMatch_woGwoGSF_EC_

private

Definition at line 214 of file AntiElectronIDMVA6.h.

◆ mva_NoEleMatch_woGwoGSF_VFEC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_NoEleMatch_woGwoGSF_VFEC_

private

Definition at line 218 of file AntiElectronIDMVA6.h.

◆ mva_wGwGSF_BL_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_wGwGSF_BL_

private

Definition at line 213 of file AntiElectronIDMVA6.h.

◆ mva_wGwGSF_EC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_wGwGSF_EC_

private

Definition at line 217 of file AntiElectronIDMVA6.h.

◆ mva_wGwGSF_VFEC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_wGwGSF_VFEC_

private

Definition at line 221 of file AntiElectronIDMVA6.h.

◆ mva_woGwGSF_BL_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_woGwGSF_BL_

private

Definition at line 212 of file AntiElectronIDMVA6.h.

◆ mva_woGwGSF_EC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_woGwGSF_EC_

private

Definition at line 216 of file AntiElectronIDMVA6.h.

◆ mva_woGwGSF_VFEC_

template<class TauType , class ElectronType >

const GBRForest* AntiElectronIDMVA6< TauType, ElectronType >::mva_woGwGSF_VFEC_

private

Definition at line 220 of file AntiElectronIDMVA6.h.

◆ mvaName_NoEleMatch_wGwoGSF_BL_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_wGwoGSF_BL_

private

Definition at line 170 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_NoEleMatch_wGwoGSF_EC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_wGwoGSF_EC_

private

Definition at line 174 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_NoEleMatch_wGwoGSF_VFEC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_wGwoGSF_VFEC_

private

Definition at line 178 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_NoEleMatch_woGwoGSF_BL_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_woGwoGSF_BL_

private

Definition at line 169 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_NoEleMatch_woGwoGSF_EC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_woGwoGSF_EC_

private

Definition at line 173 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_NoEleMatch_woGwoGSF_VFEC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_woGwoGSF_VFEC_

private

Definition at line 177 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_wGwGSF_BL_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_wGwGSF_BL_

private

Definition at line 172 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_wGwGSF_EC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_wGwGSF_EC_

private

Definition at line 176 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_wGwGSF_VFEC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_wGwGSF_VFEC_

private

Definition at line 180 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_woGwGSF_BL_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_woGwGSF_BL_

private

Definition at line 171 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_woGwGSF_EC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_woGwGSF_EC_

private

Definition at line 175 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaName_woGwGSF_VFEC_

template<class TauType , class ElectronType >

std::string AntiElectronIDMVA6< TauType, ElectronType >::mvaName_woGwGSF_VFEC_

private

Definition at line 179 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_NoEleMatch_wGwoGSF_BL_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_wGwoGSF_BL_

private

Definition at line 183 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_NoEleMatch_wGwoGSF_EC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_wGwoGSF_EC_

private

Definition at line 187 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_NoEleMatch_wGwoGSF_VFEC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_wGwoGSF_VFEC_

private

Definition at line 191 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_NoEleMatch_woGwoGSF_BL_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_woGwoGSF_BL_

private

Definition at line 182 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_NoEleMatch_woGwoGSF_EC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_woGwoGSF_EC_

private

Definition at line 186 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_NoEleMatch_woGwoGSF_VFEC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_woGwoGSF_VFEC_

private

Definition at line 190 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_wGwGSF_BL_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_wGwGSF_BL_

private

Definition at line 185 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_wGwGSF_EC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_wGwGSF_EC_

private

Definition at line 189 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_wGwGSF_VFEC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_wGwGSF_VFEC_

private

Definition at line 193 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_woGwGSF_BL_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_woGwGSF_BL_

private

Definition at line 184 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_woGwGSF_EC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_woGwGSF_EC_

private

Definition at line 188 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ mvaToken_woGwGSF_VFEC_

template<class TauType , class ElectronType >

edm::ESGetToken<GBRForest, GBRWrapperRcd> AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_woGwGSF_VFEC_

private

Definition at line 192 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ positionAtECalEntrance_

template<class TauType , class ElectronType >

PositionAtECalEntranceComputer AntiElectronIDMVA6< TauType, ElectronType >::positionAtECalEntrance_

private

Definition at line 227 of file AntiElectronIDMVA6.h.

◆ usePhiAtEcalEntranceExtrapolation_

template<class TauType , class ElectronType >

bool AntiElectronIDMVA6< TauType, ElectronType >::usePhiAtEcalEntranceExtrapolation_

private

Definition at line 195 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_NoEleMatch_wGwoGSF_Barrel_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_wGwoGSF_Barrel_

private

Definition at line 198 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_NoEleMatch_wGwoGSF_Endcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_wGwoGSF_Endcap_

private

Definition at line 202 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_NoEleMatch_wGwoGSF_VFEndcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_wGwoGSF_VFEndcap_

private

Definition at line 206 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_NoEleMatch_woGwoGSF_Barrel_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_woGwoGSF_Barrel_

private

Definition at line 197 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_NoEleMatch_woGwoGSF_Endcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_woGwoGSF_Endcap_

private

Definition at line 201 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_NoEleMatch_woGwoGSF_VFEndcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_woGwoGSF_VFEndcap_

private

Definition at line 205 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_wGwGSF_Barrel_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_wGwGSF_Barrel_

private

Definition at line 200 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_wGwGSF_Endcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_wGwGSF_Endcap_

private

Definition at line 204 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_wGwGSF_VFEndcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_wGwGSF_VFEndcap_

private

Definition at line 208 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_woGwGSF_Barrel_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_woGwGSF_Barrel_

private

Definition at line 199 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_woGwGSF_Endcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_woGwGSF_Endcap_

private

Definition at line 203 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ var_woGwGSF_VFEndcap_

template<class TauType , class ElectronType >

std::vector<float> AntiElectronIDMVA6< TauType, ElectronType >::var_woGwGSF_VFEndcap_

private

Definition at line 207 of file AntiElectronIDMVA6.h.

Referenced by AntiElectronIDMVA6< TauType, ElectronType >::AntiElectronIDMVA6().

◆ verbosity_

template<class TauType , class ElectronType >

const int AntiElectronIDMVA6< TauType, ElectronType >::verbosity_

private

Definition at line 232 of file AntiElectronIDMVA6.h.

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

template<class TauType, class ElectronType> class AntiElectronIDMVA6< TauType, ElectronType >

Member Typedef Documentation

◆ ElectronCollection

◆ ElectronRef

Constructor & Destructor Documentation

◆ AntiElectronIDMVA6()

◆ ~AntiElectronIDMVA6()

Member Function Documentation

◆ beginEvent()

◆ dCrackEta()

◆ dCrackPhi()

◆ energyWeightedEtaAndPhiAtECal()

◆ getElecVars()

◆ getElecVarsHGCalTypeSpecific() [1/2]

◆ getElecVarsHGCalTypeSpecific() [2/2]

◆ getTauGammaVecs()

◆ getTauVars()

◆ getTauVarsTypeSpecific() [1/2]

◆ getTauVarsTypeSpecific() [2/2]

◆ minimum()

◆ mvaValue() [1/4]

◆ mvaValue() [2/4]

◆ mvaValue() [3/4]

◆ mvaValue() [4/4]

◆ mvaValuePhase2()

Member Data Documentation

◆ ecalBarrelEndcapEtaBorder_

◆ ecalEndcapVFEndcapEtaBorder_

◆ electronIds_

◆ electronIds_tokens_

◆ inputFileName_

◆ inputFilesToDelete_

◆ isInitialized_

◆ isPhase2_

◆ loadMVAfromDB_

◆ mva_NoEleMatch_wGwoGSF_BL_

◆ mva_NoEleMatch_wGwoGSF_EC_

◆ mva_NoEleMatch_wGwoGSF_VFEC_

◆ mva_NoEleMatch_woGwoGSF_BL_

◆ mva_NoEleMatch_woGwoGSF_EC_

◆ mva_NoEleMatch_woGwoGSF_VFEC_

◆ mva_wGwGSF_BL_

◆ mva_wGwGSF_EC_

◆ mva_wGwGSF_VFEC_

◆ mva_woGwGSF_BL_

◆ mva_woGwGSF_EC_

◆ mva_woGwGSF_VFEC_

◆ mvaName_NoEleMatch_wGwoGSF_BL_

◆ mvaName_NoEleMatch_wGwoGSF_EC_

◆ mvaName_NoEleMatch_wGwoGSF_VFEC_

◆ mvaName_NoEleMatch_woGwoGSF_BL_

◆ mvaName_NoEleMatch_woGwoGSF_EC_

◆ mvaName_NoEleMatch_woGwoGSF_VFEC_

◆ mvaName_wGwGSF_BL_

◆ mvaName_wGwGSF_EC_

◆ mvaName_wGwGSF_VFEC_

◆ mvaName_woGwGSF_BL_

◆ mvaName_woGwGSF_EC_

◆ mvaName_woGwGSF_VFEC_

◆ mvaToken_NoEleMatch_wGwoGSF_BL_

◆ mvaToken_NoEleMatch_wGwoGSF_EC_

◆ mvaToken_NoEleMatch_wGwoGSF_VFEC_

◆ mvaToken_NoEleMatch_woGwoGSF_BL_

◆ mvaToken_NoEleMatch_woGwoGSF_EC_

◆ mvaToken_NoEleMatch_woGwoGSF_VFEC_

◆ mvaToken_wGwGSF_BL_

◆ mvaToken_wGwGSF_EC_

◆ mvaToken_wGwGSF_VFEC_

◆ mvaToken_woGwGSF_BL_

◆ mvaToken_woGwGSF_EC_

◆ mvaToken_woGwGSF_VFEC_

◆ positionAtECalEntrance_

◆ usePhiAtEcalEntranceExtrapolation_

◆ var_NoEleMatch_wGwoGSF_Barrel_

template<class TauType, class ElectronType>
class AntiElectronIDMVA6< TauType, ElectronType >