BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
Single_histo.cxx
Go to the documentation of this file.
1void Single_histo( char* filein = "Single.root", char* fileout = "Single_histo.root" ) {
2 TFile* fin = new TFile( filein );
3 TTree* mdctrack = (TTree*)fin->Get( "mdctrack" );
4 TTree* dedx = (TTree*)fin->Get( "dedx" );
5 TTree* ext = (TTree*)fin->Get( "ext" );
6 TTree* tof = (TTree*)fin->Get( "tof" );
7 TTree* emc = (TTree*)fin->Get( "emc" );
8 TTree* muc = (TTree*)fin->Get( "muc" );
9 TTree* runinfo = (TTree*)fin->Get( "runinfo" );
10 TTree* mcpart = (TTree*)fin->Get( "mcpart" );
11 TTree* pid_kal = (TTree*)fin->Get( "pid_kal" );
12
13 double mdc_vx0;
14 double mdc_vy0;
15 double mdc_vz0;
16 double mdc_pch;
17 double mdc_trackid;
18 double mdc_charge;
19 double mdc_pxy0;
20 double mdc_px0;
21 double mdc_py0;
22 double mdc_pz0;
23 double mdc_r0;
24 double mdc_phi;
25 double mdc_theta;
26 double mdc_ndof;
27 double mdc_nster;
28 double mdc_stat;
29 double mdc_angd;
30 double mdc_rvxy0;
31 double mdc_rvz0;
32 double mdc_d0;
33 double mdc_phi0;
34 double mdc_kappa;
35 double mdc_dz;
36 double mdc_tanlamda;
37 double mdc_err11;
38 double mdc_err21;
39 double mdc_err22;
40 double mdc_err31;
41 double mdc_err32;
42 double mdc_err33;
43 double mdc_err41;
44 double mdc_err42;
45 double mdc_err43;
46 double mdc_err44;
47 double mdc_err51;
48 double mdc_err52;
49 double mdc_err53;
50 double mdc_err54;
51 double mdc_err55;
52
53 double dedx_ptrk;
54 double dedx_chie;
55 double dedx_chimu;
56 double dedx_chipi;
57 double dedx_chik;
58 double dedx_chip;
59 double dedx_probPH;
60 double dedx_normPH;
61 double dedx_ghit;
62 double dedx_thit;
63
64 double ext_tof1path;
65 double ext_tof1z;
66 double ext_tof1t;
67 double ext_tof1x;
68 double ext_tof1y;
69 double ext_tof2path;
70 double ext_tof2z;
71 double ext_tof2t;
72 double ext_tof2x;
73 double ext_tof2y;
74 double ext_emctheta;
75 double ext_emcphi;
76 double ext_emcpath;
77 double ext_mucz;
78 double ext_muct;
79 double ext_mucx;
80 double ext_mucy;
81
82 double tof_path;
83 double tof_zrhit;
84 double tof_ph;
85 double tof_tof;
86 double tof_errtof;
87 double tof_beta;
88 double tof_texpe;
89 double tof_texpmu;
90 double tof_texppi;
91 double tof_texpka;
92 double tof_texppro;
93 double tof_toffsete;
94 double tof_toffsetmu;
95 double tof_toffsetpi;
96 double tof_toffsetka;
97 double tof_toffsetpro;
98 double tof_toffsetatpr;
99 double tof_sigmae;
100 double tof_sigmamu;
101 double tof_sigmapi;
102 double tof_sigmaka;
103 double tof_sigmapro;
104 double tof_sigmaatpr;
105 double tof_quality;
106 double tof_t0;
107 double tof_errt0;
108 double tof_errz;
109 double tof_phi;
110 double tof_errphi;
111 double tof_energy;
112 double tof_errenergy;
113
114 double emc_x;
115 double emc_y;
116 double emc_z;
117 double emc_theta;
118 double emc_phi;
119 double emc_dx;
120 double emc_dy;
121 double emc_dz;
122 double emc_dtheta;
123 double emc_dphi;
124 double emc_energy;
125 double emc_dE;
126 double emc_eseed;
127 double emc_e3x3;
128 double emc_e5x5;
129 double emc_secp;
130 double emc_latp;
131
132 double muc_depth;
133 double muc_chi2;
134 double muc_rms;
135 double muc_xpos;
136 double muc_ypos;
137 double muc_zpos;
138 double muc_xsigma;
139 double muc_ysigma;
140 double muc_zsigma;
141 double muc_px;
142 double muc_py;
143 double muc_pz;
144 double muc_distance;
145 double muc_deltaphi;
146
147 double runinfo_runNo;
148 double runinfo_event;
149 double runinfo_totcharged;
150 double runinfo_totneutral;
151 double runinfo_tottracks;
152 double runinfo_ngood;
153
154 double mcpart_vx;
155 double mcpart_vy;
156 double mcpart_vz;
157 double mcpart_pxy;
158 double mcpart_px;
159 double mcpart_py;
160 double mcpart_pz;
161 double mcpart_pro;
162 double mcpart_e;
163 double mcpart_p;
164 double mcpart_theta;
165 double mcpart_phi;
166 double mcpart_m;
167 double mcpart_index;
168
169 double pid_ptrk;
170 double pid_cost;
171 double pid_dedx;
172 double pid_tof1;
173 double pid_tof2;
174 double prob_pion;
175 double prob_kaon;
176 double prob_proton;
177 double kal_px;
178 double kal_py;
179 double kal_pz;
180 double kal_p;
181 double kal_pxy;
182 double kal_n;
183 double kal_d0;
184 double kal_phi0;
185 double kal_kappa;
186 double kal_dz;
187 double kal_tanlamda;
188 double kal_err11;
189 double kal_err21;
190 double kal_err22;
191 double kal_err31;
192 double kal_err32;
193 double kal_err33;
194 double kal_err41;
195 double kal_err42;
196 double kal_err43;
197 double kal_err44;
198 double kal_err51;
199 double kal_err52;
200 double kal_err53;
201 double kal_err54;
202 double kal_err55;
203
204 mdctrack->SetBranchAddress( "mdc_vx0", &mdc_vx0 );
205 mdctrack->SetBranchAddress( "mdc_vy0", &mdc_vy0 );
206 mdctrack->SetBranchAddress( "mdc_vz0", &mdc_vz0 );
207 mdctrack->SetBranchAddress( "mdc_pch", &mdc_pch );
208 mdctrack->SetBranchAddress( "mdc_trackid", &mdc_trackid );
209 mdctrack->SetBranchAddress( "mdc_charge", &mdc_charge );
210 mdctrack->SetBranchAddress( "mdc_pxy0", &mdc_pxy0 );
211 mdctrack->SetBranchAddress( "mdc_px0", &mdc_px0 );
212 mdctrack->SetBranchAddress( "mdc_py0", &mdc_py0 );
213 mdctrack->SetBranchAddress( "mdc_pz0", &mdc_pz0 );
214 mdctrack->SetBranchAddress( "mdc_r0", &mdc_r0 );
215 mdctrack->SetBranchAddress( "mdc_phi", &mdc_phi );
216 mdctrack->SetBranchAddress( "mdc_theta", &mdc_theta );
217 mdctrack->SetBranchAddress( "mdc_ndof", &mdc_ndof );
218 mdctrack->SetBranchAddress( "mdc_nster", &mdc_nster );
219 mdctrack->SetBranchAddress( "mdc_stat", &mdc_stat );
220 mdctrack->SetBranchAddress( "mdc_angd", &mdc_angd );
221 mdctrack->SetBranchAddress( "mdc_rvxy0", &mdc_rvxy0 );
222 mdctrack->SetBranchAddress( "mdc_rvz0", &mdc_rvz0 );
223 mdctrack->SetBranchAddress( "mdc_d0", &mdc_d0 );
224 mdctrack->SetBranchAddress( "mdc_phi0", &mdc_phi0 );
225 mdctrack->SetBranchAddress( "mdc_kappa", &mdc_kappa );
226 mdctrack->SetBranchAddress( "mdc_dz", &mdc_dz );
227 mdctrack->SetBranchAddress( "mdc_tanlamda", &mdc_tanlamda );
228 mdctrack->SetBranchAddress( "mdc_err11", &mdc_err11 );
229 mdctrack->SetBranchAddress( "mdc_err21", &mdc_err21 );
230 mdctrack->SetBranchAddress( "mdc_err22", &mdc_err22 );
231 mdctrack->SetBranchAddress( "mdc_err31", &mdc_err31 );
232 mdctrack->SetBranchAddress( "mdc_err32", &mdc_err32 );
233 mdctrack->SetBranchAddress( "mdc_err33", &mdc_err33 );
234 mdctrack->SetBranchAddress( "mdc_err41", &mdc_err41 );
235 mdctrack->SetBranchAddress( "mdc_err42", &mdc_err42 );
236 mdctrack->SetBranchAddress( "mdc_err43", &mdc_err43 );
237 mdctrack->SetBranchAddress( "mdc_err44", &mdc_err44 );
238 mdctrack->SetBranchAddress( "mdc_err51", &mdc_err51 );
239 mdctrack->SetBranchAddress( "mdc_err52", &mdc_err52 );
240 mdctrack->SetBranchAddress( "mdc_err53", &mdc_err53 );
241 mdctrack->SetBranchAddress( "mdc_err54", &mdc_err54 );
242 mdctrack->SetBranchAddress( "mdc_err55", &mdc_err55 );
243
244 dedx->SetBranchAddress( "dedx_ptrk", &dedx_ptrk );
245 dedx->SetBranchAddress( "dedx_chie", &dedx_chie );
246 dedx->SetBranchAddress( "dedx_chimu", &dedx_chimu );
247 dedx->SetBranchAddress( "dedx_chipi", &dedx_chipi );
248 dedx->SetBranchAddress( "dedx_chik", &dedx_chik );
249 dedx->SetBranchAddress( "dedx_chip", &dedx_chip );
250 dedx->SetBranchAddress( "dedx_probPH", &dedx_probPH );
251 dedx->SetBranchAddress( "dedx_normPH", &dedx_normPH );
252 dedx->SetBranchAddress( "dedx_ghit", &dedx_ghit );
253 dedx->SetBranchAddress( "dedx_thit", &dedx_thit );
254
255 ext->SetBranchAddress( "ext_tof1path", &ext_tof1path );
256 ext->SetBranchAddress( "ext_tof1z", &ext_tof1z );
257 ext->SetBranchAddress( "ext_tof1t", &ext_tof1t );
258 ext->SetBranchAddress( "ext_tof1x", &ext_tof1x );
259 ext->SetBranchAddress( "ext_tof1y", &ext_tof1y );
260 ext->SetBranchAddress( "ext_tof2path", &ext_tof2path );
261 ext->SetBranchAddress( "ext_tof2z", &ext_tof2z );
262 ext->SetBranchAddress( "ext_tof2t", &ext_tof2t );
263 ext->SetBranchAddress( "ext_tof2x", &ext_tof2x );
264 ext->SetBranchAddress( "ext_tof2y", &ext_tof2y );
265 ext->SetBranchAddress( "ext_emctheta", &ext_emctheta );
266 ext->SetBranchAddress( "ext_emcphi", &ext_emcphi );
267 ext->SetBranchAddress( "ext_emcpath", &ext_emcpath );
268 ext->SetBranchAddress( "ext_mucz", &ext_mucz );
269 ext->SetBranchAddress( "ext_muct", &ext_muct );
270 ext->SetBranchAddress( "ext_mucx", &ext_mucx );
271 ext->SetBranchAddress( "ext_mucy", &ext_mucy );
272
273 tof->SetBranchAddress( "tof_path", &tof_path );
274 tof->SetBranchAddress( "tof_zrhit", &tof_zrhit );
275 tof->SetBranchAddress( "tof_ph", &tof_ph );
276 tof->SetBranchAddress( "tof_tof", &tof_tof );
277 tof->SetBranchAddress( "tof_errtof", &tof_errtof );
278 tof->SetBranchAddress( "tof_beta", &tof_beta );
279 tof->SetBranchAddress( "tof_texpe", &tof_texpe );
280 tof->SetBranchAddress( "tof_texpmu", &tof_texpmu );
281 tof->SetBranchAddress( "tof_texppi", &tof_texppi );
282 tof->SetBranchAddress( "tof_texpka", &tof_texpka );
283 tof->SetBranchAddress( "tof_texppro", &tof_texppro );
284 tof->SetBranchAddress( "tof_toffsete", &tof_toffsete );
285 tof->SetBranchAddress( "tof_toffsetmu", &tof_toffsetmu );
286 tof->SetBranchAddress( "tof_toffsetpi", &tof_toffsetpi );
287 tof->SetBranchAddress( "tof_toffsetka", &tof_toffsetka );
288 tof->SetBranchAddress( "tof_toffsetpro", &tof_toffsetpro );
289 tof->SetBranchAddress( "tof_toffsetatpr", &tof_toffsetatpr );
290 tof->SetBranchAddress( "tof_sigmae", &tof_sigmae );
291 tof->SetBranchAddress( "tof_sigmamu", &tof_sigmamu );
292 tof->SetBranchAddress( "tof_sigmapi", &tof_sigmapi );
293 tof->SetBranchAddress( "tof_sigmaka", &tof_sigmaka );
294 tof->SetBranchAddress( "tof_sigmapro", &tof_sigmapro );
295 tof->SetBranchAddress( "tof_sigmaatpr", &tof_sigmaatpr );
296 tof->SetBranchAddress( "tof_quality", &tof_quality );
297 tof->SetBranchAddress( "tof_t0", &tof_t0 );
298 tof->SetBranchAddress( "tof_errt0", &tof_errt0 );
299 tof->SetBranchAddress( "tof_errz", &tof_errz );
300 tof->SetBranchAddress( "tof_phi", &tof_phi );
301 tof->SetBranchAddress( "tof_errphi", &tof_errphi );
302 tof->SetBranchAddress( "tof_energy", &tof_energy );
303 tof->SetBranchAddress( "tof_errenergy", &tof_errenergy );
304
305 emc->SetBranchAddress( "emc_x", &emc_x );
306 emc->SetBranchAddress( "emc_y", &emc_y );
307 emc->SetBranchAddress( "emc_z", &emc_z );
308 emc->SetBranchAddress( "emc_theta", &emc_theta );
309 emc->SetBranchAddress( "emc_phi", &emc_phi );
310 emc->SetBranchAddress( "emc_dx", &emc_dx );
311 emc->SetBranchAddress( "emc_dy", &emc_dy );
312 emc->SetBranchAddress( "emc_dz", &emc_dz );
313 emc->SetBranchAddress( "emc_dtheta", &emc_dtheta );
314 emc->SetBranchAddress( "emc_dphi", &emc_dphi );
315 emc->SetBranchAddress( "emc_energy", &emc_energy );
316 emc->SetBranchAddress( "emc_dE", &emc_dE );
317 emc->SetBranchAddress( "emc_eseed", &emc_eseed );
318 emc->SetBranchAddress( "emc_e3x3", &emc_e3x3 );
319 emc->SetBranchAddress( "emc_e5x5", &emc_e5x5 );
320 emc->SetBranchAddress( "emc_secp", &emc_secp );
321 emc->SetBranchAddress( "emc_latp", &emc_latp );
322
323 muc->SetBranchAddress( "muc_depth", &muc_depth );
324 muc->SetBranchAddress( "muc_chi2", &muc_chi2 );
325 muc->SetBranchAddress( "muc_rms", &muc_rms );
326 muc->SetBranchAddress( "muc_xpos", &muc_xpos );
327 muc->SetBranchAddress( "muc_ypos", &muc_ypos );
328 muc->SetBranchAddress( "muc_zpos", &muc_zpos );
329 muc->SetBranchAddress( "muc_xsigma", &muc_xsigma );
330 muc->SetBranchAddress( "muc_ysigma", &muc_ysigma );
331 muc->SetBranchAddress( "muc_zsigma", &muc_zsigma );
332 muc->SetBranchAddress( "muc_px", &muc_px );
333 muc->SetBranchAddress( "muc_py", &muc_py );
334 muc->SetBranchAddress( "muc_pz", &muc_pz );
335 muc->SetBranchAddress( "muc_distance", &muc_distance );
336 muc->SetBranchAddress( "muc_deltaphi", &muc_deltaphi );
337
338 runinfo->SetBranchAddress( "runinfo_runNo", &runinfo_runNo );
339 runinfo->SetBranchAddress( "runinfo_event", &runinfo_event );
340 runinfo->SetBranchAddress( "runinfo_totcharged", &runinfo_totcharged );
341 runinfo->SetBranchAddress( "runinfo_totneutral", &runinfo_totneutral );
342 runinfo->SetBranchAddress( "runinfo_tottracks", &runinfo_tottracks );
343 runinfo->SetBranchAddress( "runinfo_ngood", &runinfo_ngood );
344
345 mcpart->SetBranchAddress( "mcpart_vx", &mcpart_vx );
346 mcpart->SetBranchAddress( "mcpart_vy", &mcpart_vy );
347 mcpart->SetBranchAddress( "mcpart_vz", &mcpart_vz );
348 mcpart->SetBranchAddress( "mcpart_pxy", &mcpart_pxy );
349 mcpart->SetBranchAddress( "mcpart_px", &mcpart_px );
350 mcpart->SetBranchAddress( "mcpart_py", &mcpart_py );
351 mcpart->SetBranchAddress( "mcpart_pz", &mcpart_pz );
352 mcpart->SetBranchAddress( "mcpart_pro", &mcpart_pro );
353 mcpart->SetBranchAddress( "mcpart_e", &mcpart_e );
354 mcpart->SetBranchAddress( "mcpart_p", &mcpart_p );
355 mcpart->SetBranchAddress( "mcpart_theta", &mcpart_theta );
356 mcpart->SetBranchAddress( "mcpart_phi", &mcpart_phi );
357 mcpart->SetBranchAddress( "mcpart_m", &mcpart_m );
358 mcpart->SetBranchAddress( "mcpart_index", &mcpart_index );
359
360 pid_kal->SetBranchAddress( "pid_ptrk", &pid_ptrk );
361 pid_kal->SetBranchAddress( "pid_cost", &pid_cost );
362 pid_kal->SetBranchAddress( "pid_dedx", &pid_dedx );
363 pid_kal->SetBranchAddress( "pid_tof1", &pid_tof1 );
364 pid_kal->SetBranchAddress( "pid_tof2", &pid_tof2 );
365 pid_kal->SetBranchAddress( "prob_pion", &prob_pion );
366 pid_kal->SetBranchAddress( "prob_kaon", &prob_kaon );
367 pid_kal->SetBranchAddress( "prob_proton", &prob_proton );
368 pid_kal->SetBranchAddress( "kal_px", &kal_px );
369 pid_kal->SetBranchAddress( "kal_py", &kal_py );
370 pid_kal->SetBranchAddress( "kal_pz", &kal_pz );
371 pid_kal->SetBranchAddress( "kal_p", &kal_p );
372 pid_kal->SetBranchAddress( "kal_pxy", &kal_pxy );
373 pid_kal->SetBranchAddress( "kal_n", &kal_n );
374 pid_kal->SetBranchAddress( "kal_d0", &kal_d0 );
375 pid_kal->SetBranchAddress( "kal_phi0", &kal_phi0 );
376 pid_kal->SetBranchAddress( "kal_kappa", &kal_kappa );
377 pid_kal->SetBranchAddress( "kal_dz", &kal_dz );
378 pid_kal->SetBranchAddress( "kal_tanlamda", &kal_tanlamda );
379 pid_kal->SetBranchAddress( "kal_err11", &kal_err11 );
380 pid_kal->SetBranchAddress( "kal_err21", &kal_err21 );
381 pid_kal->SetBranchAddress( "kal_err22", &kal_err22 );
382 pid_kal->SetBranchAddress( "kal_err31", &kal_err31 );
383 pid_kal->SetBranchAddress( "kal_err32", &kal_err32 );
384 pid_kal->SetBranchAddress( "kal_err33", &kal_err33 );
385 pid_kal->SetBranchAddress( "kal_err41", &kal_err41 );
386 pid_kal->SetBranchAddress( "kal_err42", &kal_err42 );
387 pid_kal->SetBranchAddress( "kal_err43", &kal_err43 );
388 pid_kal->SetBranchAddress( "kal_err44", &kal_err44 );
389 pid_kal->SetBranchAddress( "kal_err51", &kal_err51 );
390 pid_kal->SetBranchAddress( "kal_err52", &kal_err52 );
391 pid_kal->SetBranchAddress( "kal_err53", &kal_err53 );
392 pid_kal->SetBranchAddress( "kal_err54", &kal_err54 );
393 pid_kal->SetBranchAddress( "kal_err55", &kal_err55 );
394
395 // TFile *fout = new TFile("/besfs2/users/yuanwl/pion/pion_test_histo.root","recreate");
396 TFile* fout = new TFile( fileout, "recreate" );
397
398 TH1F* Mdc_vx0 = new TH1F( "Mdc_vx0", "vx0", 200, -1, 1 );
399 TH1F* Mdc_vy0 = new TH1F( "Mdc_vy0", "vy0", 200, -1, 1 );
400 TH1F* Mdc_vz0 = new TH1F( "Mdc_vz0", "vz0", 200, -8, 8 );
401 TH1F* Mdc_pch = new TH1F( "Mdc_pch", "pch", 200, -0.5, 3 );
402 TH1F* Mdc_trackid = new TH1F( "Mdc_trackid", "trackid", 200, -1, 1 );
403 TH1F* Mdc_charge = new TH1F( "Mdc_charge", "charge", 200, -2, 2 );
404 TH1F* Mdc_pxy0 = new TH1F( "Mdc_pxy0", "pxy0", 200, -0.5, 5 );
405 TH1F* Mdc_px0 = new TH1F( "Mdc_px0", "px0", 200, -3, 3 );
406 TH1F* Mdc_py0 = new TH1F( "Mdc_py0", "py0", 200, -3, 3 );
407 TH1F* Mdc_pz0 = new TH1F( "Mdc_pz0", "pz0", 200, -3, 3 );
408 TH1F* Mdc_r0 = new TH1F( "Mdc_r0", "r0", 200, -0.5, 1 );
409 TH1F* Mdc_phi = new TH1F( "Mdc_phi", "phi", 200, -200, 200 );
410 TH1F* Mdc_theta = new TH1F( "Mdc_theta", "theta", 200, -180, 180 );
411 TH1F* Mdc_ndof = new TH1F( "Mdc_ndof", "ndof", 200, 10, 60 );
412 TH1F* Mdc_nster = new TH1F( "Mdc_nster", "nster", 200, 5, 35 );
413 TH1F* Mdc_stat = new TH1F( "Mdc_stat", "stat", 200, -0.5, 0.5 );
414 TH1F* Mdc_angd = new TH1F( "Mdc_angd", "angd", 200, -0.1, 3 );
415 TH1F* Mdc_rvxy0 = new TH1F( "Mdc_rvxy0", "rvxy0", 200, 0, 0.4 );
416 TH1F* Mdc_rvz0 = new TH1F( "Mdc_rvz0", "rvz0", 200, -8, 8 );
417 TH1F* Mdc_d0 = new TH1F( "Mdc_d0", "d0", 200, -1, 1 );
418 TH1F* Mdc_phi0 = new TH1F( "Mdc_phi0", "phi0", 200, 0, 8 );
419 TH1F* Mdc_kappa = new TH1F( "Mdc_kappa", "kappa", 200, -10, 8 );
420 TH1F* Mdc_dz = new TH1F( "Mdc_dz", "dz", 200, -8, 8 );
421 TH1F* Mdc_tanlamda = new TH1F( "Mdc_tanlamda", "tan(lamda)", 200, -3, 3 );
422 TH1F* Mdc_err11 = new TH1F( "Mdc_err11", "error", 200, -0.0005, 0.0005 );
423 TH1F* Mdc_err21 = new TH1F( "Mdc_err21", "error", 200, -5e-5, 5e-5 );
424 TH1F* Mdc_err22 = new TH1F( "Mdc_err22", "error", 200, -5e-4, 5e-4 );
425 TH1F* Mdc_err31 = new TH1F( "Mdc_err31", "error", 200, -5e-3, 5e-3 );
426 TH1F* Mdc_err32 = new TH1F( "Mdc_err32", "error", 200, -5e-5, 5e-5 );
427 TH1F* Mdc_err33 = new TH1F( "Mdc_err33", "error", 200, -5e-6, 5e-6 );
428 TH1F* Mdc_err41 = new TH1F( "Mdc_err41", "error", 200, -5e-6, 5e-6 );
429 TH1F* Mdc_err42 = new TH1F( "Mdc_err42", "error", 200, -5e-5, 5e-5 );
430 TH1F* Mdc_err43 = new TH1F( "Mdc_err43", "error", 200, -5e-6, 5e-6 );
431 TH1F* Mdc_err44 = new TH1F( "Mdc_err44", "error", 200, -5e-5, 5e-5 );
432 TH1F* Mdc_err51 = new TH1F( "Mdc_err51", "error", 200, -5e-4, 5e-4 );
433 TH1F* Mdc_err52 = new TH1F( "Mdc_err52", "error", 200, -8e-6, 5e-6 );
434 TH1F* Mdc_err53 = new TH1F( "Mdc_err53", "error", 200, 0, 0.04 );
435 TH1F* Mdc_err54 = new TH1F( "Mdc_err54", "error", 200, -8e-4, 8e-4 );
436 TH1F* Mdc_err55 = new TH1F( "Mdc_err55", "error", 200, -5e-5, 5e-5 );
437
438 TH1F* Dedx_ptrk = new TH1F( "Dedx_ptrk", "ptrk", 200, 0, 5 );
439 TH1F* Dedx_chie = new TH1F( "Dedx_chie", "chie", 200, -10, 150 );
440 TH1F* Dedx_chimu = new TH1F( "Dedx_chimu", "chimu", 200, -1e18, 1e18 );
441 TH1F* Dedx_chipi = new TH1F( "Dedx_chipi", "chipi", 200, -8, 150 );
442 TH1F* Dedx_chik = new TH1F( "Dedx_chik", "chik", 200, -20, 40 );
443 TH1F* Dedx_chip = new TH1F( "Dedx_chip", "chip", 200, -20, 20 );
444 TH1F* Dedx_probPH = new TH1F( "Dedx_probPH", "probPH", 200, 200, 4500 );
445 TH1F* Dedx_normPH = new TH1F( "Dedx_normPH", "normPH", 200, -80, 80 );
446 TH1F* Dedx_ghit = new TH1F( "Dedx_ghit", "ghit", 200, -10, 30 );
447 TH1F* Dedx_thit = new TH1F( "Dedx_thit", "thit", 200, 0, 100 );
448
449 TH1F* Ext_tof1path = new TH1F( "Ext_tof1path", "tof1path", 200, 80, 200 );
450 TH1F* Ext_tof1z = new TH1F( "Ext_tof1z", "tof1z", 200, 0, 4 );
451 TH1F* Ext_tof1t = new TH1F( "Ext_tof1t", "tof1t", 200, 0, 5 );
452 TH1F* Ext_tof1x = new TH1F( "Ext_tof1x", "tof1x", 200, -2, 20 );
453 TH1F* Ext_tof1y = new TH1F( "Ext_tof1y", "tof1y", 200, -1, 25 );
454 TH1F* Ext_tof2path = new TH1F( "Ext_tof2path", "tof2path", 200, 50, 200 );
455 TH1F* Ext_tof2z = new TH1F( "Ext_tof2z", "tof2z", 200, -1, 15 );
456 TH1F* Ext_tof2t = new TH1F( "Ext_tof2t", "tof2t", 200, -1, 20 );
457 TH1F* Ext_tof2x = new TH1F( "Ext_tof2x", "tof2x", 200, -1, 5 );
458 TH1F* Ext_tof2y = new TH1F( "Ext_tof2y", "tof2y", 200, -1, 5 );
459 TH1F* Ext_emctheta = new TH1F( "Ext_emctheta", "emctheta", 200, -50, 50 );
460 TH1F* Ext_emcphi = new TH1F( "Ext_emcphi", "emcphi", 200, -50, 50 );
461 TH1F* Ext_emcpath = new TH1F( "Ext_emcpath", "emcpath", 200, -5, 40 );
462 TH1F* Ext_mucz = new TH1F( "Ext_mucz", "mucz", 200, -5, 60 );
463 TH1F* Ext_muct = new TH1F( "Ext_muct", "muct", 200, -5, 40 );
464 TH1F* Ext_mucx = new TH1F( "Ext_mucx", "mucx", 200, -5, 80 );
465 TH1F* Ext_mucy = new TH1F( "Ext_mucy", "mucy", 200, -5, 80 );
466
467 TH1F* Tof_path = new TH1F( "Tof_path", "path", 200, 0, 250 );
468 TH1F* Tof_zrhit = new TH1F( "Tof_zrhit", "zrhit", 200, -300, 300 );
469 TH1F* Tof_ph = new TH1F( "Tof_ph", "ph", 200, 0, 8000 );
470 TH1F* Tof_tof = new TH1F( "Tof_tof", "tof", 200, -100, 100 );
471 TH1F* Tof_errtof = new TH1F( "Tof_errtof", "errtof", 200, -50, 300 );
472 TH1F* Tof_beta = new TH1F( "Tof_beta", "beta", 200, -1, 1 );
473 TH1F* Tof_texpe = new TH1F( "Tof_texpe", "texpe", 200, -50, 50 );
474 TH1F* Tof_texpmu = new TH1F( "Tof_texpmu", "texpmu", 200, -50, 50 );
475 TH1F* Tof_texppi = new TH1F( "Tof_texppi", "texppi", 200, -50, 50 );
476 TH1F* Tof_texpka = new TH1F( "Tof_texpka", "texpka", 200, -50, 50 );
477 TH1F* Tof_texppro = new TH1F( "Tof_texppro", "texppro", 200, -50, 50 );
478 TH1F* Tof_toffsete = new TH1F( "Tof_toffsete", "toffsete", 200, -1, 1 );
479 TH1F* Tof_toffsetmu = new TH1F( "Tof_toffsetmu", "toffsetmu", 200, -1, 1 );
480 TH1F* Tof_toffsetpi = new TH1F( "Tof_toffsetpi", "toffsetpi", 200, -1, 1 );
481 TH1F* Tof_toffsetka = new TH1F( "Tof_toffsetka", "toffsetka", 200, -1, 1 );
482 TH1F* Tof_toffsetpro = new TH1F( "Tof_toffsetpro", "toffsetpro", 200, -1, 1 );
483 TH1F* Tof_toffsetatpr = new TH1F( "Tof_toffsetatpr", "toffsetatpr", 200, -1, 1 );
484 TH1F* Tof_sigmae = new TH1F( "Tof_sigmae", "sigmae", 200, -50, 250 );
485 TH1F* Tof_sigmamu = new TH1F( "Tof_sigmamu", "sigmamu", 200, -1, 1 );
486 TH1F* Tof_sigmapi = new TH1F( "Tof_sigmapi", "sigmapi", 200, -1, 1 );
487 TH1F* Tof_sigmaka = new TH1F( "Tof_sigmaka", "sigmaka", 200, -1, 1 );
488 TH1F* Tof_sigmapro = new TH1F( "Tof_sigmapro", "sigmapro", 200, -1, 1 );
489 TH1F* Tof_sigmaatpr = new TH1F( "Tof_sigmaatpr", "sigmaatpr", 200, -1, 1 );
490 TH1F* Tof_quality = new TH1F( "Tof_quality", "quality", 200, -1, 1 );
491 TH1F* Tof_t0 = new TH1F( "Tof_t0", "t0", 200, 0, 1000 );
492 TH1F* Tof_errt0 = new TH1F( "Tof_errt0", "errt0", 200, -1, 1 );
493 TH1F* Tof_errz = new TH1F( "Tof_errz", "errz", 200, 0, 5 );
494 TH1F* Tof_phi = new TH1F( "Tof_phi", "phi", 200, 4e5, 6e5 );
495 TH1F* Tof_errphi = new TH1F( "Tof_errphi", "errphi", 200, 4e5, 6e5 );
496 TH1F* Tof_energy = new TH1F( "Tof_energy", "energy", 200, 0, 150 );
497 TH1F* Tof_errenergy = new TH1F( "Tof_errenergy", "errenergy", 200, 0, 15000 );
498
499 TH1F* Emc_x = new TH1F( "Emc_x", "x", 200, -100, 100 );
500 TH1F* Emc_y = new TH1F( "Emc_y", "y", 200, -100, 100 );
501 TH1F* Emc_z = new TH1F( "Emc_z", "z", 200, -150, 150 );
502 TH1F* Emc_theta = new TH1F( "Emc_theta", "theta", 200, 0, 180 );
503 TH1F* Emc_phi = new TH1F( "Emc_phi", "phi", 200, -200, 200 );
504 TH1F* Emc_dx = new TH1F( "Emc_dx", "dx", 200, 0, 5 );
505 TH1F* Emc_dy = new TH1F( "Emc_dy", "dy", 200, 0, 5 );
506 TH1F* Emc_dz = new TH1F( "Emc_dz", "dz", 200, 0, 3 );
507 TH1F* Emc_dtheta = new TH1F( "Emc_dtheta", "dtheta", 200, 0, 2 );
508 TH1F* Emc_dphi = new TH1F( "Emc_dphi", "dphi", 200, 0, 3 );
509 TH1F* Emc_energy = new TH1F( "Emc_energy", "", 200, -1, 2 );
510 TH1F* Emc_dE = new TH1F( "Emc_dE", "dE", 200, 0, 0.1 );
511 TH1F* Emc_eseed = new TH1F( "Emc_eseed", "eseed", 200, 0, 2 );
512 TH1F* Emc_e3x3 = new TH1F( "Emc_e3x3", "e3x3", 200, 0, 2 );
513 TH1F* Emc_e5x5 = new TH1F( "Emc_e5x5", "e5x5", 200, 0, 2 );
514 TH1F* Emc_secp = new TH1F( "Emc_secp", "secp", 200, 0, 150 );
515 TH1F* Emc_latp = new TH1F( "Emc_latp", "latp", 200, 0, 1 );
516
517 TH1F* Muc_depth = new TH1F( "Muc_depth", "depth", 200, 0, 60 );
518 TH1F* Muc_chi2 = new TH1F( "Muc_chi2", "chi2", 200, 0, 2 );
519 TH1F* Muc_rms = new TH1F( "Muc_rms", "rms", 200, -1, 1 );
520 TH1F* Muc_xpos = new TH1F( "Muc_xpos", "xpos", 200, -200, 200 );
521 TH1F* Muc_ypos = new TH1F( "Muc_ypos", "ypos", 200, -200, 200 );
522 TH1F* Muc_zpos = new TH1F( "Muc_zpos", "zpos", 200, -300, 300 );
523 TH1F* Muc_xsigma = new TH1F( "Muc_xsigma", "xsigma", 200, 0, 20000 );
524 TH1F* Muc_ysigma = new TH1F( "Muc_ysigma", "ysigma", 200, 0, 10 );
525 TH1F* Muc_zsigma = new TH1F( "Muc_zsigma", "zsigma", 200, -1, 1 );
526 TH1F* Muc_px = new TH1F( "Muc_px", "px", 200, -500, 500 );
527 TH1F* Muc_py = new TH1F( "Muc_py", "py", 200, -1000, 1000 );
528 TH1F* Muc_pz = new TH1F( "Muc_pz", "pz", 200, -10, 10 );
529 TH1F* Muc_distance = new TH1F( "Muc_distance", "distance", 200, -30, 30 );
530 TH1F* Muc_deltaphi = new TH1F( "Muc_deltaphi", "deltaphi", 200, -0.5, 0.5 );
531
532 TH1F* Runinfo_runNo = new TH1F( "Runinfo_runNo", "runNo", 200, -10000, 10000 );
533 TH1F* Runinfo_event = new TH1F( "Runinfo_event", "event", 200, 0, 2e5 );
534 TH1F* Runinfo_totcharged = new TH1F( "Runinfo_totcharged", "totcharged", 200, 0, 9 );
535 TH1F* Runinfo_totneutral = new TH1F( "Runinfo_totneutral", "totneutral", 200, 0, 9 );
536 TH1F* Runinfo_tottracks = new TH1F( "Runinfo_tottracks", "tottracks", 200, 0, 10 );
537 TH1F* Runinfo_ngood = new TH1F( "Runinfo_ngood", "ngood", 200, 0, 4 );
538
539 TH1F* Mcpart_vx = new TH1F( "Mcpart_vx", "vx", 200, -1, 1 );
540 TH1F* Mcpart_vy = new TH1F( "Mcpart_vy", "vy", 200, -1, 1 );
541 TH1F* Mcpart_vz = new TH1F( "Mcpart_vz", "vz", 200, -8, 8 );
542 TH1F* Mcpart_pxy = new TH1F( "Mcpart_pxy", "pxy", 200, -0.5, 5 );
543 TH1F* Mcpart_px = new TH1F( "Mcpart_px", "px", 200, -3, 3 );
544 TH1F* Mcpart_py = new TH1F( "Mcpart_py", "py", 200, -3, 3 );
545 TH1F* Mcpart_pz = new TH1F( "Mcpart_pz", "pz", 200, -3, 3 );
546 TH1F* Mcpart_e = new TH1F( "Mcpart_e", "e", 200, 0, 2.5 );
547 TH1F* Mcpart_p = new TH1F( "Mcpart_p", "p", 200, -0.1, 2.1 );
548 TH1F* Mcpart_theta = new TH1F( "Mcpart_theta", "theta", 200, 0, 180 );
549 TH1F* Mcpart_phi = new TH1F( "Mcpart_phi", "phi", 200, -200, 200 );
550 TH1F* Mcpart_m = new TH1F( "Mcpart_m", "m", 200, 0, 0.2 );
551 TH1F* Mcpart_pro = new TH1F( "Mcpart_pro", "pro", 200, -2500, 2500 );
552 TH1F* Mcpart_index = new TH1F( "Mcpart_index", "index", 200, -1, 1 );
553
554 TH1F* Pid_ptrk = new TH1F( "Pid_ptrk", "ptrk", 200, 0, 3 );
555 TH1F* Pid_cost = new TH1F( "Pid_cost", "cost", 200, -1, 1 );
556 TH1F* Pid_dedx = new TH1F( "Pid_dedx", "dedx", 200, -10, 100 );
557 TH1F* Pid_tof1 = new TH1F( "Pid_tof1", "tof1", 200, -10, 150 );
558 TH1F* Pid_tof2 = new TH1F( "Pid_tof2", "tof2", 200, -10, 70 );
559 TH1F* Prob_pion = new TH1F( "Prob_pion", "pion", 200, 0, 1 );
560 TH1F* Prob_kaon = new TH1F( "Prob_kaon", "kaon", 200, 0, 1 );
561 TH1F* Prob_proton = new TH1F( "Prob_proton", "proton", 200, 0, 1 );
562 TH1F* Kal_px = new TH1F( "Kal_px", "px", 200, -2.5, 2.5 );
563 TH1F* Kal_py = new TH1F( "Kal_py", "py", 200, -2.5, 2.5 );
564 TH1F* Kal_pz = new TH1F( "Kal_pz", "pz", 200, -2.5, 2.5 );
565 TH1F* Kal_p = new TH1F( "Kal_p", "p", 200, 0, 2.5 );
566 TH1F* Kal_pxy = new TH1F( "Kal_pxy", "pxy", 200, 0, 2.5 );
567 TH1F* Kal_n = new TH1F( "Kal_n", "n", 200, 0, 1.1 );
568 TH1F* Kal_d0 = new TH1F( "Kal_d0", "d0", 200, -1, 1 );
569 TH1F* Kal_phi0 = new TH1F( "Kal_phi0", "phi0", 200, 0, 8 );
570 TH1F* Kal_kappa = new TH1F( "Kal_kappa", "kappa", 200, -10, 5 );
571 TH1F* Kal_dz = new TH1F( "Kal_dz", "dz", 200, -8, 8 );
572 TH1F* Kal_tanlamda = new TH1F( "Kal_tanlamda", "tan(lamda)", 200, -3, 3 );
573 TH1F* Kal_err11 = new TH1F( "Kal_err11", "error", 200, -0.0008, 0.0009 );
574 TH1F* Kal_err21 = new TH1F( "Kal_err21", "error", 200, -5e-5, 30e-5 );
575 TH1F* Kal_err22 = new TH1F( "Kal_err22", "error", 200, -5e-4, 5e-4 );
576 TH1F* Kal_err31 = new TH1F( "Kal_err31", "error", 200, -5e-3, 5e-3 );
577 TH1F* Kal_err32 = new TH1F( "Kal_err32", "error", 200, -5e-5, 5e-5 );
578 TH1F* Kal_err33 = new TH1F( "Kal_err33", "error", 200, 0, 0.01 );
579 TH1F* Kal_err41 = new TH1F( "Kal_err41", "error", 200, -0.2, 0.2 );
580 TH1F* Kal_err42 = new TH1F( "Kal_err42", "error", 200, -5e-5, 5e-5 );
581 TH1F* Kal_err43 = new TH1F( "Kal_err43", "error", 200, -0.2, 0.2 );
582 TH1F* Kal_err44 = new TH1F( "Kal_err44", "error", 200, -1, 20 );
583 TH1F* Kal_err51 = new TH1F( "Kal_err51", "error", 200, -5e-4, 5e-4 );
584 TH1F* Kal_err52 = new TH1F( "Kal_err52", "error", 200, -5e-6, 5e-6 );
585 TH1F* Kal_err53 = new TH1F( "Kal_err53", "error", 200, 0, 1e-3 );
586 TH1F* Kal_err54 = new TH1F( "Kal_err54", "error", 200, -1e-3, 8e-4 );
587 TH1F* Kal_err55 = new TH1F( "Kal_err55", "error", 200, -5e-5, 8e-5 );
588
589 for ( int i = 0; i < mdctrack->GetEntries(); i++ )
590 {
591
592 mdctrack->GetEntry( i );
593 Mdc_vx0->Fill( mdc_vx0 );
594 Mdc_vy0->Fill( mdc_vy0 );
595 Mdc_vz0->Fill( mdc_vz0 );
596 Mdc_pch->Fill( mdc_pch );
597 Mdc_trackid->Fill( mdc_trackid );
598 Mdc_charge->Fill( mdc_charge );
599 Mdc_pxy0->Fill( mdc_pxy0 );
600 Mdc_px0->Fill( mdc_px0 );
601 Mdc_py0->Fill( mdc_py0 );
602 Mdc_pz0->Fill( mdc_pz0 );
603 Mdc_r0->Fill( mdc_r0 );
604 Mdc_phi->Fill( mdc_phi );
605 Mdc_theta->Fill( mdc_theta );
606 Mdc_ndof->Fill( mdc_ndof );
607 Mdc_nster->Fill( mdc_nster );
608 Mdc_stat->Fill( mdc_stat );
609 Mdc_angd->Fill( mdc_angd );
610 Mdc_rvxy0->Fill( mdc_rvxy0 );
611 Mdc_rvz0->Fill( mdc_rvz0 );
612 Mdc_d0->Fill( mdc_d0 );
613 Mdc_phi0->Fill( mdc_phi0 );
614 Mdc_kappa->Fill( mdc_kappa );
615 Mdc_dz->Fill( mdc_dz );
616 Mdc_tanlamda->Fill( mdc_tanlamda );
617 Mdc_err11->Fill( mdc_err11 );
618 Mdc_err21->Fill( mdc_err21 );
619 Mdc_err22->Fill( mdc_err22 );
620 Mdc_err31->Fill( mdc_err31 );
621 Mdc_err32->Fill( mdc_err32 );
622 Mdc_err33->Fill( mdc_err33 );
623 Mdc_err41->Fill( mdc_err41 );
624 Mdc_err42->Fill( mdc_err42 );
625 Mdc_err43->Fill( mdc_err43 );
626 Mdc_err44->Fill( mdc_err44 );
627 Mdc_err51->Fill( mdc_err51 );
628 Mdc_err52->Fill( mdc_err52 );
629 Mdc_err53->Fill( mdc_err53 );
630 Mdc_err54->Fill( mdc_err54 );
631 Mdc_err55->Fill( mdc_err55 );
632 }
633
634 for ( int i = 0; i < dedx->GetEntries(); i++ )
635 {
636
637 dedx->GetEntry( i );
638 Dedx_ptrk->Fill( dedx_ptrk );
639 Dedx_chie->Fill( dedx_chie );
640 Dedx_chimu->Fill( dedx_chimu );
641 Dedx_chipi->Fill( dedx_chipi );
642 Dedx_chik->Fill( dedx_chik );
643 Dedx_chip->Fill( dedx_chip );
644 Dedx_probPH->Fill( dedx_probPH );
645 Dedx_normPH->Fill( dedx_normPH );
646 Dedx_ghit->Fill( dedx_ghit );
647 Dedx_thit->Fill( dedx_thit );
648 }
649
650 for ( int i = 0; i < ext->GetEntries(); i++ )
651 {
652
653 ext->GetEntry( i );
654 Ext_tof1path->Fill( ext_tof1path );
655 Ext_tof1z->Fill( ext_tof1z );
656 Ext_tof1t->Fill( ext_tof1t );
657 Ext_tof1x->Fill( ext_tof1x );
658 Ext_tof1y->Fill( ext_tof1y );
659 Ext_tof2path->Fill( ext_tof2path );
660 Ext_tof2z->Fill( ext_tof2z );
661 Ext_tof2t->Fill( ext_tof2t );
662 Ext_tof2x->Fill( ext_tof2x );
663 Ext_tof2y->Fill( ext_tof2y );
664 Ext_emctheta->Fill( ext_emctheta );
665 Ext_emcphi->Fill( ext_emcphi );
666 Ext_emcpath->Fill( ext_emcpath );
667 Ext_mucz->Fill( ext_mucz );
668 Ext_muct->Fill( ext_muct );
669 Ext_mucx->Fill( ext_mucx );
670 Ext_mucy->Fill( ext_mucy );
671 }
672
673 for ( int i = 0; i < tof->GetEntries(); i++ )
674 {
675
676 tof->GetEntry( i );
677 Tof_path->Fill( tof_path );
678 Tof_zrhit->Fill( tof_zrhit );
679 Tof_ph->Fill( tof_ph );
680 Tof_tof->Fill( tof_tof );
681 Tof_errtof->Fill( tof_errtof );
682 Tof_beta->Fill( tof_beta );
683 Tof_texpe->Fill( tof_texpe );
684 Tof_texpmu->Fill( tof_texpmu );
685 Tof_texppi->Fill( tof_texppi );
686 Tof_texpka->Fill( tof_texpka );
687 Tof_texppro->Fill( tof_texppro );
688 Tof_toffsete->Fill( tof_toffsete );
689 Tof_toffsetmu->Fill( tof_toffsetmu );
690 Tof_toffsetpi->Fill( tof_toffsetpi );
691 Tof_toffsetka->Fill( tof_toffsetka );
692 Tof_toffsetpro->Fill( tof_toffsetpro );
693 Tof_toffsetatpr->Fill( tof_toffsetatpr );
694 Tof_sigmae->Fill( tof_sigmae );
695 Tof_sigmamu->Fill( tof_sigmamu );
696 Tof_sigmapi->Fill( tof_sigmapi );
697 Tof_sigmaka->Fill( tof_sigmaka );
698 Tof_sigmapro->Fill( tof_sigmapro );
699 Tof_sigmaatpr->Fill( tof_sigmaatpr );
700 Tof_quality->Fill( tof_quality );
701 Tof_t0->Fill( tof_t0 );
702 Tof_errt0->Fill( tof_errt0 );
703 Tof_errz->Fill( tof_errz );
704 Tof_phi->Fill( tof_phi );
705 Tof_errphi->Fill( tof_errphi );
706 Tof_energy->Fill( tof_energy );
707 Tof_errenergy->Fill( tof_errenergy );
708 }
709
710 for ( int i = 0; i < emc->GetEntries(); i++ )
711 {
712
713 emc->GetEntry( i );
714 Emc_x->Fill( emc_x );
715 Emc_y->Fill( emc_y );
716 Emc_z->Fill( emc_z );
717 Emc_theta->Fill( emc_theta );
718 Emc_phi->Fill( emc_phi );
719 Emc_dx->Fill( emc_dx );
720 Emc_dy->Fill( emc_dy );
721 Emc_dz->Fill( emc_dz );
722 Emc_dtheta->Fill( emc_dtheta );
723 Emc_dphi->Fill( emc_dphi );
724 Emc_energy->Fill( emc_energy );
725 Emc_dE->Fill( emc_dE );
726 Emc_eseed->Fill( emc_eseed );
727 Emc_e3x3->Fill( emc_e3x3 );
728 Emc_e5x5->Fill( emc_e5x5 );
729 Emc_secp->Fill( emc_secp );
730 Emc_latp->Fill( emc_latp );
731 }
732
733 for ( int i = 0; i < muc->GetEntries(); i++ )
734 {
735
736 muc->GetEntry( i );
737 Muc_depth->Fill( muc_depth );
738 Muc_chi2->Fill( muc_chi2 );
739 Muc_rms->Fill( muc_rms );
740 Muc_xpos->Fill( muc_xpos );
741 Muc_ypos->Fill( muc_ypos );
742 Muc_zpos->Fill( muc_zpos );
743 Muc_xsigma->Fill( muc_xsigma );
744 Muc_ysigma->Fill( muc_ysigma );
745 Muc_zsigma->Fill( muc_zsigma );
746 Muc_px->Fill( muc_px );
747 Muc_py->Fill( muc_py );
748 Muc_pz->Fill( muc_pz );
749 Muc_distance->Fill( muc_distance );
750 Muc_deltaphi->Fill( muc_deltaphi );
751 }
752
753 for ( int i = 0; i < runinfo->GetEntries(); i++ )
754 {
755
756 runinfo->GetEntry( i );
757 Runinfo_runNo->Fill( runinfo_runNo );
758 Runinfo_event->Fill( runinfo_event );
759 Runinfo_totcharged->Fill( runinfo_totcharged );
760 Runinfo_totneutral->Fill( runinfo_totneutral );
761 Runinfo_tottracks->Fill( runinfo_tottracks );
762 Runinfo_ngood->Fill( runinfo_ngood );
763 }
764
765 for ( int i = 0; i < mcpart->GetEntries(); i++ )
766 {
767
768 mcpart->GetEntry( i );
769 Mcpart_vx->Fill( mcpart_vx );
770 Mcpart_vy->Fill( mcpart_vy );
771 Mcpart_vz->Fill( mcpart_vz );
772 Mcpart_pxy->Fill( mcpart_pxy );
773 Mcpart_px->Fill( mcpart_px );
774 Mcpart_py->Fill( mcpart_py );
775 Mcpart_pz->Fill( mcpart_pz );
776 Mcpart_e->Fill( mcpart_e );
777 Mcpart_p->Fill( mcpart_p );
778 Mcpart_theta->Fill( mcpart_theta );
779 Mcpart_phi->Fill( mcpart_phi );
780 Mcpart_m->Fill( mcpart_m );
781 Mcpart_pro->Fill( mcpart_pro );
782 Mcpart_index->Fill( mcpart_index );
783 }
784
785 for ( int i = 0; i < pid_kal->GetEntries(); i++ )
786 {
787
788 pid_kal->GetEntry( i );
789 Pid_ptrk->Fill( pid_ptrk );
790 Pid_cost->Fill( pid_cost );
791 Pid_dedx->Fill( pid_dedx );
792 Pid_tof1->Fill( pid_tof1 );
793 Pid_tof2->Fill( pid_tof2 );
794 Prob_pion->Fill( prob_pion );
795 Prob_kaon->Fill( prob_kaon );
796 Prob_proton->Fill( prob_proton );
797 Kal_px->Fill( kal_px );
798 Kal_py->Fill( kal_py );
799 Kal_pz->Fill( kal_pz );
800 Kal_p->Fill( kal_p );
801 Kal_pxy->Fill( kal_pxy );
802 Kal_n->Fill( kal_n );
803 Kal_d0->Fill( kal_d0 );
804 Kal_phi0->Fill( kal_phi0 );
805 Kal_kappa->Fill( kal_kappa );
806 Kal_dz->Fill( kal_dz );
807 Kal_tanlamda->Fill( kal_tanlamda );
808 Kal_err11->Fill( kal_err11 );
809 Kal_err21->Fill( kal_err21 );
810 Kal_err22->Fill( kal_err22 );
811 Kal_err31->Fill( kal_err31 );
812 Kal_err32->Fill( kal_err32 );
813 Kal_err33->Fill( kal_err33 );
814 Kal_err41->Fill( kal_err41 );
815 Kal_err42->Fill( kal_err42 );
816 Kal_err43->Fill( kal_err43 );
817 Kal_err44->Fill( kal_err44 );
818 Kal_err51->Fill( kal_err51 );
819 Kal_err52->Fill( kal_err52 );
820 Kal_err53->Fill( kal_err53 );
821 Kal_err54->Fill( kal_err54 );
822 Kal_err55->Fill( kal_err55 );
823 }
824
825 Mdc_vx0->Write();
826 Mdc_vy0->Write();
827 Mdc_vz0->Write();
828 Mdc_pch->Write();
829 Mdc_trackid->Write();
830 Mdc_charge->Write();
831 Mdc_pxy0->Write();
832 Mdc_px0->Write();
833 Mdc_py0->Write();
834 Mdc_pz0->Write();
835 Mdc_r0->Write();
836 Mdc_phi->Write();
837 Mdc_theta->Write();
838 Mdc_ndof->Write();
839 Mdc_nster->Write();
840 Mdc_stat->Write();
841 Mdc_angd->Write();
842 Mdc_rvxy0->Write();
843 Mdc_rvz0->Write();
844 Mdc_d0->Write();
845 Mdc_phi0->Write();
846 Mdc_kappa->Write();
847 Mdc_dz->Write();
848 Mdc_tanlamda->Write();
849 Mdc_err11->Write();
850 Mdc_err21->Write();
851 Mdc_err22->Write();
852 Mdc_err31->Write();
853 Mdc_err32->Write();
854 Mdc_err33->Write();
855 Mdc_err41->Write();
856 Mdc_err42->Write();
857 Mdc_err43->Write();
858 Mdc_err44->Write();
859 Mdc_err51->Write();
860 Mdc_err52->Write();
861 Mdc_err53->Write();
862 Mdc_err54->Write();
863 Mdc_err55->Write();
864
865 Dedx_ptrk->Write();
866 Dedx_chie->Write();
867 Dedx_chimu->Write();
868 Dedx_chipi->Write();
869 Dedx_chik->Write();
870 Dedx_chip->Write();
871 Dedx_probPH->Write();
872 Dedx_normPH->Write();
873 Dedx_ghit->Write();
874 Dedx_thit->Write();
875
876 Ext_tof1path->Write();
877 Ext_tof1z->Write();
878 Ext_tof1t->Write();
879 Ext_tof1x->Write();
880 Ext_tof1y->Write();
881 Ext_tof2path->Write();
882 Ext_tof2z->Write();
883 Ext_tof2t->Write();
884 Ext_tof2x->Write();
885 Ext_tof2y->Write();
886 Ext_emctheta->Write();
887 Ext_emcphi->Write();
888 Ext_emcpath->Write();
889 Ext_mucz->Write();
890 Ext_muct->Write();
891 Ext_mucx->Write();
892 Ext_mucy->Write();
893
894 Tof_path->Write();
895 Tof_zrhit->Write();
896 Tof_ph->Write();
897 Tof_tof->Write();
898 Tof_errtof->Write();
899 Tof_beta->Write();
900 Tof_texpe->Write();
901 Tof_texpmu->Write();
902 Tof_texppi->Write();
903 Tof_texpka->Write();
904 Tof_texppro->Write();
905 Tof_toffsete->Write();
906 Tof_toffsetmu->Write();
907 Tof_toffsetpi->Write();
908 Tof_toffsetka->Write();
909 Tof_toffsetpro->Write();
910 Tof_toffsetatpr->Write();
911 Tof_sigmae->Write();
912 Tof_sigmamu->Write();
913 Tof_sigmapi->Write();
914 Tof_sigmaka->Write();
915 Tof_sigmapro->Write();
916 Tof_sigmaatpr->Write();
917 Tof_quality->Write();
918 Tof_t0->Write();
919 Tof_errt0->Write();
920 Tof_errz->Write();
921 Tof_phi->Write();
922 Tof_errphi->Write();
923 Tof_energy->Write();
924 Tof_errenergy->Write();
925
926 Emc_x->Write();
927 Emc_y->Write();
928 Emc_z->Write();
929 Emc_theta->Write();
930 Emc_phi->Write();
931 Emc_dx->Write();
932 Emc_dy->Write();
933 Emc_dz->Write();
934 Emc_dtheta->Write();
935 Emc_dphi->Write();
936 Emc_energy->Write();
937 Emc_dE->Write();
938 Emc_eseed->Write();
939 Emc_e3x3->Write();
940 Emc_e5x5->Write();
941 Emc_secp->Write();
942 Emc_latp->Write();
943
944 Muc_depth->Write();
945 Muc_chi2->Write();
946 Muc_rms->Write();
947 Muc_xpos->Write();
948 Muc_ypos->Write();
949 Muc_zpos->Write();
950 Muc_xsigma->Write();
951 Muc_ysigma->Write();
952 Muc_zsigma->Write();
953 Muc_px->Write();
954 Muc_py->Write();
955 Muc_pz->Write();
956 Muc_distance->Write();
957 Muc_deltaphi->Write();
958
959 Runinfo_runNo->Write();
960 Runinfo_event->Write();
961 Runinfo_totcharged->Write();
962 Runinfo_totneutral->Write();
963 Runinfo_tottracks->Write();
964 Runinfo_ngood->Write();
965
966 Mcpart_vx->Write();
967 Mcpart_vy->Write();
968 Mcpart_vz->Write();
969 Mcpart_pxy->Write();
970 Mcpart_px->Write();
971 Mcpart_py->Write();
972 Mcpart_pz->Write();
973 Mcpart_e->Write();
974 Mcpart_p->Write();
975 Mcpart_theta->Write();
976 Mcpart_phi->Write();
977 Mcpart_m->Write();
978 Mcpart_pro->Write();
979 Mcpart_index->Write();
980
981 Pid_ptrk->Write();
982 Pid_cost->Write();
983 Pid_dedx->Write();
984 Pid_tof1->Write();
985 Pid_tof2->Write();
986 Prob_pion->Write();
987 Prob_kaon->Write();
988 Prob_proton->Write();
989 Kal_px->Write();
990 Kal_py->Write();
991 Kal_pz->Write();
992 Kal_p->Write();
993 Kal_pxy->Write();
994 Kal_n->Write();
995 Kal_d0->Write();
996 Kal_phi0->Write();
997 Kal_kappa->Write();
998 Kal_dz->Write();
999 Kal_tanlamda->Write();
1000 Kal_err11->Write();
1001 Kal_err21->Write();
1002 Kal_err22->Write();
1003 Kal_err31->Write();
1004 Kal_err32->Write();
1005 Kal_err33->Write();
1006 Kal_err41->Write();
1007 Kal_err42->Write();
1008 Kal_err43->Write();
1009 Kal_err44->Write();
1010 Kal_err51->Write();
1011 Kal_err52->Write();
1012 Kal_err53->Write();
1013 Kal_err54->Write();
1014 Kal_err55->Write();
1015
1016 fin->Close();
1017 fout->Close();
1018}
Single_histo
void Single_histo(char *filein="Single.root", char *fileout="Single_histo.root")
Definition Single_histo.cxx:1