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