BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
KsKpi.h
Go to the documentation of this file.
1#ifndef Physics_Analysis_KsKpi_H
2#define Physics_Analysis_KsKpi_H
3
4#include "GaudiKernel/Algorithm.h"
5#include "GaudiKernel/ITHistSvc.h"
6#include "GaudiKernel/NTuple.h"
7#include "TH1.h"
8
9class KsKpi : public Algorithm {
10
11public:
12 KsKpi( const std::string& name, ISvcLocator* pSvcLocator );
13 StatusCode initialize();
14 StatusCode execute();
15 StatusCode finalize();
16
17private:
18 ITHistSvc* m_thistsvc;
19
20 bool m_tagKsKpi;
21
22 TH1F* m_kskpi_vx_pi1;
23 TH1F* m_kskpi_vy_pi1;
24 TH1F* m_kskpi_vz_pi1;
25 TH1F* m_kskpi_vr_pi1;
26 TH1F* m_kskpi_px_pi1;
27 TH1F* m_kskpi_py_pi1;
28 TH1F* m_kskpi_pz_pi1;
29 TH1F* m_kskpi_pp_pi1;
30 TH1F* m_kskpi_e_pi1;
31 TH1F* m_kskpi_cos_pi1;
32 TH1F* m_kskpi_emc_pi1;
33
34 TH1F* m_kskpi_vx_pi2;
35 TH1F* m_kskpi_vy_pi2;
36 TH1F* m_kskpi_vz_pi2;
37 TH1F* m_kskpi_vr_pi2;
38 TH1F* m_kskpi_px_pi2;
39 TH1F* m_kskpi_py_pi2;
40 TH1F* m_kskpi_pz_pi2;
41 TH1F* m_kskpi_pp_pi2;
42 TH1F* m_kskpi_e_pi2;
43 TH1F* m_kskpi_cos_pi2;
44 TH1F* m_kskpi_emc_pi2;
45
46 TH1F* m_kskpi_vx_pi;
47 TH1F* m_kskpi_vy_pi;
48 TH1F* m_kskpi_vz_pi;
49 TH1F* m_kskpi_vr_pi;
50 TH1F* m_kskpi_px_pi;
51 TH1F* m_kskpi_py_pi;
52 TH1F* m_kskpi_pz_pi;
53 TH1F* m_kskpi_pp_pi;
54 TH1F* m_kskpi_e_pi;
55 TH1F* m_kskpi_cos_pi;
56 TH1F* m_kskpi_emc_pi;
57
58 TH1F* m_kskpi_vx_k;
59 TH1F* m_kskpi_vy_k;
60 TH1F* m_kskpi_vz_k;
61 TH1F* m_kskpi_vr_k;
62 TH1F* m_kskpi_px_k;
63 TH1F* m_kskpi_py_k;
64 TH1F* m_kskpi_pz_k;
65 TH1F* m_kskpi_pp_k;
66 TH1F* m_kskpi_e_k;
67 TH1F* m_kskpi_cos_k;
68 TH1F* m_kskpi_emc_k;
69
70 TH1F* m_kskpi_pidchidedx_1;
71 TH1F* m_kskpi_pidchitof1_1;
72 TH1F* m_kskpi_pidchitof2_1;
73 TH1F* m_kskpi_pidchidedx_2;
74 TH1F* m_kskpi_pidchitof1_2;
75 TH1F* m_kskpi_pidchitof2_2;
76 TH1F* m_kskpi_pidchidedx_3;
77 TH1F* m_kskpi_pidchitof1_3;
78 TH1F* m_kskpi_pidchitof2_3;
79 TH1F* m_kskpi_pidchidedx_4;
80 TH1F* m_kskpi_pidchitof1_4;
81 TH1F* m_kskpi_pidchitof2_4;
82
83 TH1F* m_kskpi_vfitp_chi;
84 TH1F* m_kskpi_vfitp_vx;
85 TH1F* m_kskpi_vfitp_vy;
86 TH1F* m_kskpi_vfitp_vz;
87 TH1F* m_kskpi_vfitp_vr;
88
89 TH1F* m_kskpi_vfits_chi;
90 TH1F* m_kskpi_vfits_vx;
91 TH1F* m_kskpi_vfits_vy;
92 TH1F* m_kskpi_vfits_vz;
93 TH1F* m_kskpi_vfits_vr;
94
95 TH1F* m_kskpi_vfit2_chi;
96 TH1F* m_kskpi_vfit2_mks;
97 TH1F* m_kskpi_vfit2_ct;
98 TH1F* m_kskpi_vfit2_dl;
99 TH1F* m_kskpi_vfit2_dle;
100
101 TH1F* m_kskpi_4c_chi;
102 TH1F* m_kskpi_4c_mks;
103 TH1F* m_kskpi_4c_mksk;
104 TH1F* m_kskpi_4c_mkspi;
105 TH1F* m_kskpi_4c_mkpi;
106 TH1F* m_kskpi_4c_ks_px;
107 TH1F* m_kskpi_4c_ks_py;
108 TH1F* m_kskpi_4c_ks_pz;
109 TH1F* m_kskpi_4c_ks_p;
110 TH1F* m_kskpi_4c_ks_cos;
111
112 // Declare r0, z0 cut for charged tracks
113 double m_vr0cut;
114 double m_vz0cut;
115 double m_coscut;
116 bool m_useKalTrk;
117
118 // Declare energy, dphi, dthe cuts for fake gamma's
119 double m_energyThreshold;
120 double m_gammaPhiCut;
121 double m_gammaThetaCut;
122
123 //
124 int m_test4C;
125 int m_test5C;
126
127 //
128 int m_checkDedx;
129 int m_checkTof;
130
131 // define Ntuples here
132
133 NTuple::Tuple* m_tuple1;
134
135 NTuple::Item<long> m_run;
136 NTuple::Item<long> m_event;
137 NTuple::Item<long> m_nchrg;
138 NTuple::Item<long> m_nneu;
139 NTuple::Item<long> m_ngch;
140
141 NTuple::Array<long> m_charge;
142 NTuple::Array<double> m_vx0;
143 NTuple::Array<double> m_vy0;
144 NTuple::Array<double> m_vz0;
145 NTuple::Array<double> m_vr0;
146
147 NTuple::Array<double> m_vx;
148 NTuple::Array<double> m_vy;
149 NTuple::Array<double> m_vz;
150 NTuple::Array<double> m_vr;
151
152 NTuple::Array<double> m_px;
153 NTuple::Array<double> m_py;
154 NTuple::Array<double> m_pz;
155 NTuple::Array<double> m_p;
156 NTuple::Array<double> m_cost;
157
158 NTuple::Array<float> m_probPH;
159 NTuple::Array<float> m_normPH;
160 NTuple::Array<float> m_chie;
161 NTuple::Array<float> m_chimu;
162 NTuple::Array<float> m_chipi;
163 NTuple::Array<float> m_chik;
164 NTuple::Array<float> m_chip;
165 NTuple::Array<float> m_ghit;
166 NTuple::Array<float> m_thit;
167
168 NTuple::Array<float> m_e_emc;
169
170 NTuple::Array<double> m_qual_etof;
171 NTuple::Array<double> m_tof_etof;
172 NTuple::Array<double> m_te_etof;
173 NTuple::Array<double> m_tmu_etof;
174 NTuple::Array<double> m_tpi_etof;
175 NTuple::Array<double> m_tk_etof;
176 NTuple::Array<double> m_tp_etof;
177
178 NTuple::Array<double> m_qual_btof1;
179 NTuple::Array<double> m_tof_btof1;
180 NTuple::Array<double> m_te_btof1;
181 NTuple::Array<double> m_tmu_btof1;
182 NTuple::Array<double> m_tpi_btof1;
183 NTuple::Array<double> m_tk_btof1;
184 NTuple::Array<double> m_tp_btof1;
185
186 NTuple::Array<double> m_qual_btof2;
187 NTuple::Array<double> m_tof_btof2;
188 NTuple::Array<double> m_te_btof2;
189 NTuple::Array<double> m_tmu_btof2;
190 NTuple::Array<double> m_tpi_btof2;
191 NTuple::Array<double> m_tk_btof2;
192 NTuple::Array<double> m_tp_btof2;
193
194 NTuple::Array<long> m_pidcode;
195 NTuple::Array<double> m_pidprob;
196 NTuple::Array<double> m_pidchiDedx;
197 NTuple::Array<double> m_pidchiTof1;
198 NTuple::Array<double> m_pidchiTof2;
199
200 NTuple::Item<long> m_npip;
201 NTuple::Item<long> m_npim;
202 NTuple::Item<long> m_nkp;
203 NTuple::Item<long> m_nkm;
204 NTuple::Item<long> m_np;
205 NTuple::Item<long> m_npb;
206
207 NTuple::Item<double> m_vfits_chi;
208 NTuple::Item<double> m_vfits_vx;
209 NTuple::Item<double> m_vfits_vy;
210 NTuple::Item<double> m_vfits_vz;
211 NTuple::Item<double> m_vfits_vr;
212
213 NTuple::Item<double> m_vfitp_chi;
214 NTuple::Item<double> m_vfitp_vx;
215 NTuple::Item<double> m_vfitp_vy;
216 NTuple::Item<double> m_vfitp_vz;
217 NTuple::Item<double> m_vfitp_vr;
218
219 NTuple::Item<double> m_vfit2_chi;
220 NTuple::Item<double> m_vfit2_mks;
221 NTuple::Item<double> m_vfit2_ct;
222 NTuple::Item<double> m_vfit2_dl;
223 NTuple::Item<double> m_vfit2_dle;
224
225 NTuple::Item<double> m_chi2_fs4c;
226 NTuple::Item<double> m_mks_fs4c;
227 NTuple::Item<double> m_mkspi_fs4c;
228 NTuple::Item<double> m_mksk_fs4c;
229 NTuple::Item<double> m_mkpi_fs4c;
230
231 NTuple::Item<double> m_4c_chi2;
232 NTuple::Item<double> m_4c_mks;
233 NTuple::Item<double> m_4c_mkspi;
234 NTuple::Item<double> m_4c_mksk;
235 NTuple::Item<double> m_4c_mkpi;
236 NTuple::Item<double> m_4c_ks_px;
237 NTuple::Item<double> m_4c_ks_py;
238 NTuple::Item<double> m_4c_ks_pz;
239 NTuple::Item<double> m_4c_ks_p;
240 NTuple::Item<double> m_4c_ks_cos;
241
242 NTuple::Item<double> m_chi2_5c;
243 NTuple::Item<double> m_mks_5c;
244 NTuple::Item<double> m_mkspi_5c;
245 NTuple::Item<double> m_mksk_5c;
246 NTuple::Item<double> m_mkpi_5c;
247};
248
249#endif
KsKpi::initialize
StatusCode initialize()
Definition KsKpi.cxx:82
KsKpi::KsKpi
KsKpi(const std::string &name, ISvcLocator *pSvcLocator)
Definition KsKpi.cxx:61
KsKpi::execute
StatusCode execute()
Definition KsKpi.cxx:409
KsKpi::finalize
StatusCode finalize()
Definition KsKpi.cxx:1109