1 |
anlevin |
1.2 |
#ifndef WWAnalysis_AnalysisStep_RochCor2012_h
|
2 |
|
|
#define WWAnalysis_AnalysisStep_RochCor2012_h
|
3 |
|
|
|
4 |
|
|
//// VERSION for 2012 received from Jiyeon on 30 september 2012
|
5 |
|
|
//// moved static const float from .h to .cc to make the gcc happy
|
6 |
|
|
|
7 |
sixie |
1.1 |
#include <iostream>
|
8 |
anlevin |
1.2 |
|
9 |
sixie |
1.1 |
#include <TChain.h>
|
10 |
|
|
#include <TClonesArray.h>
|
11 |
|
|
#include <TString.h>
|
12 |
|
|
#include <map>
|
13 |
|
|
|
14 |
|
|
#include <TSystem.h>
|
15 |
|
|
#include <TROOT.h>
|
16 |
|
|
#include <TMath.h>
|
17 |
|
|
#include <TLorentzVector.h>
|
18 |
|
|
#include <TRandom3.h>
|
19 |
|
|
|
20 |
|
|
|
21 |
|
|
class rochcor2012 {
|
22 |
|
|
public:
|
23 |
|
|
rochcor2012();
|
24 |
|
|
rochcor2012(int seed);
|
25 |
|
|
~rochcor2012();
|
26 |
|
|
|
27 |
anlevin |
1.3 |
void momcor_mc(TLorentzVector&, float, float, int, bool);
|
28 |
anlevin |
1.2 |
void momcor_data(TLorentzVector&, float, float, int);
|
29 |
sixie |
1.1 |
|
30 |
|
|
void musclefit_data(TLorentzVector& , TLorentzVector&);
|
31 |
anlevin |
1.2 |
|
32 |
sixie |
1.1 |
float zptcor(float);
|
33 |
|
|
int etabin(float);
|
34 |
|
|
int phibin(float);
|
35 |
|
|
|
36 |
|
|
private:
|
37 |
|
|
|
38 |
|
|
TRandom3 eran;
|
39 |
|
|
TRandom3 sran;
|
40 |
|
|
|
41 |
anlevin |
1.2 |
// static float netabin[9] = {-2.4,-2.1,-1.4,-0.7,0.0,0.7,1.4,2.1,2.4};
|
42 |
|
|
static const float netabin[9];
|
43 |
|
|
////^^^^^------------ GP BEGIN
|
44 |
|
|
static const double pi;
|
45 |
|
|
|
46 |
|
|
static const float genm_smr;
|
47 |
|
|
static const float genm;
|
48 |
|
|
|
49 |
|
|
static const float mrecm;
|
50 |
|
|
static const float drecm;
|
51 |
|
|
static const float mgscl_stat;
|
52 |
|
|
static const float mgscl_syst;
|
53 |
|
|
static const float dgscl_stat;
|
54 |
|
|
static const float dgscl_syst;
|
55 |
|
|
|
56 |
|
|
//iteration2 after FSR : after Z Pt correction
|
57 |
|
|
static const float delta;
|
58 |
|
|
static const float delta_stat;
|
59 |
|
|
static const float delta_syst;
|
60 |
|
|
|
61 |
|
|
static const float sf;
|
62 |
|
|
static const float sf_stat;
|
63 |
|
|
static const float sf_syst;
|
64 |
|
|
|
65 |
|
|
static const float apar;
|
66 |
|
|
static const float bpar;
|
67 |
|
|
static const float cpar;
|
68 |
|
|
static const float d0par;
|
69 |
|
|
static const float e0par;
|
70 |
|
|
static const float d1par;
|
71 |
|
|
static const float e1par;
|
72 |
|
|
static const float d2par;
|
73 |
|
|
static const float e2par;
|
74 |
|
|
////^^^^^------------ GP END
|
75 |
|
|
|
76 |
sixie |
1.1 |
//---------------------------------------------------------------------------------------------
|
77 |
|
|
|
78 |
|
|
static const float dcor_bf[8][8];
|
79 |
|
|
static const float dcor_ma[8][8];
|
80 |
|
|
static const float mcor_bf[8][8];
|
81 |
|
|
static const float mcor_ma[8][8];
|
82 |
|
|
static const float dcor_bfer[8][8];
|
83 |
|
|
static const float dcor_maer[8][8];
|
84 |
|
|
static const float mcor_bfer[8][8];
|
85 |
|
|
static const float mcor_maer[8][8];
|
86 |
|
|
|
87 |
|
|
//=======================================================================================================
|
88 |
|
|
|
89 |
|
|
static const float dmavg[8][8];
|
90 |
|
|
static const float dpavg[8][8];
|
91 |
|
|
static const float mmavg[8][8];
|
92 |
|
|
static const float mpavg[8][8];
|
93 |
|
|
|
94 |
|
|
//===============================================================================================
|
95 |
|
|
//parameters for Z pt correction
|
96 |
|
|
static const int nptbins=84;
|
97 |
|
|
static const float ptlow[85];
|
98 |
|
|
|
99 |
|
|
static const float zptscl[84];
|
100 |
|
|
static const float zptscler[84];
|
101 |
|
|
|
102 |
|
|
float mptsys_mc_dm[8][8];
|
103 |
|
|
float mptsys_mc_da[8][8];
|
104 |
|
|
float mptsys_da_dm[8][8];
|
105 |
|
|
float mptsys_da_da[8][8];
|
106 |
|
|
|
107 |
|
|
float gscler_mc_dev;
|
108 |
|
|
float gscler_da_dev;
|
109 |
|
|
};
|
110 |
anlevin |
1.2 |
|
111 |
sixie |
1.1 |
#endif
|