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