StRoot  1
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
StQABookHist.cxx
1 // $Id: StQABookHist.cxx,v 2.80 2019/03/14 02:31:53 genevb Exp $
2 // $Log: StQABookHist.cxx,v $
3 // Revision 2.80 2019/03/14 02:31:53 genevb
4 // Introduce iTPC plots
5 //
6 // Revision 2.79 2018/07/03 21:33:34 genevb
7 // Introduce EPD (code provided by J. Ewigleben)
8 //
9 // Revision 2.78 2018/05/02 21:07:40 genevb
10 // Initial accomodation for iTPC
11 //
12 // Revision 2.77 2017/02/25 03:24:30 genevb
13 // Run 17: remove HFT
14 //
15 // Revision 2.76 2016/05/25 03:59:44 genevb
16 // fix uninit members
17 //
18 // Revision 2.75 2016/05/13 22:04:49 genevb
19 // Address coverity findings: uninit vars, dead code, one PMD error, and one TOF error
20 //
21 // Revision 2.74 2016/02/19 03:52:15 genevb
22 // Expand track detector ID histograms
23 //
24 // Revision 2.73 2015/07/17 19:09:03 genevb
25 // SSD copied for SST, and HFT histogams use SST now too
26 //
27 // Revision 2.72 2015/04/02 19:53:47 genevb
28 // TPC dE/dx changes: Bethe-Bloch => Bichsel, and tighter cuts against pile-up tracks
29 //
30 // Revision 2.71 2015/01/16 21:08:28 genevb
31 // Initial versions of HFT histograms
32 //
33 // Revision 2.70 2014/07/22 20:39:28 genevb
34 // Add MTD to Offline QA
35 //
36 // Revision 2.69 2013/03/12 03:06:02 genevb
37 // Add FMS/FPD histograms for Run 13+
38 //
39 // Revision 2.68 2012/03/05 03:42:32 genevb
40 // Remove TPC XY dist, add TPC RPhi charge
41 //
42 // Revision 2.67 2012/02/09 03:01:05 genevb
43 // No FTPC histograms for Run 12+
44 //
45 // Revision 2.66 2012/02/08 22:10:35 genevb
46 // Updates for Run 12
47 //
48 // Revision 2.65 2011/05/31 21:35:49 genevb
49 // TPC request: add time bucket distribution of hits
50 //
51 // Revision 2.64 2009/03/19 01:08:08 genevb
52 // Show both xy and rphi TPC hit hists
53 //
54 // Revision 2.63 2009/03/13 19:27:24 genevb
55 // Now draw TPC xy hits in polar coords
56 //
57 // Revision 2.62 2009/03/09 16:34:54 genevb
58 // Tweak a couple histogram ranges
59 //
60 // Revision 2.61 2007/12/10 19:58:20 genevb
61 // Use log10 for number of secondary vertices
62 //
63 // Revision 2.60 2007/12/03 16:51:09 genevb
64 // Out-of-order statements for pre-run8 data
65 //
66 // Revision 2.59 2007/11/30 05:38:50 genevb
67 // Changes for Run8: mostly silicon removal, TOF addition
68 //
69 // Revision 2.58 2007/04/25 18:35:57 genevb
70 // Additional SSD hists
71 //
72 // Revision 2.57 2007/04/24 00:33:58 genevb
73 // SSD hists
74 //
75 // Revision 2.56 2007/04/12 20:39:48 genevb
76 // Cleanup (removal) of CalibVtx, Nfitpnt, Chisq1, Rich, histograms
77 //
78 // Revision 2.55 2007/04/07 04:40:30 genevb
79 // Remove fit pnts/tot; retitle log as log10
80 //
81 // Revision 2.54 2007/04/02 05:36:23 genevb
82 // PMD always in now
83 //
84 // Revision 2.53 2007/03/13 18:46:28 genevb
85 // Added Svt Laser Diff
86 //
87 // Revision 2.52 2007/02/26 20:45:01 genevb
88 // SVT drift hist
89 //
90 // Revision 2.51 2006/05/18 16:37:11 genevb
91 // Change FTPC hist ranges
92 //
93 // Revision 2.50 2005/02/22 19:38:39 genevb
94 // Do PMD hists only for year 4 and later (real data)
95 //
96 // Revision 2.49 2005/02/08 17:22:46 genevb
97 // PMD histo changes, handle estGlobal/ITTF tracks
98 //
99 // Revision 2.48 2005/01/27 05:28:25 genevb
100 // PMD changes
101 //
102 // Revision 2.47 2004/12/13 15:52:36 genevb
103 // Numerous updates: PMD, primtrk, FPD, QAShift lists
104 //
105 // Revision 2.46 2004/10/04 16:40:42 genevb
106 // FTPC radial histos
107 //
108 // Revision 2.45 2004/04/23 23:15:29 genevb
109 // Added signedDCA (Impact) plots for globals
110 //
111 // Revision 2.44 2004/02/12 05:03:10 genevb
112 // Year 4 AuAu changes. New SVT histos.
113 //
114 // Revision 2.43 2004/01/27 20:08:14 genevb
115 // Modify SVT point hist
116 //
117 // Revision 2.42 2004/01/10 01:10:18 genevb
118 // Preparations for Year 5, added some svt plots
119 //
120 // Revision 2.41 2003/09/02 17:59:21 perev
121 // gcc 3.2 updates + WarnOff
122 //
123 // Revision 2.40 2003/04/19 00:17:49 genevb
124 // Updated for dAu/pp running
125 //
126 // Revision 2.39 2003/02/20 20:09:54 genevb
127 // Several changes for new trigger scheme, dAu data
128 //
129 // Revision 2.38 2003/02/19 06:38:29 genevb
130 // Rework trigger and mult/event class sections
131 //
132 // Revision 2.37 2003/02/15 21:59:08 genevb
133 // dedx lable
134 //
135 // Revision 2.36 2003/01/24 15:24:48 genevb
136 // Further dAu changes
137 //
138 // Revision 2.35 2003/01/23 20:53:11 genevb
139 // Additional dAu changes
140 //
141 // Revision 2.34 2003/01/23 04:03:18 jeromel
142 // Include fixed
143 //
144 // Revision 2.33 2003/01/13 19:09:43 genevb
145 // Minor FTPC histo changes
146 //
147 // Revision 2.32 2002/07/11 12:26:22 genevb
148 // Some FTPC histogram changes
149 //
150 // Revision 2.31 2002/05/29 13:54:30 genevb
151 // Some changes to FTPC chisq histos
152 //
153 // Revision 2.30 2002/04/23 01:59:56 genevb
154 // Addition of BBC/FPD histos
155 //
156 // Revision 2.29 2002/02/12 18:41:59 genevb
157 // Additional FTPC histograms
158 //
159 // Revision 2.28 2002/01/21 22:09:24 genevb
160 // Include some ftpc histograms from StFtpcClusterMaker
161 //
162 // Revision 2.27 2001/12/28 09:19:13 genevb
163 // Adjustments for pp running
164 //
165 // Revision 2.26 2001/11/20 21:53:45 lansdell
166 // added x-y dist of hits, tpc east&west histos
167 //
168 // Revision 2.25 2001/11/13 21:17:56 lansdell
169 // changed limits on dedx error on mean, tpc
170 //
171 // Revision 2.24 2001/11/02 20:50:03 genevb
172 // Changed histogram ranges for momenta
173 //
174 // Revision 2.23 2001/10/25 20:45:12 lansdell
175 // changed ftpc # hits range to 25k
176 //
177 // Revision 2.22 2001/10/25 20:43:43 lansdell
178 // changed 2d ftpc globtrk histo ranges
179 //
180 // Revision 2.21 2001/10/10 05:00:19 genevb
181 // Changed ftpc point/track ranges
182 //
183 // Revision 2.20 2001/08/07 07:51:27 lansdell
184 // primvtx check for different multiplicities crashed for MC data, now fixed
185 //
186 // Revision 2.19 2001/08/03 23:31:59 lansdell
187 // added missing axis labels to 2d histos
188 //
189 // Revision 2.18 2001/08/03 20:33:55 lansdell
190 // added primvtx check histos for different multiplicities; separated x-y plot of first point on track, tpc into east and west histos
191 //
192 // Revision 2.17 2001/07/31 23:21:42 lansdell
193 // added last point, hit-helix histos
194 //
195 // Revision 2.16 2001/06/06 18:44:59 lansdell
196 // changed domain on some qa_shift histograms
197 //
198 // Revision 2.15 2001/05/29 23:23:06 lansdell
199 // removed impact param plots for FTPC from qa_shift list
200 //
201 // Revision 2.14 2001/05/25 17:20:09 lansdell
202 // squashed a small bug
203 //
204 // Revision 2.13 2001/05/25 16:31:21 lansdell
205 // more updates to qa shift histograms
206 //
207 // Revision 2.12 2001/05/24 01:48:13 lansdell
208 // qa_shift histograms updated
209 //
210 // Revision 2.11 2001/05/23 00:14:52 lansdell
211 // more changes for qa_shift histograms
212 //
213 // Revision 2.10 2001/05/16 20:57:03 lansdell
214 // new histograms added for qa_shift printlist; some histogram ranges changed; StMcEvent now used in StEventQA
215 //
216 // Revision 2.9 2001/05/02 16:10:47 lansdell
217 // changed some histogram limits
218 //
219 // Revision 2.8 2001/05/01 15:17:36 genevb
220 // Execute EMC code only if EMC libs loaded
221 //
222 // Revision 2.7 2001/04/28 22:05:13 genevb
223 // Added EMC histograms
224 //
225 // Revision 2.6 2001/04/26 20:45:19 lansdell
226 // changed some histogram ranges (TPC+SVT radius at start, impact param for primary tracks)
227 //
228 // Revision 2.5 2001/04/26 16:34:28 genevb
229 // Fixed some histogram ranges
230 //
231 // Revision 2.4 2001/04/25 21:35:26 genevb
232 // Added V0 phi distributions
233 //
234 // Revision 2.3 2001/04/24 22:53:51 lansdell
235 // Removed redundant radial position of first hit histograms
236 //
237 // Revision 2.2 2000/09/08 18:55:53 lansdell
238 // turned on FTPC primary track histograms
239 //
240 // Revision 2.1 2000/09/06 15:31:13 lansdell
241 // changed FTPC radial position of first hit histogram limits
242 //
243 // Revision 2.0 2000/08/25 16:02:40 genevb
244 // New revision: new structure, multiplicity classes
245 //
246 //
248 // //
249 // StQABookHist class for multiplicity-based QA histograms //
250 // //
252 
253 #include <Stiostream.h>
254 #include <math.h>
255 #include "QAH.h"
256 #include "StQABookHist.h"
257 #include "StQAMakerBase.h"
258 #include "StEmcUtil/others/StEmcMath.h"
259 #include "StEmcUtil/others/emcDetectorName.h"
260 #include "TROOT.h"
261 
262 
263 ClassImp(StQABookHist)
264 
265 //_____________________________________________________________________________
266 StQABookHist::StQABookHist(const char* type) : QAHistType(type) {
267 
268  histsSet = StQA_Undef;
269  silHists = kFALSE;
270  ftpHists = kFALSE;
271 
272  int i = 0;
273 
274 // - Set all the histogram booking constants
275 
276  nxpT = 50;
277  nyeta = 50;
278  xminpT = 0.0;
279  xmaxpT = 5.0;
280 // ymineta = -2.0;
281 // ymaxeta = 2.0;
282 
283  nchisq = 50;
284 // nmass = 40;
285 // ntau = 40;
286  ndedx = 50;
287  npnt = 50;
288  nleng = 50;
289  npsi = 36;
290  knpsi = 42;
291  nvrt = 100;
292  knyeta = 60;
293  knid = 10;
294  cnp = 50;
295  cndedx = 50;
296 
297  kminnid = 0.0;
298  kmaxnid = 10.0;
299  minpsi = 0.0;
300  kminpsi = -60.0;
301  maxpsi = 360.0;
302  minchisq = 0.;
303  maxchisq = 10.0;
304  minmass = 0.0;
305  maxmass = 2.0;
306  minpnt = 0.0;
307  maxpnt = 50.0;
308  minleng = 0.0;
309  maxleng = 200.0;
310  mintau = 0.0;
311  maxtau = 20.0;
312  mintrk = 0.0;
313  maxtrk = 8000.0;
314  minvrt = 2000.0;
315  maxvrt = 4000.0;
316  minmpt = 0.0;
317  maxmpt = 2.0;
318  minmeta = -0.2;
319  maxmeta = 0.2;
320  kmineta = -3.0;
321  kmaxeta = 3.0;
322  minxyz = 0.0;
323  maxxyz = 50.0;
324  cminp = 0.0;
325  cmaxp = 2.0;
326  cmindedx = 0.0;
327  cmaxdedx = 0.1e-04*1e6; // change from GeV to keV per cm
328 
329 
330 // - Zero all histogram pointers
331 
332  mNullPrimVtxClass = 0;
333 
334 // for method MakeGlob - from table globtrk
335  m_globtrk_fit_prob=0;
336  m_globtrk_tot=0;
337  m_globtrk_good=0;
338  m_globtrk_good_sm=0;
339  m_globtrk_goodTTS=0;
340  m_globtrk_goodF=0;
341  m_globtrk_iflag=0;
342  m_det_id=0;
343  m_dcaToBeamXY=0;
344  m_dcaToBeamZ1=0;
345  m_dcaToBeamZ2=0;
346  m_dcaToBeamZ3=0;
347  m_zDcaTanl=0;
348  m_zDcaZf=0;
349  m_zDcaPsi=0;
350  m_zDcaPhi0=0;
351 
352  m_pointT=0;
353  m_pointF=0;
354  m_pointFE=0;
355  m_pointFW=0;
356  m_max_pointT=0;
357  m_max_pointF=0;
358  m_max_pointFE=0;
359  m_max_pointFW=0;
360  m_fit_pointT=0;
361  m_fit_pointTTS=0;
362  m_glb_sptsTS=0;
363  m_glb_ratiomTTS=0;
364  m_glb_ratiomT=0;
365  m_glb_ratiomF=0;
366  m_glb_ratiomFE=0;
367  m_glb_ratiomFW=0;
368  m_glb_chargeT=0;
369  m_glb_chargeF=0;
370  m_glb_chargeFE=0;
371  m_glb_chargeFW=0;
372  m_glb_r0T=0;
373  m_glb_phi0T=0;
374  m_glb_z0T=0;
375  m_glb_curvT=0;
376  m_glb_padfT=0;
377  m_glb_padfTEW=0;
378  m_glb_xf0=0;
379  m_glb_xfT=0;
380  m_glb_xfF=0;
381  m_glb_xfFE=0;
382  m_glb_xfFW=0;
383  m_glb_yf0=0;
384  m_glb_yfT=0;
385  m_glb_yfF=0;
386  m_glb_yfFE=0;
387  m_glb_yfFW=0;
388  m_glb_zf0=0;
389  m_glb_zfT=0;
390  m_glb_zfF=0;
391  m_glb_zfFE=0;
392  m_glb_zfFW=0;
393  m_glb_planefF=0;
394  m_glb_f0=0;
395  m_glb_rzf0=0;
396  m_glb_rzl0=0;
397  m_glb_radfT=0;
398  m_glb_radfF=0;
399  m_glb_radfFE=0;
400  m_glb_radfFW=0;
401  m_glb_phifT=0;
402  m_glb_phifTS=0;
403  m_psiT=0;
404  m_psiTTS=0;
405  m_psiF=0;
406  m_psiFE=0;
407  m_psiFW=0;
408  m_tanlT=0;
409  m_glb_thetaT=0;
410  m_etaT=0;
411  m_etaTTS=0;
412  m_etaF=0;
413  m_etaFE=0;
414  m_etaFW=0;
415  m_pTT=0;
416  m_pTTTS=0;
417  m_pTF=0;
418  m_pTFE=0;
419  m_pTFW=0;
420  m_momT=0;
421  m_momF=0;
422  m_momFE=0;
423  m_momFW=0;
424  m_lengthT=0;
425  m_lengthF=0;
426  m_lengthFE=0;
427  m_lengthFW=0;
428  m_chisq0T=0;
429  m_chisq0TTS=0;
430  m_chisq0F=0;
431  m_chisq0FE=0;
432  m_chisq0FW=0;
433  m_chisq1F=0;
434  m_chisq1FE=0;
435  m_chisq1FW=0;
436  m_glb_impactT=0;
437  m_glb_simpactT=0;
438  m_glb_impactrT=0;
439  m_glb_impactTTS=0;
440  m_glb_impactrTTS=0;
441  m_glb_impactF=0;
442  m_glb_impactrF=0;
443 
444  m_pointTS=0;
445  m_max_pointTS=0;
446  m_fit_pointTS=0;
447  m_glb_ratiomTS=0;
448  m_glb_chargeTS=0;
449  m_glb_r0TS=0;
450  m_glb_phi0TS=0;
451  m_glb_z0TS=0;
452  m_glb_curvTS=0;
453  m_glb_xfTS=0;
454  m_glb_yfTS=0;
455  m_glb_zfTS=0;
456  m_glb_xf0TS=0;
457  m_glb_yf0TS=0;
458  m_glb_zf0TS=0;
459  m_glb_f0TS=0;
460  m_glb_rzf0TS=0;
461  m_glb_rzl0TS=0;
462  m_glb_radfTS=0;
463  m_psiTS=0;
464  m_tanlTS=0;
465  m_glb_thetaTS=0;
466  m_etaTS=0;
467  m_momTS=0;
468  m_pTTS=0;
469  m_lengthTS=0;
470  m_chisq0TS=0;
471  m_glb_impactTS=0;
472  m_glb_simpactTS=0;
473  m_glb_impactrTS=0;
474 
475  m_pT_eta_recT = 0;
476  m_pT_eta_recFE = 0;
477  m_pT_eta_recFW = 0;
478  m_globtrk_xf_yfTE = 0;
479  m_globtrk_xf_yfTW = 0;
480  m_globtrk_xf_yfFE = 0;
481  m_globtrk_xf_yfFW = 0;
482  m_globtrk_padtimeFE = 0;
483  m_globtrk_padtimeFW = 0;
484  m_tanl_zfT = 0;
485  m_mom_trklengthT = 0;
486  m_eta_trklengthT = 0;
487  m_eta_trklengthFE = 0;
488  m_eta_trklengthFW = 0;
489  m_npoint_lengthT = 0;
490  m_npoint_lengthFE = 0;
491  m_npoint_lengthFW = 0;
492  m_fpoint_lengthT = 0;
493  m_fpoint_lengthTTS = 0;
494  m_chisq0_momT = 0;
495  m_chisq0_etaT = 0;
496  m_chisq0_dipT = 0;
497  m_chisq0_zfT = 0;
498  m_chisq0_phiT = 0;
499  m_psi_phiT = 0;
500 
501  m_pT_eta_recTS= 0;
502  m_globtrk_xf_yfTS = 0;
503  m_tanl_zfTS = 0;
504  m_mom_trklengthTS = 0;
505  m_eta_trklengthTS = 0;
506  m_npoint_lengthTS = 0;
507  m_fpoint_lengthTS = 0;
508  m_chisq0_momTS = 0;
509  m_chisq0_etaTS = 0;
510  m_chisq0_dipTS = 0;
511  m_chisq0_zfTS = 0;
512  m_chisq0_phiTS = 0;
513  m_psi_phiTS = 0;
514 
515 // for method MakeDE - from table dst_dedx
516  m_ndedxr=0;
517 
518  m_ndedxT=0;
519  m_dedx0T=0;
520  m_dedx1T=0;
521  m_dedxTTS=0;
522  // east and west on same plot
523  m_ndedxF=0;
524  m_dedx0F=0;
525  // east and west on separate plots
526  m_ndedxFE=0;
527  m_dedx0FE=0;
528 
529  m_ndedxFW=0;
530  m_dedx0FW=0;
531 
532 // for method MakeHistPrim - from table primtrk
533  m_primtrk_tot=0;
534  m_primtrk_tot_sm=0;
535  m_primtrk_good=0;
536  m_primtrk_good_sm=0;
537  m_primtrk_goodTTS=0;
538  m_primtrk_goodF=0;
539  m_primtrk_iflag=0;
540  m_primglob_good=0;
541  m_primglob_fit=0;
542  m_pdet_id=0;
543  m_primtrk_meanptTTS=0;
544  m_primtrk_meanptF=0;
545  m_primtrk_meanetaTTS=0;
546  m_primtrk_meanetaF=0;
547 
548  m_ppointT=0;
549  m_ppointF=0;
550  m_ppointFE=0;
551  m_ppointFW=0;
552  m_pmax_pointT=0;
553  m_pmax_pointF=0;
554  m_pmax_pointFE=0;
555  m_pmax_pointFW=0;
556  m_pfit_pointT=0;
557  m_prim_ratiomT=0;
558  m_prim_ratioF=0;
559  m_prim_ratioFE=0;
560  m_prim_ratioFW=0;
561  m_prim_ratiomF=0;
562  m_prim_ratiomFE=0;
563  m_prim_ratiomFW=0;
564  m_prim_chargeT=0;
565  m_prim_chargeF=0;
566  m_prim_chargeFE=0;
567  m_prim_chargeFW=0;
568  m_prim_r0T=0;
569  m_prim_phi0T=0;
570  m_prim_z0T=0;
571  m_prim_curvT=0;
572  m_prim_xf0=0;
573  m_prim_xfT=0;
574  m_prim_xfF=0;
575  m_prim_xfFE=0;
576  m_prim_xfFW=0;
577  m_prim_yf0=0;
578  m_prim_yfT=0;
579  m_prim_yfF=0;
580  m_prim_yfFE=0;
581  m_prim_yfFW=0;
582  m_prim_zf0=0;
583  m_prim_zfT=0;
584  m_prim_zfF=0;
585  m_prim_zfFE=0;
586  m_prim_zfFW=0;
587  m_prim_f0=0;
588  m_prim_rzf0=0;
589  m_prim_rzl0=0;
590  m_prim_radfT=0;
591  m_prim_radfF=0;
592  m_prim_radfFE=0;
593  m_prim_radfFW=0;
594  m_ppsiT=0;
595  m_ppsiTTS=0;
596  m_ppsiF=0;
597  m_ppsiFE=0;
598  m_ppsiFW=0;
599  m_ptanlT=0;
600  m_prim_thetaT=0;
601  m_petaT=0;
602  m_petaTTS=0;
603  m_petaF=0;
604  m_petaFE=0;
605  m_petaFW=0;
606  m_ppTT=0;
607  m_ppTTTS=0;
608  m_ppTF=0;
609  m_ppTFE=0;
610  m_ppTFW=0;
611  m_pmomT=0;
612  m_pmomTS=0;
613  m_pmomF=0;
614  m_pmomFE=0;
615  m_pmomFW=0;
616  m_plengthT=0;
617  m_plengthF=0;
618  m_plengthFE=0;
619  m_plengthFW=0;
620  m_pchisq0T=0;
621  m_pchisq0TTS=0;
622  m_pchisq0F=0;
623  m_pchisq0FE=0;
624  m_pchisq0FW=0;
625  m_pchisq1F=0;
626  m_pchisq1FE=0;
627  m_pchisq1FW=0;
628 
629  m_ppointTS=0;
630  m_pmax_pointTS=0;
631  m_pfit_pointTS=0;
632  m_prim_ratiomTS=0;
633  m_prim_chargeTS=0;
634  m_prim_r0TS=0;
635  m_prim_phi0TS=0;
636  m_prim_z0TS=0;
637  m_prim_curvTS=0;
638  m_prim_xfTS=0;
639  m_prim_yfTS=0;
640  m_prim_zfTS=0;
641  m_prim_xf0TS=0;
642  m_prim_yf0TS=0;
643  m_prim_zf0TS=0;
644  m_prim_f0TS=0;
645  m_prim_rzf0TS=0;
646  m_prim_rzl0TS=0;
647  m_prim_radfTS=0;
648  m_ppsiTS=0;
649  m_ptanlTS=0;
650  m_prim_thetaTS=0;
651  m_petaTS=0;
652  m_momTS=0;
653  m_ppTTS=0;
654  m_plengthTS=0;
655  m_pchisq0TS=0;
656 
657  m_ppT_eta_recT = 0;
658  m_ppT_eta_recFE = 0;
659  m_ppT_eta_recFW = 0;
660  m_primtrk_xf_yfTE = 0;
661  m_primtrk_xf_yfTW = 0;
662  m_primtrk_xf_yfFE = 0;
663  m_primtrk_xf_yfFW = 0;
664  m_ptanl_zfT = 0;
665  m_pmom_trklengthT = 0;
666  m_peta_trklengthT = 0;
667  m_peta_trklengthFE = 0;
668  m_peta_trklengthFW = 0;
669  m_pnpoint_lengthT = 0;
670  m_pnpoint_lengthFE = 0;
671  m_pnpoint_lengthFW = 0;
672  m_pfpoint_lengthT = 0;
673  m_pfpoint_lengthTTS = 0;
674  m_pchisq0_momT = 0;
675  m_pchisq0_etaT = 0;
676  m_pchisq0_dipT = 0;
677  m_pchisq0_zfT = 0;
678  m_ppsi_phiT = 0;
679 
680  m_ppT_eta_recTS= 0;
681  m_primtrk_xf_yfTS = 0;
682  m_ptanl_zfTS = 0;
683  m_pmom_trklengthTS = 0;
684  m_peta_trklengthTS = 0;
685  m_pnpoint_lengthTS = 0;
686  m_pfpoint_lengthTS = 0;
687  m_pchisq0_momTS = 0;
688  m_pchisq0_etaTS = 0;
689  m_pchisq0_dipTS = 0;
690  m_pchisq0_zfTS = 0;
691  m_ppsi_phiTS = 0;
692 
693  // for MakeHistPID - from tables primtrk & dst_dedx
694  m_p_dedx_rec=0;
695 
696  // for method MakeHistVertex - from table dst_vertex
697  m_v_num=0;
698  m_v_num_sm=0;
699  m_v_vtxid=0;
700  m_v_x=0;
701  m_v_y=0;
702  m_v_z=0;
703  m_v_pchi2=0;
704  m_v_r=0;
705 
706  m_pv_vtxid=0;
707  m_pv_x=0;
708  m_pv_y=0;
709  m_pv_z=0;
710  m_pv_xy=0;
711  m_pv_pchi2=0;
712  m_pv_r=0;
713 
714  m_vtx_phi_dist=0;
715  m_vtx_r_dist=0;
716  m_vtx_z_dist=0;
717 
718  m_v0 =0;
719  m_ev0_lama_hist=0;
720  m_ev0_k0ma_hist=0;
721  m_xi_tot=0;
722  m_xi_ma_hist=0;
723  m_kink_tot=0;
724 
725  m_vtx_FtpcEastTpc_xy=0;
726  m_vtx_FtpcEastTpc_z=0;
727  m_vtx_FtpcWestTpc_xy=0;
728  m_vtx_FtpcWestTpc_z=0;
729 
730 
731 // for method MakeHistPoint
732  m_z_hits=0;
733  m_pnt_tot=0;
734  m_pnt_tot_med=0;
735  m_pnt_tot_sm=0;
736  m_pnt_id=0;
737  m_pnt_xyS=0;
738  m_pnt_rpTW=0;
739  m_pnt_rpTE=0;
740  m_pnt_phiT=0;
741  m_pnt_padrowT=0;
742  m_pnt_timeT=0;
743  m_pnt_zS=0;
744  m_pnt_phiS=0;
745  m_pnt_barrelS=0;
746  m_pnt_planeF=0;
747  m_pnt_tpc=0;
748  m_pnt_xyFE=0;
749  m_pnt_xyFW=0;
750  m_pnt_padtimeFE=0;
751  m_pnt_padtimeFW=0;
752  // east and west on same plot
753  m_pnt_ftpc=0;
754  // east and west on separate plots
755  m_pnt_ftpcE=0;
756  m_pnt_ftpcW=0;
757 
758 // for method MakeHistEval
759  m_geant_reco_pvtx_x=0;
760  m_geant_reco_pvtx_y=0;
761  m_geant_reco_pvtx_z=0;
762  m_geant_reco_vtx_z_z=0;
763 
764 // for EMC hits
765  m_emc_nhit=0;
766  m_emc_etot=0;
767  for (i=0; i<4; i++) {
768  m_emc_hits[i]=0;
769  m_emc_energy2D[i]=0;
770  m_emc_adc[i]=0;
771  m_emc_energy[i]=0;
772  }
773 
774 // for EMC-BSMD hits
775  for (i=0; i<8; i++) {
776  m_emc_hits_per_module[i]=0;
777  m_emc_energy_per_module[i]=0;
778  m_emc_strip_hits_per_module[i]=0;
779  m_emc_strip_energy_per_module[i]=0;
780  }
781 
782 // for EMC cluster finder
783  m_emc_ncl=0;
784  m_emc_etotCl=0;
785  m_emc_sig_e=0;
786  m_emc_sig_p=0;
787  for (i=0; i<4; i++) {
788  m_emc_cl[i]=0;
789  m_emc_energyCl[i]=0;
790  m_emc_HitsInCl[i]=0;
791  m_emc_EnergyCl[i]=0;
792  m_emc_EtaInCl[i]=0;
793  m_emc_PhiInCl[i]=0;
794  }
795 
796 // for EMC points
797  for (i=0; i<4; i++) {
798  m_emc_point_energy[i]=0;
799  m_emc_point_eta[i]=0;
800  m_emc_point_phi[i]=0;
801  m_emc_point_sigeta[i]=0;
802  m_emc_point_sigphi[i]=0;
803  m_emc_point_deleta[i]=0;
804  m_emc_point_delphi[i]=0;
805  m_emc_point_trmom[i]=0;
806  m_emc_points[i]=0;
807  }
808  m_emc_point_flag=0;
809 
810 // for BBC
811  for (i=0; i<4; i++) {
812  m_bbc_adc[i] = 0;
813  m_bbc_tdc[i] = 0;
814  }
815 
816 // for FPD
817  for (i=0; i<2; i++) {
818  m_fpd_top[i] = 0;
819  m_fpd_bottom[i] = 0;
820  m_fpd_south[i] = 0;
821  m_fpd_north[i] = 0;
822  }
823  for (i=0; i<8; i++) m_fpd_sums[i] = 0;
824 
825 // for PMD
826  for (i=0; i<12; i++) {
827  m_pmd_sm_adc[i] = 0;
828  m_pmd_sm_hit[i] = 0;
829  }
830  for (i=0; i<24; i++) {
831  m_pmd_chain_adc[i] = 0;
832  m_pmd_chain_hit[i] = 0;
833  }
834  m_pmd_total_hit = 0;
835  m_pmd_total_adc = 0;
836  m_cpv_total_hit = 0;
837  m_cpv_total_adc = 0;
838 
839 // for SVT
840  m_pnt_svt=0;
841  m_pnt_svtLaser=0;
842  m_pnt_svtLaserDiff=0;
843  m_svt_loc=0;
844 
845 // for SSD
846  m_pnt_ssd=0;
847  m_pnt_phiSSD = 0;
848  m_pnt_lwSSD = 0;
849  m_glb_ssd_phi = 0;
850  m_prim_ssd_phi = 0;
851  m_pnt_sizeSSD = 0;
852  m_pnt_eSSD = 0;
853 
854 // for SST
855  m_pnt_sst=0;
856  m_pnt_phiSST = 0;
857  m_pnt_lwSST = 0;
858  m_glb_sst_phi = 0;
859  m_prim_sst_phi = 0;
860  m_pnt_sizeSST = 0;
861  m_pnt_eSST = 0;
862 
863  // HFT histograms
864  // correlation plots: hits in HFT subsystems
865  m_nhit_Pxl_Ist = 0;
866  m_nhit_Pxl_Sst = 0;
867  m_nhit_Ist_Sst = 0;
868  m_global_hft_hit = 0;
869  m_primary_hft_hit = 0;
870 
871  // Hists for PIXEL
872  m_pxl_hit_phi_z_Pxl1 = 0;
873  m_pxl_hit_phi_z_Pxl2 = 0;
874  m_pxl_hit_ladder = 0;
875  m_pxl_hit_sector_sensor_Pxl1 = 0;
876  m_pxl_hit_sector_sensor_Pxl2 = 0;
877  m_pxl_nhit_Pxl1_tpc_mult = 0;
878  m_pxl_nhit_Pxl2_tpc_mult = 0;
879  m_pxl_nhit_Pxl1_tof_mult = 0;
880  m_pxl_nhit_Pxl2_tof_mult = 0;
881  m_pxl_nhit_Pxl1_Pxl2 = 0;
882  m_global_pxl_hit = 0;
883  m_primary_pxl_hit = 0;
884 
885  // Hists for IST
886  m_ist_hit_phi_z = 0;
887  m_ist_hit_ladder = 0;
888  m_ist_hit_ladder_sensor = 0;
889  m_ist_nhit_tpc_mult = 0;
890  m_ist_nhit_tof_mult = 0;
891  m_global_ist_hit = 0;
892  m_primary_ist_hit = 0;
893 
894 // for TOF
895  m_tof_hit_tray=0;
896  m_tof_hit_module=0;
897  m_tof_match_tray=0;
898  m_tof_match_module=0;
899  m_tof_vpd_hit=0;
900  m_tof_vtx_z=0;
901  m_tof_PID=0;
902 
903 // for MTD
904  m_MtdNHits = 0;
905  m_MtdHitMap = 0;
906  m_MtdNMatchHits = 0;
907  m_MtdMatchHitMap = 0;
908 
909 // for EPD
910  for (int i=0; i<24; i++) {
911  m_epd_adc[i] = 0;
912  m_epd_tac[i] = 0;
913  }
914 
915 
916 // for iTPC
917  for(int i = 0; i < 24; i++) {
918  m_TPC_ch_time_inner[i] = 0;
919  m_TPC_ch_time_outer[i] = 0;
920  m_TPC_ch_nrow[i] = 0;
921  }
922  m_TPC_adc_sec_inner = 0;
923  m_TPC_adc_sec_outer = 0;
924 
925 
926 }
927 //_____________________________________________________________________________
928 void StQABookHist::BookHist(Int_t hSet){
929 
930  histsSet = hSet;
931 
932  QAH::preString = QAHistType;
933 //book histograms --------------
934  if (histsSet != StQA_MC) {
935  mNullPrimVtxClass = QAH::H1F("QaNullPrimVtxMult","event primary vertex check",40,-2,2);
936  mNullPrimVtxClass->SetXTitle("has primary vertex? (yes = 1, no = -1)");
937  mNullPrimVtxClass->SetYTitle("# of events");
938  }
939  if (histsSet<StQA_run8) silHists = kTRUE;
940  if (histsSet<StQA_run12all) ftpHists = kTRUE;
941  BookHistPoint();
942  BookHistEMC();
943  if (histsSet == StQA_AuAuOld) {
944  BookHistBBC();
945  BookHistFPD();
946  }
947  BookHistGlob();
948  BookHistPrim();
949  BookHistDE();
950  BookHistPID();
951  BookHistVertex();
952  if (histsSet<StQA_run12all) BookHistPMD();
953  if (histsSet==StQA_MC) BookHistEval();
954  if (histsSet>=StQA_run8) BookHistTOF();
955  if (histsSet>=StQA_run12all) BookHistMTD();
956  if (histsSet>=StQA_run14 && histsSet<StQA_run17) {
957  BookHistPXL();
958  BookHistIST();
959  BookHistHFT();
960  }
961  if (histsSet>=StQA_run18) {
962  BookHistEPD();
963  BookHistiTPC();
964  }
965 
966 
967 }
968 //_____________________________________________________________________________
969 void StQABookHist::BookHistGlob(){
970 
971 // for method MakeGlob - from table globtrk
972 
973 // general
974 
975  m_globtrk_tot = QAH::H1F("QaGtrkTot","globtrk: tot num tracks - all",
976  40,0.,12500.);
977  m_globtrk_iflag = QAH::H1F("QaGtrkFlag","globtrk: iflag - all ",200,-999.,1001.);
978  m_globtrk_good = QAH::H1F("QaGtrkGood","globtrk: tot good tracks - all",40,0.,10000.);
979  m_globtrk_good_sm = QAH::H1F("QaGtrkGoodsm","globtrk: tot good tracks - all",40,0.,1000.);
980  if (silHists)
981  m_globtrk_goodTTS = QAH::H1F("QaGtrkGoodTTS","globtrk: tot good tracks - tpc,tpc+svt",150,0.,9000.);
982  m_globtrk_goodF = QAH::H2F("QaGtrkGoodF","globtrk: tot good tracks - ftpc",150,0.,1500.,150,0.,1500.);
983  m_globtrk_goodF->SetXTitle("FTPC East");
984  m_globtrk_goodF->SetYTitle("FTPC West");
985  m_globtrk_fit_prob = QAH::H1F("QaGtrkFitProb","globtrk: prob. fit is correct",100,0,1.2);
986  m_det_id = QAH::H1F("QaGtrkDetId","globtrk: Detector ID good tracks - all",48,-0.5,47.5);
987  m_dcaToBeamXY = QAH::H2F("QaGtrkDcaBeamXY","globtrk: xy-DCA to Beam Axis (z=0)",80,-4,4,80,-4,4);
988  m_dcaToBeamXY->SetXTitle("x");
989  m_dcaToBeamXY->SetYTitle("y");
990  m_dcaToBeamZ1 = QAH::H1F("QaGtrkDcaBeamZ1","globtrk: z-DCA to Beam Axis, coarse scale",100,-200,200);
991  m_dcaToBeamZ2 = QAH::H1F("QaGtrkDcaBeamZ2","globtrk: z-DCA to Beam Axis",100,-50,50);
992  m_dcaToBeamZ3 = QAH::H1F("QaGtrkDcaBeamZ3","globtrk: z-DCA to Beam Axis, near z=0",80,-20,20);
993  m_zDcaTanl = QAH::H2F("QaGtrkZdcaTanl","globtrk: z-DCA to Beam Axis vs tanl, tpc",100,-25,25,32,-4,4);
994  m_zDcaTanl->SetXTitle("z-DCA");
995  m_zDcaTanl->SetYTitle("tanl");
996  m_zDcaZf = QAH::H2F("QaGtrkZdcaZf","globtrk: z-DCA to Beam Axis vs z-first",100,-25,25,50,-300,300);
997  m_zDcaZf->SetXTitle("z-DCA");
998  m_zDcaZf->SetYTitle("z of first point");
999  m_zDcaPsi = QAH::H2F("QaGtrkZdcaPsi","globtrk: z-DCA to Beam Axis vs psi",100,-25,25,64,0,360);
1000  m_zDcaPsi->SetXTitle("z-DCA");
1001  m_zDcaPsi->SetYTitle("psi");
1002  m_zDcaPhi0 = QAH::H2F("QaGtrkZdcaPhi0","globtrk: z-DCA to Beam Axis vs azimuth (phi0) at start",80,-20,20,64,0,360);
1003  m_zDcaPhi0->SetXTitle("z-DCA");
1004  m_zDcaPhi0->SetYTitle("phi0");
1005  if (silHists) {
1006  m_glb_sptsTS = QAH::H1F("QaGtrkSptsTS","globtrk: number of svt points, tpc+svt",4,0.,4.);
1007  m_glb_ratiomTTS = QAH::MH1F("QaGtrkRnmfTTS","globtrk: ratio Nfit/max pnt, tpc,svt", 55,0.,1.1,2);
1008  m_glb_ratiomTTS->SetMinimum(10);
1009  m_glb_ratiomTTS->Rebin(0,"TPC+SVT");
1010  m_glb_ratiomTTS->Rebin(1,"TPC");
1011  m_glb_ratiomTTS->SetStats(kFALSE);
1012  m_psiTTS = QAH::MH1F("QaGtrkPsiTTS", "globtrk: psi (deg), tpc,svt", 64, 0.,360.,2);
1013  m_psiTTS->Rebin(0,"TPC+SVT");
1014  m_psiTTS->Rebin(1,"TPC");
1015  m_psiTTS->SetStats(kFALSE);
1016  m_etaTTS = QAH::MH1F("QaGtrkEtaTTS", "globtrk: eta, tpc,svt",40,-2.,2.,2);
1017  m_etaTTS->Rebin(0,"TPC+SVT");
1018  m_etaTTS->Rebin(1,"TPC");
1019  m_etaTTS->SetStats(kFALSE);
1020  m_pTTTS = QAH::MH1F("QaGtrkPtTTS", "globtrk: log10 pT, tpc,svt",50,-1.,4.,2);
1021  m_pTTTS->Rebin(0,"TPC+SVT");
1022  m_pTTTS->Rebin(1,"TPC");
1023  m_pTTTS->SetStats(kFALSE);
1024  m_pTTTS->SetXTitle("log10 pT (MeV)");
1025  m_chisq0TTS = QAH::MH1F("QaGtrkChisq0TTS", "globtrk: chisq0, tpc,svt", 50, 0.,5.,2);
1026  m_chisq0TTS->Rebin(0,"TPC+SVT");
1027  m_chisq0TTS->Rebin(1,"TPC");
1028  m_chisq0TTS->SetStats(kFALSE);
1029  m_fit_pointTTS = QAH::MH1F("QaGtrkNPntFitTTS","globtrk: N fit pnts on trk, tpc,svt",50,0.,50.,2);
1030  m_fit_pointTTS->Rebin(0,"TPC+SVT");
1031  m_fit_pointTTS->Rebin(1,"TPC");
1032  m_fit_pointTTS->SetStats(kFALSE);
1033  m_glb_impactTTS = QAH::MH1F("QaGtrkImpactTTS", "globtrk: log10 impact param from prim vtx, tpc,svt",120,-3.,3.,2);
1034  m_glb_impactTTS->Rebin(0,"TPC+SVT");
1035  m_glb_impactTTS->Rebin(1,"TPC");
1036  m_glb_impactTTS->SetStats(kFALSE);
1037  m_glb_impactrTTS = QAH::MH1F("QaGtrkImpactrTTS", "globtrk: impact param from prim vtx, tpc,svt",100,0.,300.,2);
1038  m_glb_impactrTTS->Rebin(0,"TPC+SVT");
1039  m_glb_impactrTTS->Rebin(1,"TPC");
1040  m_glb_impactrTTS->SetStats(kFALSE);
1041  m_fpoint_lengthTTS = QAH::H2F("QaGtrkFitPntLTTS","globtrk: N fit pnts vs length, tpc,tpc+svt",
1042  25,0.,250.,25,0.,50.);
1043  m_fpoint_lengthTTS->SetXTitle("trk length");
1044  m_fpoint_lengthTTS->SetYTitle("Npoints on trk");
1045  }
1046 
1047 // 1D tpc
1048 
1049  m_pointT = QAH::H1F("QaGtrkNPntT", "globtrk: N points on trk,tpc", 81,-0.5, 80.5);
1050  m_max_pointT = QAH::H1F("QaGtrkNPntMaxT","globtrk: N max pnts on trk, tpc", 76, -1, 151);
1051  m_fit_pointT = QAH::H1F("QaGtrkNPntFitT","globtrk: N fit pnts on trk, tpc", 81, -0.5, 80.5);
1052  m_glb_ratiomT = QAH::H1F("QaGtrkRnmT", "globtrk: ratio Nfit/max pnt, tpc", 55, 0., 1.1);
1053  m_glb_chargeT = QAH::H1F("QaGtrkChrgT", "globtrk: charge, tpc ", 20,-2.,2.);
1054  m_glb_r0T = QAH::H1F("QaGtrkR0T", "globtrk: radius at start (cm), tpc ", 50,0.,200.);
1055  m_glb_r0T->SetMinimum(100);
1056  m_glb_phi0T = QAH::H1F("QaGtrkPhi0T", "globtrk: azimuth (phi) at start (deg,force 0,360), tpc ", 64,0.,360.);
1057  m_glb_z0T = QAH::H1F("QaGtrkZ0T", "globtrk: z-coord at start (cm), tpc ", 50, -300.,300.);
1058  m_glb_z0T->SetMinimum(500);
1059  m_glb_curvT = QAH::H1F("QaGtrkCurvT", "globtrk: log10 curvature (1/cm), tpc ", 80,-3.5,0.5);
1060  m_glb_padfTEW = QAH::MH1F("QaGtrkPadfTEW","globtrk: padrow of first hit on trk, tpc", 72,0.5,72.5,2);
1061  m_glb_padfTEW->Rebin(0,"East");
1062  m_glb_padfTEW->Rebin(1,"West");
1063  m_glb_padfTEW->SetStats(kFALSE);
1064  m_glb_padfT = QAH::H1F("QaGtrkPadfT", "globtrk: padrow of first hit on trk, tpc", 72,0.5,72.5);
1065  m_glb_f0 = QAH::MH1F("QaGtrkf0", "globtrk: first point: hit-helix, tpc",60,-3.,3.,3);
1066  m_glb_f0->Rebin(0,"X");
1067  m_glb_f0->Rebin(1,"Y");
1068  m_glb_f0->Rebin(2,"Z");
1069  m_glb_f0->SetStats(kFALSE);
1070  m_glb_xfT = QAH::H1F("QaGtrkXfT", "globtrk: x of first hit on trk, tpc", 50,-200.,200.);
1071  m_glb_xf0 = QAH::H1F("QaGtrkXf0", "globtrk: first point: x_hit - x_helix, tpc",60,-3.,3.);
1072  m_glb_yfT = QAH::H1F("QaGtrkYfT", "globtrk: y of first hit on trk, tpc", 50,-200.,200.);
1073  m_glb_yf0 = QAH::H1F("QaGtrkYf0", "globtrk: first point: y_hit - y_helix, tpc",60,-3.,3.);
1074  m_glb_zfT = QAH::H1F("QaGtrkZfT", "globtrk: z of first hit on trk, tpc", 50,-300.,300.);
1075  m_glb_zf0 = QAH::H1F("QaGtrkZf0", "globtrk: first point: z_hit - z_helix, tpc",60,-3.,3.);
1076  m_glb_radfT = QAH::H1F("QaGtrkRT", "globtrk: radial position of first hit, tpc", 50,0.,200.);
1077  m_glb_rzf0 = QAH::MH1F("QaGtrkRZf0", "globtrk: first point: hit - helix, tpc",60,-3.,3.,2);
1078  m_glb_rzf0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1079  m_glb_rzf0->Rebin(1,"z_{dif}");
1080  m_glb_rzf0->SetStats(kFALSE);
1081  m_glb_rzl0 = QAH::MH1F("QaGtrkRZl0", "globtrk: last point: hit - helix, tpc",60,-3.,3.,2);
1082  m_glb_rzl0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1083  m_glb_rzl0->Rebin(1,"z_{dif}");
1084  m_glb_rzl0->SetStats(kFALSE);
1085  m_glb_phifT = QAH::MH1F("QaGtrkPhifT", "globtrk: phi of first point on trk, tpc",64,0,360,2);
1086  m_glb_phifT->Rebin(0,"East");
1087  m_glb_phifT->Rebin(1,"West");
1088  m_glb_phifT->SetStats(kFALSE);
1089  m_lengthT = QAH::H1F("QaGtrkLengthT", "globtrk: track length, tpc", 50,0.,300.);
1090  m_psiT = QAH::H1F("QaGtrkPsiT", "globtrk: psi, tpc (deg)", 64, 0.,360.);
1091  m_tanlT = QAH::H1F("QaGtrkTanlT", "globtrk: tanl, tpc",32,-4.,4.);
1092  m_glb_thetaT = QAH::H1F("QaGtrkThetaT", "globtrk: theta (degrees), tpc",36,0.,180.);
1093  m_etaT = QAH::H1F("QaGtrkEtaT", "globtrk: eta, tpc",40,-2.,2.);
1094  m_pTT = QAH::H1F("QaGtrkPtT", "globtrk: pT, tpc",50,0.,10.);
1095  m_momT = QAH::H1F("QaGtrkPT", "globtrk: momentum, tpc",50,0.,10.);
1096  m_chisq0T = QAH::H1F("QaGtrkChisq0T", "globtrk: chisq0, tpc", 50, 0.,5.);
1097  m_glb_impactT = QAH::MH1F("QaGtrkImpactT", "globtrk: log10 impact param from prim vtx, tpc",
1098  120,-3.0,3.0,3);
1099  m_glb_impactT->Rebin(0,"East");
1100  m_glb_impactT->Rebin(1,"West");
1101  m_glb_impactT->Rebin(2,"All");
1102  m_glb_impactT->SetStats(kFALSE);
1103  m_glb_simpactT = QAH::MH1F("QaGtrkSImpactT", "globtrk: signed impact param from prim vtx, tpc",
1104  50,-0.8,0.8,3);
1105  m_glb_simpactT->Rebin(0,"East");
1106  m_glb_simpactT->Rebin(1,"West");
1107  m_glb_simpactT->Rebin(2,"All");
1108  m_glb_simpactT->SetStats(kFALSE);
1109  m_glb_impactrT = QAH::H1F("QaGtrkImpactrT", "globtrk: impact param from prim vtx, tpc",100,0.,300.);
1110 
1111 
1112 // 2D - tpc
1113 
1114  m_pT_eta_recT = QAH::H2F("QaGtrkPtVsEtaT","globtrk: log10 pT vs eta, tpc", 20,-2.,2.,40,1.,4.);
1115  m_pT_eta_recT->SetXTitle("eta");
1116  m_pT_eta_recT->SetYTitle(" log10 pT (MeV)");
1117 
1118  m_globtrk_xf_yfTE = QAH::H2F("QaGtrkXfYfTE", "globtrk: Y vs X of first hit on trk, tpc east", 40,-200.,200.,40,-200.,200.);
1119  m_globtrk_xf_yfTE->SetXTitle("x first");
1120  m_globtrk_xf_yfTE->SetYTitle("y first");
1121 
1122  m_globtrk_xf_yfTW = QAH::H2F("QaGtrkXfYfTW", "globtrk: Y vs X of first hit on trk, tpc west", 40,-200.,200.,40,-200.,200.);
1123  m_globtrk_xf_yfTW->SetXTitle("x first");
1124  m_globtrk_xf_yfTW->SetYTitle("y first");
1125 
1126  m_tanl_zfT = QAH::H2F("QaGtrkTanlzf","globtrk: tanl(dip) vs. (zfirst-zvtx)/arc length, tpc,tpc+svt",60,-3.,3.,60,-3.,3.);
1127  m_tanl_zfT->SetXTitle("(zfirst-zvtx)/arc length");
1128  m_tanl_zfT->SetYTitle("tanl");
1129 
1130  m_mom_trklengthT = QAH::H2F("QaGtrkPVsTrkLength","globtrk: log10 mom vs trk length, tpc",
1131  50,0.,250.,40,1.,4.);
1132  m_mom_trklengthT->SetXTitle("trk length");
1133  m_mom_trklengthT->SetYTitle("log10 P (MeV)");
1134 
1135  m_eta_trklengthT = QAH::H2F("QaGtrkLengthVEtaT","globtrk: trk length vs eta, tpc",
1136  20,-2.,2.,50,0.,250.);
1137  m_eta_trklengthT->SetXTitle("eta");
1138  m_eta_trklengthT->SetYTitle("length");
1139 
1140  m_npoint_lengthT = QAH::H2F("QaGtrkNPntLengthT","globtrk: N pnts vs length, tpc",
1141  25,0.,250.,40,1.,81.);
1142  m_npoint_lengthT->SetXTitle("trk length");
1143  m_npoint_lengthT->SetYTitle("Npoints on trk");
1144 
1145  m_fpoint_lengthT = QAH::H2F("QaGtrkFitPntLengthT","globtrk: N fit pnts vs length, tpc",
1146  25,0.,250.,40,1.,81.);
1147  m_fpoint_lengthT->SetXTitle("trk length");
1148  m_fpoint_lengthT->SetYTitle("Npoints on trk");
1149 
1150  m_chisq0_momT = QAH::H2F("QaGtrkChi0MomT","globtrk: Chisq0 vs log10 mom, tpc",40,1.,4.,50,0.,5.);
1151  m_chisq0_momT->SetXTitle("log10 P (MeV)");
1152  m_chisq0_momT->SetYTitle("chisq0") ;
1153 
1154  m_chisq0_etaT = QAH::H2F("QaGtrkChi0EtaT","globtrk: Chisq0 vs eta, tpc",20,-2.,2.,20,0.,5.);
1155  m_chisq0_etaT->SetXTitle("eta");
1156  m_chisq0_etaT->SetYTitle("chisq0");
1157 
1158  m_chisq0_dipT = QAH::H2F("QaGtrkChi0TanlT","globtrk: Chisq0 vs tanl(dip), tpc",20,-3.2,3.2,20,0.,5.);
1159  m_chisq0_dipT->SetXTitle("dip angle");
1160  m_chisq0_dipT->SetYTitle("chisq0");
1161 
1162  m_chisq0_zfT = QAH::H2F("QaGtrkChi0zfT","globtrk: Chisq0 vs zfirst, tpc",20,-250.,250.,20,0.,5.);
1163  m_chisq0_zfT->SetXTitle("zfirst");
1164  m_chisq0_zfT->SetYTitle("chisq0");
1165 
1166  m_chisq0_phiT = QAH::H2F("QaGtrkChi0PhiT","globtrk: Chisq0 vs phi, tpc",20,0.,360,20,0.,5.);
1167  m_chisq0_phiT->SetXTitle("phi");
1168  m_chisq0_phiT->SetYTitle("chisq0");
1169 
1170  m_psi_phiT = QAH::H2F("QaGtrkPsiPhiT","globtrk: psi vs phi, tpc",36, 0.,360.,36,0.,360.);
1171  m_psi_phiT->SetXTitle("phi");
1172  m_psi_phiT->SetYTitle("psi");
1173 
1174  if (silHists) {
1175 
1176 // 1D tpc + silicon (svt+ssd)
1177 
1178  m_pointTS = QAH::H1F("QaGtrkNPntTS", "globtrk: N points on trk,tpc+svt", 60, 0.,60.);
1179  m_max_pointTS = QAH::H1F("QaGtrkNPntMaxTS","globtrk: N max pnts on trk, tpc+svt", 50, 0.,100.);
1180  m_fit_pointTS = QAH::H1F("QaGtrkNPntFitTS","globtrk: N fit pnts on trk, tpc+svt", 60, 0.,60.);
1181  m_glb_ratiomTS = QAH::H1F("QaGtrkRnmTS", "globtrk: ratio Nfit/max pnt, tpc+svt", 55, 0., 1.1);
1182  m_glb_chargeTS = QAH::H1F("QaGtrkChrgTS", "globtrk: charge, tpc+svt ", 20,-2.,2.);
1183  m_glb_r0TS = QAH::H1F("QaGtrkR0TS", "globtrk: radius at start (cm), tpc+svt", 100,0.,25.);
1184  m_glb_r0TS->SetMinimum(100);
1185  m_glb_phi0TS = QAH::H1F("QaGtrkPhi0TS", "globtrk: azimuth (phi) at start (deg,force 0-360),tpc+svt", 64, 0.,360.);
1186  m_glb_z0TS = QAH::H1F("QaGtrkZ0TS", "globtrk: z-coord at start (cm), tpc+svt", 50, -50.,50.);
1187  m_glb_z0TS->SetMinimum(500);
1188  m_glb_curvTS = QAH::H1F("QaGtrkCurvTS", "globtrk: log10 curvature (1/cm), tpc+svt", 80,-3.5,0.5);
1189  m_glb_f0TS = QAH::MH1F("QaGtrkf0TS", "globtrk: first point: hit-helix, tpc+svt",50,-5.,5.,3);
1190  m_glb_f0TS->Rebin(0,"X");
1191  m_glb_f0TS->Rebin(1,"Y");
1192  m_glb_f0TS->Rebin(2,"Z");
1193  m_glb_f0TS->SetStats(kFALSE);
1194  m_glb_xfTS = QAH::H1F("QaGtrkXfTS", "globtrk: x of first hit on trk, tpc+svt", 50,-200.,200.);
1195  m_glb_xf0TS = QAH::H1F("QaGtrkXf0TS", "globtrk: x of first hit - on helix at start, tpc+svt",50,-5.,5.);
1196  m_glb_yfTS = QAH::H1F("QaGtrkYfTS", "globtrk: y of first hit on trk, tpc+svt", 50,-200.,200.);
1197  m_glb_yf0TS = QAH::H1F("QaGtrkYf0TS", "globtrk: y of first hit - on helix at start, tpc+svt",50,-5.,5.);
1198  m_glb_zfTS = QAH::H1F("QaGtrkZfTS", "globtrk: z of first hit on trk, svt", 50,-30.,30.);
1199  m_glb_zf0TS = QAH::H1F("QaGtrkZf0TS", "globtrk: z of first hit - on helix at start, tpc+svt",50,-5.,5.);
1200  m_glb_radfTS = QAH::H1F("QaGtrkRTS", "globtrk: radial position of first hit, tpc+svt",80,4.,200.);
1201  m_glb_rzf0TS = QAH::MH1F("QaGtrkRZf0TS", "globtrk: first point: hit - helix, tpc+svt",60,-3.,3.,2);
1202  m_glb_rzf0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1203  m_glb_rzf0TS->Rebin(1,"z_{dif}");
1204  m_glb_rzf0TS->SetStats(kFALSE);
1205  m_glb_rzl0TS = QAH::MH1F("QaGtrkRZl0TS", "globtrk: last point: hit - helix, tpc+svt",60,-3.,3.,2);
1206  m_glb_rzl0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1207  m_glb_rzl0TS->Rebin(1,"z_{dif}");
1208  m_glb_rzl0TS->SetStats(kFALSE);
1209  m_glb_phifTS = QAH::H1F("QaGtrkPhifTS", "globtrk: phi of first point on track, svt",64,0,360);
1210  m_glb_ssd_phi = QAH::H1F("QaGtrkPhifSSD", "globtrk: phi of ssd point (per event)",64,0,360);
1211  m_lengthTS = QAH::H1F("QaGtrkLengthTS", "globtrk: track length, tpc+svt", 50,0.,300.);
1212  m_psiTS = QAH::H1F("QaGtrkPsiTS", "globtrk: psi, tpc+svt (deg) ", 64, 0.,360.);
1213  m_tanlTS = QAH::H1F("QaGtrkTanlTS", "globtrk: tanl, tpc+svt",32,-4.,4.);
1214  m_glb_thetaTS = QAH::H1F("QaGtrkThetaTS", "globtrk: theta (degrees), tpc+svt",36,0.,180.);
1215  m_etaTS = QAH::H1F("QaGtrkEtaTS", "globtrk: eta, tpc+svt",40,-2.,2.);
1216  m_pTTS = QAH::H1F("QaGtrkPtTS", "globtrk: pT, tpc+svt",50,0.,10.);
1217  m_momTS = QAH::H1F("QaGtrkPTS", "globtrk: momentum, tpc+svt",50,0.,10.);
1218  m_chisq0TS = QAH::H1F("QaGtrkChisq0TS", "globtrk: chisq0, tpc+svt", 50, 0.,5.);
1219  m_glb_impactTS = QAH::MH1F("QaGtrkImpactTS", "globtrk: log10 impact param from prim vtx, tpc+svt",
1220  120,-3.0,3.0,3);
1221  m_glb_impactTS->Rebin(0,"East");
1222  m_glb_impactTS->Rebin(1,"West");
1223  m_glb_impactTS->Rebin(2,"All");
1224  m_glb_impactTS->SetStats(kFALSE);
1225  m_glb_simpactTS = QAH::MH1F("QaGtrkSImpactTS", "globtrk: signed impact param from prim vtx, tpc+svt",
1226  50,-0.8,0.8,3);
1227  m_glb_simpactTS->Rebin(0,"East");
1228  m_glb_simpactTS->Rebin(1,"West");
1229  m_glb_simpactTS->Rebin(2,"All");
1230  m_glb_simpactTS->SetStats(kFALSE);
1231  m_glb_impactrTS = QAH::H1F("QaGtrkImpactrTS", "globtrk: impact param from prim vtx, tpc+svt",
1232  100,0.,30.);
1233 
1234 
1235 // 2D - tpc + silicon (svt + ssd)
1236 
1237  m_pT_eta_recTS = QAH::H2F("QaGtrkPtVsEtaTS","globtrk: log10 pT vs eta, tpc+svt", 20,-2.,2.,40,1.,4.);
1238  m_pT_eta_recTS->SetXTitle("eta");
1239  m_pT_eta_recTS->SetYTitle(" log10 pT (MeV)");
1240 
1241  m_globtrk_xf_yfTS = QAH::H2F("QaGtrkXfYfTS", "globtrk: Y vs X of first hit on trk, tpc+svt", 40,-200.,200.,40,-200.,200.);
1242  m_globtrk_xf_yfTS->SetXTitle("x first");
1243  m_globtrk_xf_yfTS->SetYTitle("y first");
1244 
1245  m_tanl_zfTS = QAH::H2F("QaGtrkTanlzfTS","globtrk: tanl(dip) versus (zfirst-zvtx)/arc length, svt",60,-3.,3.,60,-3.,3.);
1246  m_tanl_zfTS->SetXTitle("(zfirst-zvtx)/arc length");
1247  m_tanl_zfTS->SetYTitle("tanl");
1248 
1249  m_mom_trklengthTS = QAH::H2F("QaGtrkPVsTrkLTS","globtrk: log10 mom vs trk length, tpc+svt",50,0.,250.,40,1.,4.);
1250  m_mom_trklengthTS->SetXTitle("trk length");
1251  m_mom_trklengthTS->SetYTitle("log10 P (MeV)");
1252 
1253  m_eta_trklengthTS = QAH::H2F("QaGtrkLVEtaTS","globtrk: trk length vs eta, tpc+svt",20,-2.,2.,50,0.,250.);
1254  m_eta_trklengthTS->SetXTitle("eta");
1255  m_eta_trklengthTS->SetYTitle("length");
1256 
1257  m_npoint_lengthTS = QAH::H2F("QaGtrkNPntLTS","globtrk: N pnts vs length, tpc+svt",25,0.,250.,25,0.,50.);
1258  m_npoint_lengthTS->SetXTitle("trk length");
1259  m_npoint_lengthTS->SetYTitle("Npoints on trk");
1260 
1261  m_fpoint_lengthTS = QAH::H2F("QaGtrkFitPntLTS","globtrk: N fit pnts vs length, tpc+svt",25,0.,250.,25,0.,50.);
1262  m_fpoint_lengthTS->SetXTitle("trk length");
1263  m_fpoint_lengthTS->SetYTitle("Npoints on trk");
1264 
1265  m_chisq0_momTS = QAH::H2F("QaGtrkChi0MomTS","globtrk: Chisq0 vs log10 mom, tpc+svt",40,1.,4.,50,0.,5.);
1266  m_chisq0_momTS->SetXTitle("log10 P (MeV)");
1267  m_chisq0_momTS->SetYTitle("chisq0") ;
1268 
1269  m_chisq0_etaTS = QAH::H2F("QaGtrkChi0EtaTS","globtrk: Chisq0 vs eta, tpc+svt",20,-2.,2.,20,0.,5.);
1270  m_chisq0_etaTS->SetXTitle("eta");
1271  m_chisq0_etaTS->SetYTitle("chisq0");
1272 
1273  m_chisq0_dipTS = QAH::H2F("QaGtrkChi0TanlTS","globtrk: Chisq0 vs tanl(dip), tpc+svt",20,-3.2,3.2,20,0.,5.);
1274  m_chisq0_dipTS->SetXTitle("dip angle");
1275  m_chisq0_dipTS->SetYTitle("chisq0");
1276 
1277  m_chisq0_zfTS = QAH::H2F("QaGtrkChi0zfTS","globtrk: Chisq0 vs zfirst, tpc+svt",20,-250.,250.,20,0.,5.);
1278  m_chisq0_zfTS->SetXTitle("zfirst");
1279  m_chisq0_zfTS->SetYTitle("chisq0");
1280 
1281  m_chisq0_phiTS = QAH::H2F("QaGtrkChi0PhiTS","globtrk: Chisq0 vs phi, tpc+svt",20,0.,360,20,0.,5.);
1282  m_chisq0_phiTS->SetXTitle("phi");
1283  m_chisq0_phiTS->SetYTitle("chisq0");
1284 
1285  m_psi_phiTS = QAH::H2F("QaGtrkPsiPhiTS","globtrk: psi vs phi, tpc+svt",36, 0.,360.,36,0.,360.);
1286  m_psi_phiTS->SetXTitle("phi");
1287  m_psi_phiTS->SetYTitle("psi");
1288 
1289  }
1290 
1291  if (ftpHists) { // FTPC hists
1292 
1293 // 1D ftpc
1294 
1295  // both east (solid) and west(dashed) on same plot
1296  m_pointF = QAH::MH1F("QaGtrkNPntF", "globtrk: N points on trk,ftpc", 8,4.,12.,2);
1297  m_pointF->Rebin(0,"East");
1298  m_pointF->Rebin(1,"West");
1299  m_pointF->SetStats(kFALSE);
1300  m_max_pointF = QAH::MH1F("QaGtrkNPntMaxF", "globtrk: N max pnts on trk, ftpc", 8,4.,12.,2);
1301  m_max_pointF->Rebin(0,"East");
1302  m_max_pointF->Rebin(1,"West");
1303  m_max_pointF->SetStats(kFALSE);
1304  m_glb_ratiomF = QAH::MH1F("QaGtrkRnmF", "globtrk: ratio Nfit/max pnt, ftpc", 11,0.,1.1,2);
1305  m_glb_ratiomF->Rebin(0,"East");
1306  m_glb_ratiomF->Rebin(1,"West");
1307  m_glb_ratiomF->SetStats(kFALSE);
1308  m_glb_chargeF = QAH::MH1F("QaGtrkChrgF", "globtrk: charge, ftpc", 20,-2.,2.,2);
1309  m_glb_chargeF->Rebin(0,"East");
1310  m_glb_chargeF->Rebin(1,"West");
1311  m_glb_chargeF->SetStats(kFALSE);
1312  m_glb_xfF = QAH::MH1F("QaGtrkXfF", "globtrk: x of first hit on trk, ftpc",70,-35.,35.,2);
1313  m_glb_xfF->Rebin(0,"East");
1314  m_glb_xfF->Rebin(1,"West");
1315  m_glb_xfF->SetStats(kFALSE);
1316  m_glb_yfF = QAH::MH1F("QaGtrkYfF", "globtrk: y of first hit on trk, ftpc",70,-35.,35.,2);
1317  m_glb_yfF->Rebin(0,"East");
1318  m_glb_yfF->Rebin(1,"West");
1319  m_glb_yfF->SetStats(kFALSE);
1320  m_glb_zfF = QAH::MH1F("QaGtrkZfF", "globtrk: z of first hit on trk, ftpc",75,-225.,225.,2);
1321  m_glb_zfF->Rebin(0,"East");
1322  m_glb_zfF->Rebin(1,"West");
1323  m_glb_zfF->SetStats(kFALSE);
1324  m_glb_planefF = QAH::MH1F("QaGtrkPlanefF", "globtrk: plane of first hit on trk, ftpc",20,0.5,20.5,2);
1325  m_glb_planefF->Rebin(0,"East");
1326  m_glb_planefF->Rebin(1,"West");
1327  m_glb_planefF->SetStats(kFALSE);
1328  m_glb_radfF = QAH::MH1F("QaGtrkRF", "globtrk: radial position of first hit, ftpc",35,0.,35.,2);
1329  m_glb_radfF->Rebin(0,"East");
1330  m_glb_radfF->Rebin(1,"West");
1331  m_glb_radfF->SetStats(kFALSE);
1332  m_lengthF = QAH::MH1F("QaGtrkLengthF", "globtrk: track length, ftpc",60,0.,120.,2);
1333  m_lengthF->Rebin(0,"East");
1334  m_lengthF->Rebin(1,"West");
1335  m_lengthF->SetStats(kFALSE);
1336  m_psiF = QAH::MH1F("QaGtrkPsiF", "globtrk: psi, ftpc",90,0.,360.,2);
1337  m_psiF->Rebin(0,"East");
1338  m_psiF->Rebin(1,"West");
1339  m_psiF->SetStats(kFALSE);
1340  m_etaF = QAH::MH1F("QaGtrkEtaF", "globtrk: |eta|, ftpc",80,2,5,2);
1341  m_etaF->Rebin(0,"East");
1342  m_etaF->Rebin(1,"West");
1343  m_etaF->SetStats(kFALSE);
1344  m_pTF = QAH::MH1F("QaGtrkPtF", "globtrk: pT, ftpc",50,0.,10.,2);
1345  m_pTF->Rebin(0,"East");
1346  m_pTF->Rebin(1,"West");
1347  m_pTF->SetStats(kFALSE);
1348  m_momF = QAH::MH1F("QaGtrkPF", "globtrk: momentum, ftpc",50,0.,10.,2);
1349  m_momF->Rebin(0,"East");
1350  m_momF->Rebin(1,"West");
1351  m_momF->SetStats(kFALSE);
1352  m_chisq0F = QAH::MH1F("QaGtrkChisq0F", "globtrk: chi2/ndf x,y fit, ftpc",50,0.,5.,2);
1353  m_chisq0F->Rebin(0,"East");
1354  m_chisq0F->Rebin(1,"West");
1355  m_chisq0F->SetStats(kFALSE);
1356  m_chisq1F = QAH::MH1F("QaGtrkChisq1F", "globtrk: chi2/ndf r,z fit, ftpc",50,0.,5.,2);
1357  m_chisq1F->Rebin(0,"East");
1358  m_chisq1F->Rebin(1,"West");
1359  m_chisq1F->SetStats(kFALSE);
1360  m_glb_impactF = QAH::MH1F("QaGtrkImpactF", "globtrk: log10 impact param from prim vtx, ftpc",120,-3.0,3.,2);
1361  m_glb_impactF->Rebin(0,"East");
1362  m_glb_impactF->Rebin(1,"West");
1363  m_glb_impactF->SetStats(kFALSE);
1364  m_glb_impactrF = QAH::MH1F("QaGtrkImpactrF", "globtrk: impact param from prim vtx, ftpc",100,0.,10.,2);
1365  m_glb_impactrF->Rebin(0,"East");
1366  m_glb_impactrF->Rebin(1,"West");
1367  m_glb_impactrF->SetStats(kFALSE);
1368  // separate east and west plots
1369  m_pointFE = QAH::H1F("QaGtrkNPntFE", "globtrk: N points on trk,ftpc east", 8, 4.,12.);
1370  m_pointFW = QAH::H1F("QaGtrkNPntFW", "globtrk: N points on trk,ftpc west", 8, 4.,12.);
1371  m_max_pointFE = QAH::H1F("QaGtrkNPntMaxFE", "globtrk: N max pnts on trk, ftpc east", 8, 4.,12.);
1372  m_max_pointFW = QAH::H1F("QaGtrkNPntMaxFW", "globtrk: N max pnts on trk, ftpc west", 8, 4.,12.);
1373  m_glb_ratiomFE = QAH::H1F("QaGtrkRnmFE", "globtrk: ratio Nfit/max pnt, ftpc east", 55, 0., 1.1);
1374  m_glb_ratiomFW = QAH::H1F("QaGtrkRnmFW", "globtrk: ratio Nfit/max pnt, ftpc west", 55, 0., 1.1);
1375  m_glb_chargeFE = QAH::H1F("QaGtrkChrgFE", "globtrk: charge, ftpc east ", 20,-2.,2.);
1376  m_glb_chargeFW = QAH::H1F("QaGtrkChrgFW", "globtrk: charge, ftpc west ", 20,-2.,2.);
1377  m_glb_xfFE = QAH::H1F("QaGtrkXfFE", "globtrk: x of first hit on trk, ftpc east", 70,-35.,35.);
1378  m_glb_xfFW = QAH::H1F("QaGtrkXfFW", "globtrk: x of first hit on trk, ftpc west", 70,-35.,35.);
1379  m_glb_yfFE = QAH::H1F("QaGtrkYfFE", "globtrk: y of first hit on trk, ftpc east", 70,-35.,35.);
1380  m_glb_yfFW = QAH::H1F("QaGtrkYfFW", "globtrk: y of first hit on trk, ftpc west", 70,-35.,35.);
1381  m_glb_zfFE = QAH::H1F("QaGtrkZfFE", "globtrk: z of first hit on trk, ftpc east", 75,-225.,-150.);
1382  m_glb_zfFW = QAH::H1F("QaGtrkZfFW", "globtrk: z of first hit on trk, ftpc west", 75,150.,225.);
1383  m_glb_radfFE = QAH::H1F("QaGtrkRFE", "globtrk: radial position of first hit, ftpc east",35,0.,35.);
1384  m_glb_radfFW = QAH::H1F("QaGtrkRFW", "globtrk: radial position of first hit, ftpc west",35,0.,35.);
1385  m_lengthFE = QAH::H1F("QaGtrkLengthFE", "globtrk: track length, ftpc east", 60,0.,120.);
1386  m_lengthFW = QAH::H1F("QaGtrkLengthFW", "globtrk: track length, ftpc west", 60,0.,120.);
1387  m_psiFE = QAH::H1F("QaGtrkPsiFE", "globtrk: psi, ftpc east", 90, 0.,360.);
1388  m_psiFW = QAH::H1F("QaGtrkPsiFW", "globtrk: psi, ftpc west", 90, 0.,360.);
1389  m_etaFE = QAH::H1F("QaGtrkEtaFE", "globtrk: eta, ftpc east",80,-4.5,-2.);
1390  m_etaFW = QAH::H1F("QaGtrkEtaFW", "globtrk: eta, ftpc west",80,2.,4.5);
1391  m_pTFE = QAH::H1F("QaGtrkPtFE", "globtrk: pT, ftpc east",50,0.,10.);
1392  m_pTFW = QAH::H1F("QaGtrkPtFW", "globtrk: pT, ftpc west",50,0.,10.);
1393  m_momFE = QAH::H1F("QaGtrkPFE", "globtrk: momentum, ftpc east ",50,0.,5.);
1394  m_momFW = QAH::H1F("QaGtrkPFW", "globtrk: momentum, ftpc west ",50,0.,5.);
1395  m_chisq0FE = QAH::H1F("QaGtrkChisq0FE", "globtrk: chi2/ndf x,y fit, ftpc east", 50, 0.,5.);
1396  m_chisq0FW = QAH::H1F("QaGtrkChisq0FW", "globtrk: chi2/ndf x,y fit, ftpc west", 50, 0.,5.);
1397  m_chisq1FE = QAH::H1F("QaGtrkChisq1FE", "globtrk: chi2/ndf r,z fit, ftpc east", 50, 0.,5.);
1398  m_chisq1FW = QAH::H1F("QaGtrkChisq1FW", "globtrk: chi2/ndf r,z fit, ftpc west", 50, 0.,5.);
1399 
1400 // 2D - ftpc
1401 
1402  m_pT_eta_recFE = QAH::H2F("QaGtrkPtVsEtaFE","globtrk: log10 pT vs eta, ftpcE",20,-4.5,-2.,40,1.,4.);
1403  m_pT_eta_recFE->SetXTitle("eta");
1404  m_pT_eta_recFE->SetYTitle("log10 pT");
1405  m_pT_eta_recFW = QAH::H2F("QaGtrkPtVsEtaFW","globtrk: log10 pT vs eta, ftpcW",20,2.,4.5,40,1.,4.);
1406  m_pT_eta_recFW->SetXTitle("eta");
1407  m_pT_eta_recFW->SetYTitle("log10 pT");
1408  m_globtrk_xf_yfFE = QAH::H2F("QaGtrkXfYfFE","globtrk: Y vs X of first hit on trk, ftpcE", 70,-35.,35.,70,-35.,35.);
1409  m_globtrk_xf_yfFE->SetXTitle("x first");
1410  m_globtrk_xf_yfFE->SetYTitle("y first");
1411  m_globtrk_xf_yfFW = QAH::H2F("QaGtrkXfYfFW","globtrk: Y vs X of first hit on trk, ftpcW", 70,-35.,35.,70,-35.,35.);
1412  m_globtrk_xf_yfFW->SetXTitle("x first");
1413  m_globtrk_xf_yfFW->SetYTitle("y first");
1414  m_globtrk_padtimeFE = QAH::H2F("QaGtrkPadTimeFtpcE","globtrk: #pads vs #timebins, ftpcE",12,0.5,12.5,10,0.5,10.5);
1415 // m_globtrk_padtimeFE = QAH::H2F("QaGtrkPadTimeFtpcE","globtrk: #pads vs #timebins of hits, ftpcE",12,0.5,12.5,10,0.5,10.5);
1416  m_globtrk_padtimeFE->SetXTitle("#timebins");
1417  m_globtrk_padtimeFE->SetYTitle("#pads");
1418  m_globtrk_padtimeFW = QAH::H2F("QaGtrkPadTimeFtpcW","globtrk: #pads vs #timebins, ftpcW",12,0.5,12.5,10,0.5,10.5);
1419 // m_globtrk_padtimeFW = QAH::H2F("QaGtrkPadTimeFtpcW","globtrk: #pads vs #timebins of hits, ftpcW",12,0.5,12.5,10,0.5,10.5);
1420  m_globtrk_padtimeFW->SetXTitle("#timebins");
1421  m_globtrk_padtimeFW->SetYTitle("#pads");
1422  m_eta_trklengthFE = QAH::H2F("QaGtrkLengthVEtaFE","globtrk: trk length vs eta, ftpcE",25,-4.5,-2.,30,0.,120.);
1423  m_eta_trklengthFE->SetXTitle("eta");
1424  m_eta_trklengthFE->SetYTitle("length");
1425  m_eta_trklengthFW = QAH::H2F("QaGtrkLengthVEtaFW","globtrk: trk length vs eta, ftpcW",25,2.,4.5,30,0.,120.);
1426  m_eta_trklengthFW->SetXTitle("eta");
1427  m_eta_trklengthFW->SetYTitle("length");
1428  m_npoint_lengthFE = QAH::H2F("QaGtrkNPntLengthFE","globtrk: N pnts vs length, ftpcE",30,0.,120.,15,0.,15.);
1429  m_npoint_lengthFE->SetXTitle("trk length");
1430  m_npoint_lengthFE->SetYTitle("Npoints on trk");
1431  m_npoint_lengthFW = QAH::H2F("QaGtrkNPntLengthFW","globtrk: N pnts vs length, ftpcW",30,0.,120.,15,0.,15.);
1432  m_npoint_lengthFW->SetXTitle("trk length");
1433  m_npoint_lengthFW->SetYTitle("Npoints on trk");
1434 
1435  } // ftpHists
1436 
1437  m_global_pxl_hit = QAH::H1F("QaGtrkPxlHit","PIXEL: Hits per global track",5,-0.5,4.5);
1438  m_global_pxl_hit->SetXTitle("hits in PIXEL");
1439  m_global_ist_hit = QAH::H1F("QaGtrkIstHit","IST: Hits per global track",5,-0.5,4.5);
1440  m_global_ist_hit->SetXTitle("hits in IST");
1441  m_global_hft_hit = QAH::H1F("QaGtrkHftHit","HFT: Hits per global track",5,-0.5,4.5);
1442  m_global_hft_hit->SetXTitle("hits in HFT subsystem");
1443 
1444  m_global_hft_hit->GetXaxis()->SetBinLabel(1,"");
1445  m_global_hft_hit->GetXaxis()->SetBinLabel(2,"PIXEL (inner)");
1446  m_global_hft_hit->GetXaxis()->SetBinLabel(3,"PIXEL (outer)");
1447  m_global_hft_hit->GetXaxis()->SetBinLabel(4,"IST");
1448  m_global_hft_hit->GetXaxis()->SetBinLabel(5,"SSD");
1449 
1450 }
1451 //____________________________________________________
1452 void StQABookHist::BookHistPrim(){
1453 
1454 // for method MakeHistPrim - from table primtrk
1455 
1456 // 1D
1457  m_primtrk_tot = QAH::H1F("QaPtrkTot", "primtrk: tot num tracks",50,0.,5000.);
1458  m_primtrk_tot_sm = QAH::H1F("QaPtrkTotsm", "primtrk: tot num tracks",50,0.,200.);
1459  m_primtrk_iflag = QAH::H1F("QaPtrkFlag", "primtrk: iflag - all",160,-799.,900.);
1460  m_primtrk_good = QAH::H1F("QaPtrkGood", "primtrk: tot num tracks iflag>0",50,0.,5000.);
1461  m_primtrk_good_sm = QAH::H1F("QaPtrkGoodsm","primtrk: tot num tracks iflag>0",50,0.,200.);
1462  if (silHists)
1463  m_primtrk_goodTTS = QAH::H1F("QaPtrkGoodTTS","primtrk: tot num tracks iflag>0, tpc,svt",150,0.,4500.);
1464  m_primtrk_goodF = QAH::H2F("QaPtrkGoodF", "primtrk: tot num tracks iflag>0, ftpc",150,0.,1500.,150,0.,1500.);
1465  m_primtrk_goodF->SetXTitle("East");
1466  m_primtrk_goodF->SetYTitle("West");
1467  m_primglob_good = QAH::H1F("QaPtrkGlob","primtrk: ratio primary/global tracks w/ iflag>0",50,0,1);
1468  m_primglob_fit = QAH::H1F("QaPtrkGlobFit","primtrk: ratio primary/global nfit points",50,0,2);
1469  m_pdet_id = QAH::H1F("QaPtrkDetId", "primtrk: Detector ID good tracks - all",48,-0.5,47.5);
1470  if (silHists) {
1471  m_primtrk_meanptTTS = QAH::MH1F("QaPtrkMeanPtTTS","primtrk: <pT>, tpc, tpc+svt",50,0.,2.,2);
1472  m_primtrk_meanptTTS->Rebin(0,"TPC+SVT");
1473  m_primtrk_meanptTTS->Rebin(1,"TPC");
1474  m_primtrk_meanptTTS->SetStats(kFALSE);
1475  }
1476  m_primtrk_meanptF = QAH::MH1F("QaPtrkMeanPtF","primtrk: <pT>, ftpc",50,0.,2.,2);
1477  m_primtrk_meanptF->Rebin(0,"East");
1478  m_primtrk_meanptF->Rebin(1,"West");
1479  m_primtrk_meanptF->SetStats(kFALSE);
1480  if (silHists) {
1481  m_primtrk_meanetaTTS = QAH::MH1F("QaPtrkMeanEtaTTS","primtrk: <eta>, tpc, tpc+svt",40,-2.,2.,2);
1482  m_primtrk_meanetaTTS->Rebin(0,"TPC+SVT");
1483  m_primtrk_meanetaTTS->Rebin(1,"TPC");
1484  m_primtrk_meanetaTTS->SetStats(kFALSE);
1485  }
1486  m_primtrk_meanetaF = QAH::MH1F("QaPtrkMeanEtaF","primtrk: |<eta>|, ftpc",40,2,5,2);
1487  m_primtrk_meanetaF->Rebin(0,"East");
1488  m_primtrk_meanetaF->Rebin(1,"West");
1489  m_primtrk_meanetaF->SetStats(kFALSE);
1490  if (silHists) {
1491  m_ppsiTTS = QAH::MH1F("QaPtrkPsiTTS","primtrk: psi (deg), tpc, svt", 36, 0.,360.,2);
1492  m_ppsiTTS->Rebin(0,"TPC+SVT");
1493  m_ppsiTTS->Rebin(1,"TPC");
1494  m_ppsiTTS->SetStats(kFALSE);
1495  m_petaTTS = QAH::MH1F("QaPtrkEtaTTS","primtrk: eta, tpc,svt",40,-2.,2.,2);
1496  m_petaTTS->Rebin(0,"TPC+SVT");
1497  m_petaTTS->Rebin(1,"TPC");
1498  m_petaTTS->SetStats(kFALSE);
1499  m_ppTTTS = QAH::MH1F("QaPtrkPtTTS", "primtrk: pT, tpc,svt",50,0.,5.,2);
1500  m_ppTTTS->Rebin(0,"TPC+SVT");
1501  m_ppTTTS->Rebin(1,"TPC");
1502  m_ppTTTS->SetStats(kFALSE);
1503  m_pchisq0TTS = QAH::MH1F("QaPtrkChisq0TTS", "primtrk: chisq0, tpc,svt", 50, 0.,5.,2);
1504  m_pchisq0TTS->Rebin(0,"TPC+SVT");
1505  m_pchisq0TTS->Rebin(1,"TPC");
1506  m_pchisq0TTS->SetStats(kFALSE);
1507  m_pfpoint_lengthTTS = QAH::H2F("QaPtrkFitPntLTTS","primtrk: N fit pnts vs length, tpc,tpc+svt",25,70.,350.,25,0.,50.);
1508  m_pfpoint_lengthTTS->SetXTitle("trk length");
1509  m_pfpoint_lengthTTS->SetYTitle("Npoints on trk");
1510  }
1511 
1512 // 1D tpc
1513  m_ppointT = QAH::H1F("QaPtrkNPntT", "primtrk: N points on trk,tpc", 81, -0.5,80.5);
1514  m_pmax_pointT = QAH::H1F("QaPtrkNPntMaxT","primtrk: N max pnts on trk, tpc", 76, -1,151);
1515  m_pfit_pointT = QAH::H1F("QaPtrkNPntFitT","primtrk: N fit pnts on trk, tpc", 81, -0.5,80.5);
1516  m_prim_ratiomT = QAH::H1F("QaPtrkRnmT", "primtrk: ratio Nfit/max pnt, tpc", 55, 0., 1.1);
1517  m_prim_chargeT = QAH::H1F("QaPtrkChrgT", "primtrk: charge, tpc ", 20,-2.,2.);
1518  m_prim_r0T = QAH::H1F("QaPtrkR0T", "primtrk: radius at start (cm), tpc ", 50,0.,.1);
1519  m_prim_phi0T = QAH::H1F("QaPtrkPhi0T", "primtrk: azimuth (phi) at start (deg,force 0,360), tpc ", 64,0.,360.);
1520  m_prim_z0T = QAH::H1F("QaPtrkZ0T", "primtrk: z-coord at start (cm), tpc ", 50, -300.,300.);
1521  m_prim_curvT = QAH::H1F("QaPtrkCurvT", "primtrk: log10 curvature (1/cm), tpc ", 80,-3.5,0.5);
1522  m_prim_f0 = QAH::MH1F("QaPtrkf0", "primtrk: first point: hit-helix, tpc",60,-3.,3.,3);
1523  m_prim_f0->Rebin(0,"X");
1524  m_prim_f0->Rebin(1,"Y");
1525  m_prim_f0->Rebin(2,"Z");
1526  m_prim_xfT = QAH::H1F("QaPtrkXfT", "primtrk: x of first hit on trk, tpc", 50,-200.,200.);
1527  m_prim_xf0 = QAH::H1F("QaPtrkXf0", "primtrk: first point: x_hit - x_helix, tpc",60,-3.,3.);
1528  m_prim_yfT = QAH::H1F("QaPtrkYfT", "primtrk: y of first hit on trk, tpc", 50,-200.,200.);
1529  m_prim_yf0 = QAH::H1F("QaPtrkYf0", "primtrk: first point: y_hit - y_helix, tpc",60,-3.,3.);
1530  m_prim_zfT = QAH::H1F("QaPtrkZfT", "primtrk: z of first hit on trk, tpc", 50,-200.,200.);
1531  m_prim_zf0 = QAH::H1F("QaPtrkZf0", "primtrk: first point: z_hit - z_helix, tpc",60,-3.,3.);
1532  m_prim_radfT = QAH::H1F("QaPtrkRT", "primtrk: radial position of first hit, tpc", 50,0.,200.);
1533  m_prim_rzf0 = QAH::MH1F("QaPtrkRZf0", "primtrk: first point: hit - helix (r,z), tpc",60,-3.,3.,2);
1534  m_prim_rzf0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1535  m_prim_rzf0->Rebin(1,"z_{dif}");
1536  m_prim_rzf0->SetStats(kFALSE);
1537  m_prim_rzl0 = QAH::MH1F("QaPtrkRZl0", "primtrk: last point: hit - helix (r,z), tpc",60,-3.,3.,2);
1538  m_prim_rzl0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1539  m_prim_rzl0->Rebin(1,"z_{dif}");
1540  m_prim_rzl0->SetStats(kFALSE);
1541  m_plengthT = QAH::H1F("QaPtrkLengthT", "primtrk: track length, tpc", 50,0.,300.);
1542  m_ppsiT = QAH::H1F("QaPtrkPsiT", "primtrk: psi, tpc (deg)", 36, 0.,360.);
1543  m_ptanlT = QAH::H1F("QaPtrkTanlT", "primtrk: tanl, tpc",32,-4.,4.);
1544  m_prim_thetaT = QAH::H1F("QaPtrkThetaT", "primtrk: theta (degrees), tpc",36,0.,180.);
1545  m_petaT = QAH::H1F("QaPtrkEtaT", "primtrk: eta, tpc",40,-2.,2.);
1546  m_ppTT = QAH::H1F("QaPtrkPtT", "primtrk: pT, tpc",50,0.,5.);
1547  m_pmomT = QAH::H1F("QaPtrkPT", "primtrk: momentum, tpc",50,0.,5.);
1548  m_pchisq0T = QAH::H1F("QaPtrkChisq0T", "primtrk: chisq0, tpc", 50, 0.,5.);
1549 
1550 // 2D - tpc
1551  m_ppT_eta_recT = QAH::H2F("QaPtrkPtVsEtaT","primtrk: log10 pT vs eta, tpc", 20,-2.,2.,40,1.,4.);
1552  m_ppT_eta_recT->SetXTitle("eta");
1553  m_ppT_eta_recT->SetYTitle(" log10 pT (MeV)");
1554  m_primtrk_xf_yfTE = QAH::H2F("QaPtrkXfYfTE", "primtrk: Y vs X of first hit on trk, tpc east", 40,-200.,200.,40,-200.,200.);
1555  m_primtrk_xf_yfTE->SetXTitle("x first");
1556  m_primtrk_xf_yfTE->SetYTitle("y first");
1557  m_primtrk_xf_yfTW = QAH::H2F("QaPtrkXfYfTW", "primtrk: Y vs X of first hit on trk, tpc west", 40,-200.,200.,40,-200.,200.);
1558  m_primtrk_xf_yfTW->SetXTitle("x first");
1559  m_primtrk_xf_yfTW->SetYTitle("y first");
1560  m_ptanl_zfT = QAH::H2F("QaPtrkTanlzf","primtrk: tanl(dip) versus zfirst-zvtx, tpc,tpc+svt",60,-3.,3.,60,-3.,3.);
1561  m_ptanl_zfT->SetXTitle("zfirst-zvtx");
1562  m_ptanl_zfT->SetYTitle("tanl");
1563  m_pmom_trklengthT = QAH::H2F("QaPtrkPVsTrkLength","primtrk: log10 mom vs trk length, tpc",50,70.,350.,40,1.,4.);
1564  m_pmom_trklengthT->SetXTitle("trk length");
1565  m_pmom_trklengthT->SetYTitle("log10 P (MeV)");
1566  m_peta_trklengthT = QAH::H2F("QaPtrkLengthVEtaT","primtrk: trk length vs eta, tpc",20,-2.,2.,50,70.,350.);
1567  m_peta_trklengthT->SetXTitle("eta");
1568  m_peta_trklengthT->SetYTitle("length");
1569  m_pnpoint_lengthT = QAH::H2F("QaPtrkNPntLengthT","primtrk: N pnts vs length, tpc",25,70.,350.,40,1.,81.);
1570  m_pnpoint_lengthT->SetXTitle("trk length");
1571  m_pnpoint_lengthT->SetYTitle("Npoints on trk");
1572  m_pfpoint_lengthT = QAH::H2F("QaPtrkFitPntLengthT","primtrk: N fit pnts vs length, tpc",25,70.,350.,40,1.,81.);
1573  m_pfpoint_lengthT->SetXTitle("trk length");
1574  m_pfpoint_lengthT->SetYTitle("Npoints on trk");
1575  m_pchisq0_momT = QAH::H2F("QaPtrkChi0MomT","primtrk: Chisq0 vs log10 mom, tpc",40,1.,4.,50,0.,5.);
1576  m_pchisq0_momT->SetXTitle("log10 P (MeV)");
1577  m_pchisq0_momT->SetYTitle("chisq0") ;
1578  m_pchisq0_etaT = QAH::H2F("QaPtrkChi0EtaT","primtrk: Chisq0 vs eta, tpc",20,-2.,2.,20,0.,5.);
1579  m_pchisq0_etaT->SetXTitle("eta");
1580  m_pchisq0_etaT->SetYTitle("chisq0");
1581  m_pchisq0_dipT = QAH::H2F("QaPtrkChi0TanlT","primtrk: Chisq0 vs tanl(dip), tpc",20,-3.2,3.2,20,0.,5.);
1582  m_pchisq0_dipT->SetXTitle("dip angle");
1583  m_pchisq0_dipT->SetYTitle("chisq0");
1584  m_pchisq0_zfT = QAH::H2F("QaPtrkChi0zfT","primtrk: Chisq0 vs zfirst, tpc",20,-200.,200.,20,0.,5.);
1585  m_pchisq0_zfT->SetXTitle("zfirst");
1586  m_pchisq0_zfT->SetYTitle("chisq0");
1587  m_ppsi_phiT = QAH::H2F("QaPtrkPsiPhiT","primtrk: psi vs phi, tpc",36, 0.,360.,36,0.,360.);
1588  m_ppsi_phiT->SetXTitle("phi");
1589  m_ppsi_phiT->SetYTitle("psi");
1590 
1591  if (silHists) {
1592 // 1D tpc + silicon (svt+ssd)
1593  m_ppointTS = QAH::H1F("QaPtrkNPntTS", "primtrk: N points on trk,tpc+svt", 60, 0.,60.);
1594  m_pmax_pointTS = QAH::H1F("QaPtrkNPntMaxTS","primtrk: N max pnts on trk, tpc+svt", 50, 0.,100.);
1595  m_pfit_pointTS = QAH::H1F("QaPtrkNPntFitTS","primtrk: N fit pnts on trk, tpc+svt", 60, 0.,60.);
1596  m_prim_ratiomTS = QAH::H1F("QaPtrkRnmTS", "primtrk: ratio Nfit/max pnt, tpc+svt", 55, 0., 1.2005);
1597  m_prim_chargeTS = QAH::H1F("QaPtrkChrgTS", "primtrk: charge, tpc+svt ", 20,-2.,2.);
1598  m_prim_r0TS = QAH::H1F("QaPtrkR0TS", "primtrk: radius at start (cm), tpc+svt", 50,0.,0.1);
1599  m_prim_phi0TS = QAH::H1F("QaPtrkPhi0TS", "primtrk: azimuth (phi) at start (deg,force 0-360),tpc+svt", 64, 0.,360.);
1600  m_prim_z0TS = QAH::H1F("QaPtrkZ0TS", "primtrk: z-coord at start (cm), tpc+svt", 50, -50.,50.);
1601  m_prim_curvTS = QAH::H1F("QaPtrkCurvTS", "primtrk: log10 curvature (1/cm), tpc+svt", 80,-3.5,0.5);
1602  m_prim_f0TS = QAH::MH1F("QaPtrkf0TS", "primtrk: first point: hit-helix, tpc+svt",50,-5.,5.,3);
1603  m_prim_f0TS->Rebin(0,"X");
1604  m_prim_f0TS->Rebin(1,"Y");
1605  m_prim_f0TS->Rebin(2,"Z");
1606  m_prim_f0TS->SetStats(kFALSE);
1607  m_prim_xfTS = QAH::H1F("QaPtrkXfTS", "primtrk: x of first hit on trk, tpc+svt", 50,-200.,200.);
1608  m_prim_xf0TS = QAH::H1F("QaPtrkXf0TS", "primtrk: first point: x_hit - x_helix, tpc+svt",50,-5.,5.);
1609  m_prim_yfTS = QAH::H1F("QaPtrkYfTS", "primtrk: y of first hit on trk, tpc+svt", 50,-200.,200.);
1610  m_prim_yf0TS = QAH::H1F("QaPtrkYf0TS", "primtrk: first point: y_hit - y_helix, tpc+svt",50,-5.,5.);
1611  m_prim_zfTS = QAH::H1F("QaPtrkZfTS", "primtrk: z of first hit on trk, tpc+svt", 50,-200.,200.);
1612  m_prim_zf0TS = QAH::H1F("QaPtrkZf0TS", "primtrk: first point: z_hit - z_helix, tpc+svt",50,-5.,5.);
1613  m_prim_radfTS = QAH::H1F("QaPtrkRTS", "primtrk: radial position of first hit, tpc+svt", 50,0.,200.);
1614  m_prim_rzf0TS = QAH::MH1F("QaPtrkRZf0TS", "primtrk: first point: hit - helix (r,z), tpc+svt",60,-3.,3.,2);
1615  m_prim_rzf0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1616  m_prim_rzf0TS->Rebin(1,"z_{dif}");
1617  m_prim_rzf0TS->SetStats(kFALSE);
1618  m_prim_rzl0TS = QAH::MH1F("QaPtrkRZl0TS", "primtrk: last point: hit - helix (r,z), tpc+svt",60,-3.,3.,2);
1619  m_prim_rzl0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1620  m_prim_rzl0TS->Rebin(1,"z_{dif}");
1621  m_prim_rzl0TS->SetStats(kFALSE);
1622  m_plengthTS = QAH::H1F("QaPtrkLengthTS", "primtrk: track length, tpc+svt", 50,0.,300.);
1623  m_ppsiTS = QAH::H1F("QaPtrkPsiTS", "primtrk: psi, tpc+svt (deg) ", 36, 0.,360.);
1624  m_ptanlTS = QAH::H1F("QaPtrkTanlTS", "primtrk: tanl, tpc+svt",32,-4.,4.);
1625  m_prim_thetaTS = QAH::H1F("QaPtrkThetaTS", "primtrk: theta (degrees), tpc+svt",36,0.,180.);
1626  m_petaTS = QAH::H1F("QaPtrkEtaTS", "primtrk: eta, tpc+svt",40,-2.,2.);
1627  m_ppTTS = QAH::H1F("QaPtrkPtTS", "primtrk: pT, tpc+svt",50,0.,5.);
1628  m_pmomTS = QAH::H1F("QaPtrkPTS", "primtrk: momentum, tpc+svt",50,0.,5.);
1629  m_pchisq0TS = QAH::H1F("QaPtrkChisq0TS", "primtrk: chisq0, tpc+svt", 50, 0.,5.);
1630  m_prim_ssd_phi = QAH::H1F("QaPtrkPhifSSD", "primtrk: phi of ssd point (per event)",64,0,360);
1631 
1632 // 2D - tpc + silicon (svt + ssd)
1633  m_ppT_eta_recTS = QAH::H2F("QaPtrkPtVsEtaTS","primtrk: log10 pT vs eta, tpc+svt", 20,-2.,2.,40,1.,4.);
1634  m_ppT_eta_recTS->SetXTitle("eta");
1635  m_ppT_eta_recTS->SetYTitle(" log10 pT (MeV)");
1636  m_primtrk_xf_yfTS = QAH::H2F("QaPtrkXfYfTS", "primtrk: Y vs X of first hit on trk, tpc+svt", 40,-200.,200.,40,-200.,200.);
1637  m_primtrk_xf_yfTS->SetXTitle("x first");
1638  m_primtrk_xf_yfTS->SetYTitle("y first");
1639  m_ptanl_zfTS = QAH::H2F("QaPtrkTanlzfTS","primtrk: tanl(dip) versus zfirst-zvtx, svt",60,-3.,3.,60,-3.,3.);
1640  m_ptanl_zfTS->SetXTitle("zfirst-zvtx");
1641  m_ptanl_zfTS->SetYTitle("tanl");
1642  m_pmom_trklengthTS = QAH::H2F("QaPtrkPVsTrkLTS","primtrk: log10 mom vs trk length, tpc+svt",50,70.,350.,40,1.,4.);
1643  m_pmom_trklengthTS->SetXTitle("trk length");
1644  m_pmom_trklengthTS->SetYTitle("log10 P (MeV)");
1645  m_peta_trklengthTS = QAH::H2F("QaPtrkLVEtaTS","primtrk: trk length vs eta, tpc+svt",20,-2.,2.,50,70.,350.);
1646  m_peta_trklengthTS->SetXTitle("eta");
1647  m_peta_trklengthTS->SetYTitle("length");
1648  m_pnpoint_lengthTS = QAH::H2F("QaPtrkNPntLTS","primtrk: N pnts vs length, tpc+svt",25,70.,350.,25,0.,50.);
1649  m_pnpoint_lengthTS->SetXTitle("trk length");
1650  m_pnpoint_lengthTS->SetYTitle("Npoints on trk");
1651  m_pfpoint_lengthTS = QAH::H2F("QaPtrkFitPntLTS","primtrk: N fit pnts vs length, tpc+svt",25,70.,350.,25,0.,50.);
1652  m_pfpoint_lengthTS->SetXTitle("trk length");
1653  m_pfpoint_lengthTS->SetYTitle("Npoints on trk");
1654  m_pchisq0_momTS = QAH::H2F("QaPtrkChi0MomTS","primtrk: Chisq0 vs log10 mom, tpc+svt",40,1.,4.,50,0.,5.);
1655  m_pchisq0_momTS->SetXTitle("log10 P (MeV)");
1656  m_pchisq0_momTS->SetYTitle("chisq0") ;
1657  m_pchisq0_etaTS = QAH::H2F("QaPtrkChi0EtaTS","primtrk: Chisq0 vs eta, tpc+svt",20,-2.,2.,20,0.,5.);
1658  m_pchisq0_etaTS->SetXTitle("eta");
1659  m_pchisq0_etaTS->SetYTitle("chisq0");
1660  m_pchisq0_dipTS = QAH::H2F("QaPtrkChi0TanlTS","primtrk: Chisq0 vs tanl(dip), tpc+svt",20,-3.2,3.2,20,0.,5.);
1661  m_pchisq0_dipTS->SetXTitle("dip angle");
1662  m_pchisq0_dipTS->SetYTitle("chisq0");
1663  m_pchisq0_zfTS = QAH::H2F("QaPtrkChi0zfTS","primtrk: Chisq0 vs zfirst, tpc+svt",20,-200.,200.,20,0.,5.);
1664  m_pchisq0_zfTS->SetXTitle("zfirst");
1665  m_pchisq0_zfTS->SetYTitle("chisq0");
1666  m_ppsi_phiTS = QAH::H2F("QaPtrkPsiPhiTS","primtrk: psi vs phi, tpc+svt",36, 0.,360.,36,0.,360.);
1667  m_ppsi_phiTS->SetXTitle("phi");
1668  m_ppsi_phiTS->SetYTitle("psi");
1669  } // silHists
1670 
1671  if (ftpHists) { // FTPC hists
1672 
1673 // 1D ftpc
1674  // east (solid) and west(dashed) on same plot
1675  m_ppointF = QAH::MH1F("QaPtrkNPntF", "primtrk: N points on trk,ftpc",15,0.,15.,2);
1676  m_ppointF->Rebin(0,"East");
1677  m_ppointF->Rebin(1,"West");
1678  m_ppointF->SetStats(kFALSE);
1679  m_pmax_pointF = QAH::MH1F("QaPtrkNPntMaxF", "primtrk: N max pnts on trk, ftpc",15,0.,15.,2);
1680  m_pmax_pointF->Rebin(0,"East");
1681  m_pmax_pointF->Rebin(1,"West");
1682  m_pmax_pointF->SetStats(kFALSE);
1683  m_prim_ratiomF = QAH::MH1F("QaPtrkRnmF", "primtrk: ratio Nfit/max pnt, ftpc",55,0.,1.1,2);
1684  m_prim_ratiomF->Rebin(0,"East");
1685  m_prim_ratiomF->Rebin(1,"West");
1686  m_prim_ratiomF->SetStats(kFALSE);
1687  m_prim_chargeF = QAH::MH1F("QaPtrkChrgF", "primtrk: charge, ftpc",20,-2.,2.,2);
1688  m_prim_chargeF->Rebin(0,"East");
1689  m_prim_chargeF->Rebin(1,"West");
1690  m_prim_chargeF->SetStats(kFALSE);
1691  m_prim_xfF = QAH::MH1F("QaPtrkXfF", "primtrk: x of first hit on trk, ftpc",70,-35.,35.,2);
1692  m_prim_xfF->Rebin(0,"East");
1693  m_prim_xfF->Rebin(1,"West");
1694  m_prim_xfF->SetStats(kFALSE);
1695  m_prim_yfF = QAH::MH1F("QaPtrkYfF", "primtrk: y of first hit on trk, ftpc",70,-35.,35.,2);
1696  m_prim_yfF->Rebin(0,"East");
1697  m_prim_yfF->Rebin(1,"West");
1698  m_prim_yfF->SetStats(kFALSE);
1699  m_prim_zfF = QAH::MH1F("QaPtrkZfF", "primtrk: z of first hit on trk, ftpc",75,-225.,225.,2);
1700  m_prim_zfF->Rebin(0,"East");
1701  m_prim_zfF->Rebin(1,"West");
1702  m_prim_zfF->SetStats(kFALSE);
1703  m_prim_radfF = QAH::MH1F("QaPtrkRF", "primtrk: radial position of first hit, ftpc",35,0.,35.,2);
1704  m_prim_radfF->Rebin(0,"East");
1705  m_prim_radfF->Rebin(1,"West");
1706  m_prim_radfF->SetStats(kFALSE);
1707  m_plengthF = QAH::MH1F("QaPtrkLengthF", "primtrk: track length, ftpc",50,150.,300.,2);
1708  m_plengthF->Rebin(0,"East");
1709  m_plengthF->Rebin(1,"West");
1710  m_plengthF->SetStats(kFALSE);
1711  m_ppsiF = QAH::MH1F("QaPtrkPsiF", "primtrk: psi, ftpc", 90, 0.,360.,2);
1712  m_ppsiF->Rebin(0,"East");
1713  m_ppsiF->Rebin(1,"West");
1714  m_ppsiF->SetStats(kFALSE);
1715  m_petaF = QAH::MH1F("QaPtrkEtaF", "primtrk: |eta|, ftpc",80,2,5,2);
1716  m_petaF->Rebin(0,"East");
1717  m_petaF->Rebin(1,"West");
1718  m_petaF->SetStats(kFALSE);
1719  m_ppTF = QAH::MH1F("QaPtrkPtF", "primtrk: pT, ftpc",50,0.,5.,2);
1720  m_ppTF->Rebin(0,"East");
1721  m_ppTF->Rebin(1,"West");
1722  m_ppTF->SetStats(kFALSE);
1723  m_pmomF = QAH::MH1F("QaPtrkPF", "primtrk: momentum, ftpc",50,0.,5.,2);
1724  m_pmomF->Rebin(0,"East");
1725  m_pmomF->Rebin(1,"West");
1726  m_pmomF->SetStats(kFALSE);
1727  m_pchisq0F = QAH::MH1F("QaPtrkChisq0F", "primtrk: chi2/ndf x,y fit, ftpc",50,0.,5.,2);
1728  m_pchisq0F->Rebin(0,"East");
1729  m_pchisq0F->Rebin(1,"West");
1730  m_pchisq0F->SetStats(kFALSE);
1731  m_pchisq1F = QAH::MH1F("QaPtrkChisq1F", "primtrk: chi2/ndf r,z fit, ftpc",50,0.,5.,2);
1732  m_pchisq1F->Rebin(0,"East");
1733  m_pchisq1F->Rebin(1,"West");
1734  m_pchisq1F->SetStats(kFALSE);
1735  // separate east and west histograms
1736  m_ppointFE = QAH::H1F("QaPtrkNPntFE", "primtrk: N points on trk,ftpc east", 15, 0.,15.);
1737  m_ppointFW = QAH::H1F("QaPtrkNPntFW", "primtrk: N points on trk,ftpc west", 15, 0.,15.);
1738  m_pmax_pointFE = QAH::H1F("QaPtrkNPntMaxFE", "primtrk: N max pnts on trk, ftpc east", 15, 0.,15.);
1739  m_pmax_pointFW = QAH::H1F("QaPtrkNPntMaxFW", "primtrk: N max pnts on trk, ftpc west", 15, 0.,15.);
1740  m_prim_ratiomFE = QAH::H1F("QaPtrkRnmFE", "primtrk: ratio Nfit/max pnt, ftpc east", 55, 0., 1.1);
1741  m_prim_ratiomFW = QAH::H1F("QaPtrkRnmFW", "primtrk: ratio Nfit/max pnt, ftpc west", 55, 0., 1.1);
1742  m_prim_chargeFE = QAH::H1F("QaPtrkChrgFE", "primtrk: charge, ftpc east ", 20,-2.,2.);
1743  m_prim_chargeFW = QAH::H1F("QaPtrkChrgFW", "primtrk: charge, ftpc west ", 20,-2.,2.);
1744  m_prim_xfFE = QAH::H1F("QaPtrkXfFE", "primtrk: x of first hit on trk, ftpc east",70,-35.,35.);
1745  m_prim_xfFW = QAH::H1F("QaPtrkXfFW", "primtrk: x of first hit on trk, ftpc west",70,-35.,35.);
1746  m_prim_yfFE = QAH::H1F("QaPtrkYfFE", "primtrk: y of first hit on trk, ftpc east",70,-35.,35.);
1747  m_prim_yfFW = QAH::H1F("QaPtrkYfFW", "primtrk: y of first hit on trk, ftpc west",70,-35.,35.);
1748  m_prim_zfFE = QAH::H1F("QaPtrkZfFE", "primtrk: z of first hit on trk, ftpc east",75,-225.,-150.);
1749  m_prim_zfFW = QAH::H1F("QaPtrkZfFW", "primtrk: z of first hit on trk, ftpc west",75,150.,225.);
1750  m_prim_radfFE = QAH::H1F("QaPtrkRFE", "primtrk: radial position of first hit, ftpc east",35,0.,35.);
1751  m_prim_radfFW = QAH::H1F("QaPtrkRFW", "primtrk: radial position of first hit, ftpc west",35,0.,35.);
1752  m_plengthFE = QAH::H1F("QaPtrkLengthFE", "primtrk: track length, ftpc east",50,150.,300.);
1753  m_plengthFW = QAH::H1F("QaPtrkLengthFW", "primtrk: track length, ftpc west",50,150.,300.);
1754  m_ppsiFE = QAH::H1F("QaPtrkPsiFE", "primtrk: psi, ftpc east", 90, 0.,360.);
1755  m_ppsiFW = QAH::H1F("QaPtrkPsiFW", "primtrk: psi, ftpc west", 90, 0.,360.);
1756  m_petaFE = QAH::H1F("QaPtrkEtaFE", "primtrk: eta, ftpc east",80,-4.5,-2.);
1757  m_petaFW = QAH::H1F("QaPtrkEtaFW", "primtrk: eta, ftpc west",80,2.,4.5);
1758  m_ppTFE = QAH::H1F("QaPtrkPtFE", "primtrk: pT, ftpc east",50,0.,5.);
1759  m_ppTFW = QAH::H1F("QaPtrkPtFW", "primtrk: pT, ftpc west",50,0.,5.);
1760  m_pmomFE = QAH::H1F("QaPtrkPFE", "primtrk: momentum, ftpc east ",50,0.,5.);
1761  m_pmomFW = QAH::H1F("QaPtrkPFW", "primtrk: momentum, ftpc west ",50,0.,5.);
1762  m_pchisq0FE = QAH::H1F("QaPtrkChisq0FE", "primtrk: chi2/ndf x,y fit, ftpc east", 50, 0.,5.);
1763  m_pchisq0FW = QAH::H1F("QaPtrkChisq0FW", "primtrk: chi2/ndf x,y fit, ftpc west", 50, 0.,5.);
1764  m_pchisq1FE = QAH::H1F("QaPtrkChisq1FE", "primtrk: chi2/ndf r,z fit, ftpc east", 50, 0.,5.);
1765  m_pchisq1FW = QAH::H1F("QaPtrkChisq1FW", "primtrk: chi2/ndf r,z ft, ftpc west", 50, 0.,5.);
1766 
1767 // 2D - ftpc
1768  m_ppT_eta_recFE = QAH::H2F("QaPtrkPtVsEtaFE","primtrk: log10 pT vs eta, ftpcE",20,-4.5,-2.,40,1.,4.);
1769  m_ppT_eta_recFW = QAH::H2F("QaPtrkPtVsEtaFW","primtrk: log10 pT vs eta, ftpcW",20,2.,4.5,40,1.,4.);
1770  m_primtrk_xf_yfFE = QAH::H2F("QaPtrkXfYfFE","primtrk: Y vs X of first hit on trk, ftpcE",70,-35.,35.,75,-35.,35.);
1771  m_primtrk_xf_yfFE->SetXTitle("x first");
1772  m_primtrk_xf_yfFE->SetYTitle("y first");
1773  m_primtrk_xf_yfFW = QAH::H2F("QaPtrkXfYfFW","primtrk: Y vs X of first hit on trk, ftpcW",70,-35.,35.,70,-35.,35.);
1774  m_primtrk_xf_yfFW->SetXTitle("x first");
1775  m_primtrk_xf_yfFW->SetYTitle("y first");
1776  m_peta_trklengthFE = QAH::H2F("QaPtrkLengthVEtaFE","primtrk: trk length vs eta, ftpcE",25,-4.5,-2.,50,0.,300.);
1777  m_peta_trklengthFE->SetXTitle("eta");
1778  m_peta_trklengthFE->SetYTitle("length");
1779  m_peta_trklengthFW = QAH::H2F("QaPtrkLengthVEtaFW","primtrk: trk length vs eta, ftpcW",25,2.,4.5,50,0.,300.);
1780  m_peta_trklengthFW->SetXTitle("eta");
1781  m_peta_trklengthFW->SetYTitle("length");
1782  m_pnpoint_lengthFE = QAH::H2F("QaPtrkNPntLengthFE","primtrk: N pnts vs length, ftpcE",50,0.,300.,15,0.,15.);
1783  m_pnpoint_lengthFE->SetXTitle("trk length");
1784  m_pnpoint_lengthFE->SetYTitle("Npoints on trk");
1785  m_pnpoint_lengthFW = QAH::H2F("QaPtrkNPntLengthFW","primtrk: N pnts vs length, ftpcW",50,0.,300.,15,0.,15.);
1786  m_pnpoint_lengthFW->SetXTitle("trk length");
1787  m_pnpoint_lengthFW->SetYTitle("Npoints on trk");
1788 
1789  } // ftpHists
1790 
1791 // 2D - SVT drift length
1792  if (silHists)
1793  m_svt_loc = QAH::H2F("QaPtrkSvtLoc","primtrk: SVT hit time bins",256,0,128,432,-0.5,431.5);
1794 
1795  m_primary_pxl_hit = QAH::H1F("QaPtrkPxlHit","PIXEL: Hits per primary track",5,-0.5,4.5);
1796  m_primary_pxl_hit->SetXTitle("hits in PIXEL");
1797  m_primary_ist_hit = QAH::H1F("QaPtrkIstHit","IST: Hits per primary track",5,-0.5,4.5);
1798  m_primary_ist_hit->SetXTitle("hits in IST");
1799  m_primary_hft_hit = QAH::H1F("QaPtrkHftHit","HFT: Hits per primary track",5,-0.5,4.5);
1800  m_primary_hft_hit->SetXTitle("hits in HFT subsystem");
1801 
1802  m_primary_hft_hit->GetXaxis()->SetBinLabel(1,"");
1803  m_primary_hft_hit->GetXaxis()->SetBinLabel(2,"PIXEL (inner)");
1804  m_primary_hft_hit->GetXaxis()->SetBinLabel(3,"PIXEL (outer)");
1805  m_primary_hft_hit->GetXaxis()->SetBinLabel(4,"IST");
1806  m_primary_hft_hit->GetXaxis()->SetBinLabel(5,"SSD");
1807 
1808 }
1809 //_____________________________________________________________________________
1810 void StQABookHist::BookHistDE(){
1811 
1812  // for method MakeDE - from table dst_dedx
1813  m_ndedxr = QAH::H1F("QaDedxNum", "dedx: number of tracks", 64,0., 16000.);
1814 
1815  m_ndedxT = QAH::H1F("QaDedxNdedxT", "dedx: number of point to define dE/dx, tpc",50,0., 50.);
1816  m_dedx0T = QAH::H1F("QaDedxDedx0T", "dedx: dE/dx mean (GeV/cm), tpc", ndedx, 0., 0.0005);
1817  m_dedx1T = QAH::H1F("QaDedxDedx1T", "dedx: dE/dx error on mean, tpc", ndedx, 0., 1);
1818  //m_dedxTTS = QAH::H1F("QaDedxBBTTS", "dedx: ratio <dE/dx> mean to Bethe-Bloch <dE/dx>, tpc,tpc+svt",
1819  // 50,0.,10.);
1820  m_dedxTTS = QAH::H1F("QaDedxBTTS", "dedx: ratio <dE/dx> mean to Bichsel <dE/dx>, tpc,tpc+svt",
1821  50,0.,10.);
1822  // east and west on same plot
1823  m_ndedxF = QAH::MH1F("QaDedxNdedxF", "dedx: number of point to define dE/dx, ftpc",50,0.,10.,2);
1824  m_ndedxF->Rebin(0,"East");
1825  m_ndedxF->Rebin(1,"West");
1826  m_ndedxF->SetStats(kFALSE);
1827  m_dedx0F = QAH::MH1F("QaDedxDedx0F", "dedx: dE/dx mean, ftpc",ndedx,0.,750.,2);
1828  m_dedx0F->Rebin(0,"East");
1829  m_dedx0F->Rebin(1,"West");
1830  m_dedx0F->SetStats(kFALSE);
1831  // east and west on separate plots
1832  m_ndedxFE = QAH::H1F("QaDedxNdedxFE", "dedx: number of point to define dE/dx, ftpcE",50,0.,10.);
1833  m_dedx0FE = QAH::H1F("QaDedxDedx0FE", "dedx: dE/dx mean, ftpcE", ndedx, 0., 750.);
1834 
1835  m_ndedxFW = QAH::H1F("QaDedxNdedxFW", "dedx: number of point to define dE/dx, ftpcW",50,0.,10.);
1836  m_dedx0FW = QAH::H1F("QaDedxDedx0FW", "dedx: dE/dx mean, ftpcW", ndedx, 0., 750.);
1837 }
1838 //_____________________________________________________________________________
1839 void StQABookHist::BookHistPID(){
1840 
1841  // for MakeHistPID - from tables globtrk & dst_dedx
1842  // Spectra/pid histograms. C.Ogilvie
1843 
1844  m_p_dedx_rec = QAH::H2F("QaPidGlobtrkDstdedxPVsDedx","PID: globtrk-dst_dedx, p vs dedx (reconstructed)",
1845  cnp,cminp,cmaxp,cndedx,cmindedx,cmaxdedx);
1846  m_p_dedx_rec->SetXTitle("p (GeV)");
1847  m_p_dedx_rec->SetYTitle("dedx (keV/cm)");
1848 
1849 }
1850 //_____________________________________________________________________________
1851 void StQABookHist::BookHistVertex(){
1852  // for MakeHistVertex - from table dst_vertex
1853 
1854 
1855  m_v_num = QAH::H1F("QaVtxNum", " vertex: num vertices ",50,0.,10000.);
1856  m_v_num_sm= QAH::H1F("QaVtxNumsm", " vertex: num vertices ",50,0.,500.);
1857 
1858  m_v_vtxid = QAH::H1F("QaVtxVtxId"," vertex,2ndary: Vertex ID ",10,0.,10.);
1859  m_v_x = QAH::H1F("QaVtxX", " vertex,2ndary: x ",50,-250.,250.);
1860  m_v_y = QAH::H1F("QaVtxY", " vertex,2ndary: y ",50,-250.,250.);
1861  m_v_z = QAH::H1F("QaVtxZ", " vertex,2ndary: z ",50,-250.,250.);
1862  m_v_pchi2 = QAH::H1F("QaVtxChisq"," vertex,2ndary: chisq/dof ",50,0.,50.);
1863  m_v_r = QAH::H1F("QaVtxR", " vertex,2ndary: r ",100,0,50);
1864 
1865  m_pv_vtxid = QAH::H1F("QaVtxPrVtxId"," vertex,prim: Vertex ID ",10,0.,10.);
1866  m_pv_x = QAH::H1F("QaVtxPrX", " vertex,prim: x ",50,-5.,5.);
1867  m_pv_y = QAH::H1F("QaVtxPrY", " vertex,prim: y ",50,-5.,5.);
1868  m_pv_z = QAH::H1F("QaVtxPrZ", " vertex,prim: z ",50,-50.,50.);
1869  m_pv_xy = QAH::H2F("QaVtxPrXY", " vertex,prim: x versus y",50,-5.,5.,50,-5.,5.);
1870  m_pv_pchi2 = QAH::H1F("QaVtxPrChisq"," vertex,prim: chisq/dof ",40,0.,20.);
1871  m_pv_r = QAH::H1F("QaVtxPrR", " vertex,prim: r ",100,0,0.1);
1872 
1873 // m_vtx_phi_dist = QAH::H2F("QaV0VtxPhiDist",
1874 // "V0 azimuthal distribution",36,0.,360.,25,1.,101.);
1875 // m_vtx_phi_dist->SetXTitle("Mean of phi(V0)");
1876 // m_vtx_phi_dist->SetYTitle("RMS of phi(V0)");
1877  m_vtx_phi_dist = QAH::H1F("QaV0VtxPhiDist",
1878  "V0 azimuthal distribution relative to primvtx",36,0.,360.);
1879  m_vtx_r_dist = QAH::H1F("QaV0VtxRDist",
1880  "V0 radial distribution relative to primvtx",50,0.,200.);
1881  m_vtx_z_dist = QAH::H1F("QaV0VtxZDist",
1882  "V0 Z distribution relative to primvtx",60,-30.,30.);
1883 
1884  m_v0 = QAH::H1F("QaV0Vtx","log10 total V0s (0 => -0.5)",45,-0.55,3.95);
1885  m_ev0_lama_hist = QAH::H1F("QaV0LambdaMass","V0: Lambda mass",25,1.05,1.15);
1886  m_ev0_k0ma_hist = QAH::H1F("QaV0K0Mass","V0: k0 mass",25,.4,.6);
1887 
1888  m_xi_tot = QAH::H1F("QaXiVtxTot", "log10 total Xis (0 => -0.5)",45,-0.55,3.95);
1889  m_xi_ma_hist = QAH::H1F("QaXiaMass", "Xi: Xi mass",25,1.2,1.4);
1890 
1891  m_kink_tot = QAH::H1F("QaKinkTot", "log10 total kinks (0 => -0.5)",35,-0.55,2.95);
1892 
1893  m_vtx_FtpcEastTpc_xy = QAH::H2F("QaVtxFtpcETpcXY",
1894  " vertex,prim: x(ftpcE)-x(tpc) vs y(ftpcE)-y(tpc)",
1895  800, -20., 20., 800, -20., 20.);
1896  m_vtx_FtpcEastTpc_z = QAH::H1F("QaVtxFtpcETpcZ",
1897  " vertex,prim: z(ftpcE)-z(tpc)",
1898  100, -10., 10.);
1899  m_vtx_FtpcWestTpc_xy = QAH::H2F("QaVtxFtpcWTpcXY",
1900  " vertex,prim: x(ftpcW)-x(tpc) vs y(ftpcW)-y(tpc)",
1901  800, -20., 20., 800, -20., 20.);
1902  m_vtx_FtpcWestTpc_z = QAH::H1F("QaVtxFtpcWTpcZ",
1903  " vertex,prim: z(ftpcW)-z(tpc)",
1904  100, -10., 10.);
1905 
1906 }
1907 //_____________________________________________________________________________
1908 void StQABookHist::BookHistPoint(){
1909 
1910  m_pnt_tot = QAH::H1F("QaPointTot", "point: # hits total ",100, 0.,400000.);
1911  m_pnt_tot_med = QAH::H1F("QaPointTotmed","point: # hits total ",100, 0.,25000.);
1912  m_pnt_tot_sm = QAH::H1F("QaPointTotsm", "point: # hits total ",100, 0.,2500.);
1913  m_pnt_id = QAH::H1F("QaPointId","point: detector ID of hit",30,0.,30.);
1914 
1915  m_pnt_tpc = QAH::H1F("QaPointTpc", "point: # hits tpc ",100, 0.,300000.);
1916  // east and west on same plot
1917  m_pnt_ftpc = QAH::MH1F("QaPointFtpc", "point: # hits ftpc",100,0.,25000.,2);
1918  m_pnt_ftpc->Rebin(0,"East");
1919  m_pnt_ftpc->Rebin(1,"West");
1920  m_pnt_ftpc->SetStats(kFALSE);
1921  // east and west on separate plots
1922  m_pnt_ftpcE = QAH::H1F("QaPointFtpcE","point: # hits ftpcE ",100, 0.,25000.);
1923  m_pnt_ftpcW = QAH::H1F("QaPointFtpcW","point: # hits ftpcW ",100, 0.,25000.);
1924  if (silHists) {
1925  m_pnt_svt = QAH::H1F("QaPointSvt", "point: # hits svt ",600, 0.,15000.);
1926  m_pnt_svtLaser= QAH::H2F("QaPointSvtLaser","point: laser spots, svt ",150,0,600,65,0.,130.);
1927  m_pnt_svtLaser->SetXTitle("event in file");
1928  m_pnt_svtLaserDiff= QAH::MH2F("QaPointSvtLaserDiff","point: diff of laser spots, svt ",150,0,600,101,9.8,50.2,2);
1929  m_pnt_svtLaserDiff->SetXTitle("event in file");
1930  m_pnt_svtLaserDiff->SetYTitle("time buckets");
1931  m_pnt_svtLaserDiff->Rebin(0,"Laser 1");
1932  m_pnt_svtLaserDiff->Rebin(1,"Laser 2");
1933  }
1934  if (silHists || (histsSet>=StQA_run14 && histsSet<StQA_run17)) {
1935  m_pnt_xyS = QAH::H2F("QaPointXYSvt","point: x-y distribution of hits, svt,ssd",125,-25,25,125,-25,25);
1936  m_pnt_xyS->SetStats(kFALSE);
1937  m_pnt_xyS->SetXTitle("x [cm]");
1938  m_pnt_xyS->SetYTitle("y [cm]");
1939 
1940  // for run >= 14, change the title of the histogram with vertex detector hit distribution in x-y
1941  // FIXME: move this part to BookHistPoint and pass histsSet there?
1942  if (histsSet>=StQA_run14 && histsSet<StQA_run17)
1943  m_pnt_xyS->SetTitle("PIXEL, IST, SSD: Distribution of hits in XY");
1944  }
1945  // Now using polar coords (lego didn't work well because of inner radius and colored zero)
1946  // Drawing a lego plot requires r-phi,r. Drawing a plain polar plot reguirs r,r-phi
1947  m_pnt_rpTW = QAH::H2F("QaPointRPTpcW","point: r-phi distribution of hits, tpcW",20,58.75,196.75,72,0,TMath::TwoPi());
1948  m_pnt_rpTE = QAH::H2F("QaPointRPTpcE","point: r-phi distribution of hits, tpcE",23,58.75,196.75,72,0,TMath::TwoPi());
1949  m_z_hits = QAH::H1F("QaPointZhits","point: z distribution of hits, tpc",100,-210,210);
1950  m_pnt_timeT = QAH::MH1F("QaPointTimeT","point: time bucket distribution of hits, tpc",45,0,450,2);
1951  m_pnt_timeT->Rebin(0,"East");
1952  m_pnt_timeT->Rebin(1,"West");
1953  m_pnt_timeT->SetStats(kFALSE);
1954  m_pnt_phiT = QAH::MH1F("QaPointPhiT","point: #phi distribution of hits, tpc",36,0,360,2);
1955  m_pnt_phiT->Rebin(0,"East");
1956  m_pnt_phiT->Rebin(1,"West");
1957  m_pnt_phiT->SetStats(kFALSE);
1958  m_pnt_padrowT = QAH::MH1F("QaPointPadrowT","point: padrow distribution of hits, tpc",72,0.5,72.5,2);
1959  m_pnt_padrowT->Rebin(0,"East");
1960  m_pnt_padrowT->Rebin(1,"West");
1961  m_pnt_padrowT->SetStats(kFALSE);
1962  m_pnt_padrowT->SetXTitle("padrow number");
1963 
1964  if (silHists) {
1965  m_pnt_zS = QAH::H1F("QaPointZhitsS","point: z distribution of hits, svt",100,-35,35);
1966  m_pnt_phiS = QAH::H1F("QaPointPhiS","point: #phi distribution of hits, svt",36,0,360);
1967  m_pnt_barrelS = QAH::H1F("QaPointBarrelS","point: barrel distribution of hits, svt",3,0.5,3.5);
1968  m_pnt_barrelS->SetXTitle("barrel number");
1969  }
1970  if (silHists) {
1971  m_pnt_ssd = QAH::H1F("QaPointSsd", "point: # hits ssd ",200, 0.,5000.);
1972  m_pnt_phiSSD = QAH::H1F("QaPointPhiSSD","SSD: #phi of hits (per event)",36,0,360);
1973  m_pnt_lwSSD = QAH::H2F("QaPointLWSSD","SSD: wafer id vs ladder id (per event)",20,0.5,20.5,16,0.5,16.5);
1974  m_pnt_lwSSD->SetXTitle("Ladder #");
1975  m_pnt_lwSSD->SetYTitle("Wafer (sensor) #");
1976  m_pnt_lwSSD->SetStats(kFALSE);
1977  m_pnt_sizeSSD = QAH::MH1F("QaPointSizeSSD","SSD: size of clusters",10,0.5,10.5,2);
1978  m_pnt_sizeSSD->Rebin(0,"P-side");
1979  m_pnt_sizeSSD->Rebin(1,"N-side");
1980  m_pnt_sizeSSD->SetStats(kFALSE);
1981  m_pnt_eSSD = QAH::H1F("QaPointESSD","SSD: log10(energy) of hits",90,-5,-2);
1982  }
1983  if (histsSet>=StQA_run14 && histsSet<StQA_run17) {
1984  m_pnt_sst = QAH::H1F("QaPointSst", "point: # hits sst ",200, 0.,5000.);
1985  m_pnt_phiSST = QAH::H1F("QaPointPhiSST","SST: #phi of hits (per event)",36,0,360);
1986  m_pnt_lwSST = QAH::H2F("QaPointLWSST","SST: wafer id vs ladder id (per event)",20,0.5,20.5,16,0.5,16.5);
1987  m_pnt_lwSST->SetXTitle("Ladder #");
1988  m_pnt_lwSST->SetYTitle("Wafer (sensor) #");
1989  m_pnt_lwSST->SetStats(kFALSE);
1990  m_pnt_sizeSST = QAH::MH1F("QaPointSizeSST","SST: size of clusters",10,0.5,10.5,2);
1991  m_pnt_sizeSST->Rebin(0,"P-side");
1992  m_pnt_sizeSST->Rebin(1,"N-side");
1993  m_pnt_sizeSST->SetStats(kFALSE);
1994  m_pnt_eSST = QAH::H1F("QaPointESST","SST: log10(energy) of hits",90,-5,-2);
1995  }
1996 
1997  m_pnt_xyFE = QAH::H2F("QaPointXYFtpcE","point: x-y distribution of hits, ftpcE",70,-35,35,70,-35,35);
1998  m_pnt_xyFW = QAH::H2F("QaPointXYFtpcW","point: x-y distribution of hits, ftpcW",70,-35,35,70,-35,35);
1999  m_pnt_planeF = QAH::MH1F("QaPointPlaneF","point: plane distribution of hits, ftpc",20,0.5,20.5,2);
2000  m_pnt_planeF->Rebin(0,"East");
2001  m_pnt_planeF->Rebin(1,"West");
2002  m_pnt_planeF->SetStats(kFALSE);
2003  m_pnt_planeF->SetXTitle("plane number");
2004  m_pnt_padtimeFE = QAH::H2F("QaPointPadTimeFtpcE","point: #pads vs #timebins of hits, ftpcE",12,0.5,12.5,10,0.5,10.5);
2005  m_pnt_padtimeFE->SetXTitle("#timebins");
2006  m_pnt_padtimeFE->SetYTitle("#pads");
2007  m_pnt_padtimeFW = QAH::H2F("QaPointPadTimeFtpcW","point: #pads vs #timebins of hits, ftpcW",12,0.5,12.5,10,0.5,10.5);
2008  m_pnt_padtimeFW->SetXTitle("#timebins");
2009  m_pnt_padtimeFW->SetYTitle("#pads");
2010 
2011 
2012 }
2013 //_____________________________________________________________________________
2014 void StQABookHist::BookHistEMC(){
2015 
2016  if (!((gROOT->GetClass("StEmcMath")) && (gROOT->GetClass("StEmcGeom"))))
2017  return;
2018 // Book the hists for SimulatorMaker
2019  m_emc_nhit=QAH::H2F("EmcNHitsVsDet","emc: Number of hit(log) .vs. Detector #",100,0.0,4.5,4,0.5,4.5);
2020  m_emc_etot=QAH::H2F("EmcEtotVsDet","emc: Total energy(log) .vs. Detector #",100,-4.0,4.5,8,0.5,4.5);
2021 
2022  const Char_t* tit={"Barrel"};
2023  const Int_t nx[4] = {40,40,300,20};
2024  const Int_t ny[4] = {120, 120, 60, 900};
2025 
2026  // for BSMD per-module histos
2027  Axis_t ModNumLo[4] = {1.,31.,61.,91.};
2028  Axis_t ModNumHi[4] = {31.,61.,91.,121.};
2029  const TString PerModuleHistName[4] = {"West1","West2","East1","East2"};
2030 
2031  Float_t rpi = M_PI + 0.00001;
2032  TString name, title;
2033  TArrayD *xb = StEmcMath::binForSmde();
2034 
2035  for(Int_t i=0; i<4; i++){
2036  name = detname[i] + "Hits";
2037  title = tit + detname[i] + " hits dist.";
2038  if(i==2) m_emc_hits[i]=QAH::H2F(name,title, xb->GetSize()-1,xb->GetArray(),ny[i],-rpi,rpi);
2039  else m_emc_hits[i] = QAH::H2F(name,title, nx[i],-1.,+1., ny[i],-rpi, rpi);
2040 
2041  name = detname[i] + "Energy2D";
2042  title = tit + detname[i] + " energy dist. in eta-phi";
2043  if(i==2) m_emc_energy2D[i] = QAH::H2F(name,title, xb->GetSize()-1,xb->GetArray(), ny[i],-rpi,rpi);
2044  else m_emc_energy2D[i] = QAH::H2F(name,title, nx[i],-1.,+1., ny[i],-rpi, rpi);
2045 
2046  name = detname[i] + "Adc";
2047  title = tit + detname[i] + " ADC dist.";
2048  m_emc_adc[i] = QAH::H1F(name,title, 5001, -0.5, 5000.5);
2049 
2050  name = detname[i] + "Energy";
2051  title = tit + detname[i] + " energy dist.";
2052  m_emc_energy[i] = QAH::H1F(name,title, 600, 0.0, 60.0);
2053  }
2054  delete xb;
2055 
2056  for(Int_t i=2; i<4; i++){ // Detector ID for BSMDE and BSMDP (BSMD eta and phi)
2057  for (Int_t j=0; j<4; j++){ // split 120 modules into 4 histos
2058 
2059  Int_t k = j;
2060  if (i>2) k = j + 4;
2061 
2062  name = detname[i] + PerModuleHistName[j] + "HitsPerModule";
2063  title = tit + detname[i] + " " + PerModuleHistName[j] + " - hits per module";
2064  m_emc_hits_per_module[k] = QAH::H1F(name,title, 30, ModNumLo[j], ModNumHi[j]);
2065 
2066  name = detname[i] + PerModuleHistName[j] + "EnergyPerModule";
2067  title = tit + detname[i] + " " + PerModuleHistName[j] + " - energy-weighted hits per module";
2068  m_emc_energy_per_module[k] = QAH::H1F(name,title, 30, ModNumLo[j], ModNumHi[j]);
2069 
2070  name = detname[i] + PerModuleHistName[j] + "StripHitsPerModule";
2071  title = tit + detname[i] + " " + PerModuleHistName[j] + " - hits in strip (within mod) vs. mod";
2072  m_emc_strip_hits_per_module[k] = QAH::H2F(name,title, 30, ModNumLo[j], ModNumHi[j], 150, 1., 151.);
2073 
2074  name = detname[i] + PerModuleHistName[j] + "StripEnergyPerModule";
2075  title = tit + detname[i] + " " + PerModuleHistName[j] + " - energy in strip (within mod) vs. mod";
2076  m_emc_strip_energy_per_module[k] = QAH::H2F(name,title, 30, ModNumLo[j], ModNumHi[j], 150, 1., 151.);
2077  }
2078  }
2079 
2080 // Book the hists for cluster finder
2081  Int_t greta[4]={40,40,300,20}; // eta bins
2082  Int_t grphi[4]={120,120,60,900}; // phi bins => 16-apr by PAI
2083  Float_t myPI = M_PI + 0.0001;
2084 
2085  m_emc_ncl=QAH::H2F("EmcNcluster","emc: Number of cluster(log) .vs. Detector #",40,0.0,4.0, 4,0.5,4.5);
2086  m_emc_etotCl=QAH::H2F("EmcEcluster" ,"emc: Total PreCluster Energy(log) .vs. Detector #", 60,-2.0,4.0, 4,0.5,4.5);
2087 
2088  Float_t rmsMax=0.026;
2089  Int_t rmsN=52;
2090  m_emc_sig_e= QAH::H2F("EmcRMSeta" ,"emc: Sigma(eta) .vs. Detector #",rmsN,0.0,rmsMax,4,0.5,4.5);
2091  m_emc_sig_p= QAH::H2F("EmcRMSphi" ,"emc: Sigma(phi) .vs. Detector #",rmsN,0.0,rmsMax,4,0.5,4.5);
2092  for (Int_t i=0; i<4; i++)
2093  {
2094  TString name_h = detname[i] + "_cluster";
2095  TString name_e = detname[i] + "_cluster_energy";
2096  TString tit_h = detname[i] + " cluster";
2097  TString tit_e = detname[i] + " energy of cluster";
2098  if(i==2) {
2099  m_emc_cl[i] = QAH::H2F(name_h,tit_h,greta[i],-1.0,1.0,grphi[i],-M_PI*1.015, M_PI*0.985);
2100  m_emc_energyCl[i] = QAH::H2F(name_e,tit_e,greta[i],-1.0,1.0,grphi[i],-M_PI*1.015, M_PI*0.985);
2101  } else {
2102  m_emc_cl[i] = QAH::H2F(name_h,tit_h,greta[i],-1.0,1.0,grphi[i],-myPI, myPI);
2103  m_emc_energyCl[i] = QAH::H2F(name_e,tit_e,greta[i],-1.0,1.0,grphi[i],-myPI, myPI);
2104  }
2105 
2106 
2107  name = detname[i] + "ClNum";
2108  title = "Number hits in cluster for " + detname[i];
2109  m_emc_HitsInCl[i] = QAH::H1F(name, title, 21, -0.5, 20.5);
2110 
2111  name = detname[i] + "ClEnergy";
2112  title = "Energy of cluster for " + detname[i];
2113  m_emc_EnergyCl[i] = QAH::H1F(name, title, 2000, 0.0, 20.0);
2114 
2115  TString name_eta = detname[i] + "Eta";
2116  TString tit_eta = "Eta of clusters for " + detname[i];
2117  TString name_phi = detname[i] + "Phi";
2118  TString tit_phi = "Phi of clusters for " + detname[i];
2119  if(i==2) {
2120  TArrayD *xb = StEmcMath::binForSmde();
2121  if(xb) {
2122  m_emc_EtaInCl[i] = QAH::H1F(name_eta, tit_eta, xb->GetSize()-1, xb->GetArray());
2123  delete xb;
2124  }
2125  m_emc_PhiInCl[i] = QAH::H1F(name_phi, tit_phi, grphi[i], -M_PI*1.015, M_PI*0.985);
2126  } else {
2127  m_emc_EtaInCl[i] = QAH::H1F(name_eta, tit_eta, greta[i], -1., 1.);
2128  m_emc_PhiInCl[i] = QAH::H1F(name_phi, tit_phi, grphi[i], -myPI, myPI);
2129  }
2130  }
2131  // Book the hists for Emc Points
2132  const TString catname[] = {"EmcCat1", "EmcCat2", "EmcCat3", "EmcCat4"};
2133 
2134  for (Int_t i=0; i<4; i++) {
2135  name = catname[i] + "_Point_Energy";
2136  title = catname[i] + " Point Energy";
2137  m_emc_point_energy[i]= QAH::H1F(name,title,100,0.,10.);
2138 
2139  name = catname[i] + "_Point_Eta";
2140  title = catname[i] + " Point Eta";
2141  m_emc_point_eta[i]= QAH::H1F(name,title,100,-1.,1.);
2142 
2143  name = catname[i] + "_Point_Phi";
2144  title = catname[i] + " Point Phi";
2145  m_emc_point_phi[i]= QAH::H1F(name,title,100,-3.14,3.14);
2146 
2147  name = catname[i] + "_Sigma_Eta";
2148  title = catname[i] + " Sigma Eta";
2149  m_emc_point_sigeta[i]= QAH::H1F(name,title,100,0.,.05);
2150 
2151  name = catname[i] + "_Sigma_Phi";
2152  title = catname[i] + " Sigma Phi";
2153  m_emc_point_sigphi[i]= QAH::H1F(name,title,100,0.,.05);
2154 
2155  name = catname[i] + "_Delta_Eta";
2156  title = catname[i] + " Delta Eta";
2157  m_emc_point_deleta[i]= QAH::H1F(name,title,100,-.2,.2);
2158 
2159  name = catname[i] + "_Delta_Phi";
2160  title = catname[i] + " Delta Phi";
2161  m_emc_point_delphi[i]= QAH::H1F(name,title,100,-.2,.2);
2162 
2163  name = catname[i] + "_Points_Multiplicity";
2164  title = catname[i] + " Points Multiplicity";
2165  m_emc_points[i]= QAH::H1F(name,title,200,0.,1000.);
2166 
2167  name = catname[i] + "_Track_Momenta";
2168  title = catname[i] + " Track Momenta ";
2169  m_emc_point_trmom[i]= QAH::H1F(name,title,100,0.,10.);
2170  }
2171 
2172  m_emc_point_flag= QAH::H1F(" Point Flag","Point Flag",5,0.5,5.5);
2173 
2174 }
2175 //_____________________________________________________________________________
2176 void StQABookHist::BookHistEval(){
2177 
2178 // these only get filled if the geant dataset is available!
2179 
2180  m_geant_reco_pvtx_x = QAH::H1F("QaGRpvtxDx"," diff geant - reco prim vtx X",
2181  100, -0.25,0.25);
2182  m_geant_reco_pvtx_y = QAH::H1F("QaGRpvtxDy"," diff geant - reco prim vtx Y",
2183  100, -0.25,0.25);
2184  m_geant_reco_pvtx_z = QAH::H1F("QaGRpvtxDz"," diff geant - reco prim vtx Z",
2185  100, -0.25,0.25);
2186  m_geant_reco_vtx_z_z = QAH::H2F("QaGRpvtxDzZ",
2187  " reco pvtx Z vs diff geant - reco Z", 100, -0.5,0.5,100,-50.,50.);
2188  m_geant_reco_vtx_z_z->SetXTitle("z vtx resolution (cm)");
2189  m_geant_reco_vtx_z_z->SetYTitle("z position of vtx (cm)");
2190 
2191 }
2192 //_____________________________________________________________________________
2193 void StQABookHist::BookHistBBC(){
2194 
2195  Char_t ID[4];
2196  Int_t i,j;
2197 
2198  m_bbc_adc[0] = QAH::MH1F("QaBbcAdcES","BBC East Small ADC",100,0.5,400.5,8);
2199  m_bbc_adc[1] = QAH::MH1F("QaBbcAdcEL","BBC East Large ADC",100,0.5,400.5,8);
2200  m_bbc_adc[2] = QAH::MH1F("QaBbcAdcWS","BBC West Small ADC",100,0.5,400.5,8);
2201  m_bbc_adc[3] = QAH::MH1F("QaBbcAdcWL","BBC West Large ADC",100,0.5,400.5,8);
2202  m_bbc_tdc[0] = QAH::MH1F("QaBbcTdcES","BBC East Small TDC",100,0.5,750.5,8);
2203  m_bbc_tdc[1] = QAH::MH1F("QaBbcTdcEL","BBC East Large TDC",100,0.5,2000.5,8);
2204  m_bbc_tdc[2] = QAH::MH1F("QaBbcTdcWS","BBC West Small TDC",100,0.5,750.5,8);
2205  m_bbc_tdc[3] = QAH::MH1F("QaBbcTdcWL","BBC West Large TDC",100,0.5,2000.5,8);
2206  for (i=0; i<8; i++) {
2207  sprintf(ID,"%d",i+1);
2208  for (j=0; j<4; j++) {
2209  m_bbc_adc[j]->Rebin(i,ID);
2210  m_bbc_tdc[j]->Rebin(i,ID);
2211  }
2212  }
2213 
2214 }
2215 //_____________________________________________________________________________
2216 void StQABookHist::BookHistFPD(){
2217 
2218  QAH::MMH1F(m_fpd_top,2,"QaFpdTop%d","FPD Top ADC %d-%d",100,0.5,1500.5,8,1);
2219  QAH::MMH1F(m_fpd_bottom,2,"QaFpdBottom%d","FPD Bottom ADC %d-%d",100,0.5,1500.5,8,1);
2220  QAH::MMH1F(m_fpd_south,2,"QaFpdSouth%d","FPD South ADC %d-%d",100,0.5,1500.5,8,1);
2221  QAH::MMH1F(m_fpd_north,2,"QaFpdNorth%d","FPD North ADC %d-%d",100,0.5,1500.5,6,1);
2222 
2223 // Book ADC histograms: FPD SUM signals
2224  m_fpd_sums[0] = QAH::H1F("QaFpdSums0","FPD SUM Top",100,0.5,2050.5);
2225  m_fpd_sums[1] = QAH::H1F("QaFpdSums1","FPD SUM Bottom",100,0.5,2050.5);
2226  m_fpd_sums[2] = QAH::H1F("QaFpdSums2","FPD SUM South",100,0.5,2050.5);
2227  m_fpd_sums[3] = QAH::H1F("QaFpdSums3","FPD SUM North",100,0.5,2050.5);
2228  m_fpd_sums[4] = QAH::H1F("QaFpdSums4","FPD SUM SmdX",100,0.5,6000.5);
2229  m_fpd_sums[5] = QAH::H1F("QaFpdSums5","FPD SUM SmdY",100,0.5,6000.5);
2230  m_fpd_sums[6] = QAH::H1F("QaFpdSums6","FPD SUM Pres1",100,0.5,1500.5);
2231  m_fpd_sums[7] = QAH::H1F("QaFpdSums7","FPD SUM Pres2",100,0.5,1500.5);
2232 
2233 }
2234 //_____________________________________________________________________________
2235 void StQABookHist::BookHistPMD(){
2236 
2237  QAH::MMH1F(m_pmd_sm_hit,12,"QaPmdSmHit%d",
2238  "PMD SM-wise Hit Multiplicity %02d-%02d",1000,0.,500.,2);
2239  QAH::MMH1F(m_pmd_sm_adc,12,"QaPmdSmAdcHit%d",
2240  "PMD SM-wise ADC/Hit Multiplicity %02d-%02d",100,0.,600.,2);
2241  QAH::MMH1F(m_pmd_chain_hit,24,"QaPmdChHitChain%d",
2242  "PMD Channel-wise Hit, Chain %02d-%02d",1728,-0.5,1727.5,2);
2243  QAH::MMH1F(m_pmd_chain_adc,24,"QaPmdChAdcChain%d",
2244  "PMD Channel-wise ADC, Chain %02d-%02d",1728,-0.5,1727.5,2);
2245 
2246  m_pmd_total_hit = QAH::H2F("QaPmdTotalHit","PMD Total Hits",100,0.,2e5,100,0.,4.);
2247  m_pmd_total_hit->SetXTitle("event id");
2248  m_pmd_total_hit->SetYTitle("log10");
2249  m_pmd_total_adc = QAH::H2F("QaPmdTotalAdc","PMD Total ADC",100,0.,2e5,100,0.,6.);
2250  m_pmd_total_adc->SetXTitle("event id");
2251  m_pmd_total_adc->SetYTitle("log10");
2252  m_cpv_total_hit = QAH::H2F("QaCpvTotalHit","CPV Total Hits",100,0.,2e5,100,0.,4.);
2253  m_cpv_total_hit->SetXTitle("event id");
2254  m_cpv_total_hit->SetYTitle("log10");
2255  m_cpv_total_adc = QAH::H2F("QaCpvTotalAdc","CPV Total ADC",100,0.,2e5,100,0.,6.);
2256  m_cpv_total_adc->SetXTitle("event id");
2257  m_cpv_total_adc->SetYTitle("log10");
2258 
2259 }
2260 //_____________________________________________________________________________
2261 void StQABookHist::BookHistTOF(){
2262 
2263  m_tof_hit_tray = QAH::H2F("QaTofHitvsTray","TOF Hits vs tray",120,0.5,120.5,100,0.,100.);
2264  m_tof_hit_module = QAH::H2F("QaTofHitvsModule","TOF Hits vs Module",65,-32.5,32.5,100,0.,100.);
2265  m_tof_match_tray = QAH::H2F("QaTofMatchvsTray","TOF Matched Hits vs tray",120,0.5,120.5,100,0.,100.);
2266  m_tof_match_module = QAH::H2F("QaTofMatchvsModule","TOF Matched Hits vs Module",65,-32.5,32.5,100,0.,100.);
2267  m_tof_vpd_hit = QAH::H2F("QaTofHitvsVpdHit","TOF Hits vs Vpd Hits",50,0.,50.,100,0.,5000.);
2268  m_tof_vtx_z = QAH::H2F("QaTofVpdZvsTpcZ","VPD vtxz vs TPC vtxz",100,-100.,100.,100,-100.,100.);
2269  m_tof_PID = QAH::H2F("QaTofPID","TOF InvBeta vs p",100,0.,5.,100,0.,4.);
2270 
2271 }
2272 //_____________________________________________________________________________
2273 void StQABookHist::BookHistHFT(){
2274 
2275  m_nhit_Pxl_Ist = QAH::H2F("QaPxlvsIstHit","PIXEL hits vs IST hits",500,0.,20000.,100,0.,5000.);
2276  m_nhit_Pxl_Ist->SetXTitle("PIXEL hits");
2277  m_nhit_Pxl_Ist->SetYTitle("IST hits");
2278 
2279  m_nhit_Pxl_Sst = QAH::H2F("QaPxlvsSstHit","PIXEL hits vs SST hits",500,0.,20000.,100,0.,5000.);
2280  m_nhit_Pxl_Sst->SetXTitle("PIXEL hits");
2281  m_nhit_Pxl_Sst->SetYTitle("SST hits");
2282 
2283  m_nhit_Ist_Sst = QAH::H2F("QaIstvsSsdHit","IST hits vs SST hits",100,0.,5000.,100,0.,5000.);
2284  m_nhit_Ist_Sst->SetXTitle("IST hits");
2285  m_nhit_Ist_Sst->SetYTitle("SST hits");
2286 
2287 }
2288 //_____________________________________________________________________________
2289 void StQABookHist::BookHistPXL(){
2290 
2291  m_pxl_hit_phi_z_Pxl1= QAH::H2F("QaPxlHitPhivsZPxl1","PIXEL: hits vs phi vs z in inner layer (per event)",20,-10.,10.,100,-TMath::Pi(),TMath::Pi());
2292  m_pxl_hit_phi_z_Pxl1->SetXTitle("z [cm]");
2293  m_pxl_hit_phi_z_Pxl1->SetYTitle("#phi");
2294  m_pxl_hit_phi_z_Pxl2= QAH::H2F("QaPxlHitPhivsZPxl2","PIXEL: hits vs phi vs z in outer layer (per event)",20,-10.,10.,100,-TMath::Pi(),TMath::Pi());
2295  m_pxl_hit_phi_z_Pxl2->SetXTitle("z [cm]");
2296  m_pxl_hit_phi_z_Pxl2->SetYTitle("#phi");
2297 
2298  m_pxl_hit_ladder = QAH::H1F("QaPxlHitvsLadder","PIXEL: hits per ladder (per event)",5,0.,5.);
2299  m_pxl_hit_ladder->SetXTitle("Ladder ID");
2300 
2301  m_pxl_hit_sector_sensor_Pxl1 = QAH::H2F("QaPxlHitvsSectorvsSensorPxl1","PIXEL: hits vs sector vs sensor in inner layer (per event)",10,0.5,10.5,10,0.5,10.5);
2302  m_pxl_hit_sector_sensor_Pxl1->SetXTitle("Sector ID");
2303  m_pxl_hit_sector_sensor_Pxl1->SetYTitle("Sensor ID");
2304  m_pxl_hit_sector_sensor_Pxl1->SetStats(kFALSE);
2305 
2306  m_pxl_hit_sector_sensor_Pxl2 = QAH::H2F("QaPxlHitvsSectorvsSensorPxl2","PIXEL: hits vs ladder vs sensor in outer layer (per event)",40,0.5,40.5,10,0.5,10.5);
2307  m_pxl_hit_sector_sensor_Pxl2->SetXTitle("Ladder ID");
2308  m_pxl_hit_sector_sensor_Pxl2->SetYTitle("Sensor ID");
2309  m_pxl_hit_sector_sensor_Pxl2->SetStats(kFALSE);
2310 
2311  m_pxl_nhit_Pxl1_tpc_mult = QAH::H2F("QaPxl1HitvsTpcMult","PIXEL: Hits in inner layer vs TPC multiplicity",100,0.,12000.,500,0.,10000.);
2312  m_pxl_nhit_Pxl1_tpc_mult->SetXTitle("PIXEL hits (inner layer)");
2313  m_pxl_nhit_Pxl1_tpc_mult->SetYTitle("TPC multiplicity");
2314 
2315  m_pxl_nhit_Pxl2_tpc_mult = QAH::H2F("QaPxl2HitvsTpcMult","PIXEL: Hits in outer layer vs TPC multiplicity",100,0.,8000.,500,0.,10000.);
2316  m_pxl_nhit_Pxl2_tpc_mult->SetXTitle("PIXEL hits (outer layer)");
2317  m_pxl_nhit_Pxl2_tpc_mult->SetYTitle("TPC multiplicity");
2318 
2319  m_pxl_nhit_Pxl1_tof_mult = QAH::H2F("QaPxl1HitvsTofMult","PIXEL: Hits in inner layer vs TOF multiplicity",100,0.,12000.,100,0.,5000.);
2320  m_pxl_nhit_Pxl1_tof_mult->SetXTitle("PIXEL hits (inner layer)");
2321  m_pxl_nhit_Pxl1_tof_mult->SetYTitle("ToF hits");
2322 
2323  m_pxl_nhit_Pxl2_tof_mult = QAH::H2F("QaPxl2HitvsTofMult","PIXEL: Hits in outer layer vs TOF multiplicity",100,0.,8000.,100,0.,5000.);
2324  m_pxl_nhit_Pxl2_tof_mult->SetXTitle("PIXEL hits (outer layer)");
2325  m_pxl_nhit_Pxl2_tof_mult->SetYTitle("ToF hits");
2326 
2327  m_pxl_nhit_Pxl1_Pxl2 = QAH::H2F("QaPxlHitLayer1vsLayer2","PIXEL: Hits in inner vs outer layer (per event)",100,0.,12000.,100,0.,8000.);
2328  m_pxl_nhit_Pxl1_Pxl2->SetXTitle("hits inner layer");
2329  m_pxl_nhit_Pxl1_Pxl2->SetYTitle("hits outer layer");
2330 
2331 }
2332 //_____________________________________________________________________________
2333 void StQABookHist::BookHistIST(){
2334 
2335 
2336  m_ist_hit_phi_z = QAH::H2F("QaIstHitPhivsZ","IST: Hits vs phi vs z (per event)",60,-30.,30.,100,-TMath::Pi(),TMath::Pi());
2337  m_ist_hit_phi_z->SetXTitle("z [cm]");
2338  m_ist_hit_phi_z->SetYTitle("#phi");
2339 
2340  m_ist_hit_ladder = QAH::H1F("QaIstHitvsLadder","IST: Hits per ladder (per event)",25,0.,25.);
2341  m_ist_hit_ladder->SetXTitle("Ladder ID");
2342 
2343  m_ist_hit_ladder_sensor = QAH::H2F("QaIstHitvsLaddervsSensor","IST: Hits vs ladder vs sensor (per event)", 24,0.5,24.5, 6, 0.5, 6.5);
2344  m_ist_hit_ladder_sensor->SetXTitle("Ladder ID");
2345  m_ist_hit_ladder_sensor->SetYTitle("Sensor ID");
2346  m_ist_hit_ladder_sensor->SetStats(kFALSE);
2347 
2348 
2349  m_ist_nhit_tpc_mult = QAH::H2F("QaIstHitvsTpcMult","IST hits vs TPC multiplicity",100,0.,5000.,500,0.,10000.);
2350  m_ist_nhit_tpc_mult->SetXTitle("IST hits");
2351  m_ist_nhit_tpc_mult->SetYTitle("TPC multiplicity");
2352 
2353  m_ist_nhit_tof_mult = QAH::H2F("QaIstHitvsTofMult","IST hits vs TOF multiplicity",100,0.,5000.,100,0.,5000.);
2354  m_ist_nhit_tof_mult->SetXTitle("IST hits");
2355  m_ist_nhit_tof_mult->SetYTitle("ToF hits");
2356 
2357 }
2358 //_____________________________________________________________________________
2359 void StQABookHist::BookHistMTD(){
2360 
2361  m_MtdNHits = QAH::H1F("QaMtdNHits","Number of MTD hits per event;N",50,0,50);
2362  m_MtdHitMap = QAH::H2F("QaMtdHitMap","MTD: channel vs backleg of hits;backleg;channel",30,0.5,30.5,60,-0.5,59.5);
2363  m_MtdNMatchHits = QAH::H1F("QaMtdNMatchHits","Number of matched MTD hits per event;N",20,0,20);
2364  m_MtdMatchHitMap = QAH::H2F("QaMtdMatchHitMap","MTD: channel vs backleg of matched hits;backleg;channel",30,0.5,30.5,60,-0.5,59.5);
2365 
2366 }
2367 //_____________________________________________________________________________
2368 void StQABookHist::BookHistEPD(){
2369 
2370  TString name, title;
2371  const TString side[2] = {"East", "West"};
2372 
2373  const int nPP = 12;
2374  for (int i=0; i<2; i++){
2375  for (int j=0; j<nPP; j++){
2376  (((name = "QaEpdAdc") += side[i]) += "PP") += (j+1);
2377  ((((title = "EPD ") += side[i]) += " PP") += (j+1)) += " ADC";
2378  m_epd_adc[(i*nPP)+j] = QAH::H2F(name,title,31,1,31,4096,0,4095);
2379  m_epd_adc[(i*nPP)+j]->SetXTitle("Tile");
2380  m_epd_adc[(i*nPP)+j]->SetYTitle("ADC");
2381 
2382  (((name = "QaEpdTac") += side[i]) += "PP") += (j+1);
2383  ((((title = "EPD ") += side[i]) += " PP") += (j+1)) += " TAC";
2384  m_epd_tac[(i*nPP)+j] = QAH::H2F(name,title,9,1,9,4096,0,4095);
2385  m_epd_tac[(i*nPP)+j]->SetXTitle("Tile");
2386  m_epd_tac[(i*nPP)+j]->SetYTitle("TAC");
2387  }
2388  }
2389 
2390 }
2391 //_____________________________________________________________________________
2392 void StQABookHist::BookHistiTPC(){
2393  char nameiTPC[60], titleiTPC[60];
2394  for(int i = 0; i < 24; i++)
2395  {
2396  sprintf(nameiTPC,"TPC_charge_time_inner_s%d", i+1);
2397  sprintf(titleiTPC,"TPC adc vs time-bin inner sector %d", i+1);
2398  m_TPC_ch_time_inner[i] = QAH::H2F(nameiTPC, titleiTPC, 50, 0, 511, 250 , 0, 500);
2399  m_TPC_ch_time_inner[i] ->SetXTitle("time bin");
2400  m_TPC_ch_time_inner[i] ->SetYTitle("adc");
2401  }
2402 
2403  for(int i = 0; i < 24; i++)
2404  {
2405  sprintf(nameiTPC,"TPC_charge_time_outer_s%d", i+1);
2406  sprintf(titleiTPC,"TPC adc vs time-bin outer sector %d", i+1);
2407  m_TPC_ch_time_outer[i] = QAH::H2F(nameiTPC, titleiTPC, 50, 0, 511, 250 , 0, 500);
2408  m_TPC_ch_time_outer[i] ->SetXTitle("time bin");
2409  m_TPC_ch_time_outer[i] ->SetYTitle("adc");
2410  }
2411 
2412  for(int i = 0; i < 24; i++)
2413  {
2414  sprintf(nameiTPC,"TPC_charge_nrow_s%d", i+1);
2415  sprintf(titleiTPC,"TPC adc vs pad-row sector %d", i+1);
2416  m_TPC_ch_nrow[i] = QAH::H2F(nameiTPC, titleiTPC, 72 ,0.5, 72.5, 250 , 0, 500);
2417  m_TPC_ch_nrow[i] ->SetXTitle("pad row");
2418  m_TPC_ch_nrow[i] ->SetYTitle("adc");
2419  }
2420 
2421  sprintf(nameiTPC,"TPC_adc_sec_inner");
2422  sprintf(titleiTPC,"TPC adc vs sector inner");
2423  m_TPC_adc_sec_inner = QAH::H2F(nameiTPC, titleiTPC, 24 ,0.5, 24.5, 250 , 0, 500);
2424  m_TPC_adc_sec_inner ->SetXTitle("sector");
2425  m_TPC_adc_sec_inner ->SetYTitle("adc");
2426 
2427  sprintf(nameiTPC,"TPC_adc_sec_outer");
2428  sprintf(titleiTPC,"TPC adc vs sector outer");
2429  m_TPC_adc_sec_outer = QAH::H2F(nameiTPC, titleiTPC, 24 ,0.5, 24.5, 250 , 0, 500);
2430  m_TPC_adc_sec_outer ->SetXTitle("sector");
2431  m_TPC_adc_sec_outer ->SetYTitle("adc");
2432 }
2433 //_____________________________________________________________________________
2434 
TH2F * m_glb_xfF
x-coord. of first hit on trk, tpc+svt
Definition: StQABookHist.h:277
TH2F * m_npoint_lengthTS
tot num points vs length, tpc
Definition: StQABookHist.h:385
TH1F * m_glb_impactrTS
log impact parameter from primary vertex, ftpc
Definition: StQABookHist.h:311
TH1F * m_lengthFE
length of track, ftpc
Definition: StQABookHist.h:346
TH1F * m_glb_yfFW
y-coord. of first hit on trk, ftpc east
Definition: StQABookHist.h:284
TH3F * m_pnt_svtLaserDiff
laser spots in svt
Definition: StQABookHist.h:714
TH1F * m_ndedxFW
dE/dx [0], ftpcE
Definition: StQABookHist.h:415
TH1F * m_plengthFE
length of track, ftpc
Definition: StQABookHist.h:527
TH2F * m_primtrk_meanptTTS
detector id of track
Definition: StQABookHist.h:429
TH1F * m_chisq0T
length of track, ftpc west
Definition: StQABookHist.h:348
TH2F * m_pT_eta_recT
overlayed hist of first point - helix point
Definition: StQABookHist.h:362
TH1F * m_primglob_good
tracks in table with iflag&gt;0, ftpc
Definition: StQABookHist.h:425
TH1F * m_glb_xf0TS
x-coord. of first hit - at start of helix
Definition: StQABookHist.h:292
TH2F * m_globtrk_padtimeFE
Y vs X of first hit on trk, ftpc west.
Definition: StQABookHist.h:371
TH1F * m_vtx_FtpcWestTpc_z
FtpcWest prim vtx - TPC prim vtx, x vs y.
Definition: StQABookHist.h:617
TH1F * m_pv_x
row1-vertex type
Definition: StQABookHist.h:590
TH2F * m_pmom_trklengthT
tanl(dip angle) vs zfirst-zvtx, tpc+svt
Definition: StQABookHist.h:552
TH1F * m_ndedxr
psi vs phi, tpc+svt
Definition: StQABookHist.h:402
TH2F * m_pchisq0_zfTS
chisq0 vs zfirst, tpc
Definition: StQABookHist.h:572
TH2F * m_tof_hit_module
of hits vs tray # (int over modules)
Definition: StQABookHist.h:769
TH1F * m_dcaToBeamZ1
xy-dca to beam axis
Definition: StQABookHist.h:232
TH1F * m_geant_reco_pvtx_z
prim vtx y, diff geant - reco
Definition: StQABookHist.h:649
TH1F * m_prim_zfT
y-coord. of first hit on trk, ftpc west
Definition: StQABookHist.h:477
TH1F * m_global_ist_hit
of hits in IST vs TOF track multiplicity
Definition: StQABookHist.h:763
TH2F * m_glb_impactF
log impact parameter from primary vertex, tpc,tpc+svt
Definition: StQABookHist.h:310
TH1F * m_kink_tot
xi Mass
Definition: StQABookHist.h:611
TH1F * m_vtx_FtpcEastTpc_z
FtpcEast prim vtx - TPC prim vtx, x vs y.
Definition: StQABookHist.h:615
TH1F * m_glb_chargeFW
particle charge in units of |e| - ftpc east
Definition: StQABookHist.h:264
TH1F * m_ppsiT
radial (xy) coordinate of first hit, ftpc west
Definition: StQABookHist.h:497
TH2F * m_prim_f0
chi square [1], ftpc west
Definition: StQABookHist.h:539
TH2F * m_primtrk_xf_yfTW
Y vs X of first hit on trk, tpc east.
Definition: StQABookHist.h:546
TH1F * m_glb_xf0
plane of first hit on trk, ftpc
Definition: StQABookHist.h:291
TH1F * m_petaTS
eta, tpc
Definition: StQABookHist.h:508
TH1F * m_dedx0FW
number of point to find dE/dx, ftpcW
Definition: StQABookHist.h:416
TH1F * m_etaTS
eta, tpc
Definition: StQABookHist.h:328
TH2F * m_pTTTS
pT, tpc+svt
Definition: StQABookHist.h:340
TH1F * m_pnt_tot_med
number of hits total
Definition: StQABookHist.h:637
TH1F * m_petaT
theta - tpc+svt
Definition: StQABookHist.h:507
TH1F * m_pchisq0T
length of track, ftpc west
Definition: StQABookHist.h:529
TH1F * m_pnt_tot_sm
number of hits total, med range
Definition: StQABookHist.h:638
TH2F * m_chisq0_zfTS
chisq0 vs zfirst, tpc
Definition: StQABookHist.h:395
TH2F * m_peta_trklengthT
mom vs. trk length, tpc+svt
Definition: StQABookHist.h:554
TH1F * m_ppointFE
number of points on the track - ftpc
Definition: StQABookHist.h:437
TH2F * m_glb_rzf0
z-coord. of first hit - at start of helix+svt
Definition: StQABookHist.h:297
TH2F * m_primtrk_meanetaTTS
&lt;pT&gt;, ftpc
Definition: StQABookHist.h:431
TH2F * m_petaF
eta, tpc,tpc+svt
Definition: StQABookHist.h:510
TH1F * m_prim_zfFE
z-coord. of first hit on trk, ftpc
Definition: StQABookHist.h:480
TH2F * m_pnt_timeT
z dist. of hits, svt
Definition: StQABookHist.h:623
TH1F * m_pTFW
pT, ftpc east
Definition: StQABookHist.h:343
TH2F * m_ist_nhit_tpc_mult
of hits vs ladder vs sensor
Definition: StQABookHist.h:761
TH1F * m_global_hft_hit
of hits in ITS vs hits in ITS
Definition: StQABookHist.h:740
TH2F * m_pfpoint_lengthTS
num fit points vs length, tpc
Definition: StQABookHist.h:559
TH1F * m_primary_hft_hit
HFT hits per global track.
Definition: StQABookHist.h:741
TH1F * m_glb_yfT
x-coord. of first hit on trk, ftpc west
Definition: StQABookHist.h:280
TH2F * m_pnt_xyFE
barrel dist. of hits, svt
Definition: StQABookHist.h:631
TH1F * m_prim_xfFW
x-coord. of first hit on trk, ftpc east
Definition: StQABookHist.h:471
TH2F * m_chisq0_phiT
chisq0 vs zfirst, tpc+svt
Definition: StQABookHist.h:396
TH2F * m_pmomF
momentum, tpc+svt
Definition: StQABookHist.h:515
TH2F * m_glb_impactT
chi square [1], ftpc west
Definition: StQABookHist.h:356
TH2F * m_chisq0_momTS
chisq0 vs momentum, tpc
Definition: StQABookHist.h:389
TH1F * m_prim_xf0TS
x-coord. of first hit - at start of helix
Definition: StQABookHist.h:483
TH2F * m_primtrk_xf_yfTE
pT versus eta, ftpcW
Definition: StQABookHist.h:545
TH1F * m_prim_radfFE
radial (xy) coordinate of first hit, ftpc
Definition: StQABookHist.h:495
TH1F * m_glb_radfT
transverse &amp; longitudinal helix center of last hit - helix hit+svt
Definition: StQABookHist.h:301
TH2F * m_pT_eta_recTS
pT versus eta, tpc
Definition: StQABookHist.h:363
TH1F * m_primtrk_good_sm
tracks in table with iflag&gt;0
Definition: StQABookHist.h:422
TH1F * m_prim_radfFW
radial (xy) coordinate of first hit, ftpc east
Definition: StQABookHist.h:496
TH1F * m_prim_xfFE
x-coord. of first hit on trk, ftpc
Definition: StQABookHist.h:470
TH2F * m_pT_eta_recFE
pT versus eta, tpc+svt
Definition: StQABookHist.h:364
TH1F * m_pTT
momentum, ftpc west
Definition: StQABookHist.h:338
TH1F * m_prim_radfTS
radial (xy) coordinate of first hit, tpc
Definition: StQABookHist.h:493
TH1F * m_v_pchi2
system
Definition: StQABookHist.h:586
TH2F * m_prim_rzf0TS
transverse &amp; longitudinal helix center of last hit - helix hit
Definition: StQABookHist.h:490
TH1F * m_xi_tot
K0 mass.
Definition: StQABookHist.h:607
TH2F * m_ppTF
pT, tpc,tpc+svt
Definition: StQABookHist.h:521
TH1F * m_global_pxl_hit
of hits in inner vs outer PIXEL layer
Definition: StQABookHist.h:754
TH1F * m_v_z
STAR reference.
Definition: StQABookHist.h:585
TH2F * m_glb_impactTTS
signed impact parameter from primary vertex, tpc+svt
Definition: StQABookHist.h:309
TH1F * m_ppointFW
number of points on the track - ftpc east
Definition: StQABookHist.h:438
TH2F * m_tof_match_module
of matched hits vs tray
Definition: StQABookHist.h:771
TH2F * m_npoint_lengthT
num fit points vs length, tpc,tpc+svt
Definition: StQABookHist.h:384
TH2F * m_tof_hit_tray
of IST hits per primary track
Definition: StQABookHist.h:768
TH1F * m_dcaToBeamZ2
z-dca to beam axis -210 to -105
Definition: StQABookHist.h:233
TH2F * m_ppT_eta_recTS
pT versus eta, tpc
Definition: StQABookHist.h:542
TH2F * m_chisq0_phiTS
chisq0 vs phi, tpc
Definition: StQABookHist.h:397
TH1F * m_prim_z0T
azimuthal angle at start (deg), tpc+svt
Definition: StQABookHist.h:463
TH1F * m_emc_point_sigeta[4]
Point Phi spectra.
Definition: StQABookHist.h:682
TH1F * m_chisq1FE
chi square [1], ftpc
Definition: StQABookHist.h:354
TH2F * m_ppsi_phiT
chisq0 vs zfirst, tpc+svt
Definition: StQABookHist.h:573
TH1F * m_glb_thetaT
tan(dip) =pz/pt at start, tpc+svt
Definition: StQABookHist.h:325
TH2F * m_prim_ratiomF
ratio of n fit pnts over max n pnts - tpc+svt
Definition: StQABookHist.h:451
TH1F * m_pchisq0FW
chi square [0], ftpc east
Definition: StQABookHist.h:534
TH1F * m_v0
z distribution of V0s relative to primVtx
Definition: StQABookHist.h:602
TH1F * m_pmax_pointTS
number of max possible track points - tpc
Definition: StQABookHist.h:440
TH1F * m_v_num_sm
number of vertices
Definition: StQABookHist.h:581
TH1F * m_prim_yf0TS
y-coord. of first hit - at start of helix
Definition: StQABookHist.h:485
TH1F * m_pTFE
pT, ftpc
Definition: StQABookHist.h:342
TH2F * m_pchisq0_dipTS
chisq0 vs dip angle, tpc
Definition: StQABookHist.h:570
TH2F * m_mom_trklengthTS
mom vs. trk length, tpc
Definition: StQABookHist.h:376
TH1F * m_chisq0FE
chi square [0], ftpc
Definition: StQABookHist.h:351
TH2F * m_ppT_eta_recFW
pT versus eta, ftpcE
Definition: StQABookHist.h:544
TH2F * m_chisq0F
chi square [0], tpc,svt
Definition: StQABookHist.h:350
TH2F * m_eta_trklengthFE
trk length vs. eta, tpc+svt
Definition: StQABookHist.h:379
TH1F * m_dedx1T
dE/dx [0], tpc
Definition: StQABookHist.h:406
TH1F * m_ppTFE
pT, ftpc
Definition: StQABookHist.h:522
TH2F * m_pxl_nhit_Pxl1_Pxl2
of hits in outer PIXEL layer vs TOF track multiplicity
Definition: StQABookHist.h:753
TH2F * m_globtrk_xf_yfTE
pT versus eta, ftpcW
Definition: StQABookHist.h:366
TH2F * m_pchisq1F
chi square [0], ftpc west
Definition: StQABookHist.h:535
TH1F * m_emc_point_delphi[4]
Point DeltaEta spectra.
Definition: StQABookHist.h:685
TH1F * m_ist_hit_ladder
of hits vs phi vs z position
Definition: StQABookHist.h:759
TH1F * m_psiTS
psi reconstructed, tpc
Definition: StQABookHist.h:318
TH1F * m_pnt_tpc
detector ID of the hit
Definition: StQABookHist.h:641
TH1F * m_vtx_phi_dist
radius to primary vertex
Definition: StQABookHist.h:597
TH2F * m_globtrk_padtimeFW
padlength vs timelength of hits on track, ftpcE
Definition: StQABookHist.h:372
TH2F * m_pnt_rpTE
rphi dist. of hits, tpcW
Definition: StQABookHist.h:626
TH1F * m_primtrk_tot_sm
tracks in table
Definition: StQABookHist.h:420
TH1F * m_v_y
vertex coordinates in
Definition: StQABookHist.h:584
TH2F * m_pfpoint_lengthTTS
num fit points vs length, tpc+svt
Definition: StQABookHist.h:560
TH2F * m_petaTTS
eta, tpc+svt
Definition: StQABookHist.h:509
TH2F * m_tof_vtx_z
of tof hits vs # of vpd hits
Definition: StQABookHist.h:773
TH2F * m_fit_pointTTS
number of track points used for fitting - tpc+svt
Definition: StQABookHist.h:252
TH2F * m_primtrk_meanptF
&lt;pT&gt;, tpc, tpc+svt
Definition: StQABookHist.h:430
TH2F * m_chisq0_dipT
chisq0 vs eta, tpc+svt
Definition: StQABookHist.h:392
TH1F * m_glb_xfFW
x-coord. of first hit on trk, ftpc east
Definition: StQABookHist.h:279
TH2F * m_ptanl_zfTS
tanl(dip angle) vs zfirst-zvtx, tpc
Definition: StQABookHist.h:551
TH2F * m_geant_reco_vtx_z_z
prim vtx z, diff geant - reco
Definition: StQABookHist.h:650
TH1F * m_etaFE
eta, ftpc
Definition: StQABookHist.h:331
TH1F * m_pchisq0TS
chi square [0], tpc
Definition: StQABookHist.h:530
TH1F * m_pv_y
row1-vertex coordinates in
Definition: StQABookHist.h:591
TH2F * m_chisq0_dipTS
chisq0 vs dip angle, tpc
Definition: StQABookHist.h:393
TH1F * m_prim_zf0
y-coord. of first hit - at start of helix+svt
Definition: StQABookHist.h:486
TH1F * m_pfit_pointTS
number of track points used for fitting - tpc
Definition: StQABookHist.h:445
TH1F * m_prim_radfT
transverse &amp; longitudinal helix center of last hit - helix hit+svt
Definition: StQABookHist.h:492
TH1F * m_glb_zfFW
z-coord. of first hit on trk, ftpc east
Definition: StQABookHist.h:289
TH2F * m_vtx_FtpcEastTpc_xy
number of kinks
Definition: StQABookHist.h:614
TH1F * m_glb_zf0TS
z-coord. of first hit - at start of helix
Definition: StQABookHist.h:296
TH1F * m_dcaToBeamZ3
z-dca to beam axis -105 to 0
Definition: StQABookHist.h:234
TH2F * m_fpoint_lengthT
trk length vs. eta, ftpc west
Definition: StQABookHist.h:381
TH1F * m_pv_pchi2
x versus y
Definition: StQABookHist.h:594
TH1F * m_glb_phi0TS
azimuthal angle at start (deg), tpc
Definition: StQABookHist.h:268
TH1F * m_pfit_pointT
number of max possible track points - ftpc west
Definition: StQABookHist.h:444
TH1F * m_max_pointFW
number of max possible track points - ftpc east
Definition: StQABookHist.h:249
TH2F * m_glb_radfF
impact parameter from primary vertex, ftpc
Definition: StQABookHist.h:314
TH1F * m_prim_chargeFW
particle charge in units of |e| - ftpc east
Definition: StQABookHist.h:458
TH1F * m_pointFW
number of points on the track - ftpc east
Definition: StQABookHist.h:244
TH1F * m_pnt_zS
z dist. of hits, tpc
Definition: StQABookHist.h:622
TH2F * m_pxl_nhit_Pxl1_tof_mult
of hits in outer PIXEL layer vs TPC track multiplicity
Definition: StQABookHist.h:751
TH2F * m_chisq0_momT
tot num points vs length, ftpc west
Definition: StQABookHist.h:388
TH1F * m_chisq1FW
chi square [1], ftpc east
Definition: StQABookHist.h:355
TH1F * m_glb_ratiomT
ratio of n fit pnts over max n pnts - tpc,tpc+svt
Definition: StQABookHist.h:255
TH2F * m_p_dedx_rec
psi vs phi, tpc+svt
Definition: StQABookHist.h:577
TH1F * m_glb_padfT
z-coord at start (cm), tpc+svt
Definition: StQABookHist.h:271
TH1F * m_ppTT
momentum, ftpc west
Definition: StQABookHist.h:518
TH2F * m_tof_PID
vertex z from vpd vs verex z from TPC
Definition: StQABookHist.h:774
TH2F * m_glb_simpactT
log impact parameter from primary vertex, tpc
Definition: StQABookHist.h:357
TH2F * m_pchisq0_dipT
chisq0 vs eta, tpc+svt
Definition: StQABookHist.h:569
TH1F * m_globtrk_fit_prob
Definition: StQABookHist.h:223
TH2F * m_ppTTTS
pT, tpc+svt
Definition: StQABookHist.h:520
TH2F * m_glb_padfTEW
padrow of first hit on trk, tpc
Definition: StQABookHist.h:272
TH2F * m_glb_impactTS
chi square [0], tpc+svt
Definition: StQABookHist.h:307
TH2F * m_prim_radfF
radial (xy) coordinate of first hit, tpc+svt
Definition: StQABookHist.h:494
TH1F * m_pnt_barrelS
padrow dist. of hits, tpc
Definition: StQABookHist.h:630
TH1F * m_ppsiFW
psi reconstructed, ftpc east
Definition: StQABookHist.h:502
TH2F * m_ist_hit_ladder_sensor
of hits vs x ladder (integrated over ladder, per event)
Definition: StQABookHist.h:760
TH1F * m_prim_yfFE
y-coord. of first hit on trk, ftpc
Definition: StQABookHist.h:475
TH2F * m_pnt_planeF
padlength vs timelength of hits, ftpcW
Definition: StQABookHist.h:635
TH1F * m_prim_curvTS
curvature (1/cm), tpc
Definition: StQABookHist.h:466
TH1F * m_prim_zfTS
z-coord. of first hit on trk, tpc
Definition: StQABookHist.h:478
TH1F * m_prim_zfFW
z-coord. of first hit on trk, ftpc east
Definition: StQABookHist.h:481
TH2F * m_glb_chargeF
particle charge in units of |e| - tpc+svt
Definition: StQABookHist.h:262
TH1F * m_ppointT
&lt;eta&gt;, ftpc
Definition: StQABookHist.h:434
TH1F * m_glb_xfT
curvature (1/cm), tpc+svt
Definition: StQABookHist.h:275
TH1F * m_prim_chargeFE
particle charge in units of |e| - ftpc
Definition: StQABookHist.h:457
TH1F * m_glb_phifTS
phi dist. of first point on trk, tpc
Definition: StQABookHist.h:304
TH2F * m_ppT_eta_recFE
pT versus eta, tpc+svt
Definition: StQABookHist.h:543
TH1F * m_petaFE
eta, ftpc
Definition: StQABookHist.h:511
TH2F * m_vtx_FtpcWestTpc_xy
FtpcEast prim vtx z - TPC prim vtx z.
Definition: StQABookHist.h:616
TH2F * m_glb_rzf0TS
transverse &amp; longitudinal helix center of last hit - helix hit
Definition: StQABookHist.h:299
TH2F * m_glb_impactrTTS
impact parameter from primary vertex, tpc+svt
Definition: StQABookHist.h:312
TH2F * m_glb_zfF
z-coord. of first hit on trk, tpc+svt
Definition: StQABookHist.h:287
TH1F * m_pointFE
number of points on the track - ftpc
Definition: StQABookHist.h:243
TH2F * m_eta_trklengthTS
trk length vs. eta, tpc
Definition: StQABookHist.h:378
TH1F * m_globtrk_goodTTS
tracks in table with iflag&gt;0,small range
Definition: StQABookHist.h:227
TH1F * m_psiFE
psi reconstructed, ftpc
Definition: StQABookHist.h:321
TH1F * m_glb_phi0T
radius at start (cm), tpc+svt
Definition: StQABookHist.h:267
TH2F * m_eta_trklengthT
mom vs. trk length, tpc+svt
Definition: StQABookHist.h:377
TH1F * m_ppTTS
pT, tpc
Definition: StQABookHist.h:519
TH1F * m_glb_xfFE
x-coord. of first hit on trk, ftpc
Definition: StQABookHist.h:278
TH2F * m_eta_trklengthFW
trk length vs. eta, ftpc east
Definition: StQABookHist.h:380
TH2F * m_pnt_padtimeFW
padlength vs timelength of hits, ftpcE
Definition: StQABookHist.h:634
TH1F * m_prim_ratiomT
ratio of n fit pnts over tot n pnts - ftpc west
Definition: StQABookHist.h:449
TH1F * m_prim_xf0
z-coord. of first hit on trk, ftpc west
Definition: StQABookHist.h:482
TH1F * m_pnt_ssd
SVT drift.
Definition: StQABookHist.h:718
TH2F * m_pchisq0F
chi square [0], tpc,tpc+svt
Definition: StQABookHist.h:532
TH2F * m_pointF
number of points on the track - tpc+svt
Definition: StQABookHist.h:242
TH1F * m_pmax_pointFE
number of max possible track points - ftpc
Definition: StQABookHist.h:442
TH1F * m_v_x
vertex type
Definition: StQABookHist.h:583
TH2F * m_glb_f0
impact parameter from primary vertex, tpc
Definition: StQABookHist.h:360
TH2F * m_globtrk_xf_yfFE
Y vs X of first hit on trk, tpc+svt.
Definition: StQABookHist.h:369
TH1F * m_glb_xfTS
x-coord. of first hit on trk, tpc
Definition: StQABookHist.h:276
TH1F * m_prim_zf0TS
z-coord. of first hit - at start of helix
Definition: StQABookHist.h:487
TH1F * m_dedx0T
number of point to find dE/dx, tpc
Definition: StQABookHist.h:405
TH1F * m_pmomTS
momentum, tpc
Definition: StQABookHist.h:514
TH2F * m_primtrk_goodF
tracks in table with iflag&gt;0, tpc,svt
Definition: StQABookHist.h:424
TH2F * m_ptanl_zfT
Y vs X of first hit on trk, ftpc west.
Definition: StQABookHist.h:550
TH1F * m_pnt_phiS
phi dist. of hits, tpc
Definition: StQABookHist.h:628
TH1F * m_glb_radfFW
radial (xy) coordinate of first hit, ftpc east
Definition: StQABookHist.h:316
TH1F * m_prim_xfTS
x-coord. of first hit on trk, tpc
Definition: StQABookHist.h:468
TH1F * m_glb_yf0
x-coord. of first hit - at start of helix+svt
Definition: StQABookHist.h:293
TH1F * m_glb_impactrT
signed impact parameter from primary vertex, tpc
Definition: StQABookHist.h:358
TH1F * m_ndedxFE
dE/dx [0], ftpc
Definition: StQABookHist.h:412
TH2F * m_prim_f0TS
overlayed hist of first point - helix point
Definition: StQABookHist.h:540
TH1F * m_momT
eta, ftpc west
Definition: StQABookHist.h:333
TH1F * m_glb_chargeTS
particle charge in units of |e| - tpc
Definition: StQABookHist.h:261
TH1F * m_ev0_k0ma_hist
Lambda mass.
Definition: StQABookHist.h:604
TH2F * m_glb_rzl0
transverse &amp; longitudinal helix center of first hit - helix hit
Definition: StQABookHist.h:298
TH2F * m_peta_trklengthTS
trk length vs. eta, tpc
Definition: StQABookHist.h:555
TH2F * m_pxl_nhit_Pxl2_tof_mult
of hits in inner PIXEL layer vs TOF track multiplicity
Definition: StQABookHist.h:752
TH1F * m_chisq0TS
length of track, tpc+svt
Definition: StQABookHist.h:306
TH2F * m_primtrk_xf_yfFW
Y vs X of first hit on trk, ftpc east.
Definition: StQABookHist.h:549
TH1F * m_globtrk_iflag
tracks in table with iflag&gt;0, ftpc
Definition: StQABookHist.h:229
TH1F * m_pmomT
eta, ftpc west
Definition: StQABookHist.h:513
TH1F * m_v_num
dedx vs p
Definition: StQABookHist.h:580
TH2F * m_zDcaZf
z-dca vs tanl
Definition: StQABookHist.h:236
TH1F * m_petaFW
eta, ftpc east
Definition: StQABookHist.h:512
TH1F * m_glb_radfTS
radial (xy) coordinate of first hit, tpc
Definition: StQABookHist.h:302
TH1F * m_emc_points[4]
Point TrMom spectra.
Definition: StQABookHist.h:687
TH2F * m_ppsiTTS
psi reconstructed, tpc+svt
Definition: StQABookHist.h:499
TH1F * m_primtrk_goodTTS
tracks in table with iflag&gt;0, small range
Definition: StQABookHist.h:423
TH1F * m_vtx_z_dist
radial distribution of V0s relative to primVtx
Definition: StQABookHist.h:599
TH1F * m_pmax_pointT
number of points on the track - ftpc west
Definition: StQABookHist.h:439
TH2F * m_tof_match_tray
of hits vs module # (int over trays, east side use -32 - -1)
Definition: StQABookHist.h:770
TH1F * m_tanlTS
tan(dip) =pz/pt at start, tpc
Definition: StQABookHist.h:324
TH2F * m_chisq0_etaTS
chisq0 vs eta, tpc
Definition: StQABookHist.h:391
TH2F * m_fpoint_lengthTS
num fit points vs length, tpc
Definition: StQABookHist.h:382
TH2F * m_prim_yfF
y-coord. of first hit on trk, tpc+svt
Definition: StQABookHist.h:474
TH2F * m_pnt_ftpc
number of hits tpc
Definition: StQABookHist.h:642
TH1F * m_prim_thetaT
tan(dip) =pz/pt at start, tpc+svt
Definition: StQABookHist.h:505
TH1F * m_max_pointFE
number of max possible track points - ftpc
Definition: StQABookHist.h:248
TH2F * m_glb_f0TS
overlayed hist of first point - helix point
Definition: StQABookHist.h:361
TH1F * m_glb_zf0
y-coord. of first hit - at start of helix+svt
Definition: StQABookHist.h:295
TH2F * m_pv_xy
system
Definition: StQABookHist.h:593
TH2F * m_glb_impactrF
impact parameter from primary vertex, tpc,tpc+svt
Definition: StQABookHist.h:313
TH1F * m_det_id
iflag value
Definition: StQABookHist.h:230
TH1F * m_pnt_phiSSD
number of hits ssd
Definition: StQABookHist.h:719
TH2F * m_pnt_xyFW
xy dist. of hits, ftpcE
Definition: StQABookHist.h:632
TH1F * m_pv_z
STAR reference.
Definition: StQABookHist.h:592
TH2F * m_prim_rzl0TS
transverse &amp; longitudinal helix center of first hit - helix hit+svt
Definition: StQABookHist.h:491
TH2F * m_dcaToBeamXY
detector id of track
Definition: StQABookHist.h:231
TH1F * m_pv_vtxid
radius to vertex
Definition: StQABookHist.h:589
TH1F * m_pointTS
number of points on the track - tpc
Definition: StQABookHist.h:241
TH2F * m_psiTTS
psi reconstructed, tpc+svt
Definition: StQABookHist.h:319
TH1F * m_glb_yfTS
y-coord. of first hit on trk, tpc
Definition: StQABookHist.h:281
TH2F * m_prim_chargeF
particle charge in units of |e| - tpc+svt
Definition: StQABookHist.h:456
TH2F * m_zDcaTanl
z-dca to beam axis 0 to 105
Definition: StQABookHist.h:235
TH1F * m_etaFW
eta, ftpc east
Definition: StQABookHist.h:332
TH2F * m_glb_phifT
radial (xy) coordinate of first hit, tpc+svt
Definition: StQABookHist.h:303
TH2F * m_chisq0_zfT
chisq0 vs dip angle, tpc+svt
Definition: StQABookHist.h:394
TH2F * m_pT_eta_recFW
pT versus eta, ftpcE
Definition: StQABookHist.h:365
TH1F * m_pmomFE
momentum, ftpc
Definition: StQABookHist.h:516
TH2F * m_prim_rzl0
transverse &amp; longitudinal helix center of first hit - helix hit
Definition: StQABookHist.h:489
TH1F * m_momTS
momentum, tpc
Definition: StQABookHist.h:334
TH1F * m_pchisq1FE
chi square [1], ftpc
Definition: StQABookHist.h:536
TH1F * m_glb_ratiomFE
ratio of n fit pnts over max n pnts - ftpc
Definition: StQABookHist.h:258
TH1F * m_glb_curvTS
curvature (1/cm), tpc
Definition: StQABookHist.h:274
TH1F * m_plengthFW
length of track, ftpc east
Definition: StQABookHist.h:528
TH1F * m_momFW
momentum, ftpc east
Definition: StQABookHist.h:337
TH1F * m_ev0_lama_hist
v0 vertices
Definition: StQABookHist.h:603
TH2F * m_pchisq0_etaTS
chisq0 vs eta, tpc
Definition: StQABookHist.h:568
TH2F * m_primtrk_xf_yfTS
Y vs X of first hit on trk, tpc west.
Definition: StQABookHist.h:547
TH2F * m_pxl_hit_phi_z_Pxl1
HFT hits per primary track.
Definition: StQABookHist.h:744
TH1F * m_plengthTS
length of track, tpc
Definition: StQABookHist.h:525
TH2F * m_pxl_hit_sector_sensor_Pxl1
of hits vs x ladder (integrated over ladder, per event)
Definition: StQABookHist.h:747
TH1F * m_lengthFW
length of track, ftpc east
Definition: StQABookHist.h:347
TH1F * m_primary_pxl_hit
of PIXEL hits per global track
Definition: StQABookHist.h:755
TH1F * m_glb_curvT
padrow of first hit on trk, tpc east and west
Definition: StQABookHist.h:273
TH1F * m_pnt_tot
plane dist. of hits, ftpc
Definition: StQABookHist.h:636
TH1F * m_prim_chargeT
ratio of n fit pnts over max n pnts - ftpc west
Definition: StQABookHist.h:454
TH2F * m_globtrk_goodF
tracks in table with iflag&gt;0, tpc, tpc+svt
Definition: StQABookHist.h:228
TH1F * m_primtrk_good
tracks in table, small range
Definition: StQABookHist.h:421
TH2F * m_pnpoint_lengthTS
tot num points vs length, tpc
Definition: StQABookHist.h:562
TH2F * m_nhit_Pxl_Sst
of hits in PIXEL vs hits in IST
Definition: StQABookHist.h:738
TH1F * m_glb_ratiomFW
ratio of n fit pnts over max n pnts - ftpc east
Definition: StQABookHist.h:259
TH1F * m_prim_yfTS
y-coord. of first hit on trk, tpc
Definition: StQABookHist.h:473
TH2F * m_glb_yfF
y-coord. of first hit on trk, tpc+svt
Definition: StQABookHist.h:282
TH2F * m_pnpoint_lengthT
num fit points vs length, tpc
Definition: StQABookHist.h:561
TH2F * m_svt_loc
diff of laser spots in svt
Definition: StQABookHist.h:715
TH2F * m_pnpoint_lengthFE
tot num points vs length, tpc+svt
Definition: StQABookHist.h:563
TH1F * m_pdet_id
iflag value
Definition: StQABookHist.h:428
TH1F * m_pTTS
pT, tpc
Definition: StQABookHist.h:339
TH2F * m_glb_rzl0TS
transverse &amp; longitudinal helix center of first hit - helix hit+svt
Definition: StQABookHist.h:300
TH1F * m_ppointTS
number of points on the track - tpc
Definition: StQABookHist.h:435
TH1F * m_glb_r0TS
radius at start (cm), tpc
Definition: StQABookHist.h:266
TH1F * m_pchisq1FW
chi square [1], ftpc east
Definition: StQABookHist.h:537
TH1F * m_primtrk_tot
dE/dx [0], ftpcW
Definition: StQABookHist.h:419
TH2F * m_nhit_Ist_Sst
of hits in PIXEL vs hits in SST
Definition: StQABookHist.h:739
TH1F * m_prim_r0TS
radius at start (cm), tpc
Definition: StQABookHist.h:460
TH2F * m_pchisq0_etaT
chisq0 vs momentum, tpc+svt
Definition: StQABookHist.h:567
TH1F * m_glb_yfFE
y-coord. of first hit on trk, ftpc
Definition: StQABookHist.h:283
TH2F * m_ppsi_phiTS
psi vs phi, tpc
Definition: StQABookHist.h:574
TH1F * m_emc_point_trmom[4]
Point DeltaPhi spectra.
Definition: StQABookHist.h:686
TH1F * m_pmomFW
momentum, ftpc east
Definition: StQABookHist.h:517
TH1F * m_pointT
z-dca vs phi0
Definition: StQABookHist.h:240
TH1F * m_ptanlTS
tan(dip) =pz/pt at start, tpc
Definition: StQABookHist.h:504
TH1F * m_glb_zfT
y-coord. of first hit on trk, ftpc west
Definition: StQABookHist.h:285
TH2F * m_prim_zfF
z-coord. of first hit on trk, tpc+svt
Definition: StQABookHist.h:479
TH2F * m_pnt_xyS
time dist. of hits, tpc
Definition: StQABookHist.h:624
TH2F * m_glb_planefF
z-coord. of first hit on trk, ftpc west
Definition: StQABookHist.h:290
TH2F * m_pxl_nhit_Pxl2_tpc_mult
of hits in inner PIXEL layer vs TPC track multiplicity
Definition: StQABookHist.h:750
TH2F * m_lengthF
length of track, tpc
Definition: StQABookHist.h:345
TH2F * m_ppsiF
psi reconstructed, tpc,tpc+svt
Definition: StQABookHist.h:500
TH1F * m_glb_thetaTS
theta - tpc
Definition: StQABookHist.h:326
TH2F * m_ppT_eta_recT
overlayed hist of first point - helix point
Definition: StQABookHist.h:541
TH2F * m_chisq0_etaT
chisq0 vs momentum, tpc+svt
Definition: StQABookHist.h:390
TH2F * m_etaF
eta, tpc,tpc+svt
Definition: StQABookHist.h:330
TH2F * m_etaTTS
eta, tpc+svt
Definition: StQABookHist.h:329
TH2F * m_chisq0TTS
chi square [0], tpc
Definition: StQABookHist.h:349
TH1F * m_xi_ma_hist
number of xi vertices
Definition: StQABookHist.h:608
TH1F * m_prim_yfFW
y-coord. of first hit on trk, ftpc east
Definition: StQABookHist.h:476
TH1F * m_ppTFW
pT, ftpc east
Definition: StQABookHist.h:523
TH2F * m_plengthF
length of track, tpc+svt
Definition: StQABookHist.h:526
TH2F * m_pmax_pointF
number of max possible track points - tpc+svt
Definition: StQABookHist.h:441
TH2F * m_tanl_zfTS
tanl(dip angle) vs zfirst-zvtx, tpc
Definition: StQABookHist.h:374
TH2F * m_pchisq0_momT
tot num points vs length, ftpc west
Definition: StQABookHist.h:565
TH2F * m_pmom_trklengthTS
mom vs. trk length, tpc
Definition: StQABookHist.h:553
TH1F * m_ptanlT
psi reconstructed, ftpc west
Definition: StQABookHist.h:503
TH2F * m_psi_phiT
chisq0 vs phi, tpc+svt
Definition: StQABookHist.h:398
TH1F * m_emc_point_sigphi[4]
Point SigmaEta spectra.
Definition: StQABookHist.h:683
TH2F * m_dedx0F
number of point to find dE/dx, ftpc
Definition: StQABookHist.h:410
TH2F * m_pchisq0_zfT
chisq0 vs dip angle, tpc+svt
Definition: StQABookHist.h:571
TH1F * m_lengthTS
phi dist. of first point on trk, tpc+svt
Definition: StQABookHist.h:305
TH1F * m_prim_phi0TS
azimuthal angle at start (deg), tpc
Definition: StQABookHist.h:462
TH1F * m_pnt_ftpcE
number of hits ftpc
Definition: StQABookHist.h:643
TH1F * m_psiT
radial (xy) coordinate of first hit, ftpc west
Definition: StQABookHist.h:317
TH1F * m_primtrk_iflag
fit points ratio prim/glob, all detectors
Definition: StQABookHist.h:427
TH2F * m_max_pointF
number of max possible track points - tpc+svt
Definition: StQABookHist.h:247
TH1F * m_prim_yf0
x-coord. of first hit - at start of helix+svt
Definition: StQABookHist.h:484
TH1F * m_prim_curvT
z-coord at start (cm), tpc+svt
Definition: StQABookHist.h:465
TH2F * m_glb_ratiomF
ratio of n fit pnts over max n pnts - tpc+svt
Definition: StQABookHist.h:257
TH1F * m_ndedxT
number of tracks with dedx info
Definition: StQABookHist.h:404
TH2F * m_primtrk_meanetaF
&lt;eta&gt;, tpc, tpc+svt
Definition: StQABookHist.h:432
TH1F * m_globtrk_good_sm
tracks in table with iflag&gt;0
Definition: StQABookHist.h:226
TH2F * m_prim_rzf0
z-coord. of first hit - at start of helix+svt
Definition: StQABookHist.h:488
TH2F * m_pnt_svtLaser
number of hits svt
Definition: StQABookHist.h:713
TH2F * m_tof_vpd_hit
of matched hits vs module
Definition: StQABookHist.h:772
TH1F * m_lengthT
pT, ftpc west
Definition: StQABookHist.h:344
TH1F * m_dedxTTS
dE/dx [1], tpc
Definition: StQABookHist.h:407
TH2F * m_fpoint_lengthTTS
num fit points vs length, tpc+svt
Definition: StQABookHist.h:383
TH1F * m_glb_chargeT
ratio of n fit pnts over max n pnts - ftpc west
Definition: StQABookHist.h:260
TH1F * m_prim_ratiomTS
ratio of n fit pnts over max n pnts - tpc
Definition: StQABookHist.h:450
TH1F * m_plengthT
pT, ftpc west
Definition: StQABookHist.h:524
TH1F * m_prim_ratiomFE
ratio of n fit pnts over max n pnts - ftpc
Definition: StQABookHist.h:452
TH1F * m_dedx0FE
number of point to find dE/dx, ftpcE
Definition: StQABookHist.h:413
TH2F * m_tanl_zfT
padlength vs timelength of hits on track, ftpcW
Definition: StQABookHist.h:373
TH2F * m_ist_hit_phi_z
of PIXEL hits per primary track
Definition: StQABookHist.h:758
TH2F * m_emc_nhit
prim vtx z, diff geant - reco vs reco z
Definition: StQABookHist.h:655
TH1F * m_glb_zfFE
z-coord. of first hit on trk, ftpc
Definition: StQABookHist.h:288
TH2F * m_bbc_adc[4]
Point Flag spectra.
Definition: StQABookHist.h:691
TH1F * m_emc_point_deleta[4]
Point SigmaPhi spectra.
Definition: StQABookHist.h:684
TH1F * m_pmax_pointFW
number of max possible track points - ftpc east
Definition: StQABookHist.h:443
TH1F * m_glb_z0T
azimuthal angle at start (deg), tpc+svt
Definition: StQABookHist.h:269
TH2F * m_glb_ratiomTTS
number of track pnts on svt - tpc+svt
Definition: StQABookHist.h:254
TH2F * m_pxl_hit_sector_sensor_Pxl2
of hits vs sector vs sensor, inner PIXEL layer
Definition: StQABookHist.h:748
TH1F * m_glb_zfTS
z-coord. of first hit on trk, tpc
Definition: StQABookHist.h:286
TH1F * m_glb_r0T
particle charge in units of |e| - ftpc west
Definition: StQABookHist.h:265
TH1F * m_z_hits
FtpcWest prim vtx z - TPC prim vtx z.
Definition: StQABookHist.h:621
TH1F * m_primglob_fit
trks w/ iflag&gt;0, ratio prim/glob, all detectors
Definition: StQABookHist.h:426
TH2F * m_globtrk_xf_yfTW
Y vs X of first hit on trk, tpc east.
Definition: StQABookHist.h:367
TH1F * m_prim_r0T
particle charge in units of |e| - ftpc west
Definition: StQABookHist.h:459
TH1F * m_glb_chargeFE
particle charge in units of |e| - ftpc
Definition: StQABookHist.h:263
TH2F * m_pnt_phiT
rphi dist. of hits, tpcE
Definition: StQABookHist.h:627
TH1F * m_emc_point_phi[4]
Point Eta spectra.
Definition: StQABookHist.h:681
TH1F * m_pnt_id
number of hits total, small range
Definition: StQABookHist.h:639
TH1F * m_max_pointT
number of points on the track - ftpc west
Definition: StQABookHist.h:245
TH1F * m_MtdNHits
TOF PID: 1/beta vs p.
Definition: StQABookHist.h:776
TH1F * m_tanlT
psi reconstructed, ftpc west
Definition: StQABookHist.h:323
TH2F * m_mom_trklengthT
tanl(dip angle) vs zfirst-zvtx, tpc+svt
Definition: StQABookHist.h:375
TH2F * m_npoint_lengthFE
tot num points vs length, tpc+svt
Definition: StQABookHist.h:386
TH2F * m_pnt_padrowT
phi dist. of hits, svt
Definition: StQABookHist.h:629
TH1F * m_geant_reco_pvtx_y
prim vtx x, diff geant - reco
Definition: StQABookHist.h:648
TH1F * m_pchisq0FE
chi square [0], ftpc
Definition: StQABookHist.h:533
TH1F * m_momFE
momentum, ftpc
Definition: StQABookHist.h:336
TH1F * m_fit_pointTS
number of track points used for fitting - tpc
Definition: StQABookHist.h:251
TH2F * m_peta_trklengthFE
trk length vs. eta, tpc+svt
Definition: StQABookHist.h:556
TH1F * m_v_vtxid
number of vertices,small range
Definition: StQABookHist.h:582
TH1F * m_glb_radfFE
radial (xy) coordinate of first hit, ftpc
Definition: StQABookHist.h:315
TH2F * m_npoint_lengthFW
tot num points vs length, ftpc east
Definition: StQABookHist.h:387
TH1F * m_fit_pointT
number of max possible track points - ftpc west
Definition: StQABookHist.h:250
TH1F * m_pnt_phiSST
number of hits sst
Definition: StQABookHist.h:728
TH2F * m_globtrk_xf_yfTS
Y vs X of first hit on trk, tpc west.
Definition: StQABookHist.h:368
TH2F * m_pTF
pT, tpc,tpc+svt
Definition: StQABookHist.h:341
TH2F * m_prim_xfF
x-coord. of first hit on trk, tpc+svt
Definition: StQABookHist.h:469
TH2F * m_ppointF
number of points on the track - tpc+svt
Definition: StQABookHist.h:436
TH1F * m_ppsiFE
psi reconstructed, ftpc
Definition: StQABookHist.h:501
TH2F * m_ist_nhit_tof_mult
of hits in IST vs TPC track multiplicity
Definition: StQABookHist.h:762
TH1F * m_ppsiTS
psi reconstructed, tpc
Definition: StQABookHist.h:498
TH1F * m_max_pointTS
number of max possible track points - tpc
Definition: StQABookHist.h:246
TH1F * m_prim_chargeTS
particle charge in units of |e| - tpc
Definition: StQABookHist.h:455
TH2F * m_zDcaPhi0
z-dca vs psi
Definition: StQABookHist.h:238
TH2F * m_pchisq0_momTS
chisq0 vs momentum, tpc
Definition: StQABookHist.h:566
TH1F * m_prim_thetaTS
theta - tpc
Definition: StQABookHist.h:506
TH2F * m_psiF
psi reconstructed, tpc,tpc+svt
Definition: StQABookHist.h:320
TH1F * m_etaT
theta - tpc+svt
Definition: StQABookHist.h:327
TH1F * m_prim_yfT
x-coord. of first hit on trk, ftpc west
Definition: StQABookHist.h:472
TH2F * m_glb_simpactTS
log impact parameter from primary vertex, tpc+svt
Definition: StQABookHist.h:308
TH2F * m_pfpoint_lengthT
trk length vs. eta, ftpc west
Definition: StQABookHist.h:558
TH1F * m_pnt_ftpcW
number of hits ftpcE
Definition: StQABookHist.h:644
TH2F * m_psi_phiTS
psi vs phi, tpc
Definition: StQABookHist.h:399
TH2F * m_zDcaPsi
z-dca vs z-first
Definition: StQABookHist.h:237
TH2F * m_pnt_padtimeFE
xy dist. of hits, ftpcW
Definition: StQABookHist.h:633
TH1F * m_glb_z0TS
z-coord at start (cm), tpc
Definition: StQABookHist.h:270
TH1F * m_emc_point_eta[4]
Point Energy spectra.
Definition: StQABookHist.h:680
TH2F * m_peta_trklengthFW
trk length vs. eta, ftpc east
Definition: StQABookHist.h:557
TH1F * m_vtx_r_dist
azimuthal distribution of V0s relative to primVtx
Definition: StQABookHist.h:598
TH1F * m_pv_r
row1-chisq per dof of vertex fit
Definition: StQABookHist.h:595
TH1F * m_prim_z0TS
z-coord at start (cm), tpc
Definition: StQABookHist.h:464
TH2F * m_pnpoint_lengthFW
tot num points vs length, ftpc east
Definition: StQABookHist.h:564
TH2F * m_globtrk_xf_yfFW
Y vs X of first hit on trk, ftpc east.
Definition: StQABookHist.h:370
TH1F * m_prim_ratiomFW
ratio of n fit pnts over max n pnts - ftpc east
Definition: StQABookHist.h:453
TH2F * m_pchisq0TTS
chi square [0], tpc+svt
Definition: StQABookHist.h:531
TH1F * m_chisq0FW
chi square [0], ftpc east
Definition: StQABookHist.h:352
TH1F * m_primary_ist_hit
of IST hits per global track
Definition: StQABookHist.h:764
TH1F * m_globtrk_good
tracks in table
Definition: StQABookHist.h:225
TH1F * m_pxl_hit_ladder
of hits vs phi vs z position, outer PIXEL layer
Definition: StQABookHist.h:746
TH1F * m_psiFW
psi reconstructed, ftpc east
Definition: StQABookHist.h:322
TH2F * m_pxl_hit_phi_z_Pxl2
of hits vs phi vs z position, inner PIXEL layer
Definition: StQABookHist.h:745
TH2F * m_ndedxF
&lt;dE/dx&gt;/(Bichsel &lt;dE/dx&gt;), tpc,tpc+svt
Definition: StQABookHist.h:409
TH2F * m_pnt_rpTW
xy dist. of hits, svt
Definition: StQABookHist.h:625
TH1F * m_emc_point_flag
Emc Point multiplicity.
Definition: StQABookHist.h:688
TH1F * m_glb_ratiomTS
ratio of n fit pnts over max n pnts - tpc
Definition: StQABookHist.h:256
TH2F * m_momF
momentum, tpc+svt
Definition: StQABookHist.h:335
TH1F * m_geant_reco_pvtx_x
number of hits ftpcW
Definition: StQABookHist.h:647
TH1F * m_glb_yf0TS
y-coord. of first hit - at start of helix
Definition: StQABookHist.h:294
TH1F * m_glb_sptsTS
number of track points used for fitting - tpc,tpc+svt
Definition: StQABookHist.h:253
TH2F * m_pxl_nhit_Pxl1_tpc_mult
of hits vs sector vs sensor, outer PIXEL layer
Definition: StQABookHist.h:749
TH1F * m_prim_phi0T
radius at start (cm), tpc+svt
Definition: StQABookHist.h:461
TH1F * m_prim_xfT
curvature (1/cm), tpc+svt
Definition: StQABookHist.h:467
TH2F * m_chisq1F
chi square [0], ftpc west
Definition: StQABookHist.h:353
TH2F * m_primtrk_xf_yfFE
Y vs X of first hit on trk, tpc+svt.
Definition: StQABookHist.h:548
TH1F * m_v_r
chisq per dof of vertex fit
Definition: StQABookHist.h:587