DPGenStudy/Systematic.h at master · sckao/DPGenStudy · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
#ifndef Systematic_H
#define Systematic_H

#include "TObject.h"
#include <vector>
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <string>
#include <TMath.h>
#include <TF1.h>
#include <TH1.h>
#include <TH2.h>
#include <THStack.h>
#include <TFile.h>
#include <TTree.h>
#include <TBranch.h>
#include <TRandom3.h>
#include <TLeaf.h>
#include <TString.h>
#include <TSystem.h>
#include <TStyle.h>
#include <TFrame.h>
#include "TLorentzVector.h"

#include <algorithm>

#include "AnaInput.h"
#include "DPSelection.h"
#include "hDraw.h"

/*
#define MAXPHO 10
#define MAXVTX 10
#define MAXMU 5
*/

class Systematic : public TObject {

public:

   Systematic( string datacardfile = "DataCard.txt");
   ~Systematic();

   void McDataCompare() ;

   void RunData( string dataName, int nSample = -1 ) ;
   void RunMC( string mcName, int nSample = -1, double weight = 1 ) ;
   void ComparePlot( TH1D* hData, TH1D* hMC, string plotName );

   void ShiftPeak( TH1D* h1, TH1D* h2 ) ;
   void WriteDataHisto() ;
   void WriteMcHisto() ;

private:

   AnaInput*     Input;
   DPSelection*  select;
   hDraw*        h_draw;

   string hfolder  ;
   string plotType ;
   int    ProcessEvents ;
   int    SkipEvents ;
   // Root File for Rootuple or histograms
   //TFile *theFile ;
   //string hfName ;

   vector<double> timeCalib ;
   vector<double> normV ;

   vector<objID> selectJets ;
   vector<objID> selectPho ;

   TH1D* h_dataTime  ;
   TH1D* h_dataTimeA  ;
   TH1D* h_dataMET       ;
   TH1D* h_mcTime  ;
   TH1D* h_mcTimeC  ;
   TH1D* h_mcTimeA  ;
   TH1D* h_mcMET   ;

   //ClassDef(Systematic, 1);
};

//#if !defined(__CINT__)
//    ClassImp(Systematic);
#endif