Class Reference for E1039 Core & Analysis Software
Common.h
Go to the documentation of this file.
1 
2 // 103: J/psi mass, BG = NIM3 e906 hits
3 // 113: J/psi mass, BG = full bkg,
4 // 119: 2.0-9.0 GeV, BG = full bkg, normal KMag pol
5 // 123: 2.5-3.5 GeV, BG = full bkg, reverse KMag pol
6 // 129: 2.0-9.0 GeV, BG = full bkg, reverse KMag pol
7 // 120: 2.0-9.0 GeV, BG = full bkg, normal KMag pol, H1X y-gap
8 // 130: 2.0-9.0 GeV, BG = full bkg, reverse KMag pol, H1X y-gap
9 // 131: 2.0-9.0 GeV, BG = full bkg, normal KMag pol, H1X y-gap, symmetric
10 // 132: 2.0-9.0 GeV, BG = full bkg, reverse KMag pol, H1X y-gap, symmetric
11 // 133: 2.0-9.0 GeV, BG = full bkg, normal KMag pol, H1X y-gap, symmetric, not hot-road removal
12 // 134: 2.0-9.0 GeV, BG = full bkg, reverse KMag pol, H1X y-gap, symmetric, not hot-road removal
13 // 135: 2.0-9.0 GeV, BG = full bkg, zero KMag, H1X y-gap, symmetric, not hot-road removal
14 // 136: 2.0-9.0 GeV, BG = full bkg, normal 50% KMag, H1X y-gap, symmetric
15 
16 void GetParams(string& rs_id, double& mass_lo, double& mass_hi, int& inte_cut, double& frac_cut, string& list_sig, string& list_bg)
17 {
18  rs_id = "136"; // Select a roadset here to analyze.
19  list_sig = "list_signal.txt";
20  list_bg = "list_bg.txt";
21 
22  if (rs_id == "102") {
23  mass_lo = 4.0;
24  mass_hi = 5.0;
25  inte_cut = 1200;
26  frac_cut = 0.010;
27  } else if (rs_id == "103") {
28  mass_lo = 2.5;
29  mass_hi = 3.5;
30  inte_cut = /*1200*/10.e3;
31  frac_cut = 0.040;
32  list_bg = "list_bg_NIM3e906.txt";
33  } else if (rs_id == "104") {
34  mass_lo = 2.5;
35  mass_hi = 3.5;
36  inte_cut = 1200;
37  frac_cut = 0.033;
38  } else if (rs_id == "105") {
39  mass_lo = 7.0;
40  mass_hi = 8.0;
41  inte_cut = 1200;
42  frac_cut = 0.010;
43  } else if (rs_id == "113") {
44  mass_lo = 2.5;
45  mass_hi = 3.5;
46  inte_cut = 42000;
47  frac_cut = 0.030;
48  list_bg = "list_bg_fullsimRun06.txt";
49  } else if (rs_id == "119") {
50  mass_lo = 2.0;
51  mass_hi = 9.0;
52  inte_cut = 42000;
53  frac_cut = 0.040;
54  list_bg = "list_bg_fullsimRun06.txt";
55  } else if (rs_id == "120") {
56  mass_lo = 2.0;
57  mass_hi = 9.0;
58  inte_cut = 42000;
59  frac_cut = 0.043;
60  list_bg = "list_bg_fullsimRun06.txt";
61  } else if (rs_id == "123") { //J/Psi roadset for reverse KMAG polarity
62  mass_lo = 2.5;
63  mass_hi = 3.5;
64  inte_cut = 43000;
65  frac_cut = 0.0925;
66  list_sig = "list_signal_reverseKMAG.txt";
67  list_bg = "list_bg_fullsimRun06_reverseKMAG.txt";
68  } else if (rs_id == "129") {
69  mass_lo = 2.0;
70  mass_hi = 9.0;
71  inte_cut = 42000;
72  frac_cut = 0.12;
73  list_sig = "list_signal_reverseKMAG.txt";
74  list_bg = "list_bg_fullsimRun06_reverseKMAG.txt";
75  } else if (rs_id == "130") {
76  mass_lo = 2.0;
77  mass_hi = 9.0;
78  inte_cut = 42000;
79  frac_cut = 0.10;
80  list_sig = "list_signal_reverseKMAG.txt";
81  list_bg = "list_bg_fullsimRun06_reverseKMAG.txt";
82  } else if (rs_id == "131") {
83  mass_lo = 2.0;
84  mass_hi = 9.0;
85  inte_cut = 42000;
86  frac_cut = 0.086; // 0.100 too larger, 0.080, too small?
87  list_bg = "list_bg_fullsimRun06.txt";
88  } else if (rs_id == "132") {
89  mass_lo = 2.0;
90  mass_hi = 9.0;
91  inte_cut = 42000;
92  frac_cut = 0.36;
93  list_sig = "list_signal_reverseKMAG.txt";
94  list_bg = "list_bg_fullsimRun06_reverseKMAG.txt";
95  } else if (rs_id == "133") {
96  mass_lo = 2.0;
97  mass_hi = 9.0;
98  inte_cut = 1000; // Does not matter
99  frac_cut = 0;
100  list_bg = "list_bg_fullsimRun06.txt";
101  } else if (rs_id == "134") {
102  mass_lo = 2.0;
103  mass_hi = 9.0;
104  inte_cut = 1000; // Does not matter
105  frac_cut = 0;
106  list_sig = "list_signal_reverseKMAG.txt";
107  list_bg = "list_bg_fullsimRun06_reverseKMAG.txt";
108  } else if (rs_id == "135") {
109  mass_lo = 2.0;
110  mass_hi = 9.0;
111  inte_cut = 0; // Does not matter
112  frac_cut = 0;
113  list_sig = "list_signal_zeroKMAG.txt";
114  //list_bg = "list_bg_fullsimRun06_reverseKMAG.txt"; // Does not matter
115  } else if (rs_id == "136") {
116  mass_lo = 2.0;
117  mass_hi = 9.0;
118  inte_cut = 80500; // with 45000
119  frac_cut = 0.140; // 0.140 too small, 0.141 too large
120  list_sig = "list_signal_KMagNormal050.txt";
121  list_bg = "list_bg_KMagNormal050.txt";
122  }
123 }
GetParams
void GetParams(string &rs_id, double &mass_lo, double &mass_hi, int &inte_cut, double &frac_cut, string &list_sig, string &list_bg)
Definition: Common.h:16