| 1 |
mzanetti |
1.1 |
// ------------------------------------------------------------
|
| 2 |
|
|
//
|
| 3 |
|
|
// QGLikelihoodCalculator - Class
|
| 4 |
|
|
// for the computation of the QG likelihood.
|
| 5 |
|
|
// Needs files provided by having run the
|
| 6 |
|
|
// Ntp1Finalizer_QG on QCD samples.
|
| 7 |
|
|
//
|
| 8 |
|
|
// ------------------------------------------------------------
|
| 9 |
|
|
|
| 10 |
|
|
#ifndef QGLikelihoodCalculator_h
|
| 11 |
|
|
#define QGLikelihoodCalculator_h
|
| 12 |
|
|
|
| 13 |
|
|
#include <string>
|
| 14 |
|
|
|
| 15 |
|
|
#include "TFile.h"
|
| 16 |
|
|
#include "TH1F.h"
|
| 17 |
|
|
#include <map>
|
| 18 |
|
|
|
| 19 |
|
|
|
| 20 |
|
|
|
| 21 |
|
|
class QGLikelihoodCalculator {
|
| 22 |
|
|
|
| 23 |
|
|
public:
|
| 24 |
|
|
// QGLikelihoodCalculator( const std::string& fileName="QG_QCD_Pt-15to3000_TuneZ2_Flat_7TeV_pythia6_Summer11-PU_S3_START42_V11-v2.root", unsigned nPtBins=21, unsigned int nRhoBins=25 );
|
| 25 |
|
|
QGLikelihoodCalculator( const TString dataDir, Bool_t chs = false);
|
| 26 |
|
|
~QGLikelihoodCalculator();
|
| 27 |
|
|
|
| 28 |
|
|
float computeQGLikelihoodPU( float pt, float rhoPF, int nCharged, int nNeutral, float ptD, float rmsCand=-1. );
|
| 29 |
|
|
float computeQGLikelihood2012( float pt, float eta, float rho, int nPFCandidates_QC, float ptD_QC, float axis2_QC ); //new
|
| 30 |
|
|
Float_t QGvalue(std::map<TString, Float_t>);
|
| 31 |
|
|
|
| 32 |
|
|
float likelihoodProduct( float nCharged, float nNeutral, float ptD, float rmsCand, TH1F* h1_nCharged, TH1F* h1_nNeutral, TH1F* h1_ptD, TH1F* h1_rmsCand);
|
| 33 |
|
|
|
| 34 |
|
|
|
| 35 |
|
|
|
| 36 |
|
|
private:
|
| 37 |
|
|
|
| 38 |
|
|
TFile* histoFile_;
|
| 39 |
|
|
std::map<std::string,TH1F*> plots_;
|
| 40 |
|
|
unsigned int nPtBins_;
|
| 41 |
|
|
unsigned int nRhoBins_;
|
| 42 |
|
|
|
| 43 |
|
|
};
|
| 44 |
|
|
|
| 45 |
|
|
|
| 46 |
|
|
#endif
|
