#include <MvaSoftMuonEstimator.h>

Public Member Functions
	MvaSoftMuonEstimator (std::string)

float	mvaValue (float, float, float, float, float)

	~MvaSoftMuonEstimator ()

Private Attributes
float	mva_deltaR

float	mva_ptRel

float	mva_ratio

float	mva_sip2d

float	mva_sip3d

TMVA::Reader *	TMVAReader

std::string	weightFile

Detailed Description

Definition at line 16 of file MvaSoftMuonEstimator.h.

Constructor & Destructor Documentation

MvaSoftMuonEstimator::MvaSoftMuonEstimator ( std::string weightFile )

Definition at line 12 of file MvaSoftMuonEstimator.cc.

References reco::details::loadTMVAWeights(), mva_deltaR, mva_ptRel, mva_ratio, mva_sip2d, mva_sip3d, and TMVAReader.

                                                                {
   TMVAReader = new TMVA::Reader("Color:Silent:Error");
   TMVAReader->SetVerbose(false);
   TMVAReader->AddVariable("TagInfo1.sip3d", &mva_sip3d);
   TMVAReader->AddVariable("TagInfo1.sip2d", &mva_sip2d);
   TMVAReader->AddVariable("TagInfo1.ptRel", &mva_ptRel);
   TMVAReader->AddVariable("TagInfo1.deltaR", &mva_deltaR);
   TMVAReader->AddVariable("TagInfo1.ratio", &mva_ratio);
   reco::details::loadTMVAWeights(TMVAReader, "BDT", weightFile.c_str()); 
 }

MvaSoftMuonEstimator::~MvaSoftMuonEstimator ( )

Definition at line 23 of file MvaSoftMuonEstimator.cc.

References TMVAReader.

                                             {
   delete TMVAReader;
 }

Member Function Documentation

float MvaSoftMuonEstimator::mvaValue	(	float	sip3d,
		float	sip2d,
		float	ptRel,
		float	deltaR,
		float	ratio
	)

Definition at line 29 of file MvaSoftMuonEstimator.cc.

References deltaR(), mva_deltaR, mva_ptRel, mva_ratio, mva_sip2d, mva_sip3d, GlobalPosition_Frontier_DevDB_cff::tag, and TMVAReader.

                                                                                                      {
   mva_sip3d = sip3d;
   mva_sip2d = sip2d;
   mva_ptRel = ptRel;
   mva_deltaR = deltaR;
   mva_ratio = ratio;
   // Evaluate tagger
   float tag = TMVAReader->EvaluateMVA("BDT");
   // Transform output between approximately 0 and 1
   tag = (tag+1.)/2.;
   
   return tag;
 }