StRoot: St_QA_Maker/StQABookHist.cxx Source File

00001 // $Id: StQABookHist.cxx,v 2.68 2012/03/05 03:42:32 genevb Exp $
00002 // $Log: StQABookHist.cxx,v $
00003 // Revision 2.68  2012/03/05 03:42:32  genevb
00004 // Remove TPC XY dist, add TPC RPhi charge
00005 //
00006 // Revision 2.67  2012/02/09 03:01:05  genevb
00007 // No FTPC histograms for Run 12+
00008 //
00009 // Revision 2.66  2012/02/08 22:10:35  genevb
00010 // Updates for Run 12
00011 //
00012 // Revision 2.65  2011/05/31 21:35:49  genevb
00013 // TPC request: add time bucket distribution of hits
00014 //
00015 // Revision 2.64  2009/03/19 01:08:08  genevb
00016 // Show both xy and rphi TPC hit hists
00017 //
00018 // Revision 2.63  2009/03/13 19:27:24  genevb
00019 // Now draw TPC xy hits in polar coords
00020 //
00021 // Revision 2.62  2009/03/09 16:34:54  genevb
00022 // Tweak a couple histogram ranges
00023 //
00024 // Revision 2.61  2007/12/10 19:58:20  genevb
00025 // Use log10 for number of secondary vertices
00026 //
00027 // Revision 2.60  2007/12/03 16:51:09  genevb
00028 // Out-of-order statements for pre-run8 data
00029 //
00030 // Revision 2.59  2007/11/30 05:38:50  genevb
00031 // Changes for Run8: mostly silicon removal, TOF addition
00032 //
00033 // Revision 2.58  2007/04/25 18:35:57  genevb
00034 // Additional SSD hists
00035 //
00036 // Revision 2.57  2007/04/24 00:33:58  genevb
00037 // SSD hists
00038 //
00039 // Revision 2.56  2007/04/12 20:39:48  genevb
00040 // Cleanup (removal) of CalibVtx, Nfitpnt, Chisq1, Rich, histograms
00041 //
00042 // Revision 2.55  2007/04/07 04:40:30  genevb
00043 // Remove fit pnts/tot; retitle log as log10
00044 //
00045 // Revision 2.54  2007/04/02 05:36:23  genevb
00046 // PMD always in now
00047 //
00048 // Revision 2.53  2007/03/13 18:46:28  genevb
00049 // Added Svt Laser Diff
00050 //
00051 // Revision 2.52  2007/02/26 20:45:01  genevb
00052 // SVT drift hist
00053 //
00054 // Revision 2.51  2006/05/18 16:37:11  genevb
00055 // Change FTPC hist ranges
00056 //
00057 // Revision 2.50  2005/02/22 19:38:39  genevb
00058 // Do PMD hists only for year 4 and later (real data)
00059 //
00060 // Revision 2.49  2005/02/08 17:22:46  genevb
00061 // PMD histo changes, handle estGlobal/ITTF tracks
00062 //
00063 // Revision 2.48  2005/01/27 05:28:25  genevb
00064 // PMD changes
00065 //
00066 // Revision 2.47  2004/12/13 15:52:36  genevb
00067 // Numerous updates: PMD, primtrk, FPD, QAShift lists
00068 //
00069 // Revision 2.46  2004/10/04 16:40:42  genevb
00070 // FTPC radial histos
00071 //
00072 // Revision 2.45  2004/04/23 23:15:29  genevb
00073 // Added signedDCA (Impact) plots for globals
00074 //
00075 // Revision 2.44  2004/02/12 05:03:10  genevb
00076 // Year 4 AuAu changes. New SVT histos.
00077 //
00078 // Revision 2.43  2004/01/27 20:08:14  genevb
00079 // Modify SVT point hist
00080 //
00081 // Revision 2.42  2004/01/10 01:10:18  genevb
00082 // Preparations for Year 5, added some svt plots
00083 //
00084 // Revision 2.41  2003/09/02 17:59:21  perev
00085 // gcc 3.2 updates + WarnOff
00086 //
00087 // Revision 2.40  2003/04/19 00:17:49  genevb
00088 // Updated for dAu/pp running
00089 //
00090 // Revision 2.39  2003/02/20 20:09:54  genevb
00091 // Several changes for new trigger scheme, dAu data
00092 //
00093 // Revision 2.38  2003/02/19 06:38:29  genevb
00094 // Rework trigger and mult/event class sections
00095 //
00096 // Revision 2.37  2003/02/15 21:59:08  genevb
00097 // dedx lable
00098 //
00099 // Revision 2.36  2003/01/24 15:24:48  genevb
00100 // Further dAu changes
00101 //
00102 // Revision 2.35  2003/01/23 20:53:11  genevb
00103 // Additional dAu changes
00104 //
00105 // Revision 2.34  2003/01/23 04:03:18  jeromel
00106 // Include fixed
00107 //
00108 // Revision 2.33  2003/01/13 19:09:43  genevb
00109 // Minor FTPC histo changes
00110 //
00111 // Revision 2.32  2002/07/11 12:26:22  genevb
00112 // Some FTPC histogram changes
00113 //
00114 // Revision 2.31  2002/05/29 13:54:30  genevb
00115 // Some changes to FTPC chisq histos
00116 //
00117 // Revision 2.30  2002/04/23 01:59:56  genevb
00118 // Addition of BBC/FPD histos
00119 //
00120 // Revision 2.29  2002/02/12 18:41:59  genevb
00121 // Additional FTPC histograms
00122 //
00123 // Revision 2.28  2002/01/21 22:09:24  genevb
00124 // Include some ftpc histograms from StFtpcClusterMaker
00125 //
00126 // Revision 2.27  2001/12/28 09:19:13  genevb
00127 // Adjustments for pp running
00128 //
00129 // Revision 2.26  2001/11/20 21:53:45  lansdell
00130 // added x-y dist of hits, tpc east&west histos
00131 //
00132 // Revision 2.25  2001/11/13 21:17:56  lansdell
00133 // changed limits on dedx error on mean, tpc
00134 //
00135 // Revision 2.24  2001/11/02 20:50:03  genevb
00136 // Changed histogram ranges for momenta
00137 //
00138 // Revision 2.23  2001/10/25 20:45:12  lansdell
00139 // changed ftpc # hits range to 25k
00140 //
00141 // Revision 2.22  2001/10/25 20:43:43  lansdell
00142 // changed 2d ftpc globtrk histo ranges
00143 //
00144 // Revision 2.21  2001/10/10 05:00:19  genevb
00145 // Changed ftpc point/track ranges
00146 //
00147 // Revision 2.20  2001/08/07 07:51:27  lansdell
00148 // primvtx check for different multiplicities crashed for MC data, now fixed
00149 //
00150 // Revision 2.19  2001/08/03 23:31:59  lansdell
00151 // added missing axis labels to 2d histos
00152 //
00153 // Revision 2.18  2001/08/03 20:33:55  lansdell
00154 // added primvtx check histos for different multiplicities; separated x-y plot of first point on track, tpc into east and west histos
00155 //
00156 // Revision 2.17  2001/07/31 23:21:42  lansdell
00157 // added last point, hit-helix histos
00158 //
00159 // Revision 2.16  2001/06/06 18:44:59  lansdell
00160 // changed domain on some qa_shift histograms
00161 //
00162 // Revision 2.15  2001/05/29 23:23:06  lansdell
00163 // removed impact param plots for FTPC from qa_shift list
00164 //
00165 // Revision 2.14  2001/05/25 17:20:09  lansdell
00166 // squashed a small bug
00167 //
00168 // Revision 2.13  2001/05/25 16:31:21  lansdell
00169 // more updates to qa shift histograms
00170 //
00171 // Revision 2.12  2001/05/24 01:48:13  lansdell
00172 // qa_shift histograms updated
00173 //
00174 // Revision 2.11  2001/05/23 00:14:52  lansdell
00175 // more changes for qa_shift histograms
00176 //
00177 // Revision 2.10  2001/05/16 20:57:03  lansdell
00178 // new histograms added for qa_shift printlist; some histogram ranges changed; StMcEvent now used in StEventQA
00179 //
00180 // Revision 2.9  2001/05/02 16:10:47  lansdell
00181 // changed some histogram limits
00182 //
00183 // Revision 2.8  2001/05/01 15:17:36  genevb
00184 // Execute EMC code only if EMC libs loaded
00185 //
00186 // Revision 2.7  2001/04/28 22:05:13  genevb
00187 // Added EMC histograms
00188 //
00189 // Revision 2.6  2001/04/26 20:45:19  lansdell
00190 // changed some histogram ranges (TPC+SVT radius at start, impact param for primary tracks)
00191 //
00192 // Revision 2.5  2001/04/26 16:34:28  genevb
00193 // Fixed some histogram ranges
00194 //
00195 // Revision 2.4  2001/04/25 21:35:26  genevb
00196 // Added V0 phi distributions
00197 //
00198 // Revision 2.3  2001/04/24 22:53:51  lansdell
00199 // Removed redundant radial position of first hit histograms
00200 //
00201 // Revision 2.2  2000/09/08 18:55:53  lansdell
00202 // turned on FTPC primary track histograms
00203 //
00204 // Revision 2.1  2000/09/06 15:31:13  lansdell
00205 // changed FTPC radial position of first hit histogram limits
00206 //
00207 // Revision 2.0  2000/08/25 16:02:40  genevb
00208 // New revision: new structure, multiplicity classes
00209 //
00210 // 
00212 //                                                                       //
00213 //  StQABookHist class for multiplicity-based QA histograms              //
00214 //                                                                       //
00216 
00217 #include <Stiostream.h>
00218 #include <math.h>
00219 #include "QAH.h"
00220 #include "StQABookHist.h"
00221 #include "StQAMakerBase.h"
00222 #include "StEmcUtil/others/StEmcMath.h"
00223 #include "StEmcUtil/others/emcDetectorName.h"
00224 #include "TROOT.h"
00225 
00226 
00227 ClassImp(StQABookHist)
00228   
00229 //_____________________________________________________________________________
00230 StQABookHist::StQABookHist(const char* type) : QAHistType(type) {
00231 
00232   silHists = kFALSE;
00233   ftpHists = kFALSE;
00234 
00235   int i = 0;
00236 
00237 //  - Set all the histogram booking constants
00238 
00239   nxpT = 50;
00240   nyeta = 50;
00241   xminpT = 0.0;
00242   xmaxpT = 5.0;
00243 //  ymineta = -2.0;
00244 //  ymaxeta =  2.0;
00245 
00246   nchisq = 50;
00247 //  nmass  = 40;
00248 //  ntau   = 40;
00249   ndedx  = 50;
00250   npnt   = 50;
00251   nleng  = 50;
00252   npsi   = 36;
00253   knpsi  = 42;
00254   nvrt   = 100;
00255   knyeta = 60;
00256   knid   = 10;
00257   cnp   = 50;
00258   cndedx = 50;    
00259 
00260   kminnid  = 0.0;
00261   kmaxnid  = 10.0;
00262   minpsi   = 0.0;
00263   kminpsi  = -60.0;
00264   maxpsi   = 360.0;
00265   minchisq = 0.;
00266   maxchisq = 10.0;
00267   minmass  = 0.0;
00268   maxmass  = 2.0;
00269   minpnt   = 0.0;
00270   maxpnt   = 50.0;
00271   minleng  = 0.0;
00272   maxleng  = 200.0;
00273   mintau   = 0.0;
00274   maxtau   = 20.0;
00275   mintrk   = 0.0;
00276   maxtrk   = 8000.0;
00277   minvrt   = 2000.0;
00278   maxvrt   = 4000.0;
00279   minmpt   = 0.0;
00280   maxmpt   = 2.0;
00281   minmeta  = -0.2;
00282   maxmeta  = 0.2;
00283   kmineta  = -3.0;
00284   kmaxeta  = 3.0;
00285   minxyz   = 0.0;
00286   maxxyz   = 50.0;
00287   cminp = 0.0;
00288   cmaxp = 2.0;
00289   cmindedx = 0.0;
00290   cmaxdedx =  0.1e-04*1e6; // change from GeV to keV per cm
00291 
00292 
00293 //  - Zero all histogram pointers
00294 
00295   mNullPrimVtxClass = 0;
00296 
00297 // for method MakeGlob - from table globtrk
00298   m_globtrk_tot=0;
00299   m_globtrk_good=0;
00300   m_globtrk_good_sm=0;
00301   m_globtrk_goodTTS=0;
00302   m_globtrk_goodF=0;
00303   m_globtrk_iflag=0;
00304   m_det_id=0;
00305   m_dcaToBeamXY=0;
00306   m_dcaToBeamZ1=0;
00307   m_dcaToBeamZ2=0;
00308   m_dcaToBeamZ3=0;
00309   m_zDcaTanl=0;
00310   m_zDcaZf=0;
00311   m_zDcaPsi=0;
00312   m_zDcaPhi0=0;
00313 
00314   m_pointT=0;
00315   m_pointF=0;
00316   m_pointFE=0;
00317   m_pointFW=0;
00318   m_max_pointT=0;
00319   m_max_pointF=0;
00320   m_max_pointFE=0;
00321   m_max_pointFW=0;
00322   m_fit_pointT=0;
00323   m_fit_pointTTS=0;
00324   m_glb_sptsTS=0;
00325   m_glb_ratiomTTS=0;
00326   m_glb_ratiomT=0;
00327   m_glb_ratiomF=0;
00328   m_glb_ratiomFE=0;
00329   m_glb_ratiomFW=0;
00330   m_glb_chargeT=0;
00331   m_glb_chargeF=0;
00332   m_glb_chargeFE=0;
00333   m_glb_chargeFW=0;
00334   m_glb_r0T=0;
00335   m_glb_phi0T=0;
00336   m_glb_z0T=0;
00337   m_glb_curvT=0;
00338   m_glb_padfT=0;  
00339   m_glb_padfTEW=0;  
00340   m_glb_xf0=0;
00341   m_glb_xfT=0;
00342   m_glb_xfF=0;
00343   m_glb_xfFE=0;
00344   m_glb_xfFW=0;
00345   m_glb_yf0=0;
00346   m_glb_yfT=0;  
00347   m_glb_yfF=0;
00348   m_glb_yfFE=0;  
00349   m_glb_yfFW=0;     
00350   m_glb_zf0=0;     
00351   m_glb_zfT=0; 
00352   m_glb_zfF=0;
00353   m_glb_zfFE=0; 
00354   m_glb_zfFW=0;
00355   m_glb_planefF=0;
00356   m_glb_f0=0;
00357   m_glb_rzf0=0;
00358   m_glb_rzl0=0;
00359   m_glb_radfT=0;
00360   m_glb_radfF=0;
00361   m_glb_radfFE=0;
00362   m_glb_radfFW=0;
00363   m_glb_phifT=0;
00364   m_glb_phifTS=0;
00365   m_psiT=0;  
00366   m_psiTTS=0;  
00367   m_psiF=0;
00368   m_psiFE=0;
00369   m_psiFW=0;
00370   m_tanlT=0;
00371   m_glb_thetaT=0;  
00372   m_etaT=0;    
00373   m_etaTTS=0;    
00374   m_etaF=0;
00375   m_etaFE=0;
00376   m_etaFW=0;
00377   m_pTT=0;
00378   m_pTTTS=0;
00379   m_pTF=0;
00380   m_pTFE=0;
00381   m_pTFW=0;
00382   m_momT=0; 
00383   m_momF=0;
00384   m_momFE=0; 
00385   m_momFW=0;        
00386   m_lengthT=0;  
00387   m_lengthF=0;
00388   m_lengthFE=0;  
00389   m_lengthFW=0;  
00390   m_chisq0T=0;     
00391   m_chisq0TTS=0;     
00392   m_chisq0F=0;
00393   m_chisq0FE=0;
00394   m_chisq0FW=0;
00395   m_chisq1F=0;
00396   m_chisq1FE=0;
00397   m_chisq1FW=0;
00398   m_glb_impactT=0; 
00399   m_glb_simpactT=0; 
00400   m_glb_impactrT=0; 
00401   m_glb_impactTTS=0; 
00402   m_glb_impactrTTS=0; 
00403   m_glb_impactF=0; 
00404   m_glb_impactrF=0; 
00405 
00406   m_pointTS=0;        
00407   m_max_pointTS=0;    
00408   m_fit_pointTS=0;    
00409   m_glb_ratiomTS=0;   
00410   m_glb_chargeTS=0;   
00411   m_glb_r0TS=0;       
00412   m_glb_phi0TS=0;     
00413   m_glb_z0TS=0;       
00414   m_glb_curvTS=0;     
00415   m_glb_xfTS=0;       
00416   m_glb_yfTS=0;       
00417   m_glb_zfTS=0;       
00418   m_glb_xf0TS=0;      
00419   m_glb_yf0TS=0;      
00420   m_glb_zf0TS=0;      
00421   m_glb_f0TS=0;      
00422   m_glb_rzf0TS=0;
00423   m_glb_rzl0TS=0;
00424   m_glb_radfTS=0;     
00425   m_psiTS=0;          
00426   m_tanlTS=0;         
00427   m_glb_thetaTS=0;    
00428   m_etaTS=0;          
00429   m_momTS=0;          
00430   m_pTTS=0;           
00431   m_lengthTS=0;       
00432   m_chisq0TS=0;       
00433   m_glb_impactTS=0;   
00434   m_glb_simpactTS=0;   
00435   m_glb_impactrTS=0;   
00436 
00437   m_pT_eta_recT = 0;
00438   m_pT_eta_recFE = 0;
00439   m_pT_eta_recFW = 0;
00440   m_globtrk_xf_yfTE = 0;
00441   m_globtrk_xf_yfTW = 0;
00442   m_globtrk_xf_yfFE = 0;
00443   m_globtrk_xf_yfFW = 0;
00444   m_globtrk_padtimeFE = 0;
00445   m_globtrk_padtimeFW = 0;
00446   m_tanl_zfT  = 0;
00447   m_mom_trklengthT = 0;
00448   m_eta_trklengthT = 0;
00449   m_eta_trklengthFE = 0;
00450   m_eta_trklengthFW = 0;
00451   m_npoint_lengthT = 0; 
00452   m_npoint_lengthFE = 0;        
00453   m_npoint_lengthFW = 0;                  
00454   m_fpoint_lengthT = 0;
00455   m_fpoint_lengthTTS = 0;
00456   m_chisq0_momT = 0;
00457   m_chisq0_etaT = 0;
00458   m_chisq0_dipT = 0;
00459   m_chisq0_zfT = 0;
00460   m_chisq0_phiT = 0;
00461 
00462   m_pT_eta_recTS= 0;
00463   m_globtrk_xf_yfTS = 0;
00464   m_tanl_zfTS  = 0;
00465   m_mom_trklengthTS = 0;
00466   m_eta_trklengthTS = 0;
00467   m_npoint_lengthTS = 0;        
00468   m_fpoint_lengthTS = 0;
00469   m_chisq0_momTS = 0;
00470   m_chisq0_etaTS = 0;
00471   m_chisq0_dipTS = 0;
00472   m_chisq0_zfTS = 0;
00473   m_chisq0_phiTS = 0;
00474 
00475 // for method MakeDE - from table dst_dedx
00476   m_ndedxr=0;        
00477 
00478   m_ndedxT=0;         
00479   m_dedx0T=0;         
00480   m_dedx1T=0;         
00481   m_dedxTTS=0;        
00482   // east and west on same plot
00483   m_ndedxF=0;         
00484   m_dedx0F=0;         
00485   // east and west on separate plots
00486   m_ndedxFE=0;         
00487   m_dedx0FE=0;         
00488   
00489   m_ndedxFW=0;         
00490   m_dedx0FW=0;         
00491   
00492 // for method MakeHistPrim - from table primtrk
00493   m_primtrk_tot=0;
00494   m_primtrk_tot_sm=0;
00495   m_primtrk_good=0;
00496   m_primtrk_good_sm=0;
00497   m_primtrk_goodTTS=0;
00498   m_primtrk_goodF=0;
00499   m_primtrk_iflag=0;
00500   m_primglob_good=0;
00501   m_pdet_id=0;
00502   m_primtrk_meanptTTS=0;
00503   m_primtrk_meanptF=0;
00504   m_primtrk_meanetaTTS=0;
00505   m_primtrk_meanetaF=0;
00506 
00507   m_ppointT=0;
00508   m_ppointF=0;
00509   m_ppointFE=0;
00510   m_ppointFW=0;
00511   m_pmax_pointT=0;
00512   m_pmax_pointF=0;
00513   m_pmax_pointFE=0;
00514   m_pmax_pointFW=0;
00515   m_pfit_pointT=0;
00516   m_prim_ratiomT=0;
00517   m_prim_ratiomF=0;
00518   m_prim_ratiomFE=0;
00519   m_prim_ratiomFW=0;
00520   m_prim_chargeT=0;
00521   m_prim_chargeF=0;
00522   m_prim_chargeFE=0;
00523   m_prim_chargeFW=0;
00524   m_prim_r0T=0;
00525   m_prim_phi0T=0;
00526   m_prim_z0T=0;
00527   m_prim_curvT=0;
00528   m_prim_xf0=0;
00529   m_prim_xfT=0;
00530   m_prim_xfF=0;
00531   m_prim_xfFE=0;
00532   m_prim_xfFW=0;
00533   m_prim_yf0=0;
00534   m_prim_yfT=0;  
00535   m_prim_yfF=0;
00536   m_prim_yfFE=0;
00537   m_prim_yfFW=0;
00538   m_prim_zf0=0;     
00539   m_prim_zfT=0; 
00540   m_prim_zfF=0;
00541   m_prim_zfFE=0;
00542   m_prim_zfFW=0;
00543   m_prim_f0=0;
00544   m_prim_rzf0=0;
00545   m_prim_rzl0=0;
00546   m_prim_radfT=0;
00547   m_prim_radfF=0;
00548   m_prim_radfFE=0;
00549   m_prim_radfFW=0;
00550   m_ppsiT=0;  
00551   m_ppsiTTS=0;  
00552   m_ppsiF=0;
00553   m_ppsiFE=0;
00554   m_ppsiFW=0;
00555   m_ptanlT=0;   
00556   m_prim_thetaT=0;  
00557   m_petaT=0;    
00558   m_petaTTS=0;    
00559   m_petaF=0;
00560   m_petaFE=0;
00561   m_petaFW=0;
00562   m_ppTT=0;
00563   m_ppTTTS=0;
00564   m_ppTF=0;
00565   m_ppTFE=0;
00566   m_ppTFW=0;
00567   m_pmomT=0; 
00568   m_pmomF=0;
00569   m_pmomFE=0;
00570   m_pmomFW=0;
00571   m_plengthT=0;  
00572   m_plengthF=0;
00573   m_plengthFE=0;
00574   m_plengthFW=0;
00575   m_pchisq0T=0;     
00576   m_pchisq0TTS=0;     
00577   m_pchisq0F=0;
00578   m_pchisq0FE=0;
00579   m_pchisq0FW=0;
00580   m_pchisq1F=0;
00581   m_pchisq1FE=0;
00582   m_pchisq1FW=0;
00583 
00584   m_ppointTS=0;        
00585   m_pmax_pointTS=0;    
00586   m_pfit_pointTS=0;    
00587   m_prim_ratiomTS=0;   
00588   m_prim_chargeTS=0;   
00589   m_prim_r0TS=0;       
00590   m_prim_phi0TS=0;     
00591   m_prim_z0TS=0;       
00592   m_prim_curvTS=0;     
00593   m_prim_xfTS=0;       
00594   m_prim_yfTS=0;       
00595   m_prim_zfTS=0;       
00596   m_prim_xf0TS=0;      
00597   m_prim_yf0TS=0;
00598   m_prim_zf0TS=0;
00599   m_prim_f0TS=0;
00600   m_prim_rzf0TS=0;
00601   m_prim_rzl0TS=0;
00602   m_prim_radfTS=0;
00603   m_ppsiTS=0;          
00604   m_ptanlTS=0;         
00605   m_prim_thetaTS=0;    
00606   m_petaTS=0;          
00607   m_momTS=0;          
00608   m_ppTTS=0;           
00609   m_plengthTS=0;       
00610   m_pchisq0TS=0;       
00611 
00612   m_ppT_eta_recT = 0;
00613   m_ppT_eta_recFE = 0;
00614   m_ppT_eta_recFW = 0;
00615   m_primtrk_xf_yfTE = 0;
00616   m_primtrk_xf_yfTW = 0;
00617   m_primtrk_xf_yfFE = 0;
00618   m_primtrk_xf_yfFW = 0;
00619   m_ptanl_zfT  = 0;
00620   m_pmom_trklengthT = 0;
00621   m_peta_trklengthT = 0;
00622   m_peta_trklengthFE = 0;
00623   m_peta_trklengthFW = 0;
00624   m_pnpoint_lengthT = 0;        
00625   m_pnpoint_lengthFE = 0;       
00626   m_pnpoint_lengthFW = 0;                 
00627   m_pfpoint_lengthT = 0;
00628   m_pfpoint_lengthTTS = 0;
00629   m_pchisq0_momT = 0;
00630   m_pchisq0_etaT = 0;
00631   m_pchisq0_dipT = 0;
00632   m_pchisq0_zfT = 0;
00633 
00634   m_ppT_eta_recTS= 0;
00635   m_primtrk_xf_yfTS = 0;
00636   m_ptanl_zfTS  = 0;
00637   m_pmom_trklengthTS = 0;
00638   m_peta_trklengthTS = 0;
00639   m_pnpoint_lengthTS = 0;       
00640   m_pfpoint_lengthTS = 0;
00641   m_pchisq0_momTS = 0;
00642   m_pchisq0_etaTS = 0;
00643   m_pchisq0_dipTS = 0;
00644   m_pchisq0_zfTS = 0;
00645 
00646   // for MakeHistPID - from tables primtrk & dst_dedx 
00647   m_p_dedx_rec=0;   
00648   
00649   // for method MakeHistVertex - from table dst_vertex
00650   m_v_num=0;   
00651   m_v_num_sm=0;
00652   m_v_vtxid=0; 
00653   m_v_x=0;     
00654   m_v_y=0;     
00655   m_v_z=0;     
00656   m_v_pchi2=0; 
00657   m_v_r=0;     
00658   
00659   m_pv_vtxid=0; 
00660   m_pv_x=0;     
00661   m_pv_y=0;     
00662   m_pv_z=0;     
00663   m_pv_xy=0;    
00664   m_pv_pchi2=0; 
00665   m_pv_r=0;     
00666 
00667   m_vtx_phi_dist=0;  
00668   m_vtx_r_dist=0;  
00669   m_vtx_z_dist=0;  
00670 
00671   m_v0           =0; 
00672   m_ev0_lama_hist=0; 
00673   m_ev0_k0ma_hist=0; 
00674   m_xi_tot=0;        
00675   m_xi_ma_hist=0;    
00676   m_kink_tot=0;      
00677   
00678   m_vtx_FtpcEastTpc_xy=0; 
00679   m_vtx_FtpcEastTpc_z=0;  
00680   m_vtx_FtpcWestTpc_xy=0; 
00681   m_vtx_FtpcWestTpc_z=0;  
00682   
00683 
00684 // for method MakeHistPoint
00685   m_z_hits=0;
00686   m_pnt_tot=0;     
00687   m_pnt_tot_med=0; 
00688   m_pnt_tot_sm=0;  
00689   m_pnt_id=0;      
00690   m_pnt_xyS=0;     
00691   m_pnt_rpTW=0;    
00692   m_pnt_rpTE=0;    
00693   m_pnt_phiT=0;    
00694   m_pnt_padrowT=0; 
00695   m_pnt_timeT=0;   
00696   m_pnt_zS=0;      
00697   m_pnt_phiS=0;    
00698   m_pnt_barrelS=0; 
00699   m_pnt_planeF=0;  
00700   m_pnt_tpc=0;     
00701   m_pnt_svt=0;     
00702   m_pnt_ssd=0;     
00703   m_pnt_xyFE=0;    
00704   m_pnt_xyFW=0;    
00705   m_pnt_padtimeFE=0;    
00706   m_pnt_padtimeFW=0;    
00707   // east and west on same plot
00708   m_pnt_ftpc=0;    
00709   // east and west on separate plots
00710   m_pnt_ftpcE=0;   
00711   m_pnt_ftpcW=0;   
00712   m_pnt_svtLaser=0;     
00713   m_pnt_svtLaserDiff=0; 
00714 
00715 // for method MakeHistEval  
00716   m_geant_reco_pvtx_x=0;  
00717   m_geant_reco_pvtx_y=0;  
00718   m_geant_reco_pvtx_z=0;  
00719   m_geant_reco_vtx_z_z=0; 
00720 
00721 // for BBC
00722   for (i=0; i<4; i++) {
00723     m_bbc_adc[i] = 0;
00724     m_bbc_tdc[i] = 0;
00725   }
00726 
00727 // for FPD
00728   for (i=0; i<2; i++) {
00729     m_fpd_top[i] = 0;
00730     m_fpd_bottom[i] = 0;
00731     m_fpd_south[i] = 0;
00732     m_fpd_north[i] = 0;
00733   }
00734   for (i=0; i<8; i++) m_fpd_sums[i] = 0;
00735 
00736 // for PMD
00737   for (i=0; i<12; i++) {
00738     m_pmd_sm_adc[i] = 0;
00739     m_pmd_sm_hit[i] = 0;
00740   }
00741   for (i=0; i<24; i++) {
00742     m_pmd_chain_adc[i] = 0;
00743     m_pmd_chain_hit[i] = 0;
00744   }
00745   m_pmd_total_hit = 0;
00746   m_pmd_total_adc = 0;
00747   m_cpv_total_hit = 0;
00748   m_cpv_total_adc = 0;
00749 
00750 // for SSD
00751   m_pnt_phiSSD = 0;
00752   m_pnt_lwSSD = 0;
00753   m_glb_ssd_phi = 0;
00754   m_prim_ssd_phi = 0;
00755   m_pnt_sizeSSD = 0;
00756   m_pnt_eSSD = 0;
00757 
00758 }
00759 //_____________________________________________________________________________
00760 void StQABookHist::BookHist(Int_t histsSet){
00761 
00762   QAH::preString = QAHistType;
00763 //book histograms --------------
00764   if (histsSet != StQA_MC) {
00765     mNullPrimVtxClass = QAH::H1F("QaNullPrimVtxMult","event primary vertex check",40,-2,2);
00766     mNullPrimVtxClass->SetXTitle("has primary vertex? (yes = 1, no = -1)");
00767     mNullPrimVtxClass->SetYTitle("# of events");
00768   }
00769   if (!(histsSet==StQA_run8     ||
00770         histsSet==StQA_run12all ||
00771         histsSet==StQA_run12      )) silHists = kTRUE;
00772   if (!(histsSet==StQA_run12all ||
00773         histsSet==StQA_run12      )) ftpHists = kTRUE;
00774   BookHistPoint();
00775   BookHistEMC();
00776   if (histsSet == StQA_AuAuOld) {
00777     BookHistBBC();
00778     BookHistFPD();
00779   }
00780   BookHistGlob();
00781   BookHistPrim();
00782   BookHistDE();
00783   BookHistPID();
00784   BookHistVertex();
00785   if (!(histsSet==StQA_run12all ||
00786         histsSet==StQA_run12      )) BookHistPMD();
00787   if (histsSet==StQA_MC) BookHistEval();
00788   if (histsSet==StQA_run8     ||
00789       histsSet==StQA_run12all ||
00790       histsSet==StQA_run12      ) BookHistTOF();
00791   
00792 }
00793 //_____________________________________________________________________________
00794 void StQABookHist::BookHistGlob(){
00795   
00796 // for method MakeGlob - from table globtrk
00797 
00798 // general
00799 
00800   m_globtrk_tot      = QAH::H1F("QaGtrkTot","globtrk: tot num tracks - all",
00801                             40,0.,12500.);
00802   m_globtrk_iflag    = QAH::H1F("QaGtrkFlag","globtrk: iflag - all ",200,-999.,1001.);
00803   m_globtrk_good     = QAH::H1F("QaGtrkGood","globtrk: tot good tracks - all",40,0.,10000.);
00804   m_globtrk_good_sm  = QAH::H1F("QaGtrkGoodsm","globtrk: tot good tracks - all",40,0.,1000.);
00805   if (silHists)
00806   m_globtrk_goodTTS  = QAH::H1F("QaGtrkGoodTTS","globtrk: tot good tracks - tpc,tpc+svt",150,0.,9000.);
00807   m_globtrk_goodF    = QAH::H2F("QaGtrkGoodF","globtrk: tot good tracks - ftpc",150,0.,1500.,150,0.,1500.);
00808   m_globtrk_goodF->SetXTitle("FTPC East");
00809   m_globtrk_goodF->SetYTitle("FTPC West");
00810   m_globtrk_fit_prob = QAH::H1F("QaGtrkFitProb","globtrk: prob. fit is correct",100,0,1.2);
00811   m_det_id           = QAH::H1F("QaGtrkDetId","globtrk: Detector ID good tracks - all",25,0.,25.);
00812   m_dcaToBeamXY      = QAH::H2F("QaGtrkDcaBeamXY","globtrk: xy-DCA to Beam Axis (z=0)",80,-4,4,80,-4,4);
00813   m_dcaToBeamXY->SetXTitle("x");
00814   m_dcaToBeamXY->SetYTitle("y");
00815   m_dcaToBeamZ1      = QAH::H1F("QaGtrkDcaBeamZ1","globtrk: z-DCA to Beam Axis, coarse scale",100,-200,200);
00816   m_dcaToBeamZ2      = QAH::H1F("QaGtrkDcaBeamZ2","globtrk: z-DCA to Beam Axis",100,-50,50);
00817   m_dcaToBeamZ3      = QAH::H1F("QaGtrkDcaBeamZ3","globtrk: z-DCA to Beam Axis, near z=0",80,-20,20);
00818   m_zDcaTanl         = QAH::H2F("QaGtrkZdcaTanl","globtrk: z-DCA to Beam Axis vs tanl, tpc",100,-25,25,32,-4,4);
00819   m_zDcaTanl->SetXTitle("z-DCA");
00820   m_zDcaTanl->SetYTitle("tanl");
00821   m_zDcaZf           = QAH::H2F("QaGtrkZdcaZf","globtrk: z-DCA to Beam Axis vs z-first",100,-25,25,50,-300,300);
00822   m_zDcaZf->SetXTitle("z-DCA");
00823   m_zDcaZf->SetYTitle("z of first point");
00824   m_zDcaPsi          = QAH::H2F("QaGtrkZdcaPsi","globtrk: z-DCA to Beam Axis vs psi",100,-25,25,64,0,360);
00825   m_zDcaPsi->SetXTitle("z-DCA");
00826   m_zDcaPsi->SetYTitle("psi");
00827   m_zDcaPhi0         = QAH::H2F("QaGtrkZdcaPhi0","globtrk: z-DCA to Beam Axis vs azimuth (phi0) at start",80,-20,20,64,0,360);
00828   m_zDcaPhi0->SetXTitle("z-DCA");
00829   m_zDcaPhi0->SetYTitle("phi0");
00830   if (silHists) {
00831   m_glb_sptsTS       = QAH::H1F("QaGtrkSptsTS","globtrk: number of svt points, tpc+svt",4,0.,4.);
00832   m_glb_ratiomTTS    = QAH::MH1F("QaGtrkRnmfTTS","globtrk: ratio Nfit/max pnt, tpc,svt", 55,0.,1.1,2);
00833   m_glb_ratiomTTS->SetMinimum(10);
00834   m_glb_ratiomTTS->Rebin(0,"TPC+SVT");
00835   m_glb_ratiomTTS->Rebin(1,"TPC");
00836   m_glb_ratiomTTS->SetStats(kFALSE);
00837   m_psiTTS           = QAH::MH1F("QaGtrkPsiTTS", "globtrk: psi (deg), tpc,svt", 64, 0.,360.,2);
00838   m_psiTTS->Rebin(0,"TPC+SVT");
00839   m_psiTTS->Rebin(1,"TPC");
00840   m_psiTTS->SetStats(kFALSE);
00841   m_etaTTS           = QAH::MH1F("QaGtrkEtaTTS", "globtrk: eta, tpc,svt",40,-2.,2.,2);
00842   m_etaTTS->Rebin(0,"TPC+SVT");
00843   m_etaTTS->Rebin(1,"TPC");
00844   m_etaTTS->SetStats(kFALSE);
00845   m_pTTTS            = QAH::MH1F("QaGtrkPtTTS",  "globtrk: log10 pT, tpc,svt",50,-1.,4.,2);
00846   m_pTTTS->Rebin(0,"TPC+SVT");
00847   m_pTTTS->Rebin(1,"TPC");
00848   m_pTTTS->SetStats(kFALSE);
00849   m_pTTTS->SetXTitle("log10 pT (MeV)");
00850   m_chisq0TTS        = QAH::MH1F("QaGtrkChisq0TTS", "globtrk: chisq0, tpc,svt", 50, 0.,5.,2);
00851   m_chisq0TTS->Rebin(0,"TPC+SVT");
00852   m_chisq0TTS->Rebin(1,"TPC");
00853   m_chisq0TTS->SetStats(kFALSE);
00854   m_fit_pointTTS     = QAH::MH1F("QaGtrkNPntFitTTS","globtrk: N fit pnts on trk, tpc,svt",50,0.,50.,2);
00855   m_fit_pointTTS->Rebin(0,"TPC+SVT");
00856   m_fit_pointTTS->Rebin(1,"TPC");
00857   m_fit_pointTTS->SetStats(kFALSE);
00858   m_glb_impactTTS    = QAH::MH1F("QaGtrkImpactTTS", "globtrk: log10 impact param from prim vtx, tpc,svt",120,-3.,3.,2);
00859   m_glb_impactTTS->Rebin(0,"TPC+SVT");
00860   m_glb_impactTTS->Rebin(1,"TPC");
00861   m_glb_impactTTS->SetStats(kFALSE);
00862   m_glb_impactrTTS   = QAH::MH1F("QaGtrkImpactrTTS", "globtrk: impact param from prim vtx, tpc,svt",100,0.,300.,2);
00863   m_glb_impactrTTS->Rebin(0,"TPC+SVT");
00864   m_glb_impactrTTS->Rebin(1,"TPC");
00865   m_glb_impactrTTS->SetStats(kFALSE);
00866   m_fpoint_lengthTTS = QAH::H2F("QaGtrkFitPntLTTS","globtrk: N fit pnts vs length, tpc,tpc+svt",
00867                                 25,0.,250.,25,0.,50.);
00868   m_fpoint_lengthTTS->SetXTitle("trk length");
00869   m_fpoint_lengthTTS->SetYTitle("Npoints on trk");
00870   }
00871 
00872 // 1D tpc
00873 
00874   m_pointT      = QAH::H1F("QaGtrkNPntT",   "globtrk: N points on trk,tpc", 50, 0.,50.);
00875   m_max_pointT  = QAH::H1F("QaGtrkNPntMaxT","globtrk: N max pnts on trk, tpc", 50, 0.,100.);
00876   m_fit_pointT  = QAH::H1F("QaGtrkNPntFitT","globtrk: N fit pnts on trk, tpc", 50, 0.,50.);
00877   m_glb_ratiomT = QAH::H1F("QaGtrkRnmT",    "globtrk: ratio Nfit/max pnt, tpc", 55, 0., 1.1);
00878   m_glb_chargeT = QAH::H1F("QaGtrkChrgT",   "globtrk: charge, tpc ", 20,-2.,2.);
00879   m_glb_r0T     = QAH::H1F("QaGtrkR0T",     "globtrk: radius at start (cm), tpc ", 50,0.,200.);
00880     m_glb_r0T->SetMinimum(100);
00881   m_glb_phi0T   = QAH::H1F("QaGtrkPhi0T",   "globtrk: azimuth (phi) at start (deg,force 0,360), tpc ", 64,0.,360.);
00882   m_glb_z0T     = QAH::H1F("QaGtrkZ0T",     "globtrk: z-coord at start (cm), tpc ", 50, -300.,300.);
00883     m_glb_z0T->SetMinimum(500);
00884   m_glb_curvT   = QAH::H1F("QaGtrkCurvT",   "globtrk: log10 curvature (1/cm), tpc ", 80,-3.5,0.5);
00885   m_glb_padfTEW = QAH::MH1F("QaGtrkPadfTEW","globtrk: padrow of first hit on trk, tpc", 45,0.5,45.5,2);
00886   m_glb_padfTEW->Rebin(0,"East");
00887   m_glb_padfTEW->Rebin(1,"West");
00888   m_glb_padfTEW->SetStats(kFALSE);
00889   m_glb_padfT   = QAH::H1F("QaGtrkPadfT",   "globtrk: padrow of first hit on trk, tpc", 45,0.5,45.5);
00890   m_glb_f0      = QAH::MH1F("QaGtrkf0",     "globtrk: first point: hit-helix, tpc",60,-3.,3.,3);
00891   m_glb_f0->Rebin(0,"X");
00892   m_glb_f0->Rebin(1,"Y");
00893   m_glb_f0->Rebin(2,"Z");
00894   m_glb_f0->SetStats(kFALSE);
00895   m_glb_xfT     = QAH::H1F("QaGtrkXfT",     "globtrk: x of first hit on trk, tpc", 50,-200.,200.);
00896   m_glb_xf0     = QAH::H1F("QaGtrkXf0",     "globtrk: first point: x_hit - x_helix, tpc",60,-3.,3.);
00897   m_glb_yfT     = QAH::H1F("QaGtrkYfT",     "globtrk: y of first hit on trk, tpc", 50,-200.,200.);
00898   m_glb_yf0     = QAH::H1F("QaGtrkYf0",     "globtrk: first point: y_hit - y_helix, tpc",60,-3.,3.);
00899   m_glb_zfT     = QAH::H1F("QaGtrkZfT",     "globtrk: z of first hit on trk, tpc", 50,-300.,300.);
00900   m_glb_zf0     = QAH::H1F("QaGtrkZf0",     "globtrk: first point: z_hit - z_helix, tpc",60,-3.,3.);
00901   m_glb_radfT   = QAH::H1F("QaGtrkRT",      "globtrk: radial position of first hit, tpc", 50,0.,200.);
00902   m_glb_rzf0    = QAH::MH1F("QaGtrkRZf0",   "globtrk: first point: hit - helix, tpc",60,-3.,3.,2);
00903   m_glb_rzf0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
00904   m_glb_rzf0->Rebin(1,"z_{dif}");
00905   m_glb_rzf0->SetStats(kFALSE);
00906   m_glb_rzl0    = QAH::MH1F("QaGtrkRZl0",   "globtrk: last point: hit - helix, tpc",60,-3.,3.,2);
00907   m_glb_rzl0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
00908   m_glb_rzl0->Rebin(1,"z_{dif}");
00909   m_glb_rzl0->SetStats(kFALSE);
00910   m_glb_phifT   = QAH::MH1F("QaGtrkPhifT",  "globtrk: phi of first point on trk, tpc",64,0,360,2);
00911   m_glb_phifT->Rebin(0,"East");
00912   m_glb_phifT->Rebin(1,"West");
00913   m_glb_phifT->SetStats(kFALSE);
00914   m_lengthT     = QAH::H1F("QaGtrkLengthT", "globtrk: track length, tpc", 50,0.,300.);
00915   m_psiT        = QAH::H1F("QaGtrkPsiT",    "globtrk: psi, tpc (deg)", 64, 0.,360.);
00916   m_tanlT       = QAH::H1F("QaGtrkTanlT",   "globtrk: tanl, tpc",32,-4.,4.);
00917   m_glb_thetaT  = QAH::H1F("QaGtrkThetaT",  "globtrk: theta (degrees), tpc",36,0.,180.);
00918   m_etaT        = QAH::H1F("QaGtrkEtaT",    "globtrk: eta, tpc",40,-2.,2.);
00919   m_pTT         = QAH::H1F("QaGtrkPtT",     "globtrk: pT, tpc",50,0.,10.);
00920   m_momT        = QAH::H1F("QaGtrkPT",      "globtrk: momentum, tpc",50,0.,10.);
00921   m_chisq0T     = QAH::H1F("QaGtrkChisq0T", "globtrk: chisq0, tpc", 50, 0.,5.);
00922   m_glb_impactT = QAH::MH1F("QaGtrkImpactT", "globtrk: log10 impact param from prim vtx, tpc",
00923                             120,-3.0,3.0,3);
00924   m_glb_impactT->Rebin(0,"East");
00925   m_glb_impactT->Rebin(1,"West");
00926   m_glb_impactT->Rebin(2,"All");
00927   m_glb_impactT->SetStats(kFALSE);
00928   m_glb_simpactT = QAH::MH1F("QaGtrkSImpactT", "globtrk: signed impact param from prim vtx, tpc",
00929                             50,-0.8,0.8,3);
00930   m_glb_simpactT->Rebin(0,"East");
00931   m_glb_simpactT->Rebin(1,"West");
00932   m_glb_simpactT->Rebin(2,"All");
00933   m_glb_simpactT->SetStats(kFALSE);
00934   m_glb_impactrT = QAH::H1F("QaGtrkImpactrT", "globtrk: impact param from prim vtx, tpc",100,0.,300.);
00935 
00936 
00937 // 2D - tpc
00938 
00939   m_pT_eta_recT = QAH::H2F("QaGtrkPtVsEtaT","globtrk: log10 pT vs eta, tpc", 20,-2.,2.,40,1.,4.);
00940     m_pT_eta_recT->SetXTitle("eta");
00941     m_pT_eta_recT->SetYTitle(" log10 pT (MeV)");
00942 
00943   m_globtrk_xf_yfTE = QAH::H2F("QaGtrkXfYfTE",  "globtrk: Y vs X of first hit on trk, tpc east", 40,-200.,200.,40,-200.,200.);
00944     m_globtrk_xf_yfTE->SetXTitle("x first");
00945     m_globtrk_xf_yfTE->SetYTitle("y first");
00946 
00947   m_globtrk_xf_yfTW = QAH::H2F("QaGtrkXfYfTW",  "globtrk: Y vs X of first hit on trk, tpc west", 40,-200.,200.,40,-200.,200.);
00948     m_globtrk_xf_yfTW->SetXTitle("x first");
00949     m_globtrk_xf_yfTW->SetYTitle("y first");
00950 
00951   m_tanl_zfT = QAH::H2F("QaGtrkTanlzf","globtrk: tanl(dip) vs. (zfirst-zvtx)/arc length, tpc,tpc+svt",60,-3.,3.,60,-3.,3.);
00952     m_tanl_zfT->SetXTitle("(zfirst-zvtx)/arc length");
00953     m_tanl_zfT->SetYTitle("tanl");
00954 
00955   m_mom_trklengthT = QAH::H2F("QaGtrkPVsTrkLength","globtrk: log10 mom vs trk length, tpc",
00956                              50,0.,250.,40,1.,4.);
00957     m_mom_trklengthT->SetXTitle("trk length");  
00958     m_mom_trklengthT->SetYTitle("log10 P (MeV)");
00959 
00960   m_eta_trklengthT = QAH::H2F("QaGtrkLengthVEtaT","globtrk: trk length vs eta, tpc",
00961                              20,-2.,2.,50,0.,250.);
00962     m_eta_trklengthT->SetXTitle("eta");
00963     m_eta_trklengthT->SetYTitle("length");
00964 
00965   m_npoint_lengthT = QAH::H2F("QaGtrkNPntLengthT","globtrk: N pnts vs length, tpc",
00966                              25,0.,250.,25,0.,50.);
00967     m_npoint_lengthT->SetXTitle("trk length");
00968     m_npoint_lengthT->SetYTitle("Npoints on trk");
00969 
00970   m_fpoint_lengthT = QAH::H2F("QaGtrkFitPntLengthT","globtrk: N fit pnts vs length, tpc",
00971                              25,0.,250.,25,0.,50.);
00972     m_fpoint_lengthT->SetXTitle("trk length");
00973     m_fpoint_lengthT->SetYTitle("Npoints on trk");
00974 
00975   m_chisq0_momT = QAH::H2F("QaGtrkChi0MomT","globtrk: Chisq0 vs log10 mom, tpc",40,1.,4.,50,0.,5.);
00976     m_chisq0_momT->SetXTitle("log10 P (MeV)");
00977     m_chisq0_momT->SetYTitle("chisq0") ;
00978 
00979   m_chisq0_etaT = QAH::H2F("QaGtrkChi0EtaT","globtrk: Chisq0 vs eta, tpc",20,-2.,2.,20,0.,5.);
00980     m_chisq0_etaT->SetXTitle("eta");
00981     m_chisq0_etaT->SetYTitle("chisq0");
00982 
00983   m_chisq0_dipT = QAH::H2F("QaGtrkChi0TanlT","globtrk: Chisq0 vs tanl(dip), tpc",20,-3.2,3.2,20,0.,5.);
00984     m_chisq0_dipT->SetXTitle("dip angle");
00985     m_chisq0_dipT->SetYTitle("chisq0");
00986 
00987   m_chisq0_zfT = QAH::H2F("QaGtrkChi0zfT","globtrk: Chisq0 vs zfirst, tpc",20,-250.,250.,20,0.,5.);
00988     m_chisq0_zfT->SetXTitle("zfirst");
00989     m_chisq0_zfT->SetYTitle("chisq0");
00990 
00991   m_chisq0_phiT = QAH::H2F("QaGtrkChi0PhiT","globtrk: Chisq0 vs phi, tpc",20,0.,360,20,0.,5.);
00992     m_chisq0_phiT->SetXTitle("phi");
00993     m_chisq0_phiT->SetYTitle("chisq0");
00994 
00995   m_psi_phiT = QAH::H2F("QaGtrkPsiPhiT","globtrk: psi vs phi, tpc",36, 0.,360.,36,0.,360.);
00996      m_psi_phiT->SetXTitle("phi");
00997      m_psi_phiT->SetYTitle("psi");
00998 
00999   if (silHists) {
01000 
01001 // 1D tpc + silicon (svt+ssd)
01002 
01003   m_pointTS      = QAH::H1F("QaGtrkNPntTS",   "globtrk: N points on trk,tpc+svt", 60, 0.,60.);
01004   m_max_pointTS  = QAH::H1F("QaGtrkNPntMaxTS","globtrk: N max pnts on trk, tpc+svt", 50, 0.,100.);
01005   m_fit_pointTS  = QAH::H1F("QaGtrkNPntFitTS","globtrk: N fit pnts on trk, tpc+svt", 60, 0.,60.);
01006   m_glb_ratiomTS = QAH::H1F("QaGtrkRnmTS",    "globtrk: ratio Nfit/max pnt, tpc+svt", 55, 0., 1.1);
01007   m_glb_chargeTS = QAH::H1F("QaGtrkChrgTS",   "globtrk: charge, tpc+svt ", 20,-2.,2.);
01008   m_glb_r0TS     = QAH::H1F("QaGtrkR0TS",     "globtrk: radius at start (cm), tpc+svt", 100,0.,25.);
01009     m_glb_r0TS->SetMinimum(100);
01010   m_glb_phi0TS   = QAH::H1F("QaGtrkPhi0TS",   "globtrk: azimuth (phi) at start (deg,force 0-360),tpc+svt", 64, 0.,360.);
01011   m_glb_z0TS     = QAH::H1F("QaGtrkZ0TS",     "globtrk: z-coord at start (cm), tpc+svt", 50, -50.,50.);
01012     m_glb_z0TS->SetMinimum(500);
01013   m_glb_curvTS   = QAH::H1F("QaGtrkCurvTS",   "globtrk: log10 curvature (1/cm), tpc+svt", 80,-3.5,0.5);
01014   m_glb_f0TS     = QAH::MH1F("QaGtrkf0TS",    "globtrk: first point: hit-helix, tpc+svt",50,-5.,5.,3);
01015   m_glb_f0TS->Rebin(0,"X");
01016   m_glb_f0TS->Rebin(1,"Y");
01017   m_glb_f0TS->Rebin(2,"Z");
01018   m_glb_f0TS->SetStats(kFALSE);
01019   m_glb_xfTS     = QAH::H1F("QaGtrkXfTS",     "globtrk: x of first hit on trk, tpc+svt", 50,-200.,200.);
01020   m_glb_xf0TS    = QAH::H1F("QaGtrkXf0TS",    "globtrk: x of first hit - on helix at start, tpc+svt",50,-5.,5.);
01021   m_glb_yfTS     = QAH::H1F("QaGtrkYfTS",     "globtrk: y of first hit on trk, tpc+svt", 50,-200.,200.);
01022   m_glb_yf0TS    = QAH::H1F("QaGtrkYf0TS",    "globtrk: y of first hit - on helix at start, tpc+svt",50,-5.,5.);
01023   m_glb_zfTS     = QAH::H1F("QaGtrkZfTS",     "globtrk: z of first hit on trk, svt", 50,-30.,30.);
01024   m_glb_zf0TS    = QAH::H1F("QaGtrkZf0TS",    "globtrk: z of first hit - on helix at start, tpc+svt",50,-5.,5.);
01025   m_glb_radfTS   = QAH::H1F("QaGtrkRTS",      "globtrk: radial position of first hit, tpc+svt",80,4.,200.);
01026   m_glb_rzf0TS   = QAH::MH1F("QaGtrkRZf0TS",  "globtrk: first point: hit - helix, tpc+svt",60,-3.,3.,2);
01027   m_glb_rzf0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
01028   m_glb_rzf0TS->Rebin(1,"z_{dif}");
01029   m_glb_rzf0TS->SetStats(kFALSE);
01030   m_glb_rzl0TS   = QAH::MH1F("QaGtrkRZl0TS",  "globtrk: last point: hit - helix, tpc+svt",60,-3.,3.,2);
01031   m_glb_rzl0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
01032   m_glb_rzl0TS->Rebin(1,"z_{dif}");
01033   m_glb_rzl0TS->SetStats(kFALSE);
01034   m_glb_phifTS   = QAH::H1F("QaGtrkPhifTS",   "globtrk: phi of first point on track, svt",64,0,360);
01035   m_glb_ssd_phi  = QAH::H1F("QaGtrkPhifSSD",  "globtrk: phi of ssd point",64,0,360);
01036   m_lengthTS     = QAH::H1F("QaGtrkLengthTS", "globtrk: track length, tpc+svt", 50,0.,300.);
01037   m_psiTS        = QAH::H1F("QaGtrkPsiTS",    "globtrk: psi, tpc+svt (deg) ", 64, 0.,360.);
01038   m_tanlTS       = QAH::H1F("QaGtrkTanlTS",   "globtrk: tanl, tpc+svt",32,-4.,4.);
01039   m_glb_thetaTS  = QAH::H1F("QaGtrkThetaTS",  "globtrk: theta (degrees), tpc+svt",36,0.,180.);
01040   m_etaTS        = QAH::H1F("QaGtrkEtaTS",    "globtrk: eta, tpc+svt",40,-2.,2.);
01041   m_pTTS         = QAH::H1F("QaGtrkPtTS",     "globtrk: pT, tpc+svt",50,0.,10.);
01042   m_momTS        = QAH::H1F("QaGtrkPTS",      "globtrk: momentum, tpc+svt",50,0.,10.);
01043   m_chisq0TS     = QAH::H1F("QaGtrkChisq0TS", "globtrk: chisq0, tpc+svt", 50, 0.,5.);
01044   m_glb_impactTS = QAH::MH1F("QaGtrkImpactTS", "globtrk: log10 impact param from prim vtx, tpc+svt",
01045                             120,-3.0,3.0,3);
01046   m_glb_impactTS->Rebin(0,"East");
01047   m_glb_impactTS->Rebin(1,"West");
01048   m_glb_impactTS->Rebin(2,"All");
01049   m_glb_impactTS->SetStats(kFALSE);
01050   m_glb_simpactTS = QAH::MH1F("QaGtrkSImpactTS", "globtrk: signed impact param from prim vtx, tpc+svt",
01051                             50,-0.8,0.8,3);
01052   m_glb_simpactTS->Rebin(0,"East");
01053   m_glb_simpactTS->Rebin(1,"West");
01054   m_glb_simpactTS->Rebin(2,"All");
01055   m_glb_simpactTS->SetStats(kFALSE);
01056   m_glb_impactrTS = QAH::H1F("QaGtrkImpactrTS", "globtrk: impact param from prim vtx, tpc+svt",
01057                             100,0.,30.);
01058 
01059 
01060 // 2D - tpc + silicon (svt + ssd)
01061 
01062   m_pT_eta_recTS = QAH::H2F("QaGtrkPtVsEtaTS","globtrk: log10 pT vs eta, tpc+svt", 20,-2.,2.,40,1.,4.);
01063     m_pT_eta_recTS->SetXTitle("eta");
01064     m_pT_eta_recTS->SetYTitle(" log10 pT (MeV)");
01065 
01066   m_globtrk_xf_yfTS = QAH::H2F("QaGtrkXfYfTS",  "globtrk: Y vs X of first hit on trk, tpc+svt", 40,-200.,200.,40,-200.,200.);
01067     m_globtrk_xf_yfTS->SetXTitle("x first");
01068     m_globtrk_xf_yfTS->SetYTitle("y first");
01069 
01070   m_tanl_zfTS = QAH::H2F("QaGtrkTanlzfTS","globtrk: tanl(dip) versus (zfirst-zvtx)/arc length, svt",60,-3.,3.,60,-3.,3.);
01071     m_tanl_zfTS->SetXTitle("(zfirst-zvtx)/arc length");
01072     m_tanl_zfTS->SetYTitle("tanl");
01073 
01074   m_mom_trklengthTS = QAH::H2F("QaGtrkPVsTrkLTS","globtrk: log10 mom vs trk length, tpc+svt",50,0.,250.,40,1.,4.);
01075     m_mom_trklengthTS->SetXTitle("trk length");  
01076     m_mom_trklengthTS->SetYTitle("log10 P (MeV)");
01077 
01078   m_eta_trklengthTS = QAH::H2F("QaGtrkLVEtaTS","globtrk: trk length vs eta, tpc+svt",20,-2.,2.,50,0.,250.);
01079     m_eta_trklengthTS->SetXTitle("eta");
01080     m_eta_trklengthTS->SetYTitle("length");
01081 
01082   m_npoint_lengthTS = QAH::H2F("QaGtrkNPntLTS","globtrk: N pnts vs length, tpc+svt",25,0.,250.,25,0.,50.);
01083     m_npoint_lengthTS->SetXTitle("trk length");
01084     m_npoint_lengthTS->SetYTitle("Npoints on trk");
01085 
01086   m_fpoint_lengthTS = QAH::H2F("QaGtrkFitPntLTS","globtrk: N fit pnts vs length, tpc+svt",25,0.,250.,25,0.,50.);
01087     m_fpoint_lengthTS->SetXTitle("trk length");
01088     m_fpoint_lengthTS->SetYTitle("Npoints on trk");
01089 
01090   m_chisq0_momTS = QAH::H2F("QaGtrkChi0MomTS","globtrk: Chisq0 vs log10 mom, tpc+svt",40,1.,4.,50,0.,5.);
01091     m_chisq0_momTS->SetXTitle("log10 P (MeV)");
01092     m_chisq0_momTS->SetYTitle("chisq0") ;
01093 
01094   m_chisq0_etaTS = QAH::H2F("QaGtrkChi0EtaTS","globtrk: Chisq0 vs eta, tpc+svt",20,-2.,2.,20,0.,5.);
01095     m_chisq0_etaTS->SetXTitle("eta");
01096     m_chisq0_etaTS->SetYTitle("chisq0");
01097 
01098   m_chisq0_dipTS = QAH::H2F("QaGtrkChi0TanlTS","globtrk: Chisq0 vs tanl(dip), tpc+svt",20,-3.2,3.2,20,0.,5.);
01099     m_chisq0_dipTS->SetXTitle("dip angle");
01100     m_chisq0_dipTS->SetYTitle("chisq0");
01101 
01102   m_chisq0_zfTS = QAH::H2F("QaGtrkChi0zfTS","globtrk: Chisq0 vs zfirst, tpc+svt",20,-250.,250.,20,0.,5.);
01103     m_chisq0_zfTS->SetXTitle("zfirst");
01104     m_chisq0_zfTS->SetYTitle("chisq0");
01105 
01106   m_chisq0_phiTS = QAH::H2F("QaGtrkChi0PhiTS","globtrk: Chisq0 vs phi, tpc+svt",20,0.,360,20,0.,5.);
01107     m_chisq0_phiTS->SetXTitle("phi");
01108     m_chisq0_phiTS->SetYTitle("chisq0");
01109 
01110   m_psi_phiTS = QAH::H2F("QaGtrkPsiPhiTS","globtrk: psi vs phi, tpc+svt",36, 0.,360.,36,0.,360.);
01111      m_psi_phiTS->SetXTitle("phi");
01112      m_psi_phiTS->SetYTitle("psi");
01113 
01114   }
01115 
01116   if (ftpHists) { // FTPC hists
01117 
01118 // 1D ftpc
01119 
01120   // both east (solid) and west(dashed) on same plot
01121   m_pointF      = QAH::MH1F("QaGtrkNPntF",    "globtrk: N points on trk,ftpc", 8,4.,12.,2);
01122   m_pointF->Rebin(0,"East");
01123   m_pointF->Rebin(1,"West");
01124   m_pointF->SetStats(kFALSE);
01125   m_max_pointF  = QAH::MH1F("QaGtrkNPntMaxF", "globtrk: N max pnts on trk, ftpc", 8,4.,12.,2);
01126   m_max_pointF->Rebin(0,"East");
01127   m_max_pointF->Rebin(1,"West");
01128   m_max_pointF->SetStats(kFALSE);
01129   m_glb_ratiomF = QAH::MH1F("QaGtrkRnmF",     "globtrk: ratio Nfit/max pnt, ftpc", 11,0.,1.1,2);
01130   m_glb_ratiomF->Rebin(0,"East");
01131   m_glb_ratiomF->Rebin(1,"West");
01132   m_glb_ratiomF->SetStats(kFALSE);
01133   m_glb_chargeF = QAH::MH1F("QaGtrkChrgF",    "globtrk: charge, ftpc", 20,-2.,2.,2);
01134   m_glb_chargeF->Rebin(0,"East");
01135   m_glb_chargeF->Rebin(1,"West");
01136   m_glb_chargeF->SetStats(kFALSE);
01137   m_glb_xfF     = QAH::MH1F("QaGtrkXfF",      "globtrk: x of first hit on trk, ftpc",70,-35.,35.,2);
01138   m_glb_xfF->Rebin(0,"East");
01139   m_glb_xfF->Rebin(1,"West");
01140   m_glb_xfF->SetStats(kFALSE);
01141   m_glb_yfF     = QAH::MH1F("QaGtrkYfF",      "globtrk: y of first hit on trk, ftpc",70,-35.,35.,2);
01142   m_glb_yfF->Rebin(0,"East");
01143   m_glb_yfF->Rebin(1,"West");
01144   m_glb_yfF->SetStats(kFALSE);
01145   m_glb_zfF     = QAH::MH1F("QaGtrkZfF",      "globtrk: z of first hit on trk, ftpc",75,-225.,225.,2);
01146   m_glb_zfF->Rebin(0,"East");
01147   m_glb_zfF->Rebin(1,"West");
01148   m_glb_zfF->SetStats(kFALSE);
01149   m_glb_planefF = QAH::MH1F("QaGtrkPlanefF",  "globtrk: plane of first hit on trk, ftpc",20,0.5,20.5,2);
01150   m_glb_planefF->Rebin(0,"East");
01151   m_glb_planefF->Rebin(1,"West");
01152   m_glb_planefF->SetStats(kFALSE);
01153   m_glb_radfF   = QAH::MH1F("QaGtrkRF",       "globtrk: radial position of first hit, ftpc",35,0.,35.,2);
01154   m_glb_radfF->Rebin(0,"East");
01155   m_glb_radfF->Rebin(1,"West");
01156   m_glb_radfF->SetStats(kFALSE);
01157   m_lengthF     = QAH::MH1F("QaGtrkLengthF",  "globtrk: track length, ftpc",60,0.,120.,2);
01158   m_lengthF->Rebin(0,"East");
01159   m_lengthF->Rebin(1,"West");
01160   m_lengthF->SetStats(kFALSE);
01161   m_psiF        = QAH::MH1F("QaGtrkPsiF",     "globtrk: psi, ftpc",90,0.,360.,2);
01162   m_psiF->Rebin(0,"East");
01163   m_psiF->Rebin(1,"West");
01164   m_psiF->SetStats(kFALSE);
01165   m_etaF        = QAH::MH1F("QaGtrkEtaF",     "globtrk: |eta|, ftpc",80,2,5,2);
01166   m_etaF->Rebin(0,"East");
01167   m_etaF->Rebin(1,"West");
01168   m_etaF->SetStats(kFALSE);
01169   m_pTF         = QAH::MH1F("QaGtrkPtF",      "globtrk: pT, ftpc",50,0.,10.,2);
01170   m_pTF->Rebin(0,"East");
01171   m_pTF->Rebin(1,"West");
01172   m_pTF->SetStats(kFALSE);
01173   m_momF        = QAH::MH1F("QaGtrkPF",       "globtrk: momentum, ftpc",50,0.,10.,2);
01174   m_momF->Rebin(0,"East");
01175   m_momF->Rebin(1,"West");
01176   m_momF->SetStats(kFALSE);
01177   m_chisq0F     = QAH::MH1F("QaGtrkChisq0F",  "globtrk: chi2/ndf x,y fit, ftpc",50,0.,5.,2);
01178   m_chisq0F->Rebin(0,"East");
01179   m_chisq0F->Rebin(1,"West");
01180   m_chisq0F->SetStats(kFALSE);
01181   m_chisq1F     = QAH::MH1F("QaGtrkChisq1F",  "globtrk: chi2/ndf r,z fit, ftpc",50,0.,5.,2);
01182   m_chisq1F->Rebin(0,"East");
01183   m_chisq1F->Rebin(1,"West");
01184   m_chisq1F->SetStats(kFALSE);
01185   m_glb_impactF = QAH::MH1F("QaGtrkImpactF", "globtrk: log10 impact param from prim vtx, ftpc",120,-3.0,3.,2);
01186   m_glb_impactF->Rebin(0,"East");
01187   m_glb_impactF->Rebin(1,"West");
01188   m_glb_impactF->SetStats(kFALSE);
01189   m_glb_impactrF = QAH::MH1F("QaGtrkImpactrF", "globtrk: impact param from prim vtx, ftpc",100,0.,10.,2);
01190   m_glb_impactrF->Rebin(0,"East");
01191   m_glb_impactrF->Rebin(1,"West");
01192   m_glb_impactrF->SetStats(kFALSE);
01193   // separate east and west plots
01194   m_pointFE      = QAH::H1F("QaGtrkNPntFE",    "globtrk: N points on trk,ftpc east", 8, 4.,12.);
01195   m_pointFW      = QAH::H1F("QaGtrkNPntFW",    "globtrk: N points on trk,ftpc west", 8, 4.,12.);
01196   m_max_pointFE  = QAH::H1F("QaGtrkNPntMaxFE", "globtrk: N max pnts on trk, ftpc east", 8, 4.,12.);
01197   m_max_pointFW  = QAH::H1F("QaGtrkNPntMaxFW", "globtrk: N max pnts on trk, ftpc west", 8, 4.,12.);
01198   m_glb_ratiomFE = QAH::H1F("QaGtrkRnmFE",     "globtrk: ratio Nfit/max pnt, ftpc east", 55, 0., 1.1);
01199   m_glb_ratiomFW = QAH::H1F("QaGtrkRnmFW",     "globtrk: ratio Nfit/max pnt, ftpc west", 55, 0., 1.1);
01200   m_glb_chargeFE = QAH::H1F("QaGtrkChrgFE",    "globtrk: charge, ftpc east ", 20,-2.,2.);
01201   m_glb_chargeFW = QAH::H1F("QaGtrkChrgFW",    "globtrk: charge, ftpc west ", 20,-2.,2.);
01202   m_glb_xfFE     = QAH::H1F("QaGtrkXfFE",      "globtrk: x of first hit on trk, ftpc east", 70,-35.,35.);
01203   m_glb_xfFW     = QAH::H1F("QaGtrkXfFW",      "globtrk: x of first hit on trk, ftpc west", 70,-35.,35.);
01204   m_glb_yfFE     = QAH::H1F("QaGtrkYfFE",      "globtrk: y of first hit on trk, ftpc east", 70,-35.,35.);
01205   m_glb_yfFW     = QAH::H1F("QaGtrkYfFW",      "globtrk: y of first hit on trk, ftpc west", 70,-35.,35.);
01206   m_glb_zfFE     = QAH::H1F("QaGtrkZfFE",      "globtrk: z of first hit on trk, ftpc east", 75,-225.,-150.);
01207   m_glb_zfFW     = QAH::H1F("QaGtrkZfFW",      "globtrk: z of first hit on trk, ftpc west", 75,150.,225.);
01208   m_glb_radfFE   = QAH::H1F("QaGtrkRFE",       "globtrk: radial position of first hit, ftpc east",35,0.,35.);
01209   m_glb_radfFW   = QAH::H1F("QaGtrkRFW",       "globtrk: radial position of first hit, ftpc west",35,0.,35.);
01210   m_lengthFE     = QAH::H1F("QaGtrkLengthFE",  "globtrk: track length, ftpc east", 60,0.,120.);
01211   m_lengthFW     = QAH::H1F("QaGtrkLengthFW",  "globtrk: track length, ftpc west", 60,0.,120.);
01212   m_psiFE        = QAH::H1F("QaGtrkPsiFE",     "globtrk: psi, ftpc east", 90, 0.,360.);
01213   m_psiFW        = QAH::H1F("QaGtrkPsiFW",     "globtrk: psi, ftpc west", 90, 0.,360.);
01214   m_etaFE        = QAH::H1F("QaGtrkEtaFE",     "globtrk: eta, ftpc east",80,-4.5,-2.);
01215   m_etaFW        = QAH::H1F("QaGtrkEtaFW",     "globtrk: eta, ftpc west",80,2.,4.5);
01216   m_pTFE         = QAH::H1F("QaGtrkPtFE",      "globtrk: pT, ftpc east",50,0.,10.);
01217   m_pTFW         = QAH::H1F("QaGtrkPtFW",      "globtrk: pT, ftpc west",50,0.,10.);
01218   m_momFE        = QAH::H1F("QaGtrkPFE",       "globtrk: momentum, ftpc east ",50,0.,5.);
01219   m_momFW        = QAH::H1F("QaGtrkPFW",       "globtrk: momentum, ftpc west ",50,0.,5.);
01220   m_chisq0FE     = QAH::H1F("QaGtrkChisq0FE",  "globtrk: chi2/ndf x,y fit, ftpc east", 50, 0.,5.);
01221   m_chisq0FW     = QAH::H1F("QaGtrkChisq0FW",  "globtrk: chi2/ndf x,y fit, ftpc west", 50, 0.,5.);
01222   m_chisq1FE     = QAH::H1F("QaGtrkChisq1FE",  "globtrk: chi2/ndf r,z fit, ftpc east", 50, 0.,5.);
01223   m_chisq1FW     = QAH::H1F("QaGtrkChisq1FW",  "globtrk: chi2/ndf r,z fit, ftpc west", 50, 0.,5.);
01224 
01225 // 2D - ftpc
01226 
01227    m_pT_eta_recFE = QAH::H2F("QaGtrkPtVsEtaFE","globtrk: log10 pT vs eta, ftpcE",20,-4.5,-2.,40,1.,4.);
01228     m_pT_eta_recFE->SetXTitle("eta");
01229     m_pT_eta_recFE->SetYTitle("log10 pT");
01230    m_pT_eta_recFW = QAH::H2F("QaGtrkPtVsEtaFW","globtrk: log10 pT vs eta, ftpcW",20,2.,4.5,40,1.,4.);
01231     m_pT_eta_recFW->SetXTitle("eta");
01232     m_pT_eta_recFW->SetYTitle("log10 pT");
01233   m_globtrk_xf_yfFE = QAH::H2F("QaGtrkXfYfFE","globtrk: Y vs X of first hit on trk, ftpcE", 70,-35.,35.,70,-35.,35.);
01234     m_globtrk_xf_yfFE->SetXTitle("x first");
01235     m_globtrk_xf_yfFE->SetYTitle("y first");
01236   m_globtrk_xf_yfFW = QAH::H2F("QaGtrkXfYfFW","globtrk: Y vs X of first hit on trk, ftpcW", 70,-35.,35.,70,-35.,35.);
01237     m_globtrk_xf_yfFW->SetXTitle("x first");
01238     m_globtrk_xf_yfFW->SetYTitle("y first");
01239   m_globtrk_padtimeFE = QAH::H2F("QaGtrkPadTimeFtpcE","globtrk: #pads vs #timebins, ftpcE",12,0.5,12.5,10,0.5,10.5);
01240 //  m_globtrk_padtimeFE = QAH::H2F("QaGtrkPadTimeFtpcE","globtrk: #pads vs #timebins of hits, ftpcE",12,0.5,12.5,10,0.5,10.5);
01241     m_globtrk_padtimeFE->SetXTitle("#timebins");
01242     m_globtrk_padtimeFE->SetYTitle("#pads");
01243   m_globtrk_padtimeFW = QAH::H2F("QaGtrkPadTimeFtpcW","globtrk: #pads vs #timebins, ftpcW",12,0.5,12.5,10,0.5,10.5);
01244 //  m_globtrk_padtimeFW = QAH::H2F("QaGtrkPadTimeFtpcW","globtrk: #pads vs #timebins of hits, ftpcW",12,0.5,12.5,10,0.5,10.5);
01245     m_globtrk_padtimeFW->SetXTitle("#timebins");
01246     m_globtrk_padtimeFW->SetYTitle("#pads");
01247   m_eta_trklengthFE = QAH::H2F("QaGtrkLengthVEtaFE","globtrk: trk length vs eta, ftpcE",25,-4.5,-2.,30,0.,120.);
01248     m_eta_trklengthFE->SetXTitle("eta");
01249     m_eta_trklengthFE->SetYTitle("length");
01250   m_eta_trklengthFW = QAH::H2F("QaGtrkLengthVEtaFW","globtrk: trk length vs eta, ftpcW",25,2.,4.5,30,0.,120.);
01251     m_eta_trklengthFW->SetXTitle("eta");
01252     m_eta_trklengthFW->SetYTitle("length");
01253   m_npoint_lengthFE = QAH::H2F("QaGtrkNPntLengthFE","globtrk: N pnts vs length, ftpcE",30,0.,120.,15,0.,15.);
01254     m_npoint_lengthFE->SetXTitle("trk length");
01255     m_npoint_lengthFE->SetYTitle("Npoints on trk");
01256   m_npoint_lengthFW = QAH::H2F("QaGtrkNPntLengthFW","globtrk: N pnts vs length, ftpcW",30,0.,120.,15,0.,15.);
01257     m_npoint_lengthFW->SetXTitle("trk length");
01258     m_npoint_lengthFW->SetYTitle("Npoints on trk");
01259 
01260   } // ftpHists
01261 }
01262 //____________________________________________________
01263 void StQABookHist::BookHistPrim(){
01264 
01265 // for method MakeHistPrim - from table primtrk
01266 
01267 // 1D
01268   m_primtrk_tot     = QAH::H1F("QaPtrkTot",   "primtrk: tot num tracks",50,0.,5000.);
01269   m_primtrk_tot_sm  = QAH::H1F("QaPtrkTotsm", "primtrk: tot num tracks",50,0.,200.);
01270   m_primtrk_iflag   = QAH::H1F("QaPtrkFlag",  "primtrk: iflag - all",160,-799.,900.);
01271   m_primtrk_good    = QAH::H1F("QaPtrkGood",  "primtrk: tot num tracks iflag>0",50,0.,5000.);
01272   m_primtrk_good_sm = QAH::H1F("QaPtrkGoodsm","primtrk: tot num tracks iflag>0",50,0.,200.);
01273   if (silHists)
01274   m_primtrk_goodTTS = QAH::H1F("QaPtrkGoodTTS","primtrk: tot num tracks iflag>0, tpc,svt",150,0.,4500.);
01275   m_primtrk_goodF   = QAH::H2F("QaPtrkGoodF",  "primtrk: tot num tracks iflag>0, ftpc",150,0.,1500.,150,0.,1500.);
01276   m_primtrk_goodF->SetXTitle("East");
01277   m_primtrk_goodF->SetYTitle("West");
01278   m_primglob_good   = QAH::H1F("QaPtrkGlob","primtrk: ratio primary/global tracks w/ iflag>0",50,0,1);
01279   m_primglob_fit    = QAH::H1F("QaPtrkGlobFit","primtrk: ratio primary/global nfit points",50,0,2);
01280   m_pdet_id         = QAH::H1F("QaPtrkDetId",  "primtrk: Detector ID good tracks - all",25,0.,25.);
01281   if (silHists) {
01282   m_primtrk_meanptTTS = QAH::MH1F("QaPtrkMeanPtTTS","primtrk: <pT>, tpc, tpc+svt",50,0.,2.,2);
01283   m_primtrk_meanptTTS->Rebin(0,"TPC+SVT");
01284   m_primtrk_meanptTTS->Rebin(1,"TPC");
01285   m_primtrk_meanptTTS->SetStats(kFALSE);
01286   }
01287   m_primtrk_meanptF = QAH::MH1F("QaPtrkMeanPtF","primtrk: <pT>, ftpc",50,0.,2.,2);
01288   m_primtrk_meanptF->Rebin(0,"East");
01289   m_primtrk_meanptF->Rebin(1,"West");
01290   m_primtrk_meanptF->SetStats(kFALSE);
01291   if (silHists) {
01292   m_primtrk_meanetaTTS = QAH::MH1F("QaPtrkMeanEtaTTS","primtrk: <eta>, tpc, tpc+svt",40,-2.,2.,2);
01293   m_primtrk_meanetaTTS->Rebin(0,"TPC+SVT");
01294   m_primtrk_meanetaTTS->Rebin(1,"TPC");
01295   m_primtrk_meanetaTTS->SetStats(kFALSE);
01296   }
01297   m_primtrk_meanetaF = QAH::MH1F("QaPtrkMeanEtaF","primtrk: |<eta>|, ftpc",40,2,5,2);
01298   m_primtrk_meanetaF->Rebin(0,"East");
01299   m_primtrk_meanetaF->Rebin(1,"West");
01300   m_primtrk_meanetaF->SetStats(kFALSE);
01301   if (silHists) {
01302   m_ppsiTTS         = QAH::MH1F("QaPtrkPsiTTS","primtrk: psi (deg), tpc, svt", 36, 0.,360.,2);
01303   m_ppsiTTS->Rebin(0,"TPC+SVT");
01304   m_ppsiTTS->Rebin(1,"TPC");
01305   m_ppsiTTS->SetStats(kFALSE);
01306   m_petaTTS         = QAH::MH1F("QaPtrkEtaTTS","primtrk: eta, tpc,svt",40,-2.,2.,2);
01307   m_petaTTS->Rebin(0,"TPC+SVT");
01308   m_petaTTS->Rebin(1,"TPC");
01309   m_petaTTS->SetStats(kFALSE);
01310   m_ppTTTS          = QAH::MH1F("QaPtrkPtTTS", "primtrk: pT, tpc,svt",50,0.,5.,2);
01311   m_ppTTTS->Rebin(0,"TPC+SVT");
01312   m_ppTTTS->Rebin(1,"TPC");
01313   m_ppTTTS->SetStats(kFALSE);
01314   m_pchisq0TTS     = QAH::MH1F("QaPtrkChisq0TTS", "primtrk: chisq0, tpc,svt", 50, 0.,5.,2);
01315   m_pchisq0TTS->Rebin(0,"TPC+SVT");
01316   m_pchisq0TTS->Rebin(1,"TPC");
01317   m_pchisq0TTS->SetStats(kFALSE);
01318   m_pfpoint_lengthTTS = QAH::H2F("QaPtrkFitPntLTTS","primtrk: N fit pnts vs length, tpc,tpc+svt",25,70.,350.,25,0.,50.);
01319   m_pfpoint_lengthTTS->SetXTitle("trk length");
01320   m_pfpoint_lengthTTS->SetYTitle("Npoints on trk");
01321   }
01322 
01323 // 1D tpc
01324   m_ppointT      = QAH::H1F("QaPtrkNPntT",   "primtrk: N points on trk,tpc", 60, 0.,60.);
01325   m_pmax_pointT  = QAH::H1F("QaPtrkNPntMaxT","primtrk: N max pnts on trk, tpc", 50, 0.,100.);
01326   m_pfit_pointT  = QAH::H1F("QaPtrkNPntFitT","primtrk: N fit pnts on trk, tpc", 60, 0.,60.);
01327   m_prim_ratiomT = QAH::H1F("QaPtrkRnmT",    "primtrk: ratio Nfit/max pnt, tpc", 55, 0., 1.1);
01328   m_prim_chargeT = QAH::H1F("QaPtrkChrgT",   "primtrk: charge, tpc ", 20,-2.,2.);
01329   m_prim_r0T     = QAH::H1F("QaPtrkR0T",     "primtrk: radius at start (cm), tpc ", 50,0.,.1);
01330   m_prim_phi0T   = QAH::H1F("QaPtrkPhi0T",   "primtrk: azimuth (phi) at start (deg,force 0,360), tpc ", 64,0.,360.);
01331   m_prim_z0T     = QAH::H1F("QaPtrkZ0T",     "primtrk: z-coord at start (cm), tpc ", 50, -300.,300.);
01332   m_prim_curvT   = QAH::H1F("QaPtrkCurvT",   "primtrk: log10 curvature (1/cm), tpc ", 80,-3.5,0.5);
01333   m_prim_f0      = QAH::MH1F("QaPtrkf0",     "primtrk: first point: hit-helix, tpc",60,-3.,3.,3);
01334   m_prim_f0->Rebin(0,"X");
01335   m_prim_f0->Rebin(1,"Y");
01336   m_prim_f0->Rebin(2,"Z");
01337   m_prim_xfT     = QAH::H1F("QaPtrkXfT",     "primtrk: x of first hit on trk, tpc", 50,-200.,200.);
01338   m_prim_xf0     = QAH::H1F("QaPtrkXf0",     "primtrk: first point: x_hit - x_helix, tpc",60,-3.,3.);
01339   m_prim_yfT     = QAH::H1F("QaPtrkYfT",     "primtrk: y of first hit on trk, tpc", 50,-200.,200.);
01340   m_prim_yf0     = QAH::H1F("QaPtrkYf0",     "primtrk: first point: y_hit - y_helix, tpc",60,-3.,3.);
01341   m_prim_zfT     = QAH::H1F("QaPtrkZfT",     "primtrk: z of first hit on trk, tpc", 50,-200.,200.);
01342   m_prim_zf0     = QAH::H1F("QaPtrkZf0",     "primtrk: first point: z_hit - z_helix, tpc",60,-3.,3.);
01343   m_prim_radfT   = QAH::H1F("QaPtrkRT",      "primtrk: radial position of first hit, tpc", 50,0.,200.);
01344   m_prim_rzf0    = QAH::MH1F("QaPtrkRZf0",   "primtrk: first point: hit - helix (r,z), tpc",60,-3.,3.,2);
01345   m_prim_rzf0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
01346   m_prim_rzf0->Rebin(1,"z_{dif}");
01347   m_prim_rzf0->SetStats(kFALSE);
01348   m_prim_rzl0    = QAH::MH1F("QaPtrkRZl0",   "primtrk: last point: hit - helix (r,z), tpc",60,-3.,3.,2);
01349   m_prim_rzl0->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
01350   m_prim_rzl0->Rebin(1,"z_{dif}");
01351   m_prim_rzl0->SetStats(kFALSE);
01352   m_plengthT     = QAH::H1F("QaPtrkLengthT", "primtrk: track length, tpc", 50,0.,300.);
01353   m_ppsiT        = QAH::H1F("QaPtrkPsiT",    "primtrk: psi, tpc (deg)", 36, 0.,360.);
01354   m_ptanlT       = QAH::H1F("QaPtrkTanlT",   "primtrk: tanl, tpc",32,-4.,4.);
01355   m_prim_thetaT  = QAH::H1F("QaPtrkThetaT",  "primtrk: theta (degrees), tpc",36,0.,180.);
01356   m_petaT        = QAH::H1F("QaPtrkEtaT",    "primtrk: eta, tpc",40,-2.,2.);
01357   m_ppTT         = QAH::H1F("QaPtrkPtT",     "primtrk: pT, tpc",50,0.,5.);
01358   m_pmomT        = QAH::H1F("QaPtrkPT",      "primtrk: momentum, tpc",50,0.,5.);
01359   m_pchisq0T     = QAH::H1F("QaPtrkChisq0T", "primtrk: chisq0, tpc", 50, 0.,5.);
01360 
01361 // 2D - tpc
01362   m_ppT_eta_recT = QAH::H2F("QaPtrkPtVsEtaT","primtrk: log10 pT vs eta, tpc", 20,-2.,2.,40,1.,4.);
01363     m_ppT_eta_recT->SetXTitle("eta");
01364     m_ppT_eta_recT->SetYTitle(" log10 pT (MeV)");
01365   m_primtrk_xf_yfTE = QAH::H2F("QaPtrkXfYfTE",  "primtrk: Y vs X of first hit on trk, tpc east", 40,-200.,200.,40,-200.,200.);
01366     m_primtrk_xf_yfTE->SetXTitle("x first");
01367     m_primtrk_xf_yfTE->SetYTitle("y first");
01368   m_primtrk_xf_yfTW = QAH::H2F("QaPtrkXfYfTW",  "primtrk: Y vs X of first hit on trk, tpc west", 40,-200.,200.,40,-200.,200.);
01369     m_primtrk_xf_yfTW->SetXTitle("x first");
01370     m_primtrk_xf_yfTW->SetYTitle("y first");
01371   m_ptanl_zfT = QAH::H2F("QaPtrkTanlzf","primtrk: tanl(dip) versus zfirst-zvtx, tpc,tpc+svt",60,-3.,3.,60,-3.,3.);
01372     m_ptanl_zfT->SetXTitle("zfirst-zvtx");
01373     m_ptanl_zfT->SetYTitle("tanl");
01374   m_pmom_trklengthT = QAH::H2F("QaPtrkPVsTrkLength","primtrk: log10 mom vs trk length, tpc",50,70.,350.,40,1.,4.);
01375     m_pmom_trklengthT->SetXTitle("trk length");  
01376     m_pmom_trklengthT->SetYTitle("log10 P (MeV)");
01377   m_peta_trklengthT = QAH::H2F("QaPtrkLengthVEtaT","primtrk: trk length vs eta, tpc",20,-2.,2.,50,70.,350.);
01378     m_peta_trklengthT->SetXTitle("eta");
01379     m_peta_trklengthT->SetYTitle("length");
01380   m_pnpoint_lengthT = QAH::H2F("QaPtrkNPntLengthT","primtrk: N pnts vs length, tpc",25,70.,350.,25,0.,50.);
01381     m_pnpoint_lengthT->SetXTitle("trk length");
01382     m_pnpoint_lengthT->SetYTitle("Npoints on trk");
01383   m_pfpoint_lengthT = QAH::H2F("QaPtrkFitPntLengthT","primtrk: N fit pnts vs length, tpc",25,70.,350.,25,0.,50.);
01384     m_pfpoint_lengthT->SetXTitle("trk length");
01385     m_pfpoint_lengthT->SetYTitle("Npoints on trk");
01386   m_pchisq0_momT = QAH::H2F("QaPtrkChi0MomT","primtrk: Chisq0 vs log10 mom, tpc",40,1.,4.,50,0.,5.);
01387     m_pchisq0_momT->SetXTitle("log10 P (MeV)");
01388     m_pchisq0_momT->SetYTitle("chisq0") ;
01389   m_pchisq0_etaT = QAH::H2F("QaPtrkChi0EtaT","primtrk: Chisq0 vs eta, tpc",20,-2.,2.,20,0.,5.);
01390     m_pchisq0_etaT->SetXTitle("eta");
01391     m_pchisq0_etaT->SetYTitle("chisq0");
01392   m_pchisq0_dipT = QAH::H2F("QaPtrkChi0TanlT","primtrk: Chisq0 vs tanl(dip), tpc",20,-3.2,3.2,20,0.,5.);
01393     m_pchisq0_dipT->SetXTitle("dip angle");
01394     m_pchisq0_dipT->SetYTitle("chisq0");
01395   m_pchisq0_zfT = QAH::H2F("QaPtrkChi0zfT","primtrk: Chisq0 vs zfirst, tpc",20,-200.,200.,20,0.,5.);
01396     m_pchisq0_zfT->SetXTitle("zfirst");
01397     m_pchisq0_zfT->SetYTitle("chisq0");
01398   m_ppsi_phiT = QAH::H2F("QaPtrkPsiPhiT","primtrk: psi vs phi, tpc",36, 0.,360.,36,0.,360.);
01399      m_ppsi_phiT->SetXTitle("phi");
01400      m_ppsi_phiT->SetYTitle("psi");
01401 
01402   if (silHists) {
01403 // 1D tpc + silicon (svt+ssd)
01404   m_ppointTS      = QAH::H1F("QaPtrkNPntTS",   "primtrk: N points on trk,tpc+svt", 60, 0.,60.);
01405   m_pmax_pointTS  = QAH::H1F("QaPtrkNPntMaxTS","primtrk: N max pnts on trk, tpc+svt", 50, 0.,100.);
01406   m_pfit_pointTS  = QAH::H1F("QaPtrkNPntFitTS","primtrk: N fit pnts on trk, tpc+svt", 60, 0.,60.);
01407   m_prim_ratiomTS = QAH::H1F("QaPtrkRnmTS",    "primtrk: ratio Nfit/max pnt, tpc+svt", 55, 0., 1.2005);
01408   m_prim_chargeTS = QAH::H1F("QaPtrkChrgTS",   "primtrk: charge, tpc+svt ", 20,-2.,2.);
01409   m_prim_r0TS     = QAH::H1F("QaPtrkR0TS",     "primtrk: radius at start (cm), tpc+svt", 50,0.,0.1);
01410   m_prim_phi0TS   = QAH::H1F("QaPtrkPhi0TS",   "primtrk: azimuth (phi) at start (deg,force 0-360),tpc+svt", 64, 0.,360.);
01411   m_prim_z0TS     = QAH::H1F("QaPtrkZ0TS",     "primtrk: z-coord at start (cm), tpc+svt", 50, -50.,50.);
01412   m_prim_curvTS   = QAH::H1F("QaPtrkCurvTS",   "primtrk: log10 curvature (1/cm), tpc+svt", 80,-3.5,0.5);
01413   m_prim_f0TS     = QAH::MH1F("QaPtrkf0TS",    "primtrk: first point: hit-helix, tpc+svt",50,-5.,5.,3);
01414   m_prim_f0TS->Rebin(0,"X");
01415   m_prim_f0TS->Rebin(1,"Y");
01416   m_prim_f0TS->Rebin(2,"Z");
01417   m_prim_f0TS->SetStats(kFALSE);
01418   m_prim_xfTS     = QAH::H1F("QaPtrkXfTS",     "primtrk: x of first hit on trk, tpc+svt", 50,-200.,200.);
01419   m_prim_xf0TS    = QAH::H1F("QaPtrkXf0TS",    "primtrk: first point: x_hit - x_helix, tpc+svt",50,-5.,5.);
01420   m_prim_yfTS     = QAH::H1F("QaPtrkYfTS",     "primtrk: y of first hit on trk, tpc+svt", 50,-200.,200.);
01421   m_prim_yf0TS    = QAH::H1F("QaPtrkYf0TS",    "primtrk: first point: y_hit - y_helix, tpc+svt",50,-5.,5.);
01422   m_prim_zfTS     = QAH::H1F("QaPtrkZfTS",     "primtrk: z of first hit on trk, tpc+svt", 50,-200.,200.);
01423   m_prim_zf0TS    = QAH::H1F("QaPtrkZf0TS",    "primtrk: first point: z_hit - z_helix, tpc+svt",50,-5.,5.);
01424   m_prim_radfTS   = QAH::H1F("QaPtrkRTS",      "primtrk: radial position of first hit, tpc+svt", 50,0.,200.);
01425   m_prim_rzf0TS   = QAH::MH1F("QaPtrkRZf0TS",  "primtrk: first point: hit - helix (r,z), tpc+svt",60,-3.,3.,2);
01426   m_prim_rzf0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
01427   m_prim_rzf0TS->Rebin(1,"z_{dif}");
01428   m_prim_rzf0TS->SetStats(kFALSE);
01429   m_prim_rzl0TS   = QAH::MH1F("QaPtrkRZl0TS",  "primtrk: last point: hit - helix (r,z), tpc+svt",60,-3.,3.,2);
01430   m_prim_rzl0TS->Rebin(0,"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
01431   m_prim_rzl0TS->Rebin(1,"z_{dif}");
01432   m_prim_rzl0TS->SetStats(kFALSE);
01433   m_plengthTS     = QAH::H1F("QaPtrkLengthTS", "primtrk: track length, tpc+svt", 50,0.,300.);
01434   m_ppsiTS        = QAH::H1F("QaPtrkPsiTS",    "primtrk: psi, tpc+svt (deg) ", 36, 0.,360.);
01435   m_ptanlTS       = QAH::H1F("QaPtrkTanlTS",   "primtrk: tanl, tpc+svt",32,-4.,4.);
01436   m_prim_thetaTS  = QAH::H1F("QaPtrkThetaTS",  "primtrk: theta (degrees), tpc+svt",36,0.,180.);
01437   m_petaTS        = QAH::H1F("QaPtrkEtaTS",    "primtrk: eta, tpc+svt",40,-2.,2.);
01438   m_ppTTS         = QAH::H1F("QaPtrkPtTS",     "primtrk: pT, tpc+svt",50,0.,5.);
01439   m_pmomTS        = QAH::H1F("QaPtrkPTS",      "primtrk: momentum, tpc+svt",50,0.,5.);
01440   m_pchisq0TS     = QAH::H1F("QaPtrkChisq0TS", "primtrk: chisq0, tpc+svt", 50, 0.,5.);
01441   m_prim_ssd_phi  = QAH::H1F("QaPtrkPhifSSD",  "primtrk: phi of ssd point",64,0,360);
01442 
01443 // 2D - tpc + silicon (svt + ssd)
01444   m_ppT_eta_recTS = QAH::H2F("QaPtrkPtVsEtaTS","primtrk: log10 pT vs eta, tpc+svt", 20,-2.,2.,40,1.,4.);
01445     m_ppT_eta_recTS->SetXTitle("eta");
01446     m_ppT_eta_recTS->SetYTitle(" log10 pT (MeV)");
01447   m_primtrk_xf_yfTS = QAH::H2F("QaPtrkXfYfTS",  "primtrk: Y vs X of first hit on trk, tpc+svt", 40,-200.,200.,40,-200.,200.);
01448     m_primtrk_xf_yfTS->SetXTitle("x first");
01449     m_primtrk_xf_yfTS->SetYTitle("y first");
01450   m_ptanl_zfTS = QAH::H2F("QaPtrkTanlzfTS","primtrk: tanl(dip) versus zfirst-zvtx, svt",60,-3.,3.,60,-3.,3.);
01451     m_ptanl_zfTS->SetXTitle("zfirst-zvtx");
01452     m_ptanl_zfTS->SetYTitle("tanl");
01453   m_pmom_trklengthTS = QAH::H2F("QaPtrkPVsTrkLTS","primtrk: log10 mom vs trk length, tpc+svt",50,70.,350.,40,1.,4.);
01454     m_pmom_trklengthTS->SetXTitle("trk length");  
01455     m_pmom_trklengthTS->SetYTitle("log10 P (MeV)");
01456   m_peta_trklengthTS = QAH::H2F("QaPtrkLVEtaTS","primtrk: trk length vs eta, tpc+svt",20,-2.,2.,50,70.,350.);
01457     m_peta_trklengthTS->SetXTitle("eta");
01458     m_peta_trklengthTS->SetYTitle("length");
01459   m_pnpoint_lengthTS = QAH::H2F("QaPtrkNPntLTS","primtrk: N pnts vs length, tpc+svt",25,70.,350.,25,0.,50.);
01460     m_pnpoint_lengthTS->SetXTitle("trk length");
01461     m_pnpoint_lengthTS->SetYTitle("Npoints on trk");
01462   m_pfpoint_lengthTS = QAH::H2F("QaPtrkFitPntLTS","primtrk: N fit pnts vs length, tpc+svt",25,70.,350.,25,0.,50.);
01463     m_pfpoint_lengthTS->SetXTitle("trk length");
01464     m_pfpoint_lengthTS->SetYTitle("Npoints on trk");
01465   m_pchisq0_momTS = QAH::H2F("QaPtrkChi0MomTS","primtrk: Chisq0 vs log10 mom, tpc+svt",40,1.,4.,50,0.,5.);
01466     m_pchisq0_momTS->SetXTitle("log10 P (MeV)");
01467     m_pchisq0_momTS->SetYTitle("chisq0") ;
01468   m_pchisq0_etaTS = QAH::H2F("QaPtrkChi0EtaTS","primtrk: Chisq0 vs eta, tpc+svt",20,-2.,2.,20,0.,5.);
01469     m_pchisq0_etaTS->SetXTitle("eta");
01470     m_pchisq0_etaTS->SetYTitle("chisq0");
01471   m_pchisq0_dipTS = QAH::H2F("QaPtrkChi0TanlTS","primtrk: Chisq0 vs tanl(dip), tpc+svt",20,-3.2,3.2,20,0.,5.);
01472     m_pchisq0_dipTS->SetXTitle("dip angle");
01473     m_pchisq0_dipTS->SetYTitle("chisq0");
01474   m_pchisq0_zfTS = QAH::H2F("QaPtrkChi0zfTS","primtrk: Chisq0 vs zfirst, tpc+svt",20,-200.,200.,20,0.,5.);
01475     m_pchisq0_zfTS->SetXTitle("zfirst");
01476     m_pchisq0_zfTS->SetYTitle("chisq0");
01477   m_ppsi_phiTS = QAH::H2F("QaPtrkPsiPhiTS","primtrk: psi vs phi, tpc+svt",36, 0.,360.,36,0.,360.);
01478     m_ppsi_phiTS->SetXTitle("phi");
01479     m_ppsi_phiTS->SetYTitle("psi");
01480   } // silHists
01481 
01482   if (ftpHists) { // FTPC hists
01483 
01484 // 1D ftpc
01485   // east (solid) and west(dashed) on same plot
01486   m_ppointF      = QAH::MH1F("QaPtrkNPntF",    "primtrk: N points on trk,ftpc",15,0.,15.,2);
01487   m_ppointF->Rebin(0,"East");
01488   m_ppointF->Rebin(1,"West");
01489   m_ppointF->SetStats(kFALSE);
01490   m_pmax_pointF  = QAH::MH1F("QaPtrkNPntMaxF", "primtrk: N max pnts on trk, ftpc",15,0.,15.,2);
01491   m_pmax_pointF->Rebin(0,"East");
01492   m_pmax_pointF->Rebin(1,"West");
01493   m_pmax_pointF->SetStats(kFALSE);
01494   m_prim_ratiomF = QAH::MH1F("QaPtrkRnmF",     "primtrk: ratio Nfit/max pnt, ftpc",55,0.,1.1,2);
01495   m_prim_ratiomF->Rebin(0,"East");
01496   m_prim_ratiomF->Rebin(1,"West");
01497   m_prim_ratiomF->SetStats(kFALSE);
01498   m_prim_chargeF = QAH::MH1F("QaPtrkChrgF",    "primtrk: charge, ftpc",20,-2.,2.,2);
01499   m_prim_chargeF->Rebin(0,"East");
01500   m_prim_chargeF->Rebin(1,"West");
01501   m_prim_chargeF->SetStats(kFALSE);
01502   m_prim_xfF     = QAH::MH1F("QaPtrkXfF",      "primtrk: x of first hit on trk, ftpc",70,-35.,35.,2);
01503   m_prim_xfF->Rebin(0,"East");
01504   m_prim_xfF->Rebin(1,"West");
01505   m_prim_xfF->SetStats(kFALSE);
01506   m_prim_yfF     = QAH::MH1F("QaPtrkYfF",      "primtrk: y of first hit on trk, ftpc",70,-35.,35.,2);
01507   m_prim_yfF->Rebin(0,"East");
01508   m_prim_yfF->Rebin(1,"West");
01509   m_prim_yfF->SetStats(kFALSE);
01510   m_prim_zfF     = QAH::MH1F("QaPtrkZfF",      "primtrk: z of first hit on trk, ftpc",75,-225.,225.,2);
01511   m_prim_zfF->Rebin(0,"East");
01512   m_prim_zfF->Rebin(1,"West");
01513   m_prim_zfF->SetStats(kFALSE);
01514   m_prim_radfF   = QAH::MH1F("QaPtrkRF",       "primtrk: radial position of first hit, ftpc",35,0.,35.,2);
01515   m_prim_radfF->Rebin(0,"East");
01516   m_prim_radfF->Rebin(1,"West");
01517   m_prim_radfF->SetStats(kFALSE);
01518   m_plengthF     = QAH::MH1F("QaPtrkLengthF",  "primtrk: track length, ftpc",50,150.,300.,2);
01519   m_plengthF->Rebin(0,"East");
01520   m_plengthF->Rebin(1,"West");
01521   m_plengthF->SetStats(kFALSE);
01522   m_ppsiF        = QAH::MH1F("QaPtrkPsiF",     "primtrk: psi, ftpc", 90, 0.,360.,2);
01523   m_ppsiF->Rebin(0,"East");
01524   m_ppsiF->Rebin(1,"West");
01525   m_ppsiF->SetStats(kFALSE);
01526   m_petaF        = QAH::MH1F("QaPtrkEtaF",     "primtrk: |eta|, ftpc",80,2,5,2);
01527   m_petaF->Rebin(0,"East");
01528   m_petaF->Rebin(1,"West");
01529   m_petaF->SetStats(kFALSE);
01530   m_ppTF         = QAH::MH1F("QaPtrkPtF",      "primtrk: pT, ftpc",50,0.,5.,2);
01531   m_ppTF->Rebin(0,"East");
01532   m_ppTF->Rebin(1,"West");
01533   m_ppTF->SetStats(kFALSE);
01534   m_pmomF        = QAH::MH1F("QaPtrkPF",       "primtrk: momentum, ftpc",50,0.,5.,2);
01535   m_pmomF->Rebin(0,"East");
01536   m_pmomF->Rebin(1,"West");
01537   m_pmomF->SetStats(kFALSE);
01538   m_pchisq0F     = QAH::MH1F("QaPtrkChisq0F",  "primtrk: chi2/ndf x,y fit, ftpc",50,0.,5.,2);
01539   m_pchisq0F->Rebin(0,"East");
01540   m_pchisq0F->Rebin(1,"West");
01541   m_pchisq0F->SetStats(kFALSE);
01542   m_pchisq1F     = QAH::MH1F("QaPtrkChisq1F",  "primtrk: chi2/ndf r,z fit, ftpc",50,0.,5.,2);
01543   m_pchisq1F->Rebin(0,"East");
01544   m_pchisq1F->Rebin(1,"West");
01545   m_pchisq1F->SetStats(kFALSE);
01546   // separate east and west histograms
01547   m_ppointFE      = QAH::H1F("QaPtrkNPntFE",    "primtrk: N points on trk,ftpc east", 15, 0.,15.);
01548   m_ppointFW      = QAH::H1F("QaPtrkNPntFW",    "primtrk: N points on trk,ftpc west", 15, 0.,15.);
01549   m_pmax_pointFE  = QAH::H1F("QaPtrkNPntMaxFE", "primtrk: N max pnts on trk, ftpc east", 15, 0.,15.);
01550   m_pmax_pointFW  = QAH::H1F("QaPtrkNPntMaxFW", "primtrk: N max pnts on trk, ftpc west", 15, 0.,15.);
01551   m_prim_ratiomFE = QAH::H1F("QaPtrkRnmFE",     "primtrk: ratio Nfit/max pnt, ftpc east", 55, 0., 1.1);
01552   m_prim_ratiomFW = QAH::H1F("QaPtrkRnmFW",     "primtrk: ratio Nfit/max pnt, ftpc west", 55, 0., 1.1);
01553   m_prim_chargeFE = QAH::H1F("QaPtrkChrgFE",    "primtrk: charge, ftpc east ", 20,-2.,2.);
01554   m_prim_chargeFW = QAH::H1F("QaPtrkChrgFW",    "primtrk: charge, ftpc west ", 20,-2.,2.);
01555   m_prim_xfFE     = QAH::H1F("QaPtrkXfFE",      "primtrk: x of first hit on trk, ftpc east",70,-35.,35.);
01556   m_prim_xfFW     = QAH::H1F("QaPtrkXfFW",      "primtrk: x of first hit on trk, ftpc west",70,-35.,35.);
01557   m_prim_yfFE     = QAH::H1F("QaPtrkYfFE",      "primtrk: y of first hit on trk, ftpc east",70,-35.,35.);
01558   m_prim_yfFW     = QAH::H1F("QaPtrkYfFW",      "primtrk: y of first hit on trk, ftpc west",70,-35.,35.);
01559   m_prim_zfFE     = QAH::H1F("QaPtrkZfFE",      "primtrk: z of first hit on trk, ftpc east",75,-225.,-150.);
01560   m_prim_zfFW     = QAH::H1F("QaPtrkZfFW",      "primtrk: z of first hit on trk, ftpc west",75,150.,225.);
01561   m_prim_radfFE   = QAH::H1F("QaPtrkRFE",       "primtrk: radial position of first hit, ftpc east",35,0.,35.);
01562   m_prim_radfFW   = QAH::H1F("QaPtrkRFW",       "primtrk: radial position of first hit, ftpc west",35,0.,35.);
01563   m_plengthFE     = QAH::H1F("QaPtrkLengthFE",  "primtrk: track length, ftpc east",50,150.,300.);
01564   m_plengthFW     = QAH::H1F("QaPtrkLengthFW",  "primtrk: track length, ftpc west",50,150.,300.);
01565   m_ppsiFE        = QAH::H1F("QaPtrkPsiFE",     "primtrk: psi, ftpc east", 90, 0.,360.);
01566   m_ppsiFW        = QAH::H1F("QaPtrkPsiFW",     "primtrk: psi, ftpc west", 90, 0.,360.);
01567   m_petaFE        = QAH::H1F("QaPtrkEtaFE",     "primtrk: eta, ftpc east",80,-4.5,-2.);
01568   m_petaFW        = QAH::H1F("QaPtrkEtaFW",     "primtrk: eta, ftpc west",80,2.,4.5);
01569   m_ppTFE         = QAH::H1F("QaPtrkPtFE",      "primtrk: pT, ftpc east",50,0.,5.);
01570   m_ppTFW         = QAH::H1F("QaPtrkPtFW",      "primtrk: pT, ftpc west",50,0.,5.);
01571   m_pmomFE        = QAH::H1F("QaPtrkPFE",       "primtrk: momentum, ftpc east ",50,0.,5.);
01572   m_pmomFW        = QAH::H1F("QaPtrkPFW",       "primtrk: momentum, ftpc west ",50,0.,5.);
01573   m_pchisq0FE     = QAH::H1F("QaPtrkChisq0FE",  "primtrk: chi2/ndf x,y fit, ftpc east", 50, 0.,5.);
01574   m_pchisq0FW     = QAH::H1F("QaPtrkChisq0FW",  "primtrk: chi2/ndf x,y fit, ftpc west", 50, 0.,5.);
01575   m_pchisq1FE     = QAH::H1F("QaPtrkChisq1FE",  "primtrk: chi2/ndf r,z fit, ftpc east", 50, 0.,5.);
01576   m_pchisq1FW     = QAH::H1F("QaPtrkChisq1FW",  "primtrk: chi2/ndf r,z ft, ftpc west", 50, 0.,5.);
01577 
01578 // 2D - ftpc
01579   m_ppT_eta_recFE = QAH::H2F("QaPtrkPtVsEtaFE","primtrk: log10 pT vs eta, ftpcE",20,-4.5,-2.,40,1.,4.);
01580   m_ppT_eta_recFW = QAH::H2F("QaPtrkPtVsEtaFW","primtrk: log10 pT vs eta, ftpcW",20,2.,4.5,40,1.,4.);
01581   m_primtrk_xf_yfFE = QAH::H2F("QaPtrkXfYfFE","primtrk: Y vs X of first hit on trk, ftpcE",70,-35.,35.,75,-35.,35.);
01582     m_primtrk_xf_yfFE->SetXTitle("x first");
01583     m_primtrk_xf_yfFE->SetYTitle("y first");
01584   m_primtrk_xf_yfFW = QAH::H2F("QaPtrkXfYfFW","primtrk: Y vs X of first hit on trk, ftpcW",70,-35.,35.,70,-35.,35.);
01585     m_primtrk_xf_yfFW->SetXTitle("x first");
01586     m_primtrk_xf_yfFW->SetYTitle("y first");
01587   m_peta_trklengthFE = QAH::H2F("QaPtrkLengthVEtaFE","primtrk: trk length vs eta, ftpcE",25,-4.5,-2.,50,0.,300.);
01588     m_peta_trklengthFE->SetXTitle("eta");
01589     m_peta_trklengthFE->SetYTitle("length");
01590   m_peta_trklengthFW = QAH::H2F("QaPtrkLengthVEtaFW","primtrk: trk length vs eta, ftpcW",25,2.,4.5,50,0.,300.);
01591     m_peta_trklengthFW->SetXTitle("eta");
01592     m_peta_trklengthFW->SetYTitle("length");
01593   m_pnpoint_lengthFE = QAH::H2F("QaPtrkNPntLengthFE","primtrk: N pnts vs length, ftpcE",50,0.,300.,15,0.,15.);
01594     m_pnpoint_lengthFE->SetXTitle("trk length");
01595     m_pnpoint_lengthFE->SetYTitle("Npoints on trk");
01596   m_pnpoint_lengthFW = QAH::H2F("QaPtrkNPntLengthFW","primtrk: N pnts vs length, ftpcW",50,0.,300.,15,0.,15.);
01597     m_pnpoint_lengthFW->SetXTitle("trk length");
01598     m_pnpoint_lengthFW->SetYTitle("Npoints on trk");
01599 
01600   } // ftpHists
01601 
01602 // 2D - SVT drift length
01603   if (silHists)
01604   m_svt_loc = QAH::H2F("QaPtrkSvtLoc","primtrk: SVT hit time bins",256,0,128,432,-0.5,431.5);
01605 }
01606 //_____________________________________________________________________________
01607 void StQABookHist::BookHistDE(){
01608   
01609   // for method MakeDE - from table dst_dedx
01610   m_ndedxr  = QAH::H1F("QaDedxNum",     "dedx: number of tracks", 64,0., 16000.); 
01611 
01612   m_ndedxT   = QAH::H1F("QaDedxNdedxT", "dedx: number of point to define dE/dx, tpc",50,0., 50.);  
01613   m_dedx0T   = QAH::H1F("QaDedxDedx0T", "dedx: dE/dx mean (GeV/cm), tpc", ndedx, 0., 0.0005);
01614   m_dedx1T   = QAH::H1F("QaDedxDedx1T", "dedx: dE/dx error on mean, tpc", ndedx, 0., 1);
01615   m_dedxTTS  = QAH::H1F("QaDedxBBTTS",  "dedx: ratio <dE/dx> mean to Bethe Bloch <dE/dx>, tpc,tpc+svt",
01616                         50,0.,10.);
01617   // east and west on same plot
01618   m_ndedxF   = QAH::MH1F("QaDedxNdedxF", "dedx: number of point to define dE/dx, ftpc",50,0.,10.,2);
01619   m_ndedxF->Rebin(0,"East");
01620   m_ndedxF->Rebin(1,"West");
01621   m_ndedxF->SetStats(kFALSE);
01622   m_dedx0F   = QAH::MH1F("QaDedxDedx0F", "dedx: dE/dx mean, ftpc",ndedx,0.,750.,2);
01623   m_dedx0F->Rebin(0,"East");
01624   m_dedx0F->Rebin(1,"West");
01625   m_dedx0F->SetStats(kFALSE);
01626   // east and west on separate plots
01627   m_ndedxFE   = QAH::H1F("QaDedxNdedxFE", "dedx: number of point to define dE/dx, ftpcE",50,0.,10.);
01628   m_dedx0FE   = QAH::H1F("QaDedxDedx0FE", "dedx: dE/dx mean, ftpcE", ndedx,  0., 750.);
01629   
01630   m_ndedxFW   = QAH::H1F("QaDedxNdedxFW", "dedx: number of point to define dE/dx, ftpcW",50,0.,10.);  
01631   m_dedx0FW   = QAH::H1F("QaDedxDedx0FW", "dedx: dE/dx mean, ftpcW", ndedx,  0., 750.);
01632 }
01633 //_____________________________________________________________________________
01634 void StQABookHist::BookHistPID(){
01635   
01636   // for MakeHistPID - from tables globtrk & dst_dedx 
01637   // Spectra/pid histograms. C.Ogilvie
01638   
01639   m_p_dedx_rec = QAH::H2F("QaPidGlobtrkDstdedxPVsDedx","PID: globtrk-dst_dedx,  p vs dedx (reconstructed)",
01640                           cnp,cminp,cmaxp,cndedx,cmindedx,cmaxdedx);
01641   m_p_dedx_rec->SetXTitle("p (GeV)");
01642   m_p_dedx_rec->SetYTitle("dedx (keV/cm)");
01643   
01644 }
01645 //_____________________________________________________________________________
01646 void StQABookHist::BookHistVertex(){
01647   // for MakeHistVertex - from table dst_vertex
01648   
01649   
01650   m_v_num   = QAH::H1F("QaVtxNum",  " vertex: num vertices ",50,0.,10000.);
01651   m_v_num_sm= QAH::H1F("QaVtxNumsm",  " vertex: num vertices ",50,0.,500.);
01652 
01653   m_v_vtxid = QAH::H1F("QaVtxVtxId"," vertex,2ndary: Vertex ID ",10,0.,10.);
01654   m_v_x     = QAH::H1F("QaVtxX",    " vertex,2ndary: x ",50,-250.,250.);
01655   m_v_y     = QAH::H1F("QaVtxY",    " vertex,2ndary: y ",50,-250.,250.);
01656   m_v_z     = QAH::H1F("QaVtxZ",    " vertex,2ndary: z ",50,-250.,250.);
01657   m_v_pchi2 = QAH::H1F("QaVtxChisq"," vertex,2ndary: chisq/dof ",50,0.,50.);
01658   m_v_r     = QAH::H1F("QaVtxR",    " vertex,2ndary: r ",100,0,50);
01659   
01660   m_pv_vtxid = QAH::H1F("QaVtxPrVtxId"," vertex,prim: Vertex ID ",10,0.,10.);
01661   m_pv_x     = QAH::H1F("QaVtxPrX",    " vertex,prim: x ",50,-5.,5.);
01662   m_pv_y     = QAH::H1F("QaVtxPrY",    " vertex,prim: y ",50,-5.,5.);
01663   m_pv_z     = QAH::H1F("QaVtxPrZ",    " vertex,prim: z ",50,-50.,50.);
01664   m_pv_xy    = QAH::H2F("QaVtxPrXY",   " vertex,prim: x versus y",50,-5.,5.,50,-5.,5.);
01665   m_pv_pchi2 = QAH::H1F("QaVtxPrChisq"," vertex,prim: chisq/dof ",40,0.,20.);
01666   m_pv_r     = QAH::H1F("QaVtxPrR",    " vertex,prim: r ",100,0,0.1);
01667 
01668 //  m_vtx_phi_dist  = QAH::H2F("QaV0VtxPhiDist",
01669 //            "V0 azimuthal distribution",36,0.,360.,25,1.,101.);
01670 //  m_vtx_phi_dist->SetXTitle("Mean of phi(V0)");
01671 //  m_vtx_phi_dist->SetYTitle("RMS of phi(V0)");
01672   m_vtx_phi_dist  = QAH::H1F("QaV0VtxPhiDist",
01673             "V0 azimuthal distribution relative to primvtx",36,0.,360.);
01674   m_vtx_r_dist  = QAH::H1F("QaV0VtxRDist",
01675             "V0 radial distribution relative to primvtx",50,0.,200.);
01676   m_vtx_z_dist  = QAH::H1F("QaV0VtxZDist",
01677             "V0 Z distribution relative to primvtx",60,-30.,30.);
01678 
01679   m_v0             = QAH::H1F("QaV0Vtx","log10 total V0s (0 => -0.5)",45,-0.55,3.95);
01680   m_ev0_lama_hist  = QAH::H1F("QaV0LambdaMass","V0: Lambda mass",25,1.05,1.15);
01681   m_ev0_k0ma_hist  = QAH::H1F("QaV0K0Mass","V0: k0 mass",25,.4,.6);
01682 
01683   m_xi_tot     = QAH::H1F("QaXiVtxTot", "log10 total Xis (0 => -0.5)",45,-0.55,3.95);
01684   m_xi_ma_hist = QAH::H1F("QaXiaMass",  "Xi: Xi mass",25,1.2,1.4);
01685 
01686   m_kink_tot   = QAH::H1F("QaKinkTot",  "log10 total kinks (0 => -0.5)",35,-0.55,2.95);
01687 
01688   m_vtx_FtpcEastTpc_xy = QAH::H2F("QaVtxFtpcETpcXY",
01689                                   " vertex,prim: x(ftpcE)-x(tpc) vs y(ftpcE)-y(tpc)",
01690                                   800, -20., 20., 800, -20., 20.);
01691   m_vtx_FtpcEastTpc_z  = QAH::H1F("QaVtxFtpcETpcZ",
01692                                   " vertex,prim: z(ftpcE)-z(tpc)",
01693                                   100, -10., 10.);
01694   m_vtx_FtpcWestTpc_xy = QAH::H2F("QaVtxFtpcWTpcXY",
01695                                   " vertex,prim: x(ftpcW)-x(tpc) vs y(ftpcW)-y(tpc)",
01696                                   800, -20., 20., 800, -20., 20.);
01697   m_vtx_FtpcWestTpc_z  = QAH::H1F("QaVtxFtpcWTpcZ",
01698                                   " vertex,prim: z(ftpcW)-z(tpc)",
01699                                   100, -10., 10.);
01700   
01701 }
01702 //_____________________________________________________________________________
01703 void StQABookHist::BookHistPoint(){
01704 
01705   m_pnt_tot     = QAH::H1F("QaPointTot", "point: # hits total ",100, 0.,400000.);
01706   m_pnt_tot_med = QAH::H1F("QaPointTotmed","point: # hits total ",100, 0.,25000.);
01707   m_pnt_tot_sm  = QAH::H1F("QaPointTotsm", "point: # hits total ",100, 0.,2500.);
01708   m_pnt_id      = QAH::H1F("QaPointId","point: detector ID of hit",30,0.,30.);
01709 
01710   m_pnt_tpc     = QAH::H1F("QaPointTpc",  "point: # hits tpc ",100, 0.,300000.);
01711   if (silHists) {
01712   m_pnt_svt     = QAH::H1F("QaPointSvt",  "point: # hits svt ",600, 0.,15000.);
01713   m_pnt_ssd     = QAH::H1F("QaPointSsd",  "point: # hits ssd ",200, 0.,5000.);
01714   }
01715   // east and west on same plot
01716   m_pnt_ftpc   = QAH::MH1F("QaPointFtpc", "point: # hits ftpc",100,0.,25000.,2);
01717   m_pnt_ftpc->Rebin(0,"East");
01718   m_pnt_ftpc->Rebin(1,"West");
01719   m_pnt_ftpc->SetStats(kFALSE);
01720   // east and west on separate plots
01721   m_pnt_ftpcE   = QAH::H1F("QaPointFtpcE","point: # hits ftpcE ",100, 0.,25000.);
01722   m_pnt_ftpcW   = QAH::H1F("QaPointFtpcW","point: # hits ftpcW ",100, 0.,25000.);
01723   if (silHists) {
01724   m_pnt_svtLaser= QAH::H2F("QaPointSvtLaser","point: laser spots, svt ",150,0,600,65,0.,130.);
01725   m_pnt_svtLaser->SetXTitle("event in file");
01726   m_pnt_svtLaserDiff= QAH::MH2F("QaPointSvtLaserDiff","point: diff of laser spots, svt ",150,0,600,101,9.8,50.2,2);
01727   m_pnt_svtLaserDiff->SetXTitle("event in file");
01728   m_pnt_svtLaserDiff->SetYTitle("time buckets");
01729   m_pnt_svtLaserDiff->Rebin(0,"Laser 1");
01730   m_pnt_svtLaserDiff->Rebin(1,"Laser 2");
01731   m_pnt_xyS     = QAH::H2F("QaPointXYSvt","point: x-y distribution of hits, svt,ssd",125,-25,25,125,-25,25);
01732   }
01733   // Now using polar coords (lego didn't work well because of inner radius and colored zero)
01734   // Drawing a lego plot requires r-phi,r. Drawing a plain polar plot reguirs r,r-phi
01735   m_pnt_rpTW    = QAH::H2F("QaPointRPTpcW","point: r-phi distribution of hits, tpcW",20,58.75,196.75,72,0,TMath::TwoPi());
01736   m_pnt_rpTE    = QAH::H2F("QaPointRPTpcE","point: r-phi distribution of hits, tpcE",23,58.75,196.75,72,0,TMath::TwoPi());
01737   m_z_hits      = QAH::H1F("QaPointZhits","point: z distribution of hits, tpc",100,-210,210);
01738   m_pnt_timeT   = QAH::MH1F("QaPointTimeT","point: time bucket distribution of hits, tpc",45,0,450,2);
01739   m_pnt_timeT->Rebin(0,"East");
01740   m_pnt_timeT->Rebin(1,"West");
01741   m_pnt_timeT->SetStats(kFALSE);
01742   m_pnt_phiT    = QAH::MH1F("QaPointPhiT","point: #phi distribution of hits, tpc",36,0,360,2);
01743   m_pnt_phiT->Rebin(0,"East");
01744   m_pnt_phiT->Rebin(1,"West");
01745   m_pnt_phiT->SetStats(kFALSE);
01746   m_pnt_padrowT = QAH::MH1F("QaPointPadrowT","point: padrow distribution of hits, tpc",45,0.5,45.5,2);
01747   m_pnt_padrowT->Rebin(0,"East");
01748   m_pnt_padrowT->Rebin(1,"West");
01749   m_pnt_padrowT->SetStats(kFALSE);
01750   m_pnt_padrowT->SetXTitle("padrow number");
01751 
01752   if (silHists) {
01753   m_pnt_zS      = QAH::H1F("QaPointZhitsS","point: z distribution of hits, svt",100,-35,35);
01754   m_pnt_phiS    = QAH::H1F("QaPointPhiS","point: #phi distribution of hits, svt",36,0,360);
01755   m_pnt_barrelS = QAH::H1F("QaPointBarrelS","point: barrel distribution of hits, svt",3,0.5,3.5);
01756   m_pnt_barrelS->SetXTitle("barrel number");
01757 
01758   m_pnt_phiSSD  = QAH::H1F("QaPointPhiSSD","point: #phi of hits, ssd",36,0,360);
01759   m_pnt_lwSSD   = QAH::H2F("QaPointLWSSD","point: wafer id vs ladder id, ssd",20,0.5,20.5,16,0.5,16.5);
01760   m_pnt_lwSSD->SetXTitle("Ladder #");
01761   m_pnt_lwSSD->SetYTitle("Wafer #");
01762   m_pnt_sizeSSD = QAH::MH1F("QaPointSizeSSD","point: size of clusters, ssd",10,0.5,10.5,2);
01763   m_pnt_sizeSSD->Rebin(0,"P-side");
01764   m_pnt_sizeSSD->Rebin(1,"N-side");
01765   m_pnt_sizeSSD->SetStats(kFALSE);
01766   m_pnt_eSSD = QAH::H1F("QaPointESSD","point: log10(energy) of hits, ssd",90,-5,-2);
01767   }
01768 
01769   m_pnt_xyFE    = QAH::H2F("QaPointXYFtpcE","point: x-y distribution of hits, ftpcE",70,-35,35,70,-35,35);
01770   m_pnt_xyFW    = QAH::H2F("QaPointXYFtpcW","point: x-y distribution of hits, ftpcW",70,-35,35,70,-35,35);
01771   m_pnt_planeF  = QAH::MH1F("QaPointPlaneF","point: plane distribution of hits, ftpc",20,0.5,20.5,2);
01772   m_pnt_planeF->Rebin(0,"East");
01773   m_pnt_planeF->Rebin(1,"West");
01774   m_pnt_planeF->SetStats(kFALSE);
01775   m_pnt_planeF->SetXTitle("plane number");
01776   m_pnt_padtimeFE    = QAH::H2F("QaPointPadTimeFtpcE","point: #pads vs #timebins of hits, ftpcE",12,0.5,12.5,10,0.5,10.5);
01777   m_pnt_padtimeFE->SetXTitle("#timebins");
01778   m_pnt_padtimeFE->SetYTitle("#pads");
01779   m_pnt_padtimeFW    = QAH::H2F("QaPointPadTimeFtpcW","point: #pads vs #timebins of hits, ftpcW",12,0.5,12.5,10,0.5,10.5);
01780   m_pnt_padtimeFW->SetXTitle("#timebins");
01781   m_pnt_padtimeFW->SetYTitle("#pads");
01782 
01783 
01784 }
01785 //_____________________________________________________________________________
01786 void StQABookHist::BookHistEMC(){
01787 
01788   if (!((gROOT->GetClass("StEmcMath")) && (gROOT->GetClass("StEmcGeom"))))
01789     return;
01790 // Book the hists for SimulatorMaker
01791   m_emc_nhit=QAH::H2F("EmcNHitsVsDet","emc: Number of hit(log) .vs. Detector #",100,0.0,4.5,4,0.5,4.5);
01792   m_emc_etot=QAH::H2F("EmcEtotVsDet","emc: Total energy(log) .vs. Detector #",100,-4.0,4.5,8,0.5,4.5);
01793 
01794   const Char_t* tit={"Barrel"};
01795   const Int_t   nx[4] = {40,40,300,20};
01796   const Int_t   ny[4]  = {120, 120, 60, 900};
01797   Float_t rpi = M_PI + 0.00001;
01798   TString name, title;
01799   TArrayD *xb = StEmcMath::binForSmde();
01800 
01801   for(Int_t i=0; i<4; i++){
01802     name  = detname[i] + "Hits";
01803     title = tit  + detname[i] + " hits dist.";
01804     if(i==2) m_emc_hits[i]=QAH::H2F(name,title, xb->GetSize()-1,xb->GetArray(),ny[i],-rpi,rpi);
01805     else m_emc_hits[i] = QAH::H2F(name,title, nx[i],-1.,+1., ny[i],-rpi, rpi);
01806 
01807     name   = detname[i] + "Energy2D";
01808     title  = tit + detname[i] + " energy dist. in eta&phi";
01809     if(i==2) m_emc_energy2D[i] = QAH::H2F(name,title, xb->GetSize()-1,xb->GetArray(), ny[i],-rpi,rpi);
01810     else m_emc_energy2D[i] = QAH::H2F(name,title, nx[i],-1.,+1., ny[i],-rpi, rpi);
01811 
01812     name  = detname[i] + "Adc";
01813     title = tit + detname[i] + " ADC dist.";
01814     m_emc_adc[i]  = QAH::H1F(name,title, 5001, -0.5, 5000.5);
01815 
01816     name   = detname[i] + "Energy";
01817     title  = tit + detname[i] + " energy dist.";
01818     m_emc_energy[i] = QAH::H1F(name,title, 600, 0.0, 60.0);
01819   }
01820   delete xb;
01821 
01822 // Book the hists for cluster finder
01823   Int_t greta[4]={40,40,300,20};   // eta bins
01824   Int_t grphi[4]={120,120,60,900}; // phi bins  => 16-apr by PAI
01825   Float_t myPI = M_PI + 0.0001;
01826   
01827   m_emc_ncl=QAH::H2F("EmcNcluster","emc: Number of cluster(log) .vs. Detector #",40,0.0,4.0, 4,0.5,4.5);
01828   m_emc_etotCl=QAH::H2F("EmcEcluster" ,"emc: Total PreCluster Energy(log) .vs. Detector #", 60,-2.0,4.0, 4,0.5,4.5);
01829 
01830   Float_t rmsMax=0.026;
01831   Int_t   rmsN=52;
01832   m_emc_sig_e= QAH::H2F("EmcRMSeta" ,"emc: Sigma(eta) .vs. Detector #",rmsN,0.0,rmsMax,4,0.5,4.5);
01833   m_emc_sig_p= QAH::H2F("EmcRMSphi" ,"emc: Sigma(phi) .vs. Detector #",rmsN,0.0,rmsMax,4,0.5,4.5);
01834   for (Int_t i=0; i<4; i++)
01835   {
01836     TString name_h = detname[i] + "_cluster";
01837     TString name_e = detname[i] + "_cluster_energy";
01838     TString tit_h  = detname[i] + " cluster";
01839     TString tit_e  = detname[i] + " energy of cluster";
01840     if(i==2) {
01841       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);
01842       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);
01843     } else {
01844       m_emc_cl[i]     = QAH::H2F(name_h,tit_h,greta[i],-1.0,1.0,grphi[i],-myPI, myPI);
01845       m_emc_energyCl[i] = QAH::H2F(name_e,tit_e,greta[i],-1.0,1.0,grphi[i],-myPI, myPI);
01846     }
01847     
01848 
01849     name  = detname[i] + "ClNum";
01850     title   = "Number hits in cluster for " + detname[i];
01851     m_emc_HitsInCl[i]   = QAH::H1F(name, title, 21, -0.5, 20.5);
01852 
01853     name  = detname[i] + "ClEnergy";
01854     title   = "Energy of cluster for " + detname[i];
01855     m_emc_EnergyCl[i]    = QAH::H1F(name, title, 2000, 0.0, 20.0);
01856 
01857     TString name_eta  = detname[i] + "Eta";
01858     TString tit_eta   = "Eta of clusters for " + detname[i];
01859     TString name_phi  = detname[i] + "Phi";
01860     TString tit_phi   = "Phi of clusters for " + detname[i];
01861     if(i==2) {
01862       TArrayD *xb  = StEmcMath::binForSmde();
01863       if(xb) {
01864         m_emc_EtaInCl[i] = QAH::H1F(name_eta, tit_eta, xb->GetSize()-1, xb->GetArray());
01865         delete xb;
01866       }
01867       m_emc_PhiInCl[i]   = QAH::H1F(name_phi, tit_phi, grphi[i], -M_PI*1.015, M_PI*0.985);
01868     } else { 
01869       m_emc_EtaInCl[i]   = QAH::H1F(name_eta, tit_eta, greta[i], -1., 1.);
01870       m_emc_PhiInCl[i]   = QAH::H1F(name_phi, tit_phi, grphi[i], -myPI, myPI);
01871     }
01872   }
01873   // Book the hists for Emc Points
01874   const TString catname[] = {"EmcCat1", "EmcCat2", "EmcCat3", "EmcCat4"};
01875 
01876   for (Int_t i=0; i<4; i++) {
01877     name = catname[i] + "_Point_Energy";
01878     title = catname[i] + " Point Energy";
01879     m_emc_point_energy[i]= QAH::H1F(name,title,100,0.,10.);
01880 
01881     name = catname[i] + "_Point_Eta";
01882     title = catname[i] + " Point Eta";
01883     m_emc_point_eta[i]= QAH::H1F(name,title,100,-1.,1.);
01884 
01885     name = catname[i] + "_Point_Phi";
01886     title = catname[i] + " Point Phi";
01887     m_emc_point_phi[i]= QAH::H1F(name,title,100,-3.14,3.14);
01888    
01889     name = catname[i] + "_Sigma_Eta";
01890     title = catname[i] + " Sigma Eta";
01891     m_emc_point_sigeta[i]= QAH::H1F(name,title,100,0.,.05);
01892   
01893     name = catname[i] + "_Sigma_Phi";
01894     title = catname[i] + " Sigma Phi";
01895     m_emc_point_sigphi[i]= QAH::H1F(name,title,100,0.,.05);
01896 
01897     name = catname[i] + "_Delta_Eta";
01898     title = catname[i] + " Delta Eta";
01899     m_emc_point_deleta[i]= QAH::H1F(name,title,100,-.2,.2);
01900 
01901     name = catname[i] + "_Delta_Phi";
01902     title = catname[i] + " Delta Phi";
01903     m_emc_point_delphi[i]= QAH::H1F(name,title,100,-.2,.2);
01904 
01905     name = catname[i] + "_Points_Multiplicity";
01906     title = catname[i] + " Points Multiplicity";
01907     m_emc_points[i]= QAH::H1F(name,title,200,0.,1000.);
01908 
01909     name = catname[i] + "_Track_Momenta";
01910     title = catname[i] + " Track Momenta ";
01911     m_emc_point_trmom[i]= QAH::H1F(name,title,100,0.,10.);
01912   }
01913 
01914   m_emc_point_flag= QAH::H1F(" Point Flag","Point Flag",5,0.5,5.5);
01915 
01916 }
01917 //_____________________________________________________________________________
01918 void StQABookHist::BookHistEval(){
01919 
01920 // these only get filled if the geant dataset is available!
01921 
01922    m_geant_reco_pvtx_x  = QAH::H1F("QaGRpvtxDx"," diff geant - reco prim vtx X",
01923                               100, -0.25,0.25);
01924    m_geant_reco_pvtx_y  = QAH::H1F("QaGRpvtxDy"," diff geant - reco prim vtx Y",
01925                               100, -0.25,0.25);
01926    m_geant_reco_pvtx_z  = QAH::H1F("QaGRpvtxDz"," diff geant - reco prim vtx Z",
01927                               100, -0.25,0.25);
01928    m_geant_reco_vtx_z_z = QAH::H2F("QaGRpvtxDzZ",
01929         " reco pvtx Z vs diff geant - reco Z", 100, -0.5,0.5,100,-50.,50.);
01930      m_geant_reco_vtx_z_z->SetXTitle("z vtx resolution (cm)");
01931      m_geant_reco_vtx_z_z->SetYTitle("z position of vtx (cm)");
01932 
01933 }
01934 //_____________________________________________________________________________
01935 void StQABookHist::BookHistBBC(){
01936 
01937   Char_t ID[4];
01938   Int_t i,j;
01939 
01940   m_bbc_adc[0] = QAH::MH1F("QaBbcAdcES","BBC East Small ADC",100,0.5,400.5,8);
01941   m_bbc_adc[1] = QAH::MH1F("QaBbcAdcEL","BBC East Large ADC",100,0.5,400.5,8);
01942   m_bbc_adc[2] = QAH::MH1F("QaBbcAdcWS","BBC West Small ADC",100,0.5,400.5,8);
01943   m_bbc_adc[3] = QAH::MH1F("QaBbcAdcWL","BBC West Large ADC",100,0.5,400.5,8);
01944   m_bbc_tdc[0] = QAH::MH1F("QaBbcTdcES","BBC East Small TDC",100,0.5,750.5,8);
01945   m_bbc_tdc[1] = QAH::MH1F("QaBbcTdcEL","BBC East Large TDC",100,0.5,2000.5,8);
01946   m_bbc_tdc[2] = QAH::MH1F("QaBbcTdcWS","BBC West Small TDC",100,0.5,750.5,8);
01947   m_bbc_tdc[3] = QAH::MH1F("QaBbcTdcWL","BBC West Large TDC",100,0.5,2000.5,8);
01948   for (i=0; i<8; i++) {
01949     sprintf(ID,"%d",i+1);
01950     for (j=0; j<4; j++) {
01951       m_bbc_adc[j]->Rebin(i,ID);
01952       m_bbc_tdc[j]->Rebin(i,ID);
01953     } 
01954   } 
01955 
01956 }
01957 //_____________________________________________________________________________
01958 void StQABookHist::BookHistFPD(){
01959 
01960   QAH::MMH1F(m_fpd_top,2,"QaFpdTop%d","FPD Top ADC %d-%d",100,0.5,1500.5,8,1);
01961   QAH::MMH1F(m_fpd_bottom,2,"QaFpdBottom%d","FPD Bottom ADC %d-%d",100,0.5,1500.5,8,1);
01962   QAH::MMH1F(m_fpd_south,2,"QaFpdSouth%d","FPD South ADC %d-%d",100,0.5,1500.5,8,1);
01963   QAH::MMH1F(m_fpd_north,2,"QaFpdNorth%d","FPD North ADC %d-%d",100,0.5,1500.5,6,1);
01964   
01965 // Book ADC histograms: FPD SUM signals
01966   m_fpd_sums[0] = QAH::H1F("QaFpdSums0","FPD SUM Top",100,0.5,2050.5);
01967   m_fpd_sums[1] = QAH::H1F("QaFpdSums1","FPD SUM Bottom",100,0.5,2050.5);
01968   m_fpd_sums[2] = QAH::H1F("QaFpdSums2","FPD SUM South",100,0.5,2050.5);
01969   m_fpd_sums[3] = QAH::H1F("QaFpdSums3","FPD SUM North",100,0.5,2050.5);
01970   m_fpd_sums[4] = QAH::H1F("QaFpdSums4","FPD SUM SmdX",100,0.5,6000.5);
01971   m_fpd_sums[5] = QAH::H1F("QaFpdSums5","FPD SUM SmdY",100,0.5,6000.5);
01972   m_fpd_sums[6] = QAH::H1F("QaFpdSums6","FPD SUM Pres1",100,0.5,1500.5);
01973   m_fpd_sums[7] = QAH::H1F("QaFpdSums7","FPD SUM Pres2",100,0.5,1500.5);
01974 
01975 }
01976 //_____________________________________________________________________________
01977 void StQABookHist::BookHistPMD(){
01978 
01979   QAH::MMH1F(m_pmd_sm_hit,12,"QaPmdSmHit%d",
01980              "PMD SM-wise Hit Multiplicity %02d-%02d",1000,0.,500.,2);
01981   QAH::MMH1F(m_pmd_sm_adc,12,"QaPmdSmAdcHit%d",
01982              "PMD SM-wise ADC/Hit Multiplicity %02d-%02d",100,0.,600.,2);
01983   QAH::MMH1F(m_pmd_chain_hit,24,"QaPmdChHitChain%d",
01984              "PMD Channel-wise Hit, Chain %02d-%02d",1728,-0.5,1727.5,2);
01985   QAH::MMH1F(m_pmd_chain_adc,24,"QaPmdChAdcChain%d",
01986              "PMD Channel-wise ADC, Chain %02d-%02d",1728,-0.5,1727.5,2);
01987 
01988   m_pmd_total_hit = QAH::H2F("QaPmdTotalHit","PMD Total Hits",100,0.,2e5,100,0.,4.);
01989   m_pmd_total_hit->SetXTitle("event id");
01990   m_pmd_total_hit->SetYTitle("log10");
01991   m_pmd_total_adc = QAH::H2F("QaPmdTotalAdc","PMD Total ADC",100,0.,2e5,100,0.,6.);
01992   m_pmd_total_adc->SetXTitle("event id");
01993   m_pmd_total_adc->SetYTitle("log10");
01994   m_cpv_total_hit = QAH::H2F("QaCpvTotalHit","CPV Total Hits",100,0.,2e5,100,0.,4.);
01995   m_cpv_total_hit->SetXTitle("event id");
01996   m_cpv_total_hit->SetYTitle("log10");
01997   m_cpv_total_adc = QAH::H2F("QaCpvTotalAdc","CPV Total ADC",100,0.,2e5,100,0.,6.);
01998   m_cpv_total_adc->SetXTitle("event id");
01999   m_cpv_total_adc->SetYTitle("log10");
02000   
02001 }
02002 //_____________________________________________________________________________
02003 void StQABookHist::BookHistTOF(){
02004 
02005   m_tof_hit_tray = QAH::H2F("QaTofHitvsTray","TOF Hits vs tray",120,0.5,120.5,100,0.,100.);
02006   m_tof_hit_module = QAH::H2F("QaTofHitvsModule","TOF Hits vs Module",65,-32.5,32.5,100,0.,100.);
02007   m_tof_match_tray = QAH::H2F("QaTofMatchvsTray","TOF Matched Hits vs tray",120,0.5,120.5,100,0.,100.);
02008   m_tof_match_module = QAH::H2F("QaTofMatchvsModule","TOF Matched Hits vs Module",65,-32.5,32.5,100,0.,100.);
02009   m_tof_vpd_hit =  QAH::H2F("QaTofHitvsVpdHit","TOF Hits vs Vpd Hits",50,0.,50.,100,0.,5000.);
02010   m_tof_vtx_z =  QAH::H2F("QaTofVpdZvsTpcZ","VPD vtxz vs TPC vtxz",100,-100.,100.,100,-100.,100.);
02011   m_tof_PID =  QAH::H2F("QaTofPID","TOF InvBeta vs p",100,0.,5.,100,0.,4.);
02012 
02013 }
02014 //_____________________________________________________________________________
02015