|
|
|
00001 // $Id: TtHadLRSignalSelCalc.cc,v 1.2 2008/02/17 11:18:10 rwolf Exp $ 00002 // copied TtSemiLRSignalSelCalc.cc,v 1.2 2007/06/18 14:12:18 heyninck Exp 00003 // for fully hadronic channel 00004 00005 #include "TopQuarkAnalysis/TopEventSelection/interface/TtHadLRSignalSelCalc.h" 00006 00007 // constructor with path; default should not be used 00008 TtHadLRSignalSelCalc::TtHadLRSignalSelCalc(TString fitInputPath, std::vector<int> observables) 00009 { 00010 std::cout << "=== Constructing a TtHadLRSignalSelCalc... " << std::endl; 00011 myLR = new LRHelpFunctions(); 00012 addPurity = false; 00013 if(observables[0] == -1) addPurity = true; 00014 myLR -> readObsHistsAndFits(fitInputPath, observables, addPurity); 00015 std::cout << "=== done." << std::endl; 00016 } 00017 00018 TtHadLRSignalSelCalc::~TtHadLRSignalSelCalc() 00019 { 00020 delete myLR; 00021 } 00022 00023 void TtHadLRSignalSelCalc::operator()(TtHadEvtSolution & sol) 00024 { 00025 // find the used observables 00026 std::vector<double> obsVals; 00027 for(unsigned int o = 0; o<100; o++){ 00028 if( myLR->obsFitIncluded(o) ) obsVals.push_back(sol.getLRSignalEvtObsVal(o)); 00029 } 00030 00031 // calculate the logLR and the purity 00032 double logLR = myLR->calcLRval(obsVals); 00033 double prob = -999.; 00034 if(addPurity) prob = myLR->calcProb(logLR); 00035 00036 // fill these values to the members in the TtHadEvtSolution 00037 sol.setLRSignalEvtLRval(logLR); 00038 sol.setLRSignalEvtProb(prob); 00039 }
1.7.3