28 minMVA1prongBL_ = iConfig.
getParameter<
double>(
"minMVA1prongBL");
29 minMVA1prongStripsWgsfBL_ = iConfig.
getParameter<
double>(
"minMVA1prongStripsWgsfBL");
30 minMVA1prongStripsWOgsfBL_ = iConfig.
getParameter<
double>(
"minMVA1prongStripsWOgsfBL");
31 minMVA1prongEC_ = iConfig.
getParameter<
double>(
"minMVA1prongEC");
32 minMVA1prongStripsWgsfEC_ = iConfig.
getParameter<
double>(
"minMVA1prongStripsWgsfEC");
33 minMVA1prongStripsWOgsfEC_ = iConfig.
getParameter<
double>(
"minMVA1prongStripsWOgsfEC");
36 mva_->Initialize(method_,
37 inputFileName1prongBL_.fullPath().data(),
38 inputFileName1prongStripsWgsfBL_.
fullPath().data(),
39 inputFileName1prongStripsWOgsfBL_.
fullPath().data(),
40 inputFileName1prongEC_.
fullPath().data(),
41 inputFileName1prongStripsWgsfEC_.
fullPath().data(),
42 inputFileName1prongStripsWOgsfEC_.
fullPath().data()
48 double discriminate(
const PFTauRef& thePFTauRef);
74 double workingPoint = 0.0;
76 if( (*thePFTauRef).leadPFChargedHadrCand().isNonnull() ) {
78 mva = mva_->MVAValue( thePFTauRef );
81 ((*thePFTauRef).signalPFChargedHadrCands().size()==3 ||
82 (fabs((*thePFTauRef).eta())<1.5 && ((*thePFTauRef).signalPFGammaCands().size())<=0 && mva > minMVA1prongBL_) ||
83 (fabs((*thePFTauRef).eta())<1.5 && ((*thePFTauRef).signalPFGammaCands().size())>0 && (((*thePFTauRef).leadPFChargedHadrCand())->gsfTrackRef()).isNonnull()>0.5 && mva > minMVA1prongStripsWgsfBL_) ||
84 (fabs((*thePFTauRef).eta())<1.5 && ((*thePFTauRef).signalPFGammaCands().size())>0 && (((*thePFTauRef).leadPFChargedHadrCand())->gsfTrackRef()).isNonnull()<0.5 && mva > minMVA1prongStripsWOgsfBL_) ||
85 (fabs((*thePFTauRef).eta())>1.5 && ((*thePFTauRef).signalPFGammaCands().size())<=0 && mva > minMVA1prongEC_) ||
86 (fabs((*thePFTauRef).eta())>1.5 && ((*thePFTauRef).signalPFGammaCands().size())>0 && (((*thePFTauRef).leadPFChargedHadrCand())->gsfTrackRef()).isNonnull()>0.5 && mva > minMVA1prongStripsWgsfEC_) ||
87 (fabs((*thePFTauRef).eta())>1.5 && ((*thePFTauRef).signalPFGammaCands().size())>0 && (((*thePFTauRef).leadPFChargedHadrCand())->gsfTrackRef()).isNonnull()<0.5 && mva > minMVA1prongStripsWOgsfEC_)) ? 1.0 : 0.0;
91 return ( returnMVA_ ? mva : workingPoint);
T getParameter(std::string const &) const
PFRecoTauDiscriminationAgainstElectronMVA(const edm::ParameterSet &iConfig)
double minMVA1prongStripsWgsfBL_
~PFRecoTauDiscriminationAgainstElectronMVA()
double discriminate(const PFTauRef &thePFTauRef)
edm::FileInPath inputFileName1prongStripsWOgsfBL_
edm::FileInPath inputFileName1prongBL_
edm::FileInPath inputFileName1prongEC_
double minMVA1prongStripsWOgsfBL_
DEFINE_FWK_MODULE(CosmicTrackingParticleSelector)
edm::FileInPath inputFileName1prongStripsWgsfBL_
edm::FileInPath inputFileName1prongStripsWgsfEC_
double minMVA1prongStripsWgsfEC_
double minMVA1prongStripsWOgsfEC_
std::string fullPath() const
edm::FileInPath inputFileName1prongStripsWOgsfEC_