-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcompRDRD.h
148 lines (122 loc) · 2.92 KB
/
compRDRD.h
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
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
#include <vector>
#include <string>
#include <set>
#include <map>
#include <utility>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <TFile.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TCanvas.h>
#include <TMatrixD.h>
#include <TMatrixTUtils.h>
#include <TTree.h>
#include <TBranch.h>
#include <TLeaf.h>
#include <TLine.h>
#include <TGaxis.h>
#include <TLatex.h>
// COLORS
#define RST "\x1B[0m"
#define KRED "\x1B[31m"
#define KGRN "\x1B[32m"
#define KYEL "\x1B[33m"
#define KBLU "\x1B[34m"
#define KMAG "\x1B[35m"
#define KCYN "\x1B[36m"
#define KWHT "\x1B[37m"
#define FRED(x) KRED x RST
#define FGRN(x) KGRN x RST
#define FYEL(x) KYEL x RST
#define FBLU(x) KBLU x RST
#define FMAG(x) KMAG x RST
#define FCYN(x) KCYN x RST
#define FWHT(x) KWHT x RST
#define BOLD(x) "\x1B[1m" x RST
#define UNDL(x) "\x1B[4m" x RST
using namespace std;
//Qty containers
vector<Double_t> fQ2;
vector<Double_t> fQ2kin[5];
vector<Double_t> fQ2local;
vector<Double_t> fXBj;
vector<Double_t> fXBjkin[5];
vector<Double_t> fYBj;
vector<Double_t> fYBjkin[5];
vector<Double_t> fWBj;
vector<Double_t> fWBjkin[5];
vector<Double_t> fNu;
vector<Double_t> fNukin[5];
vector<Double_t> fMu[5];
vector<Double_t> fMup[5];
vector<Double_t> fThetaMu[3];
vector<Double_t> fTheta[5];
vector<Double_t> fPhi[5];
vector<Double_t> fVertex[5];
vector<Double_t> fQ2RD2;
vector<Double_t> fQ2kinRD2[5];
vector<Double_t> fXBjRD2;
vector<Double_t> fXBjkinRD2[5];
vector<Double_t> fYBjRD2;
vector<Double_t> fYBjkinRD2[5];
vector<Double_t> fWBjRD2;
vector<Double_t> fWBjkinRD2[5];
vector<Double_t> fNuRD2;
vector<Double_t> fNukinRD2[5];
vector<Double_t> fMuRD2[5];
vector<Double_t> fMupRD2[5];
vector<Double_t> fThetaRD2Mu[3];
vector<Double_t> fThetaRD2[5];
vector<Double_t> fPhiRD2[5];
vector<Double_t> fVertexRD2[5];
vector<Double_t> fX;
vector<Double_t> fY;
vector<Double_t> fZ;
vector<Double_t> fXRD2;
vector<Double_t> fYRD2;
vector<Double_t> fZRD2;
//Misc
set<Double_t> fLHsec_set;
Double_t* fLHsec_tab;
Double_t fLHsec;
int fId;
int fId_loose;
int fId_severe;
Double_t fNu_max[3][12];
Double_t fNu_min[3][12];
vector<Double_t> fXv;
vector<Double_t> fYv;
vector<Double_t> fZv;
vector<Double_t> fRv;
bool fAllDISflag;
bool fAllDISflag_RD2;
int xbin, ybin, zbin;
// Cuts
Double_t fXmin;
Double_t fXmax;
Double_t fYmin;
Double_t fYmax;
Double_t fWmin;
Double_t fWmax;
Double_t fPmin;
Double_t fPmax;
//Draw
TH1F* fKinematicsRD[5][13];
TH1F* fKinematicsRD2[5][13];
TH1F* fKinematicsRatio[5][13];
Int_t fCountingRD[5][13];
Int_t fCountingRD2[5][13];
vector<double> fError, fErrorRD, fErrorRD2;
Int_t fLineStyle[7] = {3,3,3,1,3,3,3};
TLine* l1[13][7];
//Graphic Style
Int_t fMarkerColor[2] = {4,2};
Int_t fMarkerStyle[2] = {24,20};
//Constants
static const Double_t fM_p = 938.272046/(1e3);
static const Double_t fM_mu = 105.6583715/(1e3);
static const Double_t fM_K = 493.677/(1e3);
static const Double_t fM_pi = 139.57018/(1e3);
string trigname[5] = {"MT","LT","OT","LAST",""};