-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsettings.h
executable file
·93 lines (72 loc) · 3.8 KB
/
settings.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
#ifndef settings_h
#define settings_h
#include <TMath.h>
#include <TLine.h>
//char hist_mass_name[50];
double PI = TMath::Pi();
//==============================================================================
// VARIABLES pT
//==============================================================================
const int N_pt_bin = 3;
int min_pt_bin[N_pt_bin] = {0,2,6};
int max_pt_bin[N_pt_bin] = {2,6,12};
//==============================================================================
// VARIABLES cost
//==============================================================================
//------------------------------------------------------------------------------
// RANGE -PI < phi < PI ; -1 < cost < 1
//------------------------------------------------------------------------------
//const int N_cost_bins = 18;
const int N_cost_bins = 22; // binning 1
const int N_TH3_cost_bins = 100;
//int min_cost_bin[N_cost_bins] = {1,26,31,36,41,43,45,47,49,51,53,55,57,59,61,66,71,76};
//int max_cost_bin[N_cost_bins] = {25,30,35,40,42,44,46,48,50,52,54,56,58,60,65,70,75,100};
//int min_cost_bin[N_cost_bins] = {1,19,26,31,36,41,43,45,47,49,51,53,55,57,59,61,66,71,76,81};
//int max_cost_bin[N_cost_bins] = {18,25,30,35,40,42,44,46,48,50,52,54,56,58,60,65,70,75,82,100};
int min_cost_bin[N_cost_bins] = {1,11,21,26,31,36,41,43,45,47,49,51,53,55,57,59,61,66,71,76,81,91}; // binning 1
int max_cost_bin[N_cost_bins] = {10,20,25,30,35,40,42,44,46,48,50,52,54,56,58,60,65,70,75,80,90,100}; // binning 1
const int dim_cost = N_cost_bins + 1;
//double bins_cost[dim_cost] = {0,25,30,35,40,42,44,46,48,50,52,54,56,58,60,65,70,75,100};
//double bins_cost[dim_cost] = {0,18,25,30,35,40,42,44,46,48,50,52,54,56,58,60,65,70,75,82,100};
double bins_cost[dim_cost] = {0,10,20,25,30,35,40,42,44,46,48,50,52,54,56,58,60,65,70,75,80,90,100}; // binning 1
double width_cost[dim_cost];
double min_cost = -1.;
double max_cost = 1.;
double value_cost[dim_cost];
double Dcost[N_cost_bins];
//==============================================================================
// VARIABLES phi
//==============================================================================
//------------------------------------------------------------------------------
// RANGE -PI < phi < PI ; -1 < cost < 1
//------------------------------------------------------------------------------
const int N_phi_bins = 10;
const int N_TH3_phi_bins = 50;
//const int N_TH3_phi_bins = 100; // double binning
//int min_phi_bin[N_phi_bins] = {1,16,19,21,23,26,29,31,33,36}; // binning 1
//int max_phi_bin[N_phi_bins] = {15,18,20,22,25,28,30,32,35,50}; // binning 1
int min_phi_bin[N_phi_bins] = {1,9,17,21,24,26,28,31,35,43}; // binning 2
int max_phi_bin[N_phi_bins] = {8,16,20,23,25,27,30,34,42,50}; // binning 2
//int min_phi_bin[N_phi_bins] = {1,31,37,41,45,51,57,61,65,71}; // double binning
//int max_phi_bin[N_phi_bins] = {30,36,40,44,50,56,60,64,70,100}; // double binning
const int dim_phi = N_phi_bins + 1;
//double bins_phi[dim_phi] = {0,15,18,20,22,25,28,30,32,35,50}; // binning 1
double bins_phi[dim_phi] = {0,8,16,20,23,25,27,30,34,42,50}; // binning 2
//double bins_phi[dim_phi] = {0,30,36,40,44,50,56,60,64,70,100}; // double binning
double width_phi[dim_phi];
double min_phi = 0;
double max_phi = PI;
double value_phi[dim_phi];
double Dphi[N_phi_bins];
//==============================================================================
// COMPUTE THE BIN AREA
//==============================================================================
double bin_area[N_cost_bins][N_phi_bins];
//==============================================================================
// PLOT THE LINES
//==============================================================================
const int N_line_cost = N_cost_bins - 1;
TLine *line_cost[N_line_cost];
const int N_line_phi = N_phi_bins - 1;
TLine *line_phi[N_line_phi];
#endif