PFRecoTauDiscriminationByMVAIsolationRun2 Class Reference

Inheritance diagram for PFRecoTauDiscriminationByMVAIsolationRun2:

Public Member Functions
void	beginEvent (const edm::Event &, const edm::EventSetup &)
double	discriminate (const PFTauRef &) const
void	endEvent (edm::Event &)
	PFRecoTauDiscriminationByMVAIsolationRun2 (const edm::ParameterSet &cfg)
	~PFRecoTauDiscriminationByMVAIsolationRun2 ()
Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
virtual double	discriminate (const TauRef &tau) const =0
void	produce (edm::Event &, const edm::EventSetup &)
	TauDiscriminationProducerBase (const edm::ParameterSet &iConfig)
	TauDiscriminationProducerBase ()
virtual	~TauDiscriminationProducerBase ()
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
enum	{   kOldDMwoLT, kOldDMwLT, kNewDMwoLT, kNewDMwLT,   kDBoldDMwLT, kDBnewDMwLT, kPWoldDMwLT, kPWnewDMwLT }
typedef edm::AssociationVector < reco::PFTauRefProd, std::vector < reco::PFTauTransverseImpactParameterRef > >	PFTauTIPAssociationByRef

Private Attributes
std::auto_ptr< PFTauDiscriminator >	category_output_
edm::EDGetTokenT < reco::PFTauDiscriminator >	ChargedIsoPtSum_token
edm::Handle < reco::PFTauDiscriminator >	chargedIsoPtSums_
edm::Handle < reco::PFTauDiscriminator >	footprintCorrection_
edm::EDGetTokenT < reco::PFTauDiscriminator >	FootprintCorrection_token
edm::FileInPath	inputFileName_
std::vector< TFile * >	inputFilesToDelete_
bool	loadMVAfromDB_
std::string	moduleLabel_
float *	mvaInput_
std::string	mvaName_
int	mvaOpt_
const GBRForest *	mvaReader_
edm::EDGetTokenT < reco::PFTauDiscriminator >	NeutralIsoPtSum_token
edm::Handle < reco::PFTauDiscriminator >	neutralIsoPtSums_
edm::Handle < reco::PFTauDiscriminator >	photonPtSumOutsideSignalCone_
edm::EDGetTokenT < reco::PFTauDiscriminator >	PhotonPtSumOutsideSignalCone_token
edm::EDGetTokenT < reco::PFTauDiscriminator >	PUcorrPtSum_token
edm::Handle < reco::PFTauDiscriminator >	puCorrPtSums_
edm::Handle < PFTauTIPAssociationByRef >	tauLifetimeInfos
edm::Handle< TauCollection >	taus_
edm::EDGetTokenT < PFTauTIPAssociationByRef >	TauTransverseImpactParameters_token
int	verbosity_

Additional Inherited Members
Public Types inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
typedef std::vector< TauType >	TauCollection
typedef edm::Ref< TauCollection >	TauRef
typedef edm::RefProd < TauCollection >	TauRefProd
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
Protected Attributes inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
std::string	moduleLabel_
double	prediscriminantFailValue_
edm::EDGetTokenT< TauCollection >	Tau_token
size_t	tauIndex_
edm::InputTag	TauProducer_

Detailed Description

Definition at line 66 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

Member Typedef Documentation

typedef edm::AssociationVector<reco::PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef> > PFRecoTauDiscriminationByMVAIsolationRun2::PFTauTIPAssociationByRef

private

Definition at line 143 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

Member Enumeration Documentation

anonymous enum

private

Enumerator
kOldDMwoLT
kOldDMwLT
kNewDMwoLT
kNewDMwLT
kDBoldDMwLT
kDBnewDMwLT
kPWoldDMwLT
kPWnewDMwLT

Definition at line 139 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

{ kOldDMwoLT, kOldDMwLT, kNewDMwoLT, kNewDMwLT, kDBoldDMwLT, kDBnewDMwLT, kPWoldDMwLT, kPWnewDMwLT };

Constructor & Destructor Documentation

PFRecoTauDiscriminationByMVAIsolationRun2::PFRecoTauDiscriminationByMVAIsolationRun2 ( const edm::ParameterSet &  cfg )

inlineexplicit

Definition at line 69 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

References assert(), Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

    : PFTauDiscriminationProducerBase(cfg),
      moduleLabel_(cfg.getParameter<std::string>("@module_label")),
      mvaReader_(0),
      mvaInput_(0),
      category_output_(0)
  {
    mvaName_ = cfg.getParameter<std::string>("mvaName");
    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";
    }    
    std::string mvaOpt_string = cfg.getParameter<std::string>("mvaOpt");
    if      ( mvaOpt_string == "oldDMwoLT" ) mvaOpt_ = kOldDMwoLT;
    else if ( mvaOpt_string == "oldDMwLT"  ) mvaOpt_ = kOldDMwLT;
    else if ( mvaOpt_string == "newDMwoLT" ) mvaOpt_ = kNewDMwoLT;
    else if ( mvaOpt_string == "newDMwLT"  ) mvaOpt_ = kNewDMwLT;
    else if ( mvaOpt_string == "DBoldDMwLT"  ) mvaOpt_ = kDBoldDMwLT;
    else if ( mvaOpt_string == "DBnewDMwLT"  ) mvaOpt_ = kDBnewDMwLT;
    else if ( mvaOpt_string == "PWoldDMwLT"  ) mvaOpt_ = kPWoldDMwLT;
    else if ( mvaOpt_string == "PWnewDMwLT"  ) mvaOpt_ = kPWnewDMwLT;
    else throw cms::Exception("PFRecoTauDiscriminationByMVAIsolationRun2")
      << " Invalid Configuration Parameter 'mvaOpt' = " << mvaOpt_string << " !!\n";
    
    if      ( mvaOpt_ == kOldDMwoLT || mvaOpt_ == kNewDMwoLT ) mvaInput_ = new float[6];
    else if ( mvaOpt_ == kOldDMwLT  || mvaOpt_ == kNewDMwLT  ) mvaInput_ = new float[12];
    else if ( mvaOpt_ == kDBoldDMwLT || mvaOpt_ == kDBnewDMwLT ||
          mvaOpt_ == kPWoldDMwLT || mvaOpt_ == kPWnewDMwLT) mvaInput_ = new float[23];
    else assert(0);

    TauTransverseImpactParameters_token = consumes<PFTauTIPAssociationByRef>(cfg.getParameter<edm::InputTag>("srcTauTransverseImpactParameters"));
    
    ChargedIsoPtSum_token = consumes<reco::PFTauDiscriminator>(cfg.getParameter<edm::InputTag>("srcChargedIsoPtSum"));
    NeutralIsoPtSum_token = consumes<reco::PFTauDiscriminator>(cfg.getParameter<edm::InputTag>("srcNeutralIsoPtSum"));
    PUcorrPtSum_token = consumes<reco::PFTauDiscriminator>(cfg.getParameter<edm::InputTag>("srcPUcorrPtSum"));
    PhotonPtSumOutsideSignalCone_token = consumes<reco::PFTauDiscriminator>(cfg.getParameter<edm::InputTag>("srcPhotonPtSumOutsideSignalCone"));
    FootprintCorrection_token = consumes<reco::PFTauDiscriminator>(cfg.getParameter<edm::InputTag>("srcFootprintCorrection"));
  
    verbosity_ = ( cfg.exists("verbosity") ) ?
      cfg.getParameter<int>("verbosity") : 0;

    produces<PFTauDiscriminator>("category");
  }

PFRecoTauDiscriminationByMVAIsolationRun2::~PFRecoTauDiscriminationByMVAIsolationRun2 ( )

inline

Definition at line 121 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

  {
    if(!loadMVAfromDB_) delete mvaReader_;
    delete[] mvaInput_;
    for ( std::vector<TFile*>::iterator it = inputFilesToDelete_.begin();
      it != inputFilesToDelete_.end(); ++it ) {
      delete (*it);
    }
  }

Member Function Documentation

void PFRecoTauDiscriminationByMVAIsolationRun2::beginEvent ( const edm::Event &  evt, const edm::EventSetup &  es  )

virtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 166 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

References edm::Event::getByToken(), and mvaPFMET_cff::loadMVAfromDB.

{
  if ( !mvaReader_ ) {
    if ( loadMVAfromDB_ ) {
      mvaReader_ = loadMVAfromDB(es, mvaName_);
    } else {
      mvaReader_ = loadMVAfromFile(inputFileName_, mvaName_, inputFilesToDelete_);
    }
  }

  evt.getByToken(TauTransverseImpactParameters_token, tauLifetimeInfos);

  evt.getByToken(ChargedIsoPtSum_token, chargedIsoPtSums_);
  evt.getByToken(NeutralIsoPtSum_token, neutralIsoPtSums_);
  evt.getByToken(PUcorrPtSum_token, puCorrPtSums_);
  evt.getByToken(PhotonPtSumOutsideSignalCone_token, photonPtSumOutsideSignalCone_);
  evt.getByToken(FootprintCorrection_token, footprintCorrection_);
  
  evt.getByToken(Tau_token, taus_);
  category_output_.reset(new PFTauDiscriminator(TauRefProd(taus_)));
}

double PFRecoTauDiscriminationByMVAIsolationRun2::discriminate

(

const PFTauRef &

tau

)

const

Definition at line 188 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

References python.rootplot.argparse::category, PFRecoTauDiscriminationByMVAIsolation2_cff::chargedIsoPtSum, reco::PFTauTransverseImpactParameter::dxy(), reco::PFTauTransverseImpactParameter::dxy_Sig(), alignCSCRings::e, reco::PFTauTransverseImpactParameter::flightLength(), reco::PFTauTransverseImpactParameter::flightLengthSig(), PFRecoTauDiscriminationByMVAIsolationRun2_cff::footprintCorrection, reco::PFTauTransverseImpactParameter::hasSecondaryVertex(), reco::PFTauTransverseImpactParameter::ip3d(), reco::PFTauTransverseImpactParameter::ip3d_Sig(), edm::Ref< C, T, F >::isNull(), dqm-mbProfile::log, bookConverter::max, min(), PFRecoTauDiscriminationByMVAIsolation2_cff::neutralIsoPtSum, PFRecoTauDiscriminationByMVAIsolationRun2_cff::photonPtSumOutsideSignalCone, PFRecoTauDiscriminationByMVAIsolation2_cff::puCorrPtSum, Sign(), mathSSE::sqrt(), and metsig::tau.

{
  // CV: define dummy category index in order to use RecoTauDiscriminantCutMultiplexer module to appy WP cuts
  double category = 0.; 
  category_output_->setValue(tauIndex_, category);

  // CV: computation of MVA value requires presence of leading charged hadron
  if ( tau->leadPFChargedHadrCand().isNull() ) return 0.;

  int tauDecayMode = tau->decayMode();

  if ( ((mvaOpt_ == kOldDMwoLT || mvaOpt_ == kOldDMwLT || mvaOpt_ == kDBoldDMwLT || mvaOpt_ == kPWoldDMwLT) && (tauDecayMode == 0 || tauDecayMode == 1 || tauDecayMode == 2 || tauDecayMode == 10)) ||
       ((mvaOpt_ == kNewDMwoLT || mvaOpt_ == kNewDMwLT || mvaOpt_ == kDBnewDMwLT || mvaOpt_ == kPWnewDMwLT) && (tauDecayMode == 0 || tauDecayMode == 1 || tauDecayMode == 2 || tauDecayMode == 5 || tauDecayMode == 6 || tauDecayMode == 10)) ) {

    double chargedIsoPtSum = (*chargedIsoPtSums_)[tau];
    double neutralIsoPtSum = (*neutralIsoPtSums_)[tau];
    double puCorrPtSum     = (*puCorrPtSums_)[tau];
    double photonPtSumOutsideSignalCone = (*photonPtSumOutsideSignalCone_)[tau];
    double footprintCorrection = (*footprintCorrection_)[tau];
    
    const reco::PFTauTransverseImpactParameter& tauLifetimeInfo = *(*tauLifetimeInfos)[tau];
    
    double decayDistX = tauLifetimeInfo.flightLength().x();
    double decayDistY = tauLifetimeInfo.flightLength().y();
    double decayDistZ = tauLifetimeInfo.flightLength().z();
    double decayDistMag = std::sqrt(decayDistX*decayDistX + decayDistY*decayDistY + decayDistZ*decayDistZ);

    double nPhoton = double(tau_n_photons_total(*tau));
    double ptWeightedDetaStrip = tau_pt_weighted_deta_strip(*tau, tauDecayMode);
    double ptWeightedDphiStrip = tau_pt_weighted_dphi_strip(*tau, tauDecayMode);
    double ptWeightedDrSignal = tau_pt_weighted_dr_signal(*tau, tauDecayMode);
    double ptWeightedDrIsolation = tau_pt_weighted_dr_iso(*tau, tauDecayMode);
    double leadingTrackChi2 = tau_leadTrackChi2(*tau);
    double eRatio = tau_Eratio(*tau);

    if ( mvaOpt_ == kOldDMwoLT || mvaOpt_ == kNewDMwoLT ) {
      mvaInput_[0]  = std::log(std::max(1., Double_t(tau->pt())));
      mvaInput_[1]  = std::fabs(tau->eta());
      mvaInput_[2]  = std::log(std::max(1.e-2, chargedIsoPtSum));
      mvaInput_[3]  = std::log(std::max(1.e-2, neutralIsoPtSum - 0.125*puCorrPtSum));
      mvaInput_[4]  = std::log(std::max(1.e-2, puCorrPtSum));
      mvaInput_[5]  = tauDecayMode;
    } else if ( mvaOpt_ == kOldDMwLT || mvaOpt_ == kNewDMwLT  ) {
      mvaInput_[0]  = std::log(std::max(1., Double_t(tau->pt())));
      mvaInput_[1]  = std::fabs(tau->eta());
      mvaInput_[2]  = std::log(std::max(1.e-2, chargedIsoPtSum));
      mvaInput_[3]  = std::log(std::max(1.e-2, neutralIsoPtSum - 0.125*puCorrPtSum));
      mvaInput_[4]  = std::log(std::max(1.e-2, puCorrPtSum));
      mvaInput_[5]  = tauDecayMode;
      mvaInput_[6]  = TMath::Sign(+1., tauLifetimeInfo.dxy());
      mvaInput_[7]  = std::sqrt(std::fabs(std::min(1., tauLifetimeInfo.dxy())));
      mvaInput_[8]  = std::min(10., std::fabs(tauLifetimeInfo.dxy_Sig()));
      mvaInput_[9]  = ( tauLifetimeInfo.hasSecondaryVertex() ) ? 1. : 0.;
      mvaInput_[10] = std::sqrt(decayDistMag);
      mvaInput_[11] = std::min(10., tauLifetimeInfo.flightLengthSig());
    } else if ( mvaOpt_ == kDBoldDMwLT || mvaOpt_ == kDBnewDMwLT ) {
      mvaInput_[0]  = std::log(std::max(1., Double_t(tau->pt())));
      mvaInput_[1]  = std::fabs(tau->eta());
      mvaInput_[2]  = std::log(std::max(1.e-2, chargedIsoPtSum));
      mvaInput_[3]  = std::log(std::max(1.e-2, neutralIsoPtSum));
      mvaInput_[4]  = std::log(std::max(1.e-2, puCorrPtSum));
      mvaInput_[5]  = std::log(std::max(1.e-2, photonPtSumOutsideSignalCone));
      mvaInput_[6]  = tauDecayMode;
      mvaInput_[7]  = std::min(30., nPhoton);
      mvaInput_[8]  = std::min(0.5, ptWeightedDetaStrip);
      mvaInput_[9]  = std::min(0.5, ptWeightedDphiStrip);
      mvaInput_[10] = std::min(0.5, ptWeightedDrSignal);
      mvaInput_[11] = std::min(0.5, ptWeightedDrIsolation);
      mvaInput_[12] = std::min(100., leadingTrackChi2);
      mvaInput_[13] = std::min(1., eRatio);
      mvaInput_[14]  = TMath::Sign(+1., tauLifetimeInfo.dxy());
      mvaInput_[15]  = std::sqrt(std::fabs(std::min(1., std::fabs(tauLifetimeInfo.dxy()))));
      mvaInput_[16]  = std::min(10., std::fabs(tauLifetimeInfo.dxy_Sig()));
      mvaInput_[17]  = TMath::Sign(+1., tauLifetimeInfo.ip3d());
      mvaInput_[18]  = std::sqrt(std::fabs(std::min(1., std::fabs(tauLifetimeInfo.ip3d()))));
      mvaInput_[19]  = std::min(10., std::fabs(tauLifetimeInfo.ip3d_Sig()));
      mvaInput_[20]  = ( tauLifetimeInfo.hasSecondaryVertex() ) ? 1. : 0.;
      mvaInput_[21] = std::sqrt(decayDistMag);
      mvaInput_[22] = std::min(10., tauLifetimeInfo.flightLengthSig());
    } else if ( mvaOpt_ == kPWoldDMwLT || mvaOpt_ == kPWnewDMwLT ) {
      mvaInput_[0]  = std::log(std::max(1., Double_t(tau->pt())));
      mvaInput_[1]  = std::fabs(tau->eta());
      mvaInput_[2]  = std::log(std::max(1.e-2, chargedIsoPtSum));
      mvaInput_[3]  = std::log(std::max(1.e-2, neutralIsoPtSum));
      mvaInput_[4]  = std::log(std::max(1.e-2, footprintCorrection));
      mvaInput_[5]  = std::log(std::max(1.e-2, photonPtSumOutsideSignalCone));
      mvaInput_[6]  = tauDecayMode;
      mvaInput_[7]  = std::min(30., nPhoton);
      mvaInput_[8]  = std::min(0.5, ptWeightedDetaStrip);
      mvaInput_[9]  = std::min(0.5, ptWeightedDphiStrip);
      mvaInput_[10] = std::min(0.5, ptWeightedDrSignal);
      mvaInput_[11] = std::min(0.5, ptWeightedDrIsolation);
      mvaInput_[12] = std::min(100., leadingTrackChi2);
      mvaInput_[13] = std::min(1., eRatio);
      mvaInput_[14]  = TMath::Sign(+1., tauLifetimeInfo.dxy());
      mvaInput_[15]  = std::sqrt(std::fabs(std::min(1., std::fabs(tauLifetimeInfo.dxy()))));
      mvaInput_[16]  = std::min(10., std::fabs(tauLifetimeInfo.dxy_Sig()));
      mvaInput_[17]  = TMath::Sign(+1., tauLifetimeInfo.ip3d());
      mvaInput_[18]  = std::sqrt(std::fabs(std::min(1., std::fabs(tauLifetimeInfo.ip3d()))));
      mvaInput_[19]  = std::min(10., std::fabs(tauLifetimeInfo.ip3d_Sig()));
      mvaInput_[20]  = ( tauLifetimeInfo.hasSecondaryVertex() ) ? 1. : 0.;
      mvaInput_[21] = std::sqrt(decayDistMag);
      mvaInput_[22] = std::min(10., tauLifetimeInfo.flightLengthSig());
    }

    double mvaValue = mvaReader_->GetClassifier(mvaInput_);
    if ( verbosity_ ) {
      edm::LogPrint("PFTauDiscByMVAIsol2") << "<PFRecoTauDiscriminationByMVAIsolationRun2::discriminate>:";
      edm::LogPrint("PFTauDiscByMVAIsol2") << " tau: Pt = " << tau->pt() << ", eta = " << tau->eta();
      edm::LogPrint("PFTauDiscByMVAIsol2") << " isolation: charged = " << chargedIsoPtSum << ", neutral = " << neutralIsoPtSum << ", PUcorr = " << puCorrPtSum;
      edm::LogPrint("PFTauDiscByMVAIsol2") << " decay mode = " << tauDecayMode;
      edm::LogPrint("PFTauDiscByMVAIsol2") << " impact parameter: distance = " << tauLifetimeInfo.dxy() << ", significance = " << tauLifetimeInfo.dxy_Sig();
      edm::LogPrint("PFTauDiscByMVAIsol2") << " has decay vertex = " << tauLifetimeInfo.hasSecondaryVertex() << ":"
                       << " distance = " << decayDistMag << ", significance = " << tauLifetimeInfo.flightLengthSig();
      edm::LogPrint("PFTauDiscByMVAIsol2") << "--> mvaValue = " << mvaValue;
    }
    return mvaValue;
  } else {
    return -1.;
  }
}

void PFRecoTauDiscriminationByMVAIsolationRun2::endEvent ( edm::Event &  evt )

virtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 310 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

References edm::Event::put().

{
  // add all category indices to event
  evt.put(category_output_, "category");
}

Member Data Documentation

std::auto_ptr<PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::category_output_

private

Definition at line 159 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::EDGetTokenT<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::ChargedIsoPtSum_token

private

Definition at line 147 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::Handle<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::chargedIsoPtSums_

private

Definition at line 148 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::Handle<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::footprintCorrection_

private

Definition at line 156 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::EDGetTokenT<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::FootprintCorrection_token

private

Definition at line 155 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::FileInPath PFRecoTauDiscriminationByMVAIsolationRun2::inputFileName_

private

Definition at line 137 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

std::vector<TFile*> PFRecoTauDiscriminationByMVAIsolationRun2::inputFilesToDelete_

private

Definition at line 161 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

bool PFRecoTauDiscriminationByMVAIsolationRun2::loadMVAfromDB_

private

Definition at line 136 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

std::string PFRecoTauDiscriminationByMVAIsolationRun2::moduleLabel_

private

Definition at line 133 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

Referenced by Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), and Modules.ESProducer::nameInProcessDesc_().

float* PFRecoTauDiscriminationByMVAIsolationRun2::mvaInput_

private

Definition at line 141 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

std::string PFRecoTauDiscriminationByMVAIsolationRun2::mvaName_

private

Definition at line 135 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

int PFRecoTauDiscriminationByMVAIsolationRun2::mvaOpt_

private

Definition at line 140 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

const GBRForest* PFRecoTauDiscriminationByMVAIsolationRun2::mvaReader_

private

Definition at line 138 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::EDGetTokenT<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::NeutralIsoPtSum_token

private

Definition at line 149 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::Handle<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::neutralIsoPtSums_

private

Definition at line 150 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::Handle<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::photonPtSumOutsideSignalCone_

private

Definition at line 154 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::EDGetTokenT<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::PhotonPtSumOutsideSignalCone_token

private

Definition at line 153 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::EDGetTokenT<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::PUcorrPtSum_token

private

Definition at line 151 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::Handle<reco::PFTauDiscriminator> PFRecoTauDiscriminationByMVAIsolationRun2::puCorrPtSums_

private

Definition at line 152 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::Handle<PFTauTIPAssociationByRef> PFRecoTauDiscriminationByMVAIsolationRun2::tauLifetimeInfos

private

Definition at line 145 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::Handle<TauCollection> PFRecoTauDiscriminationByMVAIsolationRun2::taus_

private

Definition at line 158 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

edm::EDGetTokenT<PFTauTIPAssociationByRef> PFRecoTauDiscriminationByMVAIsolationRun2::TauTransverseImpactParameters_token

private

Definition at line 144 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.

int PFRecoTauDiscriminationByMVAIsolationRun2::verbosity_

private

Definition at line 163 of file PFRecoTauDiscriminationByMVAIsolationRun2.cc.