AntiElectronIDMVA6< TauType, ElectronType > Class Template Reference

#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

template<class TauType , class ElectronType >

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

Definition at line 118 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

Definition at line 119 of file AntiElectronIDMVA6.h.

Constructor & Destructor Documentation

template<class TauType , class ElectronType >

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

Definition at line 36 of file AntiElectronIDMVA6.cc.

References AntiElectronIDMVA6< TauType, ElectronType >::electronIds_, AntiElectronIDMVA6< TauType, ElectronType >::electronIds_tokens_, Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), AntiElectronIDMVA6< TauType, ElectronType >::inputFileName_, beam_dqm_sourceclient-live_cfg::inputTag, AntiElectronIDMVA6< TauType, ElectronType >::isPhase2_, AntiElectronIDMVA6< TauType, ElectronType >::loadMVAfromDB_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_wGwoGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_wGwoGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_wGwoGSF_VFEC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_woGwoGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_woGwoGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_NoEleMatch_woGwoGSF_VFEC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_wGwGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_wGwGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_wGwGSF_VFEC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_woGwGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_woGwGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaName_woGwGSF_VFEC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_wGwoGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_wGwoGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_wGwoGSF_VFEC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_woGwoGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_woGwoGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_NoEleMatch_woGwoGSF_VFEC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_wGwGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_wGwGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_wGwGSF_VFEC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_woGwGSF_BL_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_woGwGSF_EC_, AntiElectronIDMVA6< TauType, ElectronType >::mvaToken_woGwGSF_VFEC_, AlCaHLTBitMon_QueryRunRegistry::string, AntiElectronIDMVA6< TauType, ElectronType >::usePhiAtEcalEntranceExtrapolation_, relativeConstraints::value, AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_wGwoGSF_Barrel_, AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_wGwoGSF_Endcap_, AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_wGwoGSF_VFEndcap_, AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_woGwoGSF_Barrel_, AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_woGwoGSF_Endcap_, AntiElectronIDMVA6< TauType, ElectronType >::var_NoEleMatch_woGwoGSF_VFEndcap_, AntiElectronIDMVA6< TauType, ElectronType >::var_wGwGSF_Barrel_, AntiElectronIDMVA6< TauType, ElectronType >::var_wGwGSF_Endcap_, AntiElectronIDMVA6< TauType, ElectronType >::var_wGwGSF_VFEndcap_, AntiElectronIDMVA6< TauType, ElectronType >::var_woGwGSF_Barrel_, AntiElectronIDMVA6< TauType, ElectronType >::var_woGwGSF_Endcap_, and AntiElectronIDMVA6< TauType, ElectronType >::var_woGwGSF_VFEndcap_.

    : 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>>();
      }
    }
  }
}

template<class TauType , class ElectronType >

AntiElectronIDMVA6< TauType, ElectronType >::~AntiElectronIDMVA6

(

)

Definition at line 130 of file AntiElectronIDMVA6.cc.

                                                               {
  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

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(), 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);
    }
  }
}

template<class TauType , class ElectronType >

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

private

Definition at line 1082 of file AntiElectronIDMVA6.cc.

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

                                                                      {
  //--- 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);
}

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(), M_PI, 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);
}

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 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.);
}

template<class TauType , class ElectronType >

ElecVars AntiElectronIDMVA6< TauType, ElectronType >::getElecVars

(

const ElectronRef &

theEleRef

)

Definition at line 924 of file AntiElectronIDMVA6.cc.

References antiElecIDMVA6_blocks::ElecVars::chi2NormGSF, antiElecIDMVA6_blocks::ElecVars::chi2NormKF, antiElecIDMVA6_blocks::ElecVars::deltaEta, antiElecIDMVA6_blocks::ElecVars::deltaPhi, antiElecIDMVA6_blocks::ElecVars::eEcal, antiElecIDMVA6_blocks::ElecVars::eSeedClusterOverPout, antiElecIDMVA6_blocks::ElecVars::eta, antiElecIDMVA6_blocks::ElecVars::eTotOverPin, Exception, antiElecIDMVA6_blocks::ElecVars::gsfNumHits, antiElecIDMVA6_blocks::ElecVars::gsfTracklnPt, antiElecIDMVA6_blocks::ElecVars::gsfTrackResol, edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), antiElecIDMVA6_blocks::ElecVars::kfNumHits, log, antiElecIDMVA6_blocks::ElecVars::mvaInDeltaEta, antiElecIDMVA6_blocks::ElecVars::mvaInHadEnergy, antiElecIDMVA6_blocks::ElecVars::mvaInSigmaEtaEta, antiElecIDMVA6_blocks::ElecVars::phi, antiElecIDMVA6_blocks::ElecVars::pIn, antiElecIDMVA6_blocks::ElecVars::pOut, antiElecIDMVA6_blocks::ElecVars::r9, antiElecIDMVA6_blocks::ElecVars::showerCircularity, antiElecIDMVA6_blocks::ElecVars::sigmaIEtaIEta5x5, antiElecIDMVA6_blocks::ElecVars::sigmaIPhiIPhi5x5, mathSSE::sqrt(), antiElecIDMVA6_blocks::ElecVars::superClusterEtaWidth, antiElecIDMVA6_blocks::ElecVars::superClusterPhiWidth, and relativeConstraints::value.

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;
}

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.

References antiElecIDMVA6_blocks::ElecVars::hgcalDepthCompatibility, antiElecIDMVA6_blocks::ElecVars::hgcalEcEnergyEE, antiElecIDMVA6_blocks::ElecVars::hgcalEcEnergyFH, antiElecIDMVA6_blocks::ElecVars::hgcalExpectedDepth, antiElecIDMVA6_blocks::ElecVars::hgcalExpectedSigma, antiElecIDMVA6_blocks::ElecVars::hgcalFirstLayer, antiElecIDMVA6_blocks::ElecVars::hgcalLastLayer, antiElecIDMVA6_blocks::ElecVars::hgcalLayerEfrac10, antiElecIDMVA6_blocks::ElecVars::hgcalLayerEfrac90, antiElecIDMVA6_blocks::ElecVars::hgcalMeasuredDepth, antiElecIDMVA6_blocks::ElecVars::hgcalNLayers, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaEE, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaPP, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaUU, and antiElecIDMVA6_blocks::ElecVars::hgcalSigmaVV.

                                                                                  {
  //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];
}

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.

References antiElecIDMVA6_blocks::ElecVars::hgcalDepthCompatibility, antiElecIDMVA6_blocks::ElecVars::hgcalEcEnergyEE, antiElecIDMVA6_blocks::ElecVars::hgcalEcEnergyFH, antiElecIDMVA6_blocks::ElecVars::hgcalExpectedDepth, antiElecIDMVA6_blocks::ElecVars::hgcalExpectedSigma, antiElecIDMVA6_blocks::ElecVars::hgcalFirstLayer, antiElecIDMVA6_blocks::ElecVars::hgcalLastLayer, antiElecIDMVA6_blocks::ElecVars::hgcalLayerEfrac10, antiElecIDMVA6_blocks::ElecVars::hgcalLayerEfrac90, antiElecIDMVA6_blocks::ElecVars::hgcalMeasuredDepth, antiElecIDMVA6_blocks::ElecVars::hgcalNLayers, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaEE, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaPP, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaUU, and antiElecIDMVA6_blocks::ElecVars::hgcalSigmaVV.

                                                                              {
  //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");
}

template<class TauType , class ElectronType >

TauGammaVecs AntiElectronIDMVA6< TauType, ElectronType >::getTauGammaVecs

(

const TauType &

theTau

)

Definition at line 901 of file AntiElectronIDMVA6.cc.

References HLT_FULL_cff::deltaR, 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;
}

template<class TauType , class ElectronType >

TauVars AntiElectronIDMVA6< TauType, ElectronType >::getTauVars

(

const TauType &

theTau

)

Definition at line 851 of file AntiElectronIDMVA6.cc.

References funct::abs(), antiElecIDMVA6_blocks::TauVars::dCrackEta, antiElecIDMVA6_blocks::TauVars::dCrackPhi, HLT_FULL_cff::deltaR, antiElecIDMVA6_blocks::TauVars::etaAtEcalEntrance, Exception, reco::PFCandidate::gsfTrackRef(), antiElecIDMVA6_blocks::TauVars::hasGsf, edm::Ref< C, T, F >::isNonnull(), ResonanceBuilder::mass, SiStripPI::max, pat::PackedCandidate::pdgId(), antiElecIDMVA6_blocks::TauVars::phi, antiElecIDMVA6_blocks::TauVars::pt, relativeConstraints::value, and antiElecIDMVA6_blocks::TauVars::visMassIn.

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;
}

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(), antiElecIDMVA6_blocks::TauVars::emFraction, reco::PFTau::emFraction(), pat::PackedCandidate::energy(), reco::LeafCandidate::eta(), antiElecIDMVA6_blocks::TauVars::etaAtEcalEntrance, edm::Ptr< T >::get(), reco::PFCandidate::hcalEnergy(), pat::PackedCandidate::hcalFraction(), edm::Ptr< T >::isNonnull(), reco::PFTau::leadCand(), reco::PFTau::leadChargedHadrCand(), antiElecIDMVA6_blocks::TauVars::leadChargedPFCandEtaAtEcalEntrance, antiElecIDMVA6_blocks::TauVars::leadChargedPFCandPt, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrEoP, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrHoP, SiStripPI::max, reco::LeafCandidate::p(), pat::PackedCandidate::p(), antiElecIDMVA6_blocks::TauVars::phi, reco::LeafCandidate::phi(), reco::PFCandidate::positionAtECALEntrance(), reco::TrackBase::pt(), reco::PFTau::signalCands(), and summarizeEdmComparisonLogfiles::success.

                                                                                                 {
  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;
}

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(), antiElecIDMVA6_blocks::TauVars::emFraction, pat::Tau::emFraction_MVA(), reco::LeafCandidate::eta(), antiElecIDMVA6_blocks::TauVars::etaAtEcalEntrance, pat::Tau::etaAtEcalEntrance(), pat::Tau::etaAtEcalEntranceLeadChargedCand(), validate-o2o-wbm::f, edm::Ptr< T >::get(), pat::Tau::hcalEnergyLeadChargedHadrCand(), pat::Tau::leadChargedHadrCand(), antiElecIDMVA6_blocks::TauVars::leadChargedPFCandEtaAtEcalEntrance, antiElecIDMVA6_blocks::TauVars::leadChargedPFCandPt, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrEoP, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrHoP, SiStripPI::max, antiElecIDMVA6_blocks::TauVars::phi, reco::LeafCandidate::phi(), pat::Tau::phiAtEcalEntrance(), pat::Tau::ptLeadChargedCand(), and summarizeEdmComparisonLogfiles::success.

                                                                                              {
  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;
}

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;
}

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.

References antiElecIDMVA6_blocks::TauGammaMoms::gammaEnFracIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaEnFracOut, antiElecIDMVA6_blocks::TauGammaMoms::gammaEtaMomIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaEtaMomOut, antiElecIDMVA6_blocks::TauGammaMoms::gammaPhiMomIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaPhiMomOut, antiElecIDMVA6_blocks::TauGammaVecs::gammasdEtaInSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasdEtaOutSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasdPhiInSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasdPhiOutSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasPtInSigCone, antiElecIDMVA6_blocks::TauGammaVecs::gammasPtOutSigCone, mps_fire::i, M_PI, antiElecIDMVA6_blocks::TauVars::pt, antiElecIDMVA6_blocks::TauGammaMoms::signalPFGammaCandsIn, antiElecIDMVA6_blocks::TauGammaMoms::signalPFGammaCandsOut, and mathSSE::sqrt().

                                                                                     {
  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);
  }
}

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.

References funct::abs(), antiElecIDMVA6_blocks::ElecVars::chi2NormGSF, antiElecIDMVA6_blocks::ElecVars::chi2NormKF, antiElecIDMVA6_blocks::TauVars::dCrackEta, antiElecIDMVA6_blocks::TauVars::dCrackPhi, antiElecIDMVA6_blocks::ElecVars::deltaEta, antiElecIDMVA6_blocks::ElecVars::deltaPhi, HLT_FULL_cff::deltaR, antiElecIDMVA6_blocks::ElecVars::eEcal, antiElecIDMVA6_blocks::TauVars::emFraction, antiElecIDMVA6_blocks::ElecVars::eta, antiElecIDMVA6_blocks::TauVars::etaAtEcalEntrance, antiElecIDMVA6_blocks::ElecVars::eTotOverPin, Exception, validate-o2o-wbm::f, antiElecIDMVA6_blocks::TauGammaMoms::gammaEnFracIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaEnFracOut, antiElecIDMVA6_blocks::TauGammaMoms::gammaEtaMomIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaEtaMomOut, antiElecIDMVA6_blocks::TauGammaMoms::gammaPhiMomIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaPhiMomOut, antiElecIDMVA6_blocks::ElecVars::gsfNumHits, antiElecIDMVA6_blocks::ElecVars::gsfTracklnPt, antiElecIDMVA6_blocks::ElecVars::gsfTrackResol, antiElecIDMVA6_blocks::TauVars::hasGsf, antiElecIDMVA6_blocks::ElecVars::kfNumHits, antiElecIDMVA6_blocks::TauVars::leadChargedPFCandEtaAtEcalEntrance, antiElecIDMVA6_blocks::TauVars::leadChargedPFCandPt, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrEoP, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrHoP, log, SiStripPI::max, min(), antiElecIDMVA6_blocks::ElecVars::mvaInDeltaEta, antiElecIDMVA6_blocks::ElecVars::mvaInHadEnergy, antiElecIDMVA6_blocks::ElecVars::mvaInSigmaEtaEta, antiElecIDMVA6_blocks::TauVars::phi, antiElecIDMVA6_blocks::ElecVars::phi, antiElecIDMVA6_blocks::ElecVars::pIn, antiElecIDMVA6_blocks::ElecVars::pOut, antiElecIDMVA6_blocks::TauVars::pt, antiElecIDMVA6_blocks::TauGammaMoms::signalPFGammaCandsIn, antiElecIDMVA6_blocks::TauGammaMoms::signalPFGammaCandsOut, and antiElecIDMVA6_blocks::TauVars::visMassIn.

                                                                                     {
  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;
}

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(), and AntiElectronIDMVA6< TauType, ElectronType >::getTauVars().

{
  // === 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);
}

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(), and AntiElectronIDMVA6< TauType, ElectronType >::getTauVars().

                                                                                {
  // === 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);
}

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.

References funct::abs(), antiElecIDMVA6_blocks::ElecVars::chi2NormGSF, antiElecIDMVA6_blocks::ElecVars::chi2NormKF, antiElecIDMVA6_blocks::TauVars::dCrackEta, antiElecIDMVA6_blocks::ElecVars::deltaEta, antiElecIDMVA6_blocks::ElecVars::deltaPhi, antiElecIDMVA6_blocks::ElecVars::eEcal, antiElecIDMVA6_blocks::TauVars::emFraction, antiElecIDMVA6_blocks::ElecVars::eSeedClusterOverPout, antiElecIDMVA6_blocks::TauVars::etaAtEcalEntrance, Exception, validate-o2o-wbm::f, antiElecIDMVA6_blocks::TauGammaMoms::gammaEnFracIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaEnFracOut, antiElecIDMVA6_blocks::TauGammaMoms::gammaEtaMomIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaEtaMomOut, antiElecIDMVA6_blocks::TauGammaMoms::gammaPhiMomIn, antiElecIDMVA6_blocks::TauGammaMoms::gammaPhiMomOut, antiElecIDMVA6_blocks::ElecVars::gsfNumHits, antiElecIDMVA6_blocks::ElecVars::gsfTracklnPt, antiElecIDMVA6_blocks::ElecVars::gsfTrackResol, antiElecIDMVA6_blocks::TauVars::hasGsf, antiElecIDMVA6_blocks::ElecVars::hgcalDepthCompatibility, antiElecIDMVA6_blocks::ElecVars::hgcalEcEnergyEE, antiElecIDMVA6_blocks::ElecVars::hgcalEcEnergyFH, antiElecIDMVA6_blocks::ElecVars::hgcalExpectedDepth, antiElecIDMVA6_blocks::ElecVars::hgcalExpectedSigma, antiElecIDMVA6_blocks::ElecVars::hgcalFirstLayer, antiElecIDMVA6_blocks::ElecVars::hgcalLastLayer, antiElecIDMVA6_blocks::ElecVars::hgcalLayerEfrac10, antiElecIDMVA6_blocks::ElecVars::hgcalLayerEfrac90, antiElecIDMVA6_blocks::ElecVars::hgcalMeasuredDepth, antiElecIDMVA6_blocks::ElecVars::hgcalNLayers, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaEE, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaPP, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaUU, antiElecIDMVA6_blocks::ElecVars::hgcalSigmaVV, antiElecIDMVA6_blocks::ElecVars::kfNumHits, antiElecIDMVA6_blocks::TauVars::leadChargedPFCandEtaAtEcalEntrance, antiElecIDMVA6_blocks::TauVars::leadChargedPFCandPt, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrEoP, antiElecIDMVA6_blocks::TauVars::leadPFChargedHadrHoP, log, SiStripPI::max, min(), antiElecIDMVA6_blocks::ElecVars::mvaInDeltaEta, antiElecIDMVA6_blocks::ElecVars::mvaInHadEnergy, antiElecIDMVA6_blocks::ElecVars::mvaInSigmaEtaEta, antiElecIDMVA6_blocks::ElecVars::pIn, antiElecIDMVA6_blocks::ElecVars::pOut, antiElecIDMVA6_blocks::TauVars::pt, antiElecIDMVA6_blocks::ElecVars::r9, antiElecIDMVA6_blocks::ElecVars::showerCircularity, antiElecIDMVA6_blocks::ElecVars::sigmaIEtaIEta5x5, antiElecIDMVA6_blocks::ElecVars::sigmaIPhiIPhi5x5, antiElecIDMVA6_blocks::TauGammaMoms::signalPFGammaCandsIn, antiElecIDMVA6_blocks::TauGammaMoms::signalPFGammaCandsOut, antiElecIDMVA6_blocks::ElecVars::superClusterEtaWidth, antiElecIDMVA6_blocks::ElecVars::superClusterPhiWidth, and antiElecIDMVA6_blocks::TauVars::visMassIn.

                                                                                           {
  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

template<class TauType , class ElectronType >

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

staticprivate

Definition at line 162 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

staticprivate

Definition at line 163 of file AntiElectronIDMVA6.h.

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().

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().

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().

template<class TauType , class ElectronType >

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

private

Definition at line 223 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

bool AntiElectronIDMVA6< TauType, ElectronType >::isInitialized_

private

Definition at line 165 of file AntiElectronIDMVA6.h.

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().

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().

template<class TauType , class ElectronType >

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

private

Definition at line 211 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 215 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 219 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 210 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 214 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 218 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 213 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 217 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 221 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 212 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 216 of file AntiElectronIDMVA6.h.

template<class TauType , class ElectronType >

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

private

Definition at line 220 of file AntiElectronIDMVA6.h.

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

template<class TauType , class ElectronType >

PositionAtECalEntranceComputer AntiElectronIDMVA6< TauType, ElectronType >::positionAtECalEntrance_

private

Definition at line 227 of file AntiElectronIDMVA6.h.

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

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().

template<class TauType , class ElectronType >

const int AntiElectronIDMVA6< TauType, ElectronType >::verbosity_

private

Definition at line 232 of file AntiElectronIDMVA6.h.