256 #include "StQABookHist.h"
257 #include "StQAMakerBase.h"
258 #include "StEmcUtil/others/StEmcMath.h"
259 #include "StEmcUtil/others/emcDetectorName.h"
268 histsSet = StQA_Undef;
327 cmaxdedx = 0.1e-04*1e6;
332 mNullPrimVtxClass = 0;
335 m_globtrk_fit_prob=0;
478 m_globtrk_xf_yfTE = 0;
479 m_globtrk_xf_yfTW = 0;
480 m_globtrk_xf_yfFE = 0;
481 m_globtrk_xf_yfFW = 0;
482 m_globtrk_padtimeFE = 0;
483 m_globtrk_padtimeFW = 0;
485 m_mom_trklengthT = 0;
486 m_eta_trklengthT = 0;
487 m_eta_trklengthFE = 0;
488 m_eta_trklengthFW = 0;
489 m_npoint_lengthT = 0;
490 m_npoint_lengthFE = 0;
491 m_npoint_lengthFW = 0;
492 m_fpoint_lengthT = 0;
493 m_fpoint_lengthTTS = 0;
502 m_globtrk_xf_yfTS = 0;
504 m_mom_trklengthTS = 0;
505 m_eta_trklengthTS = 0;
506 m_npoint_lengthTS = 0;
507 m_fpoint_lengthTS = 0;
543 m_primtrk_meanptTTS=0;
545 m_primtrk_meanetaTTS=0;
546 m_primtrk_meanetaF=0;
660 m_primtrk_xf_yfTE = 0;
661 m_primtrk_xf_yfTW = 0;
662 m_primtrk_xf_yfFE = 0;
663 m_primtrk_xf_yfFW = 0;
665 m_pmom_trklengthT = 0;
666 m_peta_trklengthT = 0;
667 m_peta_trklengthFE = 0;
668 m_peta_trklengthFW = 0;
669 m_pnpoint_lengthT = 0;
670 m_pnpoint_lengthFE = 0;
671 m_pnpoint_lengthFW = 0;
672 m_pfpoint_lengthT = 0;
673 m_pfpoint_lengthTTS = 0;
681 m_primtrk_xf_yfTS = 0;
683 m_pmom_trklengthTS = 0;
684 m_peta_trklengthTS = 0;
685 m_pnpoint_lengthTS = 0;
686 m_pfpoint_lengthTS = 0;
725 m_vtx_FtpcEastTpc_xy=0;
726 m_vtx_FtpcEastTpc_z=0;
727 m_vtx_FtpcWestTpc_xy=0;
728 m_vtx_FtpcWestTpc_z=0;
759 m_geant_reco_pvtx_x=0;
760 m_geant_reco_pvtx_y=0;
761 m_geant_reco_pvtx_z=0;
762 m_geant_reco_vtx_z_z=0;
767 for (i=0; i<4; i++) {
775 for (i=0; i<8; i++) {
776 m_emc_hits_per_module[i]=0;
777 m_emc_energy_per_module[i]=0;
778 m_emc_strip_hits_per_module[i]=0;
779 m_emc_strip_energy_per_module[i]=0;
787 for (i=0; i<4; i++) {
797 for (i=0; i<4; i++) {
798 m_emc_point_energy[i]=0;
799 m_emc_point_eta[i]=0;
800 m_emc_point_phi[i]=0;
801 m_emc_point_sigeta[i]=0;
802 m_emc_point_sigphi[i]=0;
803 m_emc_point_deleta[i]=0;
804 m_emc_point_delphi[i]=0;
805 m_emc_point_trmom[i]=0;
811 for (i=0; i<4; i++) {
817 for (i=0; i<2; i++) {
823 for (i=0; i<8; i++) m_fpd_sums[i] = 0;
826 for (i=0; i<12; i++) {
830 for (i=0; i<24; i++) {
831 m_pmd_chain_adc[i] = 0;
832 m_pmd_chain_hit[i] = 0;
842 m_pnt_svtLaserDiff=0;
868 m_global_hft_hit = 0;
869 m_primary_hft_hit = 0;
872 m_pxl_hit_phi_z_Pxl1 = 0;
873 m_pxl_hit_phi_z_Pxl2 = 0;
874 m_pxl_hit_ladder = 0;
875 m_pxl_hit_sector_sensor_Pxl1 = 0;
876 m_pxl_hit_sector_sensor_Pxl2 = 0;
877 m_pxl_nhit_Pxl1_tpc_mult = 0;
878 m_pxl_nhit_Pxl2_tpc_mult = 0;
879 m_pxl_nhit_Pxl1_tof_mult = 0;
880 m_pxl_nhit_Pxl2_tof_mult = 0;
881 m_pxl_nhit_Pxl1_Pxl2 = 0;
882 m_global_pxl_hit = 0;
883 m_primary_pxl_hit = 0;
887 m_ist_hit_ladder = 0;
888 m_ist_hit_ladder_sensor = 0;
889 m_ist_nhit_tpc_mult = 0;
890 m_ist_nhit_tof_mult = 0;
891 m_global_ist_hit = 0;
892 m_primary_ist_hit = 0;
898 m_tof_match_module=0;
907 m_MtdMatchHitMap = 0;
910 for (
int i=0; i<24; i++) {
917 for(
int i = 0; i < 24; i++) {
918 m_TPC_ch_time_inner[i] = 0;
919 m_TPC_ch_time_outer[i] = 0;
920 m_TPC_ch_nrow[i] = 0;
922 m_TPC_adc_sec_inner = 0;
923 m_TPC_adc_sec_outer = 0;
928 void StQABookHist::BookHist(Int_t hSet){
932 QAH::preString = QAHistType;
934 if (histsSet != StQA_MC) {
935 mNullPrimVtxClass = QAH::H1F(
"QaNullPrimVtxMult",
"event primary vertex check",40,-2,2);
936 mNullPrimVtxClass->SetXTitle(
"has primary vertex? (yes = 1, no = -1)");
937 mNullPrimVtxClass->SetYTitle(
"# of events");
939 if (histsSet<StQA_run8) silHists = kTRUE;
940 if (histsSet<StQA_run12all) ftpHists = kTRUE;
943 if (histsSet == StQA_AuAuOld) {
952 if (histsSet<StQA_run12all) BookHistPMD();
953 if (histsSet==StQA_MC) BookHistEval();
954 if (histsSet>=StQA_run8) BookHistTOF();
955 if (histsSet>=StQA_run12all) BookHistMTD();
956 if (histsSet>=StQA_run14 && histsSet<StQA_run17) {
961 if (histsSet>=StQA_run18) {
969 void StQABookHist::BookHistGlob(){
975 m_globtrk_tot = QAH::H1F(
"QaGtrkTot",
"globtrk: tot num tracks - all",
977 m_globtrk_iflag = QAH::H1F(
"QaGtrkFlag",
"globtrk: iflag - all ",200,-999.,1001.);
978 m_globtrk_good = QAH::H1F(
"QaGtrkGood",
"globtrk: tot good tracks - all",40,0.,10000.);
979 m_globtrk_good_sm = QAH::H1F(
"QaGtrkGoodsm",
"globtrk: tot good tracks - all",40,0.,1000.);
981 m_globtrk_goodTTS = QAH::H1F(
"QaGtrkGoodTTS",
"globtrk: tot good tracks - tpc,tpc+svt",150,0.,9000.);
982 m_globtrk_goodF = QAH::H2F(
"QaGtrkGoodF",
"globtrk: tot good tracks - ftpc",150,0.,1500.,150,0.,1500.);
985 m_globtrk_fit_prob = QAH::H1F(
"QaGtrkFitProb",
"globtrk: prob. fit is correct",100,0,1.2);
986 m_det_id = QAH::H1F(
"QaGtrkDetId",
"globtrk: Detector ID good tracks - all",48,-0.5,47.5);
987 m_dcaToBeamXY = QAH::H2F(
"QaGtrkDcaBeamXY",
"globtrk: xy-DCA to Beam Axis (z=0)",80,-4,4,80,-4,4);
990 m_dcaToBeamZ1 = QAH::H1F(
"QaGtrkDcaBeamZ1",
"globtrk: z-DCA to Beam Axis, coarse scale",100,-200,200);
991 m_dcaToBeamZ2 = QAH::H1F(
"QaGtrkDcaBeamZ2",
"globtrk: z-DCA to Beam Axis",100,-50,50);
992 m_dcaToBeamZ3 = QAH::H1F(
"QaGtrkDcaBeamZ3",
"globtrk: z-DCA to Beam Axis, near z=0",80,-20,20);
993 m_zDcaTanl = QAH::H2F(
"QaGtrkZdcaTanl",
"globtrk: z-DCA to Beam Axis vs tanl, tpc",100,-25,25,32,-4,4);
996 m_zDcaZf = QAH::H2F(
"QaGtrkZdcaZf",
"globtrk: z-DCA to Beam Axis vs z-first",100,-25,25,50,-300,300);
998 m_zDcaZf->SetYTitle(
"z of first point");
999 m_zDcaPsi = QAH::H2F(
"QaGtrkZdcaPsi",
"globtrk: z-DCA to Beam Axis vs psi",100,-25,25,64,0,360);
1002 m_zDcaPhi0 = QAH::H2F(
"QaGtrkZdcaPhi0",
"globtrk: z-DCA to Beam Axis vs azimuth (phi0) at start",80,-20,20,64,0,360);
1006 m_glb_sptsTS = QAH::H1F(
"QaGtrkSptsTS",
"globtrk: number of svt points, tpc+svt",4,0.,4.);
1007 m_glb_ratiomTTS = QAH::MH1F(
"QaGtrkRnmfTTS",
"globtrk: ratio Nfit/max pnt, tpc,svt", 55,0.,1.1,2);
1012 m_psiTTS = QAH::MH1F(
"QaGtrkPsiTTS",
"globtrk: psi (deg), tpc,svt", 64, 0.,360.,2);
1016 m_etaTTS = QAH::MH1F(
"QaGtrkEtaTTS",
"globtrk: eta, tpc,svt",40,-2.,2.,2);
1020 m_pTTTS = QAH::MH1F(
"QaGtrkPtTTS",
"globtrk: log10 pT, tpc,svt",50,-1.,4.,2);
1024 m_pTTTS->SetXTitle(
"log10 pT (MeV)");
1025 m_chisq0TTS = QAH::MH1F(
"QaGtrkChisq0TTS",
"globtrk: chisq0, tpc,svt", 50, 0.,5.,2);
1029 m_fit_pointTTS = QAH::MH1F(
"QaGtrkNPntFitTTS",
"globtrk: N fit pnts on trk, tpc,svt",50,0.,50.,2);
1033 m_glb_impactTTS = QAH::MH1F(
"QaGtrkImpactTTS",
"globtrk: log10 impact param from prim vtx, tpc,svt",120,-3.,3.,2);
1037 m_glb_impactrTTS = QAH::MH1F(
"QaGtrkImpactrTTS",
"globtrk: impact param from prim vtx, tpc,svt",100,0.,300.,2);
1041 m_fpoint_lengthTTS = QAH::H2F(
"QaGtrkFitPntLTTS",
"globtrk: N fit pnts vs length, tpc,tpc+svt",
1042 25,0.,250.,25,0.,50.);
1049 m_pointT = QAH::H1F(
"QaGtrkNPntT",
"globtrk: N points on trk,tpc", 81,-0.5, 80.5);
1050 m_max_pointT = QAH::H1F(
"QaGtrkNPntMaxT",
"globtrk: N max pnts on trk, tpc", 76, -1, 151);
1051 m_fit_pointT = QAH::H1F(
"QaGtrkNPntFitT",
"globtrk: N fit pnts on trk, tpc", 81, -0.5, 80.5);
1052 m_glb_ratiomT = QAH::H1F(
"QaGtrkRnmT",
"globtrk: ratio Nfit/max pnt, tpc", 55, 0., 1.1);
1053 m_glb_chargeT = QAH::H1F(
"QaGtrkChrgT",
"globtrk: charge, tpc ", 20,-2.,2.);
1054 m_glb_r0T = QAH::H1F(
"QaGtrkR0T",
"globtrk: radius at start (cm), tpc ", 50,0.,200.);
1056 m_glb_phi0T = QAH::H1F(
"QaGtrkPhi0T",
"globtrk: azimuth (phi) at start (deg,force 0,360), tpc ", 64,0.,360.);
1057 m_glb_z0T = QAH::H1F(
"QaGtrkZ0T",
"globtrk: z-coord at start (cm), tpc ", 50, -300.,300.);
1059 m_glb_curvT = QAH::H1F(
"QaGtrkCurvT",
"globtrk: log10 curvature (1/cm), tpc ", 80,-3.5,0.5);
1060 m_glb_padfTEW = QAH::MH1F(
"QaGtrkPadfTEW",
"globtrk: padrow of first hit on trk, tpc", 72,0.5,72.5,2);
1064 m_glb_padfT = QAH::H1F(
"QaGtrkPadfT",
"globtrk: padrow of first hit on trk, tpc", 72,0.5,72.5);
1065 m_glb_f0 = QAH::MH1F(
"QaGtrkf0",
"globtrk: first point: hit-helix, tpc",60,-3.,3.,3);
1070 m_glb_xfT = QAH::H1F(
"QaGtrkXfT",
"globtrk: x of first hit on trk, tpc", 50,-200.,200.);
1071 m_glb_xf0 = QAH::H1F(
"QaGtrkXf0",
"globtrk: first point: x_hit - x_helix, tpc",60,-3.,3.);
1072 m_glb_yfT = QAH::H1F(
"QaGtrkYfT",
"globtrk: y of first hit on trk, tpc", 50,-200.,200.);
1073 m_glb_yf0 = QAH::H1F(
"QaGtrkYf0",
"globtrk: first point: y_hit - y_helix, tpc",60,-3.,3.);
1074 m_glb_zfT = QAH::H1F(
"QaGtrkZfT",
"globtrk: z of first hit on trk, tpc", 50,-300.,300.);
1075 m_glb_zf0 = QAH::H1F(
"QaGtrkZf0",
"globtrk: first point: z_hit - z_helix, tpc",60,-3.,3.);
1076 m_glb_radfT = QAH::H1F(
"QaGtrkRT",
"globtrk: radial position of first hit, tpc", 50,0.,200.);
1077 m_glb_rzf0 = QAH::MH1F(
"QaGtrkRZf0",
"globtrk: first point: hit - helix, tpc",60,-3.,3.,2);
1078 m_glb_rzf0->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1081 m_glb_rzl0 = QAH::MH1F(
"QaGtrkRZl0",
"globtrk: last point: hit - helix, tpc",60,-3.,3.,2);
1082 m_glb_rzl0->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1085 m_glb_phifT = QAH::MH1F(
"QaGtrkPhifT",
"globtrk: phi of first point on trk, tpc",64,0,360,2);
1089 m_lengthT = QAH::H1F(
"QaGtrkLengthT",
"globtrk: track length, tpc", 50,0.,300.);
1090 m_psiT = QAH::H1F(
"QaGtrkPsiT",
"globtrk: psi, tpc (deg)", 64, 0.,360.);
1091 m_tanlT = QAH::H1F(
"QaGtrkTanlT",
"globtrk: tanl, tpc",32,-4.,4.);
1092 m_glb_thetaT = QAH::H1F(
"QaGtrkThetaT",
"globtrk: theta (degrees), tpc",36,0.,180.);
1093 m_etaT = QAH::H1F(
"QaGtrkEtaT",
"globtrk: eta, tpc",40,-2.,2.);
1094 m_pTT = QAH::H1F(
"QaGtrkPtT",
"globtrk: pT, tpc",50,0.,10.);
1095 m_momT = QAH::H1F(
"QaGtrkPT",
"globtrk: momentum, tpc",50,0.,10.);
1096 m_chisq0T = QAH::H1F(
"QaGtrkChisq0T",
"globtrk: chisq0, tpc", 50, 0.,5.);
1097 m_glb_impactT = QAH::MH1F(
"QaGtrkImpactT",
"globtrk: log10 impact param from prim vtx, tpc",
1103 m_glb_simpactT = QAH::MH1F(
"QaGtrkSImpactT",
"globtrk: signed impact param from prim vtx, tpc",
1109 m_glb_impactrT = QAH::H1F(
"QaGtrkImpactrT",
"globtrk: impact param from prim vtx, tpc",100,0.,300.);
1114 m_pT_eta_recT = QAH::H2F(
"QaGtrkPtVsEtaT",
"globtrk: log10 pT vs eta, tpc", 20,-2.,2.,40,1.,4.);
1118 m_globtrk_xf_yfTE = QAH::H2F(
"QaGtrkXfYfTE",
"globtrk: Y vs X of first hit on trk, tpc east", 40,-200.,200.,40,-200.,200.);
1122 m_globtrk_xf_yfTW = QAH::H2F(
"QaGtrkXfYfTW",
"globtrk: Y vs X of first hit on trk, tpc west", 40,-200.,200.,40,-200.,200.);
1126 m_tanl_zfT = QAH::H2F(
"QaGtrkTanlzf",
"globtrk: tanl(dip) vs. (zfirst-zvtx)/arc length, tpc,tpc+svt",60,-3.,3.,60,-3.,3.);
1127 m_tanl_zfT->SetXTitle(
"(zfirst-zvtx)/arc length");
1130 m_mom_trklengthT = QAH::H2F(
"QaGtrkPVsTrkLength",
"globtrk: log10 mom vs trk length, tpc",
1131 50,0.,250.,40,1.,4.);
1135 m_eta_trklengthT = QAH::H2F(
"QaGtrkLengthVEtaT",
"globtrk: trk length vs eta, tpc",
1136 20,-2.,2.,50,0.,250.);
1140 m_npoint_lengthT = QAH::H2F(
"QaGtrkNPntLengthT",
"globtrk: N pnts vs length, tpc",
1141 25,0.,250.,40,1.,81.);
1145 m_fpoint_lengthT = QAH::H2F(
"QaGtrkFitPntLengthT",
"globtrk: N fit pnts vs length, tpc",
1146 25,0.,250.,40,1.,81.);
1150 m_chisq0_momT = QAH::H2F(
"QaGtrkChi0MomT",
"globtrk: Chisq0 vs log10 mom, tpc",40,1.,4.,50,0.,5.);
1154 m_chisq0_etaT = QAH::H2F(
"QaGtrkChi0EtaT",
"globtrk: Chisq0 vs eta, tpc",20,-2.,2.,20,0.,5.);
1158 m_chisq0_dipT = QAH::H2F(
"QaGtrkChi0TanlT",
"globtrk: Chisq0 vs tanl(dip), tpc",20,-3.2,3.2,20,0.,5.);
1162 m_chisq0_zfT = QAH::H2F(
"QaGtrkChi0zfT",
"globtrk: Chisq0 vs zfirst, tpc",20,-250.,250.,20,0.,5.);
1166 m_chisq0_phiT = QAH::H2F(
"QaGtrkChi0PhiT",
"globtrk: Chisq0 vs phi, tpc",20,0.,360,20,0.,5.);
1170 m_psi_phiT = QAH::H2F(
"QaGtrkPsiPhiT",
"globtrk: psi vs phi, tpc",36, 0.,360.,36,0.,360.);
1178 m_pointTS = QAH::H1F(
"QaGtrkNPntTS",
"globtrk: N points on trk,tpc+svt", 60, 0.,60.);
1179 m_max_pointTS = QAH::H1F(
"QaGtrkNPntMaxTS",
"globtrk: N max pnts on trk, tpc+svt", 50, 0.,100.);
1180 m_fit_pointTS = QAH::H1F(
"QaGtrkNPntFitTS",
"globtrk: N fit pnts on trk, tpc+svt", 60, 0.,60.);
1181 m_glb_ratiomTS = QAH::H1F(
"QaGtrkRnmTS",
"globtrk: ratio Nfit/max pnt, tpc+svt", 55, 0., 1.1);
1182 m_glb_chargeTS = QAH::H1F(
"QaGtrkChrgTS",
"globtrk: charge, tpc+svt ", 20,-2.,2.);
1183 m_glb_r0TS = QAH::H1F(
"QaGtrkR0TS",
"globtrk: radius at start (cm), tpc+svt", 100,0.,25.);
1185 m_glb_phi0TS = QAH::H1F(
"QaGtrkPhi0TS",
"globtrk: azimuth (phi) at start (deg,force 0-360),tpc+svt", 64, 0.,360.);
1186 m_glb_z0TS = QAH::H1F(
"QaGtrkZ0TS",
"globtrk: z-coord at start (cm), tpc+svt", 50, -50.,50.);
1188 m_glb_curvTS = QAH::H1F(
"QaGtrkCurvTS",
"globtrk: log10 curvature (1/cm), tpc+svt", 80,-3.5,0.5);
1189 m_glb_f0TS = QAH::MH1F(
"QaGtrkf0TS",
"globtrk: first point: hit-helix, tpc+svt",50,-5.,5.,3);
1194 m_glb_xfTS = QAH::H1F(
"QaGtrkXfTS",
"globtrk: x of first hit on trk, tpc+svt", 50,-200.,200.);
1195 m_glb_xf0TS = QAH::H1F(
"QaGtrkXf0TS",
"globtrk: x of first hit - on helix at start, tpc+svt",50,-5.,5.);
1196 m_glb_yfTS = QAH::H1F(
"QaGtrkYfTS",
"globtrk: y of first hit on trk, tpc+svt", 50,-200.,200.);
1197 m_glb_yf0TS = QAH::H1F(
"QaGtrkYf0TS",
"globtrk: y of first hit - on helix at start, tpc+svt",50,-5.,5.);
1198 m_glb_zfTS = QAH::H1F(
"QaGtrkZfTS",
"globtrk: z of first hit on trk, svt", 50,-30.,30.);
1199 m_glb_zf0TS = QAH::H1F(
"QaGtrkZf0TS",
"globtrk: z of first hit - on helix at start, tpc+svt",50,-5.,5.);
1200 m_glb_radfTS = QAH::H1F(
"QaGtrkRTS",
"globtrk: radial position of first hit, tpc+svt",80,4.,200.);
1201 m_glb_rzf0TS = QAH::MH1F(
"QaGtrkRZf0TS",
"globtrk: first point: hit - helix, tpc+svt",60,-3.,3.,2);
1202 m_glb_rzf0TS->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1205 m_glb_rzl0TS = QAH::MH1F(
"QaGtrkRZl0TS",
"globtrk: last point: hit - helix, tpc+svt",60,-3.,3.,2);
1206 m_glb_rzl0TS->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1209 m_glb_phifTS = QAH::H1F(
"QaGtrkPhifTS",
"globtrk: phi of first point on track, svt",64,0,360);
1210 m_glb_ssd_phi = QAH::H1F(
"QaGtrkPhifSSD",
"globtrk: phi of ssd point (per event)",64,0,360);
1211 m_lengthTS = QAH::H1F(
"QaGtrkLengthTS",
"globtrk: track length, tpc+svt", 50,0.,300.);
1212 m_psiTS = QAH::H1F(
"QaGtrkPsiTS",
"globtrk: psi, tpc+svt (deg) ", 64, 0.,360.);
1213 m_tanlTS = QAH::H1F(
"QaGtrkTanlTS",
"globtrk: tanl, tpc+svt",32,-4.,4.);
1214 m_glb_thetaTS = QAH::H1F(
"QaGtrkThetaTS",
"globtrk: theta (degrees), tpc+svt",36,0.,180.);
1215 m_etaTS = QAH::H1F(
"QaGtrkEtaTS",
"globtrk: eta, tpc+svt",40,-2.,2.);
1216 m_pTTS = QAH::H1F(
"QaGtrkPtTS",
"globtrk: pT, tpc+svt",50,0.,10.);
1217 m_momTS = QAH::H1F(
"QaGtrkPTS",
"globtrk: momentum, tpc+svt",50,0.,10.);
1218 m_chisq0TS = QAH::H1F(
"QaGtrkChisq0TS",
"globtrk: chisq0, tpc+svt", 50, 0.,5.);
1219 m_glb_impactTS = QAH::MH1F(
"QaGtrkImpactTS",
"globtrk: log10 impact param from prim vtx, tpc+svt",
1225 m_glb_simpactTS = QAH::MH1F(
"QaGtrkSImpactTS",
"globtrk: signed impact param from prim vtx, tpc+svt",
1231 m_glb_impactrTS = QAH::H1F(
"QaGtrkImpactrTS",
"globtrk: impact param from prim vtx, tpc+svt",
1237 m_pT_eta_recTS = QAH::H2F(
"QaGtrkPtVsEtaTS",
"globtrk: log10 pT vs eta, tpc+svt", 20,-2.,2.,40,1.,4.);
1241 m_globtrk_xf_yfTS = QAH::H2F(
"QaGtrkXfYfTS",
"globtrk: Y vs X of first hit on trk, tpc+svt", 40,-200.,200.,40,-200.,200.);
1245 m_tanl_zfTS = QAH::H2F(
"QaGtrkTanlzfTS",
"globtrk: tanl(dip) versus (zfirst-zvtx)/arc length, svt",60,-3.,3.,60,-3.,3.);
1246 m_tanl_zfTS->SetXTitle(
"(zfirst-zvtx)/arc length");
1249 m_mom_trklengthTS = QAH::H2F(
"QaGtrkPVsTrkLTS",
"globtrk: log10 mom vs trk length, tpc+svt",50,0.,250.,40,1.,4.);
1253 m_eta_trklengthTS = QAH::H2F(
"QaGtrkLVEtaTS",
"globtrk: trk length vs eta, tpc+svt",20,-2.,2.,50,0.,250.);
1257 m_npoint_lengthTS = QAH::H2F(
"QaGtrkNPntLTS",
"globtrk: N pnts vs length, tpc+svt",25,0.,250.,25,0.,50.);
1261 m_fpoint_lengthTS = QAH::H2F(
"QaGtrkFitPntLTS",
"globtrk: N fit pnts vs length, tpc+svt",25,0.,250.,25,0.,50.);
1265 m_chisq0_momTS = QAH::H2F(
"QaGtrkChi0MomTS",
"globtrk: Chisq0 vs log10 mom, tpc+svt",40,1.,4.,50,0.,5.);
1269 m_chisq0_etaTS = QAH::H2F(
"QaGtrkChi0EtaTS",
"globtrk: Chisq0 vs eta, tpc+svt",20,-2.,2.,20,0.,5.);
1273 m_chisq0_dipTS = QAH::H2F(
"QaGtrkChi0TanlTS",
"globtrk: Chisq0 vs tanl(dip), tpc+svt",20,-3.2,3.2,20,0.,5.);
1277 m_chisq0_zfTS = QAH::H2F(
"QaGtrkChi0zfTS",
"globtrk: Chisq0 vs zfirst, tpc+svt",20,-250.,250.,20,0.,5.);
1281 m_chisq0_phiTS = QAH::H2F(
"QaGtrkChi0PhiTS",
"globtrk: Chisq0 vs phi, tpc+svt",20,0.,360,20,0.,5.);
1285 m_psi_phiTS = QAH::H2F(
"QaGtrkPsiPhiTS",
"globtrk: psi vs phi, tpc+svt",36, 0.,360.,36,0.,360.);
1296 m_pointF = QAH::MH1F(
"QaGtrkNPntF",
"globtrk: N points on trk,ftpc", 8,4.,12.,2);
1300 m_max_pointF = QAH::MH1F(
"QaGtrkNPntMaxF",
"globtrk: N max pnts on trk, ftpc", 8,4.,12.,2);
1304 m_glb_ratiomF = QAH::MH1F(
"QaGtrkRnmF",
"globtrk: ratio Nfit/max pnt, ftpc", 11,0.,1.1,2);
1308 m_glb_chargeF = QAH::MH1F(
"QaGtrkChrgF",
"globtrk: charge, ftpc", 20,-2.,2.,2);
1312 m_glb_xfF = QAH::MH1F(
"QaGtrkXfF",
"globtrk: x of first hit on trk, ftpc",70,-35.,35.,2);
1316 m_glb_yfF = QAH::MH1F(
"QaGtrkYfF",
"globtrk: y of first hit on trk, ftpc",70,-35.,35.,2);
1320 m_glb_zfF = QAH::MH1F(
"QaGtrkZfF",
"globtrk: z of first hit on trk, ftpc",75,-225.,225.,2);
1324 m_glb_planefF = QAH::MH1F(
"QaGtrkPlanefF",
"globtrk: plane of first hit on trk, ftpc",20,0.5,20.5,2);
1328 m_glb_radfF = QAH::MH1F(
"QaGtrkRF",
"globtrk: radial position of first hit, ftpc",35,0.,35.,2);
1332 m_lengthF = QAH::MH1F(
"QaGtrkLengthF",
"globtrk: track length, ftpc",60,0.,120.,2);
1336 m_psiF = QAH::MH1F(
"QaGtrkPsiF",
"globtrk: psi, ftpc",90,0.,360.,2);
1339 m_psiF->SetStats(kFALSE);
1340 m_etaF = QAH::MH1F(
"QaGtrkEtaF",
"globtrk: |eta|, ftpc",80,2,5,2);
1343 m_etaF->SetStats(kFALSE);
1344 m_pTF = QAH::MH1F(
"QaGtrkPtF",
"globtrk: pT, ftpc",50,0.,10.,2);
1345 m_pTF->Rebin(0,
"East");
1346 m_pTF->Rebin(1,
"West");
1347 m_pTF->SetStats(kFALSE);
1348 m_momF = QAH::MH1F(
"QaGtrkPF",
"globtrk: momentum, ftpc",50,0.,10.,2);
1351 m_momF->SetStats(kFALSE);
1352 m_chisq0F = QAH::MH1F(
"QaGtrkChisq0F",
"globtrk: chi2/ndf x,y fit, ftpc",50,0.,5.,2);
1356 m_chisq1F = QAH::MH1F(
"QaGtrkChisq1F",
"globtrk: chi2/ndf r,z fit, ftpc",50,0.,5.,2);
1360 m_glb_impactF = QAH::MH1F(
"QaGtrkImpactF",
"globtrk: log10 impact param from prim vtx, ftpc",120,-3.0,3.,2);
1364 m_glb_impactrF = QAH::MH1F(
"QaGtrkImpactrF",
"globtrk: impact param from prim vtx, ftpc",100,0.,10.,2);
1369 m_pointFE = QAH::H1F(
"QaGtrkNPntFE",
"globtrk: N points on trk,ftpc east", 8, 4.,12.);
1370 m_pointFW = QAH::H1F(
"QaGtrkNPntFW",
"globtrk: N points on trk,ftpc west", 8, 4.,12.);
1371 m_max_pointFE = QAH::H1F(
"QaGtrkNPntMaxFE",
"globtrk: N max pnts on trk, ftpc east", 8, 4.,12.);
1372 m_max_pointFW = QAH::H1F(
"QaGtrkNPntMaxFW",
"globtrk: N max pnts on trk, ftpc west", 8, 4.,12.);
1373 m_glb_ratiomFE = QAH::H1F(
"QaGtrkRnmFE",
"globtrk: ratio Nfit/max pnt, ftpc east", 55, 0., 1.1);
1374 m_glb_ratiomFW = QAH::H1F(
"QaGtrkRnmFW",
"globtrk: ratio Nfit/max pnt, ftpc west", 55, 0., 1.1);
1375 m_glb_chargeFE = QAH::H1F(
"QaGtrkChrgFE",
"globtrk: charge, ftpc east ", 20,-2.,2.);
1376 m_glb_chargeFW = QAH::H1F(
"QaGtrkChrgFW",
"globtrk: charge, ftpc west ", 20,-2.,2.);
1377 m_glb_xfFE = QAH::H1F(
"QaGtrkXfFE",
"globtrk: x of first hit on trk, ftpc east", 70,-35.,35.);
1378 m_glb_xfFW = QAH::H1F(
"QaGtrkXfFW",
"globtrk: x of first hit on trk, ftpc west", 70,-35.,35.);
1379 m_glb_yfFE = QAH::H1F(
"QaGtrkYfFE",
"globtrk: y of first hit on trk, ftpc east", 70,-35.,35.);
1380 m_glb_yfFW = QAH::H1F(
"QaGtrkYfFW",
"globtrk: y of first hit on trk, ftpc west", 70,-35.,35.);
1381 m_glb_zfFE = QAH::H1F(
"QaGtrkZfFE",
"globtrk: z of first hit on trk, ftpc east", 75,-225.,-150.);
1382 m_glb_zfFW = QAH::H1F(
"QaGtrkZfFW",
"globtrk: z of first hit on trk, ftpc west", 75,150.,225.);
1383 m_glb_radfFE = QAH::H1F(
"QaGtrkRFE",
"globtrk: radial position of first hit, ftpc east",35,0.,35.);
1384 m_glb_radfFW = QAH::H1F(
"QaGtrkRFW",
"globtrk: radial position of first hit, ftpc west",35,0.,35.);
1385 m_lengthFE = QAH::H1F(
"QaGtrkLengthFE",
"globtrk: track length, ftpc east", 60,0.,120.);
1386 m_lengthFW = QAH::H1F(
"QaGtrkLengthFW",
"globtrk: track length, ftpc west", 60,0.,120.);
1387 m_psiFE = QAH::H1F(
"QaGtrkPsiFE",
"globtrk: psi, ftpc east", 90, 0.,360.);
1388 m_psiFW = QAH::H1F(
"QaGtrkPsiFW",
"globtrk: psi, ftpc west", 90, 0.,360.);
1389 m_etaFE = QAH::H1F(
"QaGtrkEtaFE",
"globtrk: eta, ftpc east",80,-4.5,-2.);
1390 m_etaFW = QAH::H1F(
"QaGtrkEtaFW",
"globtrk: eta, ftpc west",80,2.,4.5);
1391 m_pTFE = QAH::H1F(
"QaGtrkPtFE",
"globtrk: pT, ftpc east",50,0.,10.);
1392 m_pTFW = QAH::H1F(
"QaGtrkPtFW",
"globtrk: pT, ftpc west",50,0.,10.);
1393 m_momFE = QAH::H1F(
"QaGtrkPFE",
"globtrk: momentum, ftpc east ",50,0.,5.);
1394 m_momFW = QAH::H1F(
"QaGtrkPFW",
"globtrk: momentum, ftpc west ",50,0.,5.);
1395 m_chisq0FE = QAH::H1F(
"QaGtrkChisq0FE",
"globtrk: chi2/ndf x,y fit, ftpc east", 50, 0.,5.);
1396 m_chisq0FW = QAH::H1F(
"QaGtrkChisq0FW",
"globtrk: chi2/ndf x,y fit, ftpc west", 50, 0.,5.);
1397 m_chisq1FE = QAH::H1F(
"QaGtrkChisq1FE",
"globtrk: chi2/ndf r,z fit, ftpc east", 50, 0.,5.);
1398 m_chisq1FW = QAH::H1F(
"QaGtrkChisq1FW",
"globtrk: chi2/ndf r,z fit, ftpc west", 50, 0.,5.);
1402 m_pT_eta_recFE = QAH::H2F(
"QaGtrkPtVsEtaFE",
"globtrk: log10 pT vs eta, ftpcE",20,-4.5,-2.,40,1.,4.);
1405 m_pT_eta_recFW = QAH::H2F(
"QaGtrkPtVsEtaFW",
"globtrk: log10 pT vs eta, ftpcW",20,2.,4.5,40,1.,4.);
1408 m_globtrk_xf_yfFE = QAH::H2F(
"QaGtrkXfYfFE",
"globtrk: Y vs X of first hit on trk, ftpcE", 70,-35.,35.,70,-35.,35.);
1411 m_globtrk_xf_yfFW = QAH::H2F(
"QaGtrkXfYfFW",
"globtrk: Y vs X of first hit on trk, ftpcW", 70,-35.,35.,70,-35.,35.);
1414 m_globtrk_padtimeFE = QAH::H2F(
"QaGtrkPadTimeFtpcE",
"globtrk: #pads vs #timebins, ftpcE",12,0.5,12.5,10,0.5,10.5);
1418 m_globtrk_padtimeFW = QAH::H2F(
"QaGtrkPadTimeFtpcW",
"globtrk: #pads vs #timebins, ftpcW",12,0.5,12.5,10,0.5,10.5);
1422 m_eta_trklengthFE = QAH::H2F(
"QaGtrkLengthVEtaFE",
"globtrk: trk length vs eta, ftpcE",25,-4.5,-2.,30,0.,120.);
1425 m_eta_trklengthFW = QAH::H2F(
"QaGtrkLengthVEtaFW",
"globtrk: trk length vs eta, ftpcW",25,2.,4.5,30,0.,120.);
1428 m_npoint_lengthFE = QAH::H2F(
"QaGtrkNPntLengthFE",
"globtrk: N pnts vs length, ftpcE",30,0.,120.,15,0.,15.);
1431 m_npoint_lengthFW = QAH::H2F(
"QaGtrkNPntLengthFW",
"globtrk: N pnts vs length, ftpcW",30,0.,120.,15,0.,15.);
1437 m_global_pxl_hit = QAH::H1F(
"QaGtrkPxlHit",
"PIXEL: Hits per global track",5,-0.5,4.5);
1439 m_global_ist_hit = QAH::H1F(
"QaGtrkIstHit",
"IST: Hits per global track",5,-0.5,4.5);
1441 m_global_hft_hit = QAH::H1F(
"QaGtrkHftHit",
"HFT: Hits per global track",5,-0.5,4.5);
1452 void StQABookHist::BookHistPrim(){
1457 m_primtrk_tot = QAH::H1F(
"QaPtrkTot",
"primtrk: tot num tracks",50,0.,5000.);
1458 m_primtrk_tot_sm = QAH::H1F(
"QaPtrkTotsm",
"primtrk: tot num tracks",50,0.,200.);
1459 m_primtrk_iflag = QAH::H1F(
"QaPtrkFlag",
"primtrk: iflag - all",160,-799.,900.);
1460 m_primtrk_good = QAH::H1F(
"QaPtrkGood",
"primtrk: tot num tracks iflag>0",50,0.,5000.);
1461 m_primtrk_good_sm = QAH::H1F(
"QaPtrkGoodsm",
"primtrk: tot num tracks iflag>0",50,0.,200.);
1463 m_primtrk_goodTTS = QAH::H1F(
"QaPtrkGoodTTS",
"primtrk: tot num tracks iflag>0, tpc,svt",150,0.,4500.);
1464 m_primtrk_goodF = QAH::H2F(
"QaPtrkGoodF",
"primtrk: tot num tracks iflag>0, ftpc",150,0.,1500.,150,0.,1500.);
1467 m_primglob_good = QAH::H1F(
"QaPtrkGlob",
"primtrk: ratio primary/global tracks w/ iflag>0",50,0,1);
1468 m_primglob_fit = QAH::H1F(
"QaPtrkGlobFit",
"primtrk: ratio primary/global nfit points",50,0,2);
1469 m_pdet_id = QAH::H1F(
"QaPtrkDetId",
"primtrk: Detector ID good tracks - all",48,-0.5,47.5);
1471 m_primtrk_meanptTTS = QAH::MH1F(
"QaPtrkMeanPtTTS",
"primtrk: <pT>, tpc, tpc+svt",50,0.,2.,2);
1476 m_primtrk_meanptF = QAH::MH1F(
"QaPtrkMeanPtF",
"primtrk: <pT>, ftpc",50,0.,2.,2);
1481 m_primtrk_meanetaTTS = QAH::MH1F(
"QaPtrkMeanEtaTTS",
"primtrk: <eta>, tpc, tpc+svt",40,-2.,2.,2);
1491 m_ppsiTTS = QAH::MH1F(
"QaPtrkPsiTTS",
"primtrk: psi (deg), tpc, svt", 36, 0.,360.,2);
1495 m_petaTTS = QAH::MH1F(
"QaPtrkEtaTTS",
"primtrk: eta, tpc,svt",40,-2.,2.,2);
1499 m_ppTTTS = QAH::MH1F(
"QaPtrkPtTTS",
"primtrk: pT, tpc,svt",50,0.,5.,2);
1503 m_pchisq0TTS = QAH::MH1F(
"QaPtrkChisq0TTS",
"primtrk: chisq0, tpc,svt", 50, 0.,5.,2);
1507 m_pfpoint_lengthTTS = QAH::H2F(
"QaPtrkFitPntLTTS",
"primtrk: N fit pnts vs length, tpc,tpc+svt",25,70.,350.,25,0.,50.);
1513 m_ppointT = QAH::H1F(
"QaPtrkNPntT",
"primtrk: N points on trk,tpc", 81, -0.5,80.5);
1514 m_pmax_pointT = QAH::H1F(
"QaPtrkNPntMaxT",
"primtrk: N max pnts on trk, tpc", 76, -1,151);
1515 m_pfit_pointT = QAH::H1F(
"QaPtrkNPntFitT",
"primtrk: N fit pnts on trk, tpc", 81, -0.5,80.5);
1516 m_prim_ratiomT = QAH::H1F(
"QaPtrkRnmT",
"primtrk: ratio Nfit/max pnt, tpc", 55, 0., 1.1);
1517 m_prim_chargeT = QAH::H1F(
"QaPtrkChrgT",
"primtrk: charge, tpc ", 20,-2.,2.);
1518 m_prim_r0T = QAH::H1F(
"QaPtrkR0T",
"primtrk: radius at start (cm), tpc ", 50,0.,.1);
1519 m_prim_phi0T = QAH::H1F(
"QaPtrkPhi0T",
"primtrk: azimuth (phi) at start (deg,force 0,360), tpc ", 64,0.,360.);
1520 m_prim_z0T = QAH::H1F(
"QaPtrkZ0T",
"primtrk: z-coord at start (cm), tpc ", 50, -300.,300.);
1521 m_prim_curvT = QAH::H1F(
"QaPtrkCurvT",
"primtrk: log10 curvature (1/cm), tpc ", 80,-3.5,0.5);
1522 m_prim_f0 = QAH::MH1F(
"QaPtrkf0",
"primtrk: first point: hit-helix, tpc",60,-3.,3.,3);
1526 m_prim_xfT = QAH::H1F(
"QaPtrkXfT",
"primtrk: x of first hit on trk, tpc", 50,-200.,200.);
1527 m_prim_xf0 = QAH::H1F(
"QaPtrkXf0",
"primtrk: first point: x_hit - x_helix, tpc",60,-3.,3.);
1528 m_prim_yfT = QAH::H1F(
"QaPtrkYfT",
"primtrk: y of first hit on trk, tpc", 50,-200.,200.);
1529 m_prim_yf0 = QAH::H1F(
"QaPtrkYf0",
"primtrk: first point: y_hit - y_helix, tpc",60,-3.,3.);
1530 m_prim_zfT = QAH::H1F(
"QaPtrkZfT",
"primtrk: z of first hit on trk, tpc", 50,-200.,200.);
1531 m_prim_zf0 = QAH::H1F(
"QaPtrkZf0",
"primtrk: first point: z_hit - z_helix, tpc",60,-3.,3.);
1532 m_prim_radfT = QAH::H1F(
"QaPtrkRT",
"primtrk: radial position of first hit, tpc", 50,0.,200.);
1533 m_prim_rzf0 = QAH::MH1F(
"QaPtrkRZf0",
"primtrk: first point: hit - helix (r,z), tpc",60,-3.,3.,2);
1534 m_prim_rzf0->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1537 m_prim_rzl0 = QAH::MH1F(
"QaPtrkRZl0",
"primtrk: last point: hit - helix (r,z), tpc",60,-3.,3.,2);
1538 m_prim_rzl0->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1541 m_plengthT = QAH::H1F(
"QaPtrkLengthT",
"primtrk: track length, tpc", 50,0.,300.);
1542 m_ppsiT = QAH::H1F(
"QaPtrkPsiT",
"primtrk: psi, tpc (deg)", 36, 0.,360.);
1543 m_ptanlT = QAH::H1F(
"QaPtrkTanlT",
"primtrk: tanl, tpc",32,-4.,4.);
1544 m_prim_thetaT = QAH::H1F(
"QaPtrkThetaT",
"primtrk: theta (degrees), tpc",36,0.,180.);
1545 m_petaT = QAH::H1F(
"QaPtrkEtaT",
"primtrk: eta, tpc",40,-2.,2.);
1546 m_ppTT = QAH::H1F(
"QaPtrkPtT",
"primtrk: pT, tpc",50,0.,5.);
1547 m_pmomT = QAH::H1F(
"QaPtrkPT",
"primtrk: momentum, tpc",50,0.,5.);
1548 m_pchisq0T = QAH::H1F(
"QaPtrkChisq0T",
"primtrk: chisq0, tpc", 50, 0.,5.);
1551 m_ppT_eta_recT = QAH::H2F(
"QaPtrkPtVsEtaT",
"primtrk: log10 pT vs eta, tpc", 20,-2.,2.,40,1.,4.);
1554 m_primtrk_xf_yfTE = QAH::H2F(
"QaPtrkXfYfTE",
"primtrk: Y vs X of first hit on trk, tpc east", 40,-200.,200.,40,-200.,200.);
1557 m_primtrk_xf_yfTW = QAH::H2F(
"QaPtrkXfYfTW",
"primtrk: Y vs X of first hit on trk, tpc west", 40,-200.,200.,40,-200.,200.);
1560 m_ptanl_zfT = QAH::H2F(
"QaPtrkTanlzf",
"primtrk: tanl(dip) versus zfirst-zvtx, tpc,tpc+svt",60,-3.,3.,60,-3.,3.);
1563 m_pmom_trklengthT = QAH::H2F(
"QaPtrkPVsTrkLength",
"primtrk: log10 mom vs trk length, tpc",50,70.,350.,40,1.,4.);
1566 m_peta_trklengthT = QAH::H2F(
"QaPtrkLengthVEtaT",
"primtrk: trk length vs eta, tpc",20,-2.,2.,50,70.,350.);
1569 m_pnpoint_lengthT = QAH::H2F(
"QaPtrkNPntLengthT",
"primtrk: N pnts vs length, tpc",25,70.,350.,40,1.,81.);
1572 m_pfpoint_lengthT = QAH::H2F(
"QaPtrkFitPntLengthT",
"primtrk: N fit pnts vs length, tpc",25,70.,350.,40,1.,81.);
1575 m_pchisq0_momT = QAH::H2F(
"QaPtrkChi0MomT",
"primtrk: Chisq0 vs log10 mom, tpc",40,1.,4.,50,0.,5.);
1578 m_pchisq0_etaT = QAH::H2F(
"QaPtrkChi0EtaT",
"primtrk: Chisq0 vs eta, tpc",20,-2.,2.,20,0.,5.);
1581 m_pchisq0_dipT = QAH::H2F(
"QaPtrkChi0TanlT",
"primtrk: Chisq0 vs tanl(dip), tpc",20,-3.2,3.2,20,0.,5.);
1584 m_pchisq0_zfT = QAH::H2F(
"QaPtrkChi0zfT",
"primtrk: Chisq0 vs zfirst, tpc",20,-200.,200.,20,0.,5.);
1587 m_ppsi_phiT = QAH::H2F(
"QaPtrkPsiPhiT",
"primtrk: psi vs phi, tpc",36, 0.,360.,36,0.,360.);
1593 m_ppointTS = QAH::H1F(
"QaPtrkNPntTS",
"primtrk: N points on trk,tpc+svt", 60, 0.,60.);
1594 m_pmax_pointTS = QAH::H1F(
"QaPtrkNPntMaxTS",
"primtrk: N max pnts on trk, tpc+svt", 50, 0.,100.);
1595 m_pfit_pointTS = QAH::H1F(
"QaPtrkNPntFitTS",
"primtrk: N fit pnts on trk, tpc+svt", 60, 0.,60.);
1596 m_prim_ratiomTS = QAH::H1F(
"QaPtrkRnmTS",
"primtrk: ratio Nfit/max pnt, tpc+svt", 55, 0., 1.2005);
1597 m_prim_chargeTS = QAH::H1F(
"QaPtrkChrgTS",
"primtrk: charge, tpc+svt ", 20,-2.,2.);
1598 m_prim_r0TS = QAH::H1F(
"QaPtrkR0TS",
"primtrk: radius at start (cm), tpc+svt", 50,0.,0.1);
1599 m_prim_phi0TS = QAH::H1F(
"QaPtrkPhi0TS",
"primtrk: azimuth (phi) at start (deg,force 0-360),tpc+svt", 64, 0.,360.);
1600 m_prim_z0TS = QAH::H1F(
"QaPtrkZ0TS",
"primtrk: z-coord at start (cm), tpc+svt", 50, -50.,50.);
1601 m_prim_curvTS = QAH::H1F(
"QaPtrkCurvTS",
"primtrk: log10 curvature (1/cm), tpc+svt", 80,-3.5,0.5);
1602 m_prim_f0TS = QAH::MH1F(
"QaPtrkf0TS",
"primtrk: first point: hit-helix, tpc+svt",50,-5.,5.,3);
1607 m_prim_xfTS = QAH::H1F(
"QaPtrkXfTS",
"primtrk: x of first hit on trk, tpc+svt", 50,-200.,200.);
1608 m_prim_xf0TS = QAH::H1F(
"QaPtrkXf0TS",
"primtrk: first point: x_hit - x_helix, tpc+svt",50,-5.,5.);
1609 m_prim_yfTS = QAH::H1F(
"QaPtrkYfTS",
"primtrk: y of first hit on trk, tpc+svt", 50,-200.,200.);
1610 m_prim_yf0TS = QAH::H1F(
"QaPtrkYf0TS",
"primtrk: first point: y_hit - y_helix, tpc+svt",50,-5.,5.);
1611 m_prim_zfTS = QAH::H1F(
"QaPtrkZfTS",
"primtrk: z of first hit on trk, tpc+svt", 50,-200.,200.);
1612 m_prim_zf0TS = QAH::H1F(
"QaPtrkZf0TS",
"primtrk: first point: z_hit - z_helix, tpc+svt",50,-5.,5.);
1613 m_prim_radfTS = QAH::H1F(
"QaPtrkRTS",
"primtrk: radial position of first hit, tpc+svt", 50,0.,200.);
1614 m_prim_rzf0TS = QAH::MH1F(
"QaPtrkRZf0TS",
"primtrk: first point: hit - helix (r,z), tpc+svt",60,-3.,3.,2);
1615 m_prim_rzf0TS->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1618 m_prim_rzl0TS = QAH::MH1F(
"QaPtrkRZl0TS",
"primtrk: last point: hit - helix (r,z), tpc+svt",60,-3.,3.,2);
1619 m_prim_rzl0TS->Rebin(0,
"#pm #sqrt{ x_{dif}^{2}+y_{dif}^{2}}");
1622 m_plengthTS = QAH::H1F(
"QaPtrkLengthTS",
"primtrk: track length, tpc+svt", 50,0.,300.);
1623 m_ppsiTS = QAH::H1F(
"QaPtrkPsiTS",
"primtrk: psi, tpc+svt (deg) ", 36, 0.,360.);
1624 m_ptanlTS = QAH::H1F(
"QaPtrkTanlTS",
"primtrk: tanl, tpc+svt",32,-4.,4.);
1625 m_prim_thetaTS = QAH::H1F(
"QaPtrkThetaTS",
"primtrk: theta (degrees), tpc+svt",36,0.,180.);
1626 m_petaTS = QAH::H1F(
"QaPtrkEtaTS",
"primtrk: eta, tpc+svt",40,-2.,2.);
1627 m_ppTTS = QAH::H1F(
"QaPtrkPtTS",
"primtrk: pT, tpc+svt",50,0.,5.);
1628 m_pmomTS = QAH::H1F(
"QaPtrkPTS",
"primtrk: momentum, tpc+svt",50,0.,5.);
1629 m_pchisq0TS = QAH::H1F(
"QaPtrkChisq0TS",
"primtrk: chisq0, tpc+svt", 50, 0.,5.);
1630 m_prim_ssd_phi = QAH::H1F(
"QaPtrkPhifSSD",
"primtrk: phi of ssd point (per event)",64,0,360);
1633 m_ppT_eta_recTS = QAH::H2F(
"QaPtrkPtVsEtaTS",
"primtrk: log10 pT vs eta, tpc+svt", 20,-2.,2.,40,1.,4.);
1636 m_primtrk_xf_yfTS = QAH::H2F(
"QaPtrkXfYfTS",
"primtrk: Y vs X of first hit on trk, tpc+svt", 40,-200.,200.,40,-200.,200.);
1639 m_ptanl_zfTS = QAH::H2F(
"QaPtrkTanlzfTS",
"primtrk: tanl(dip) versus zfirst-zvtx, svt",60,-3.,3.,60,-3.,3.);
1642 m_pmom_trklengthTS = QAH::H2F(
"QaPtrkPVsTrkLTS",
"primtrk: log10 mom vs trk length, tpc+svt",50,70.,350.,40,1.,4.);
1645 m_peta_trklengthTS = QAH::H2F(
"QaPtrkLVEtaTS",
"primtrk: trk length vs eta, tpc+svt",20,-2.,2.,50,70.,350.);
1648 m_pnpoint_lengthTS = QAH::H2F(
"QaPtrkNPntLTS",
"primtrk: N pnts vs length, tpc+svt",25,70.,350.,25,0.,50.);
1651 m_pfpoint_lengthTS = QAH::H2F(
"QaPtrkFitPntLTS",
"primtrk: N fit pnts vs length, tpc+svt",25,70.,350.,25,0.,50.);
1654 m_pchisq0_momTS = QAH::H2F(
"QaPtrkChi0MomTS",
"primtrk: Chisq0 vs log10 mom, tpc+svt",40,1.,4.,50,0.,5.);
1657 m_pchisq0_etaTS = QAH::H2F(
"QaPtrkChi0EtaTS",
"primtrk: Chisq0 vs eta, tpc+svt",20,-2.,2.,20,0.,5.);
1660 m_pchisq0_dipTS = QAH::H2F(
"QaPtrkChi0TanlTS",
"primtrk: Chisq0 vs tanl(dip), tpc+svt",20,-3.2,3.2,20,0.,5.);
1663 m_pchisq0_zfTS = QAH::H2F(
"QaPtrkChi0zfTS",
"primtrk: Chisq0 vs zfirst, tpc+svt",20,-200.,200.,20,0.,5.);
1666 m_ppsi_phiTS = QAH::H2F(
"QaPtrkPsiPhiTS",
"primtrk: psi vs phi, tpc+svt",36, 0.,360.,36,0.,360.);
1675 m_ppointF = QAH::MH1F(
"QaPtrkNPntF",
"primtrk: N points on trk,ftpc",15,0.,15.,2);
1679 m_pmax_pointF = QAH::MH1F(
"QaPtrkNPntMaxF",
"primtrk: N max pnts on trk, ftpc",15,0.,15.,2);
1683 m_prim_ratiomF = QAH::MH1F(
"QaPtrkRnmF",
"primtrk: ratio Nfit/max pnt, ftpc",55,0.,1.1,2);
1687 m_prim_chargeF = QAH::MH1F(
"QaPtrkChrgF",
"primtrk: charge, ftpc",20,-2.,2.,2);
1691 m_prim_xfF = QAH::MH1F(
"QaPtrkXfF",
"primtrk: x of first hit on trk, ftpc",70,-35.,35.,2);
1695 m_prim_yfF = QAH::MH1F(
"QaPtrkYfF",
"primtrk: y of first hit on trk, ftpc",70,-35.,35.,2);
1699 m_prim_zfF = QAH::MH1F(
"QaPtrkZfF",
"primtrk: z of first hit on trk, ftpc",75,-225.,225.,2);
1703 m_prim_radfF = QAH::MH1F(
"QaPtrkRF",
"primtrk: radial position of first hit, ftpc",35,0.,35.,2);
1707 m_plengthF = QAH::MH1F(
"QaPtrkLengthF",
"primtrk: track length, ftpc",50,150.,300.,2);
1711 m_ppsiF = QAH::MH1F(
"QaPtrkPsiF",
"primtrk: psi, ftpc", 90, 0.,360.,2);
1715 m_petaF = QAH::MH1F(
"QaPtrkEtaF",
"primtrk: |eta|, ftpc",80,2,5,2);
1719 m_ppTF = QAH::MH1F(
"QaPtrkPtF",
"primtrk: pT, ftpc",50,0.,5.,2);
1722 m_ppTF->SetStats(kFALSE);
1723 m_pmomF = QAH::MH1F(
"QaPtrkPF",
"primtrk: momentum, ftpc",50,0.,5.,2);
1727 m_pchisq0F = QAH::MH1F(
"QaPtrkChisq0F",
"primtrk: chi2/ndf x,y fit, ftpc",50,0.,5.,2);
1731 m_pchisq1F = QAH::MH1F(
"QaPtrkChisq1F",
"primtrk: chi2/ndf r,z fit, ftpc",50,0.,5.,2);
1736 m_ppointFE = QAH::H1F(
"QaPtrkNPntFE",
"primtrk: N points on trk,ftpc east", 15, 0.,15.);
1737 m_ppointFW = QAH::H1F(
"QaPtrkNPntFW",
"primtrk: N points on trk,ftpc west", 15, 0.,15.);
1738 m_pmax_pointFE = QAH::H1F(
"QaPtrkNPntMaxFE",
"primtrk: N max pnts on trk, ftpc east", 15, 0.,15.);
1739 m_pmax_pointFW = QAH::H1F(
"QaPtrkNPntMaxFW",
"primtrk: N max pnts on trk, ftpc west", 15, 0.,15.);
1740 m_prim_ratiomFE = QAH::H1F(
"QaPtrkRnmFE",
"primtrk: ratio Nfit/max pnt, ftpc east", 55, 0., 1.1);
1741 m_prim_ratiomFW = QAH::H1F(
"QaPtrkRnmFW",
"primtrk: ratio Nfit/max pnt, ftpc west", 55, 0., 1.1);
1742 m_prim_chargeFE = QAH::H1F(
"QaPtrkChrgFE",
"primtrk: charge, ftpc east ", 20,-2.,2.);
1743 m_prim_chargeFW = QAH::H1F(
"QaPtrkChrgFW",
"primtrk: charge, ftpc west ", 20,-2.,2.);
1744 m_prim_xfFE = QAH::H1F(
"QaPtrkXfFE",
"primtrk: x of first hit on trk, ftpc east",70,-35.,35.);
1745 m_prim_xfFW = QAH::H1F(
"QaPtrkXfFW",
"primtrk: x of first hit on trk, ftpc west",70,-35.,35.);
1746 m_prim_yfFE = QAH::H1F(
"QaPtrkYfFE",
"primtrk: y of first hit on trk, ftpc east",70,-35.,35.);
1747 m_prim_yfFW = QAH::H1F(
"QaPtrkYfFW",
"primtrk: y of first hit on trk, ftpc west",70,-35.,35.);
1748 m_prim_zfFE = QAH::H1F(
"QaPtrkZfFE",
"primtrk: z of first hit on trk, ftpc east",75,-225.,-150.);
1749 m_prim_zfFW = QAH::H1F(
"QaPtrkZfFW",
"primtrk: z of first hit on trk, ftpc west",75,150.,225.);
1750 m_prim_radfFE = QAH::H1F(
"QaPtrkRFE",
"primtrk: radial position of first hit, ftpc east",35,0.,35.);
1751 m_prim_radfFW = QAH::H1F(
"QaPtrkRFW",
"primtrk: radial position of first hit, ftpc west",35,0.,35.);
1752 m_plengthFE = QAH::H1F(
"QaPtrkLengthFE",
"primtrk: track length, ftpc east",50,150.,300.);
1753 m_plengthFW = QAH::H1F(
"QaPtrkLengthFW",
"primtrk: track length, ftpc west",50,150.,300.);
1754 m_ppsiFE = QAH::H1F(
"QaPtrkPsiFE",
"primtrk: psi, ftpc east", 90, 0.,360.);
1755 m_ppsiFW = QAH::H1F(
"QaPtrkPsiFW",
"primtrk: psi, ftpc west", 90, 0.,360.);
1756 m_petaFE = QAH::H1F(
"QaPtrkEtaFE",
"primtrk: eta, ftpc east",80,-4.5,-2.);
1757 m_petaFW = QAH::H1F(
"QaPtrkEtaFW",
"primtrk: eta, ftpc west",80,2.,4.5);
1758 m_ppTFE = QAH::H1F(
"QaPtrkPtFE",
"primtrk: pT, ftpc east",50,0.,5.);
1759 m_ppTFW = QAH::H1F(
"QaPtrkPtFW",
"primtrk: pT, ftpc west",50,0.,5.);
1760 m_pmomFE = QAH::H1F(
"QaPtrkPFE",
"primtrk: momentum, ftpc east ",50,0.,5.);
1761 m_pmomFW = QAH::H1F(
"QaPtrkPFW",
"primtrk: momentum, ftpc west ",50,0.,5.);
1762 m_pchisq0FE = QAH::H1F(
"QaPtrkChisq0FE",
"primtrk: chi2/ndf x,y fit, ftpc east", 50, 0.,5.);
1763 m_pchisq0FW = QAH::H1F(
"QaPtrkChisq0FW",
"primtrk: chi2/ndf x,y fit, ftpc west", 50, 0.,5.);
1764 m_pchisq1FE = QAH::H1F(
"QaPtrkChisq1FE",
"primtrk: chi2/ndf r,z fit, ftpc east", 50, 0.,5.);
1765 m_pchisq1FW = QAH::H1F(
"QaPtrkChisq1FW",
"primtrk: chi2/ndf r,z ft, ftpc west", 50, 0.,5.);
1768 m_ppT_eta_recFE = QAH::H2F(
"QaPtrkPtVsEtaFE",
"primtrk: log10 pT vs eta, ftpcE",20,-4.5,-2.,40,1.,4.);
1769 m_ppT_eta_recFW = QAH::H2F(
"QaPtrkPtVsEtaFW",
"primtrk: log10 pT vs eta, ftpcW",20,2.,4.5,40,1.,4.);
1770 m_primtrk_xf_yfFE = QAH::H2F(
"QaPtrkXfYfFE",
"primtrk: Y vs X of first hit on trk, ftpcE",70,-35.,35.,75,-35.,35.);
1773 m_primtrk_xf_yfFW = QAH::H2F(
"QaPtrkXfYfFW",
"primtrk: Y vs X of first hit on trk, ftpcW",70,-35.,35.,70,-35.,35.);
1776 m_peta_trklengthFE = QAH::H2F(
"QaPtrkLengthVEtaFE",
"primtrk: trk length vs eta, ftpcE",25,-4.5,-2.,50,0.,300.);
1779 m_peta_trklengthFW = QAH::H2F(
"QaPtrkLengthVEtaFW",
"primtrk: trk length vs eta, ftpcW",25,2.,4.5,50,0.,300.);
1782 m_pnpoint_lengthFE = QAH::H2F(
"QaPtrkNPntLengthFE",
"primtrk: N pnts vs length, ftpcE",50,0.,300.,15,0.,15.);
1785 m_pnpoint_lengthFW = QAH::H2F(
"QaPtrkNPntLengthFW",
"primtrk: N pnts vs length, ftpcW",50,0.,300.,15,0.,15.);
1793 m_svt_loc = QAH::H2F(
"QaPtrkSvtLoc",
"primtrk: SVT hit time bins",256,0,128,432,-0.5,431.5);
1795 m_primary_pxl_hit = QAH::H1F(
"QaPtrkPxlHit",
"PIXEL: Hits per primary track",5,-0.5,4.5);
1797 m_primary_ist_hit = QAH::H1F(
"QaPtrkIstHit",
"IST: Hits per primary track",5,-0.5,4.5);
1799 m_primary_hft_hit = QAH::H1F(
"QaPtrkHftHit",
"HFT: Hits per primary track",5,-0.5,4.5);
1810 void StQABookHist::BookHistDE(){
1813 m_ndedxr = QAH::H1F(
"QaDedxNum",
"dedx: number of tracks", 64,0., 16000.);
1815 m_ndedxT = QAH::H1F(
"QaDedxNdedxT",
"dedx: number of point to define dE/dx, tpc",50,0., 50.);
1816 m_dedx0T = QAH::H1F(
"QaDedxDedx0T",
"dedx: dE/dx mean (GeV/cm), tpc", ndedx, 0., 0.0005);
1817 m_dedx1T = QAH::H1F(
"QaDedxDedx1T",
"dedx: dE/dx error on mean, tpc", ndedx, 0., 1);
1820 m_dedxTTS = QAH::H1F(
"QaDedxBTTS",
"dedx: ratio <dE/dx> mean to Bichsel <dE/dx>, tpc,tpc+svt",
1823 m_ndedxF = QAH::MH1F(
"QaDedxNdedxF",
"dedx: number of point to define dE/dx, ftpc",50,0.,10.,2);
1827 m_dedx0F = QAH::MH1F(
"QaDedxDedx0F",
"dedx: dE/dx mean, ftpc",ndedx,0.,750.,2);
1832 m_ndedxFE = QAH::H1F(
"QaDedxNdedxFE",
"dedx: number of point to define dE/dx, ftpcE",50,0.,10.);
1833 m_dedx0FE = QAH::H1F(
"QaDedxDedx0FE",
"dedx: dE/dx mean, ftpcE", ndedx, 0., 750.);
1835 m_ndedxFW = QAH::H1F(
"QaDedxNdedxFW",
"dedx: number of point to define dE/dx, ftpcW",50,0.,10.);
1836 m_dedx0FW = QAH::H1F(
"QaDedxDedx0FW",
"dedx: dE/dx mean, ftpcW", ndedx, 0., 750.);
1839 void StQABookHist::BookHistPID(){
1844 m_p_dedx_rec = QAH::H2F(
"QaPidGlobtrkDstdedxPVsDedx",
"PID: globtrk-dst_dedx, p vs dedx (reconstructed)",
1845 cnp,cminp,cmaxp,cndedx,cmindedx,cmaxdedx);
1851 void StQABookHist::BookHistVertex(){
1855 m_v_num = QAH::H1F(
"QaVtxNum",
" vertex: num vertices ",50,0.,10000.);
1856 m_v_num_sm= QAH::H1F(
"QaVtxNumsm",
" vertex: num vertices ",50,0.,500.);
1858 m_v_vtxid = QAH::H1F(
"QaVtxVtxId",
" vertex,2ndary: Vertex ID ",10,0.,10.);
1859 m_v_x = QAH::H1F(
"QaVtxX",
" vertex,2ndary: x ",50,-250.,250.);
1860 m_v_y = QAH::H1F(
"QaVtxY",
" vertex,2ndary: y ",50,-250.,250.);
1861 m_v_z = QAH::H1F(
"QaVtxZ",
" vertex,2ndary: z ",50,-250.,250.);
1862 m_v_pchi2 = QAH::H1F(
"QaVtxChisq",
" vertex,2ndary: chisq/dof ",50,0.,50.);
1863 m_v_r = QAH::H1F(
"QaVtxR",
" vertex,2ndary: r ",100,0,50);
1865 m_pv_vtxid = QAH::H1F(
"QaVtxPrVtxId",
" vertex,prim: Vertex ID ",10,0.,10.);
1866 m_pv_x = QAH::H1F(
"QaVtxPrX",
" vertex,prim: x ",50,-5.,5.);
1867 m_pv_y = QAH::H1F(
"QaVtxPrY",
" vertex,prim: y ",50,-5.,5.);
1868 m_pv_z = QAH::H1F(
"QaVtxPrZ",
" vertex,prim: z ",50,-50.,50.);
1869 m_pv_xy = QAH::H2F(
"QaVtxPrXY",
" vertex,prim: x versus y",50,-5.,5.,50,-5.,5.);
1870 m_pv_pchi2 = QAH::H1F(
"QaVtxPrChisq",
" vertex,prim: chisq/dof ",40,0.,20.);
1871 m_pv_r = QAH::H1F(
"QaVtxPrR",
" vertex,prim: r ",100,0,0.1);
1878 "V0 azimuthal distribution relative to primvtx",36,0.,360.);
1880 "V0 radial distribution relative to primvtx",50,0.,200.);
1882 "V0 Z distribution relative to primvtx",60,-30.,30.);
1884 m_v0 = QAH::H1F(
"QaV0Vtx",
"log10 total V0s (0 => -0.5)",45,-0.55,3.95);
1885 m_ev0_lama_hist = QAH::H1F(
"QaV0LambdaMass",
"V0: Lambda mass",25,1.05,1.15);
1888 m_xi_tot = QAH::H1F(
"QaXiVtxTot",
"log10 total Xis (0 => -0.5)",45,-0.55,3.95);
1889 m_xi_ma_hist = QAH::H1F(
"QaXiaMass",
"Xi: Xi mass",25,1.2,1.4);
1891 m_kink_tot = QAH::H1F(
"QaKinkTot",
"log10 total kinks (0 => -0.5)",35,-0.55,2.95);
1894 " vertex,prim: x(ftpcE)-x(tpc) vs y(ftpcE)-y(tpc)",
1895 800, -20., 20., 800, -20., 20.);
1897 " vertex,prim: z(ftpcE)-z(tpc)",
1900 " vertex,prim: x(ftpcW)-x(tpc) vs y(ftpcW)-y(tpc)",
1901 800, -20., 20., 800, -20., 20.);
1903 " vertex,prim: z(ftpcW)-z(tpc)",
1908 void StQABookHist::BookHistPoint(){
1910 m_pnt_tot = QAH::H1F(
"QaPointTot",
"point: # hits total ",100, 0.,400000.);
1911 m_pnt_tot_med = QAH::H1F(
"QaPointTotmed",
"point: # hits total ",100, 0.,25000.);
1912 m_pnt_tot_sm = QAH::H1F(
"QaPointTotsm",
"point: # hits total ",100, 0.,2500.);
1913 m_pnt_id = QAH::H1F(
"QaPointId",
"point: detector ID of hit",30,0.,30.);
1915 m_pnt_tpc = QAH::H1F(
"QaPointTpc",
"point: # hits tpc ",100, 0.,300000.);
1917 m_pnt_ftpc = QAH::MH1F(
"QaPointFtpc",
"point: # hits ftpc",100,0.,25000.,2);
1922 m_pnt_ftpcE = QAH::H1F(
"QaPointFtpcE",
"point: # hits ftpcE ",100, 0.,25000.);
1923 m_pnt_ftpcW = QAH::H1F(
"QaPointFtpcW",
"point: # hits ftpcW ",100, 0.,25000.);
1925 m_pnt_svt = QAH::H1F(
"QaPointSvt",
"point: # hits svt ",600, 0.,15000.);
1926 m_pnt_svtLaser= QAH::H2F(
"QaPointSvtLaser",
"point: laser spots, svt ",150,0,600,65,0.,130.);
1928 m_pnt_svtLaserDiff= QAH::MH2F(
"QaPointSvtLaserDiff",
"point: diff of laser spots, svt ",150,0,600,101,9.8,50.2,2);
1934 if (silHists || (histsSet>=StQA_run14 && histsSet<StQA_run17)) {
1935 m_pnt_xyS = QAH::H2F(
"QaPointXYSvt",
"point: x-y distribution of hits, svt,ssd",125,-25,25,125,-25,25);
1942 if (histsSet>=StQA_run14 && histsSet<StQA_run17)
1943 m_pnt_xyS->SetTitle(
"PIXEL, IST, SSD: Distribution of hits in XY");
1947 m_pnt_rpTW = QAH::H2F(
"QaPointRPTpcW",
"point: r-phi distribution of hits, tpcW",20,58.75,196.75,72,0,TMath::TwoPi());
1948 m_pnt_rpTE = QAH::H2F(
"QaPointRPTpcE",
"point: r-phi distribution of hits, tpcE",23,58.75,196.75,72,0,TMath::TwoPi());
1949 m_z_hits = QAH::H1F(
"QaPointZhits",
"point: z distribution of hits, tpc",100,-210,210);
1950 m_pnt_timeT = QAH::MH1F(
"QaPointTimeT",
"point: time bucket distribution of hits, tpc",45,0,450,2);
1954 m_pnt_phiT = QAH::MH1F(
"QaPointPhiT",
"point: #phi distribution of hits, tpc",36,0,360,2);
1958 m_pnt_padrowT = QAH::MH1F(
"QaPointPadrowT",
"point: padrow distribution of hits, tpc",72,0.5,72.5,2);
1965 m_pnt_zS = QAH::H1F(
"QaPointZhitsS",
"point: z distribution of hits, svt",100,-35,35);
1966 m_pnt_phiS = QAH::H1F(
"QaPointPhiS",
"point: #phi distribution of hits, svt",36,0,360);
1967 m_pnt_barrelS = QAH::H1F(
"QaPointBarrelS",
"point: barrel distribution of hits, svt",3,0.5,3.5);
1971 m_pnt_ssd = QAH::H1F(
"QaPointSsd",
"point: # hits ssd ",200, 0.,5000.);
1972 m_pnt_phiSSD = QAH::H1F(
"QaPointPhiSSD",
"SSD: #phi of hits (per event)",36,0,360);
1973 m_pnt_lwSSD = QAH::H2F(
"QaPointLWSSD",
"SSD: wafer id vs ladder id (per event)",20,0.5,20.5,16,0.5,16.5);
1974 m_pnt_lwSSD->SetXTitle(
"Ladder #");
1975 m_pnt_lwSSD->SetYTitle(
"Wafer (sensor) #");
1976 m_pnt_lwSSD->SetStats(kFALSE);
1977 m_pnt_sizeSSD = QAH::MH1F(
"QaPointSizeSSD",
"SSD: size of clusters",10,0.5,10.5,2);
1978 m_pnt_sizeSSD->Rebin(0,
"P-side");
1979 m_pnt_sizeSSD->Rebin(1,
"N-side");
1980 m_pnt_sizeSSD->SetStats(kFALSE);
1981 m_pnt_eSSD = QAH::H1F(
"QaPointESSD",
"SSD: log10(energy) of hits",90,-5,-2);
1983 if (histsSet>=StQA_run14 && histsSet<StQA_run17) {
1984 m_pnt_sst = QAH::H1F(
"QaPointSst",
"point: # hits sst ",200, 0.,5000.);
1985 m_pnt_phiSST = QAH::H1F(
"QaPointPhiSST",
"SST: #phi of hits (per event)",36,0,360);
1986 m_pnt_lwSST = QAH::H2F(
"QaPointLWSST",
"SST: wafer id vs ladder id (per event)",20,0.5,20.5,16,0.5,16.5);
1987 m_pnt_lwSST->SetXTitle(
"Ladder #");
1988 m_pnt_lwSST->SetYTitle(
"Wafer (sensor) #");
1989 m_pnt_lwSST->SetStats(kFALSE);
1990 m_pnt_sizeSST = QAH::MH1F(
"QaPointSizeSST",
"SST: size of clusters",10,0.5,10.5,2);
1991 m_pnt_sizeSST->Rebin(0,
"P-side");
1992 m_pnt_sizeSST->Rebin(1,
"N-side");
1993 m_pnt_sizeSST->SetStats(kFALSE);
1994 m_pnt_eSST = QAH::H1F(
"QaPointESST",
"SST: log10(energy) of hits",90,-5,-2);
1997 m_pnt_xyFE = QAH::H2F(
"QaPointXYFtpcE",
"point: x-y distribution of hits, ftpcE",70,-35,35,70,-35,35);
1998 m_pnt_xyFW = QAH::H2F(
"QaPointXYFtpcW",
"point: x-y distribution of hits, ftpcW",70,-35,35,70,-35,35);
1999 m_pnt_planeF = QAH::MH1F(
"QaPointPlaneF",
"point: plane distribution of hits, ftpc",20,0.5,20.5,2);
2004 m_pnt_padtimeFE = QAH::H2F(
"QaPointPadTimeFtpcE",
"point: #pads vs #timebins of hits, ftpcE",12,0.5,12.5,10,0.5,10.5);
2007 m_pnt_padtimeFW = QAH::H2F(
"QaPointPadTimeFtpcW",
"point: #pads vs #timebins of hits, ftpcW",12,0.5,12.5,10,0.5,10.5);
2014 void StQABookHist::BookHistEMC(){
2016 if (!((gROOT->GetClass(
"StEmcMath")) && (gROOT->GetClass(
"StEmcGeom"))))
2019 m_emc_nhit=QAH::H2F(
"EmcNHitsVsDet",
"emc: Number of hit(log) .vs. Detector #",100,0.0,4.5,4,0.5,4.5);
2020 m_emc_etot=QAH::H2F(
"EmcEtotVsDet",
"emc: Total energy(log) .vs. Detector #",100,-4.0,4.5,8,0.5,4.5);
2022 const Char_t* tit={
"Barrel"};
2023 const Int_t nx[4] = {40,40,300,20};
2024 const Int_t ny[4] = {120, 120, 60, 900};
2027 Axis_t ModNumLo[4] = {1.,31.,61.,91.};
2028 Axis_t ModNumHi[4] = {31.,61.,91.,121.};
2029 const TString PerModuleHistName[4] = {
"West1",
"West2",
"East1",
"East2"};
2031 Float_t rpi = M_PI + 0.00001;
2032 TString name, title;
2033 TArrayD *xb = StEmcMath::binForSmde();
2035 for(Int_t i=0; i<4; i++){
2036 name = detname[i] +
"Hits";
2037 title = tit + detname[i] +
" hits dist.";
2038 if(i==2) m_emc_hits[i]=QAH::H2F(name,title, xb->GetSize()-1,xb->GetArray(),ny[i],-rpi,rpi);
2039 else m_emc_hits[i] = QAH::H2F(name,title, nx[i],-1.,+1., ny[i],-rpi, rpi);
2041 name = detname[i] +
"Energy2D";
2042 title = tit + detname[i] +
" energy dist. in eta-phi";
2043 if(i==2) m_emc_energy2D[i] = QAH::H2F(name,title, xb->GetSize()-1,xb->GetArray(), ny[i],-rpi,rpi);
2044 else m_emc_energy2D[i] = QAH::H2F(name,title, nx[i],-1.,+1., ny[i],-rpi, rpi);
2046 name = detname[i] +
"Adc";
2047 title = tit + detname[i] +
" ADC dist.";
2048 m_emc_adc[i] = QAH::H1F(name,title, 5001, -0.5, 5000.5);
2050 name = detname[i] +
"Energy";
2051 title = tit + detname[i] +
" energy dist.";
2052 m_emc_energy[i] = QAH::H1F(name,title, 600, 0.0, 60.0);
2056 for(Int_t i=2; i<4; i++){
2057 for (Int_t j=0; j<4; j++){
2062 name = detname[i] + PerModuleHistName[j] +
"HitsPerModule";
2063 title = tit + detname[i] +
" " + PerModuleHistName[j] +
" - hits per module";
2064 m_emc_hits_per_module[k] = QAH::H1F(name,title, 30, ModNumLo[j], ModNumHi[j]);
2066 name = detname[i] + PerModuleHistName[j] +
"EnergyPerModule";
2067 title = tit + detname[i] +
" " + PerModuleHistName[j] +
" - energy-weighted hits per module";
2068 m_emc_energy_per_module[k] = QAH::H1F(name,title, 30, ModNumLo[j], ModNumHi[j]);
2070 name = detname[i] + PerModuleHistName[j] +
"StripHitsPerModule";
2071 title = tit + detname[i] +
" " + PerModuleHistName[j] +
" - hits in strip (within mod) vs. mod";
2072 m_emc_strip_hits_per_module[k] = QAH::H2F(name,title, 30, ModNumLo[j], ModNumHi[j], 150, 1., 151.);
2074 name = detname[i] + PerModuleHistName[j] +
"StripEnergyPerModule";
2075 title = tit + detname[i] +
" " + PerModuleHistName[j] +
" - energy in strip (within mod) vs. mod";
2076 m_emc_strip_energy_per_module[k] = QAH::H2F(name,title, 30, ModNumLo[j], ModNumHi[j], 150, 1., 151.);
2081 Int_t greta[4]={40,40,300,20};
2082 Int_t grphi[4]={120,120,60,900};
2083 Float_t myPI = M_PI + 0.0001;
2085 m_emc_ncl=QAH::H2F(
"EmcNcluster",
"emc: Number of cluster(log) .vs. Detector #",40,0.0,4.0, 4,0.5,4.5);
2086 m_emc_etotCl=QAH::H2F(
"EmcEcluster" ,
"emc: Total PreCluster Energy(log) .vs. Detector #", 60,-2.0,4.0, 4,0.5,4.5);
2088 Float_t rmsMax=0.026;
2090 m_emc_sig_e= QAH::H2F(
"EmcRMSeta" ,
"emc: Sigma(eta) .vs. Detector #",rmsN,0.0,rmsMax,4,0.5,4.5);
2091 m_emc_sig_p= QAH::H2F(
"EmcRMSphi" ,
"emc: Sigma(phi) .vs. Detector #",rmsN,0.0,rmsMax,4,0.5,4.5);
2092 for (Int_t i=0; i<4; i++)
2094 TString name_h = detname[i] +
"_cluster";
2095 TString name_e = detname[i] +
"_cluster_energy";
2096 TString tit_h = detname[i] +
" cluster";
2097 TString tit_e = detname[i] +
" energy of cluster";
2099 m_emc_cl[i] = QAH::H2F(name_h,tit_h,greta[i],-1.0,1.0,grphi[i],-M_PI*1.015, M_PI*0.985);
2100 m_emc_energyCl[i] = QAH::H2F(name_e,tit_e,greta[i],-1.0,1.0,grphi[i],-M_PI*1.015, M_PI*0.985);
2102 m_emc_cl[i] = QAH::H2F(name_h,tit_h,greta[i],-1.0,1.0,grphi[i],-myPI, myPI);
2103 m_emc_energyCl[i] = QAH::H2F(name_e,tit_e,greta[i],-1.0,1.0,grphi[i],-myPI, myPI);
2107 name = detname[i] +
"ClNum";
2108 title =
"Number hits in cluster for " + detname[i];
2109 m_emc_HitsInCl[i] = QAH::H1F(name, title, 21, -0.5, 20.5);
2111 name = detname[i] +
"ClEnergy";
2112 title =
"Energy of cluster for " + detname[i];
2113 m_emc_EnergyCl[i] = QAH::H1F(name, title, 2000, 0.0, 20.0);
2115 TString name_eta = detname[i] +
"Eta";
2116 TString tit_eta =
"Eta of clusters for " + detname[i];
2117 TString name_phi = detname[i] +
"Phi";
2118 TString tit_phi =
"Phi of clusters for " + detname[i];
2120 TArrayD *xb = StEmcMath::binForSmde();
2122 m_emc_EtaInCl[i] = QAH::H1F(name_eta, tit_eta, xb->GetSize()-1, xb->GetArray());
2125 m_emc_PhiInCl[i] = QAH::H1F(name_phi, tit_phi, grphi[i], -M_PI*1.015, M_PI*0.985);
2127 m_emc_EtaInCl[i] = QAH::H1F(name_eta, tit_eta, greta[i], -1., 1.);
2128 m_emc_PhiInCl[i] = QAH::H1F(name_phi, tit_phi, grphi[i], -myPI, myPI);
2132 const TString catname[] = {
"EmcCat1",
"EmcCat2",
"EmcCat3",
"EmcCat4"};
2134 for (Int_t i=0; i<4; i++) {
2135 name = catname[i] +
"_Point_Energy";
2136 title = catname[i] +
" Point Energy";
2137 m_emc_point_energy[i]= QAH::H1F(name,title,100,0.,10.);
2139 name = catname[i] +
"_Point_Eta";
2140 title = catname[i] +
" Point Eta";
2143 name = catname[i] +
"_Point_Phi";
2144 title = catname[i] +
" Point Phi";
2147 name = catname[i] +
"_Sigma_Eta";
2148 title = catname[i] +
" Sigma Eta";
2151 name = catname[i] +
"_Sigma_Phi";
2152 title = catname[i] +
" Sigma Phi";
2155 name = catname[i] +
"_Delta_Eta";
2156 title = catname[i] +
" Delta Eta";
2159 name = catname[i] +
"_Delta_Phi";
2160 title = catname[i] +
" Delta Phi";
2163 name = catname[i] +
"_Points_Multiplicity";
2164 title = catname[i] +
" Points Multiplicity";
2167 name = catname[i] +
"_Track_Momenta";
2168 title = catname[i] +
" Track Momenta ";
2176 void StQABookHist::BookHistEval(){
2187 " reco pvtx Z vs diff geant - reco Z", 100, -0.5,0.5,100,-50.,50.);
2193 void StQABookHist::BookHistBBC(){
2198 m_bbc_adc[0] = QAH::MH1F(
"QaBbcAdcES",
"BBC East Small ADC",100,0.5,400.5,8);
2199 m_bbc_adc[1] = QAH::MH1F(
"QaBbcAdcEL",
"BBC East Large ADC",100,0.5,400.5,8);
2200 m_bbc_adc[2] = QAH::MH1F(
"QaBbcAdcWS",
"BBC West Small ADC",100,0.5,400.5,8);
2201 m_bbc_adc[3] = QAH::MH1F(
"QaBbcAdcWL",
"BBC West Large ADC",100,0.5,400.5,8);
2202 m_bbc_tdc[0] = QAH::MH1F(
"QaBbcTdcES",
"BBC East Small TDC",100,0.5,750.5,8);
2203 m_bbc_tdc[1] = QAH::MH1F(
"QaBbcTdcEL",
"BBC East Large TDC",100,0.5,2000.5,8);
2204 m_bbc_tdc[2] = QAH::MH1F(
"QaBbcTdcWS",
"BBC West Small TDC",100,0.5,750.5,8);
2205 m_bbc_tdc[3] = QAH::MH1F(
"QaBbcTdcWL",
"BBC West Large TDC",100,0.5,2000.5,8);
2206 for (i=0; i<8; i++) {
2207 sprintf(ID,
"%d",i+1);
2208 for (j=0; j<4; j++) {
2210 m_bbc_tdc[j]->Rebin(i,ID);
2216 void StQABookHist::BookHistFPD(){
2218 QAH::MMH1F(m_fpd_top,2,
"QaFpdTop%d",
"FPD Top ADC %d-%d",100,0.5,1500.5,8,1);
2219 QAH::MMH1F(m_fpd_bottom,2,
"QaFpdBottom%d",
"FPD Bottom ADC %d-%d",100,0.5,1500.5,8,1);
2220 QAH::MMH1F(m_fpd_south,2,
"QaFpdSouth%d",
"FPD South ADC %d-%d",100,0.5,1500.5,8,1);
2221 QAH::MMH1F(m_fpd_north,2,
"QaFpdNorth%d",
"FPD North ADC %d-%d",100,0.5,1500.5,6,1);
2224 m_fpd_sums[0] = QAH::H1F(
"QaFpdSums0",
"FPD SUM Top",100,0.5,2050.5);
2225 m_fpd_sums[1] = QAH::H1F(
"QaFpdSums1",
"FPD SUM Bottom",100,0.5,2050.5);
2226 m_fpd_sums[2] = QAH::H1F(
"QaFpdSums2",
"FPD SUM South",100,0.5,2050.5);
2227 m_fpd_sums[3] = QAH::H1F(
"QaFpdSums3",
"FPD SUM North",100,0.5,2050.5);
2228 m_fpd_sums[4] = QAH::H1F(
"QaFpdSums4",
"FPD SUM SmdX",100,0.5,6000.5);
2229 m_fpd_sums[5] = QAH::H1F(
"QaFpdSums5",
"FPD SUM SmdY",100,0.5,6000.5);
2230 m_fpd_sums[6] = QAH::H1F(
"QaFpdSums6",
"FPD SUM Pres1",100,0.5,1500.5);
2231 m_fpd_sums[7] = QAH::H1F(
"QaFpdSums7",
"FPD SUM Pres2",100,0.5,1500.5);
2235 void StQABookHist::BookHistPMD(){
2237 QAH::MMH1F(m_pmd_sm_hit,12,
"QaPmdSmHit%d",
2238 "PMD SM-wise Hit Multiplicity %02d-%02d",1000,0.,500.,2);
2239 QAH::MMH1F(m_pmd_sm_adc,12,
"QaPmdSmAdcHit%d",
2240 "PMD SM-wise ADC/Hit Multiplicity %02d-%02d",100,0.,600.,2);
2241 QAH::MMH1F(m_pmd_chain_hit,24,
"QaPmdChHitChain%d",
2242 "PMD Channel-wise Hit, Chain %02d-%02d",1728,-0.5,1727.5,2);
2243 QAH::MMH1F(m_pmd_chain_adc,24,
"QaPmdChAdcChain%d",
2244 "PMD Channel-wise ADC, Chain %02d-%02d",1728,-0.5,1727.5,2);
2246 m_pmd_total_hit = QAH::H2F(
"QaPmdTotalHit",
"PMD Total Hits",100,0.,2e5,100,0.,4.);
2247 m_pmd_total_hit->SetXTitle(
"event id");
2248 m_pmd_total_hit->SetYTitle(
"log10");
2249 m_pmd_total_adc = QAH::H2F(
"QaPmdTotalAdc",
"PMD Total ADC",100,0.,2e5,100,0.,6.);
2250 m_pmd_total_adc->SetXTitle(
"event id");
2251 m_pmd_total_adc->SetYTitle(
"log10");
2252 m_cpv_total_hit = QAH::H2F(
"QaCpvTotalHit",
"CPV Total Hits",100,0.,2e5,100,0.,4.);
2253 m_cpv_total_hit->SetXTitle(
"event id");
2254 m_cpv_total_hit->SetYTitle(
"log10");
2255 m_cpv_total_adc = QAH::H2F(
"QaCpvTotalAdc",
"CPV Total ADC",100,0.,2e5,100,0.,6.);
2256 m_cpv_total_adc->SetXTitle(
"event id");
2257 m_cpv_total_adc->SetYTitle(
"log10");
2261 void StQABookHist::BookHistTOF(){
2263 m_tof_hit_tray = QAH::H2F(
"QaTofHitvsTray",
"TOF Hits vs tray",120,0.5,120.5,100,0.,100.);
2264 m_tof_hit_module = QAH::H2F(
"QaTofHitvsModule",
"TOF Hits vs Module",65,-32.5,32.5,100,0.,100.);
2265 m_tof_match_tray = QAH::H2F(
"QaTofMatchvsTray",
"TOF Matched Hits vs tray",120,0.5,120.5,100,0.,100.);
2266 m_tof_match_module = QAH::H2F(
"QaTofMatchvsModule",
"TOF Matched Hits vs Module",65,-32.5,32.5,100,0.,100.);
2267 m_tof_vpd_hit = QAH::H2F(
"QaTofHitvsVpdHit",
"TOF Hits vs Vpd Hits",50,0.,50.,100,0.,5000.);
2268 m_tof_vtx_z = QAH::H2F(
"QaTofVpdZvsTpcZ",
"VPD vtxz vs TPC vtxz",100,-100.,100.,100,-100.,100.);
2269 m_tof_PID = QAH::H2F(
"QaTofPID",
"TOF InvBeta vs p",100,0.,5.,100,0.,4.);
2273 void StQABookHist::BookHistHFT(){
2275 m_nhit_Pxl_Ist = QAH::H2F(
"QaPxlvsIstHit",
"PIXEL hits vs IST hits",500,0.,20000.,100,0.,5000.);
2276 m_nhit_Pxl_Ist->SetXTitle(
"PIXEL hits");
2277 m_nhit_Pxl_Ist->SetYTitle(
"IST hits");
2279 m_nhit_Pxl_Sst = QAH::H2F(
"QaPxlvsSstHit",
"PIXEL hits vs SST hits",500,0.,20000.,100,0.,5000.);
2283 m_nhit_Ist_Sst = QAH::H2F(
"QaIstvsSsdHit",
"IST hits vs SST hits",100,0.,5000.,100,0.,5000.);
2289 void StQABookHist::BookHistPXL(){
2291 m_pxl_hit_phi_z_Pxl1= QAH::H2F(
"QaPxlHitPhivsZPxl1",
"PIXEL: hits vs phi vs z in inner layer (per event)",20,-10.,10.,100,-TMath::Pi(),TMath::Pi());
2294 m_pxl_hit_phi_z_Pxl2= QAH::H2F(
"QaPxlHitPhivsZPxl2",
"PIXEL: hits vs phi vs z in outer layer (per event)",20,-10.,10.,100,-TMath::Pi(),TMath::Pi());
2298 m_pxl_hit_ladder = QAH::H1F(
"QaPxlHitvsLadder",
"PIXEL: hits per ladder (per event)",5,0.,5.);
2301 m_pxl_hit_sector_sensor_Pxl1 = QAH::H2F(
"QaPxlHitvsSectorvsSensorPxl1",
"PIXEL: hits vs sector vs sensor in inner layer (per event)",10,0.5,10.5,10,0.5,10.5);
2306 m_pxl_hit_sector_sensor_Pxl2 = QAH::H2F(
"QaPxlHitvsSectorvsSensorPxl2",
"PIXEL: hits vs ladder vs sensor in outer layer (per event)",40,0.5,40.5,10,0.5,10.5);
2311 m_pxl_nhit_Pxl1_tpc_mult = QAH::H2F(
"QaPxl1HitvsTpcMult",
"PIXEL: Hits in inner layer vs TPC multiplicity",100,0.,12000.,500,0.,10000.);
2315 m_pxl_nhit_Pxl2_tpc_mult = QAH::H2F(
"QaPxl2HitvsTpcMult",
"PIXEL: Hits in outer layer vs TPC multiplicity",100,0.,8000.,500,0.,10000.);
2319 m_pxl_nhit_Pxl1_tof_mult = QAH::H2F(
"QaPxl1HitvsTofMult",
"PIXEL: Hits in inner layer vs TOF multiplicity",100,0.,12000.,100,0.,5000.);
2323 m_pxl_nhit_Pxl2_tof_mult = QAH::H2F(
"QaPxl2HitvsTofMult",
"PIXEL: Hits in outer layer vs TOF multiplicity",100,0.,8000.,100,0.,5000.);
2327 m_pxl_nhit_Pxl1_Pxl2 = QAH::H2F(
"QaPxlHitLayer1vsLayer2",
"PIXEL: Hits in inner vs outer layer (per event)",100,0.,12000.,100,0.,8000.);
2333 void StQABookHist::BookHistIST(){
2336 m_ist_hit_phi_z = QAH::H2F(
"QaIstHitPhivsZ",
"IST: Hits vs phi vs z (per event)",60,-30.,30.,100,-TMath::Pi(),TMath::Pi());
2340 m_ist_hit_ladder = QAH::H1F(
"QaIstHitvsLadder",
"IST: Hits per ladder (per event)",25,0.,25.);
2343 m_ist_hit_ladder_sensor = QAH::H2F(
"QaIstHitvsLaddervsSensor",
"IST: Hits vs ladder vs sensor (per event)", 24,0.5,24.5, 6, 0.5, 6.5);
2349 m_ist_nhit_tpc_mult = QAH::H2F(
"QaIstHitvsTpcMult",
"IST hits vs TPC multiplicity",100,0.,5000.,500,0.,10000.);
2353 m_ist_nhit_tof_mult = QAH::H2F(
"QaIstHitvsTofMult",
"IST hits vs TOF multiplicity",100,0.,5000.,100,0.,5000.);
2359 void StQABookHist::BookHistMTD(){
2361 m_MtdNHits = QAH::H1F(
"QaMtdNHits",
"Number of MTD hits per event;N",50,0,50);
2362 m_MtdHitMap = QAH::H2F(
"QaMtdHitMap",
"MTD: channel vs backleg of hits;backleg;channel",30,0.5,30.5,60,-0.5,59.5);
2363 m_MtdNMatchHits = QAH::H1F(
"QaMtdNMatchHits",
"Number of matched MTD hits per event;N",20,0,20);
2364 m_MtdMatchHitMap = QAH::H2F(
"QaMtdMatchHitMap",
"MTD: channel vs backleg of matched hits;backleg;channel",30,0.5,30.5,60,-0.5,59.5);
2368 void StQABookHist::BookHistEPD(){
2370 TString name, title;
2371 const TString side[2] = {
"East",
"West"};
2374 for (
int i=0; i<2; i++){
2375 for (
int j=0; j<nPP; j++){
2376 (((name =
"QaEpdAdc") += side[i]) +=
"PP") += (j+1);
2377 ((((title =
"EPD ") += side[i]) +=
" PP") += (j+1)) +=
" ADC";
2378 m_epd_adc[(i*nPP)+j] = QAH::H2F(name,title,31,1,31,4096,0,4095);
2379 m_epd_adc[(i*nPP)+j]->SetXTitle(
"Tile");
2380 m_epd_adc[(i*nPP)+j]->SetYTitle(
"ADC");
2382 (((name =
"QaEpdTac") += side[i]) +=
"PP") += (j+1);
2383 ((((title =
"EPD ") += side[i]) +=
" PP") += (j+1)) +=
" TAC";
2384 m_epd_tac[(i*nPP)+j] = QAH::H2F(name,title,9,1,9,4096,0,4095);
2385 m_epd_tac[(i*nPP)+j]->SetXTitle(
"Tile");
2386 m_epd_tac[(i*nPP)+j]->SetYTitle(
"TAC");
2392 void StQABookHist::BookHistiTPC(){
2393 char nameiTPC[60], titleiTPC[60];
2394 for(
int i = 0; i < 24; i++)
2396 sprintf(nameiTPC,
"TPC_charge_time_inner_s%d", i+1);
2397 sprintf(titleiTPC,
"TPC adc vs time-bin inner sector %d", i+1);
2398 m_TPC_ch_time_inner[i] = QAH::H2F(nameiTPC, titleiTPC, 50, 0, 511, 250 , 0, 500);
2399 m_TPC_ch_time_inner[i] ->SetXTitle(
"time bin");
2400 m_TPC_ch_time_inner[i] ->SetYTitle(
"adc");
2403 for(
int i = 0; i < 24; i++)
2405 sprintf(nameiTPC,
"TPC_charge_time_outer_s%d", i+1);
2406 sprintf(titleiTPC,
"TPC adc vs time-bin outer sector %d", i+1);
2407 m_TPC_ch_time_outer[i] = QAH::H2F(nameiTPC, titleiTPC, 50, 0, 511, 250 , 0, 500);
2408 m_TPC_ch_time_outer[i] ->SetXTitle(
"time bin");
2409 m_TPC_ch_time_outer[i] ->SetYTitle(
"adc");
2412 for(
int i = 0; i < 24; i++)
2414 sprintf(nameiTPC,
"TPC_charge_nrow_s%d", i+1);
2415 sprintf(titleiTPC,
"TPC adc vs pad-row sector %d", i+1);
2416 m_TPC_ch_nrow[i] = QAH::H2F(nameiTPC, titleiTPC, 72 ,0.5, 72.5, 250 , 0, 500);
2417 m_TPC_ch_nrow[i] ->SetXTitle(
"pad row");
2418 m_TPC_ch_nrow[i] ->SetYTitle(
"adc");
2421 sprintf(nameiTPC,
"TPC_adc_sec_inner");
2422 sprintf(titleiTPC,
"TPC adc vs sector inner");
2423 m_TPC_adc_sec_inner = QAH::H2F(nameiTPC, titleiTPC, 24 ,0.5, 24.5, 250 , 0, 500);
2424 m_TPC_adc_sec_inner ->SetXTitle(
"sector");
2425 m_TPC_adc_sec_inner ->SetYTitle(
"adc");
2427 sprintf(nameiTPC,
"TPC_adc_sec_outer");
2428 sprintf(titleiTPC,
"TPC adc vs sector outer");
2429 m_TPC_adc_sec_outer = QAH::H2F(nameiTPC, titleiTPC, 24 ,0.5, 24.5, 250 , 0, 500);
2430 m_TPC_adc_sec_outer ->SetXTitle(
"sector");
2431 m_TPC_adc_sec_outer ->SetYTitle(
"adc");
TH2F * m_glb_xfF
x-coord. of first hit on trk, tpc+svt
TH2F * m_npoint_lengthTS
tot num points vs length, tpc
TH1F * m_glb_impactrTS
log impact parameter from primary vertex, ftpc
TH1F * m_lengthFE
length of track, ftpc
TH1F * m_glb_yfFW
y-coord. of first hit on trk, ftpc east
TH3F * m_pnt_svtLaserDiff
laser spots in svt
TH1F * m_ndedxFW
dE/dx [0], ftpcE
TH1F * m_plengthFE
length of track, ftpc
TH2F * m_primtrk_meanptTTS
detector id of track
TH1F * m_chisq0T
length of track, ftpc west
TH2F * m_pT_eta_recT
overlayed hist of first point - helix point
TH1F * m_primglob_good
tracks in table with iflag>0, ftpc
TH1F * m_glb_xf0TS
x-coord. of first hit - at start of helix
TH2F * m_globtrk_padtimeFE
Y vs X of first hit on trk, ftpc west.
TH1F * m_vtx_FtpcWestTpc_z
FtpcWest prim vtx - TPC prim vtx, x vs y.
TH1F * m_pv_x
row1-vertex type
TH2F * m_pmom_trklengthT
tanl(dip angle) vs zfirst-zvtx, tpc+svt
TH1F * m_ndedxr
psi vs phi, tpc+svt
TH2F * m_pchisq0_zfTS
chisq0 vs zfirst, tpc
TH2F * m_tof_hit_module
of hits vs tray # (int over modules)
TH1F * m_dcaToBeamZ1
xy-dca to beam axis
TH1F * m_geant_reco_pvtx_z
prim vtx y, diff geant - reco
TH1F * m_prim_zfT
y-coord. of first hit on trk, ftpc west
TH1F * m_global_ist_hit
of hits in IST vs TOF track multiplicity
TH2F * m_glb_impactF
log impact parameter from primary vertex, tpc,tpc+svt
TH1F * m_vtx_FtpcEastTpc_z
FtpcEast prim vtx - TPC prim vtx, x vs y.
TH1F * m_glb_chargeFW
particle charge in units of |e| - ftpc east
TH1F * m_ppsiT
radial (xy) coordinate of first hit, ftpc west
TH2F * m_prim_f0
chi square [1], ftpc west
TH2F * m_primtrk_xf_yfTW
Y vs X of first hit on trk, tpc east.
TH1F * m_glb_xf0
plane of first hit on trk, ftpc
TH1F * m_dedx0FW
number of point to find dE/dx, ftpcW
TH2F * m_pTTTS
pT, tpc+svt
TH1F * m_pnt_tot_med
number of hits total
TH1F * m_petaT
theta - tpc+svt
TH1F * m_pchisq0T
length of track, ftpc west
TH1F * m_pnt_tot_sm
number of hits total, med range
TH2F * m_chisq0_zfTS
chisq0 vs zfirst, tpc
TH2F * m_peta_trklengthT
mom vs. trk length, tpc+svt
TH1F * m_ppointFE
number of points on the track - ftpc
TH2F * m_glb_rzf0
z-coord. of first hit - at start of helix+svt
TH2F * m_primtrk_meanetaTTS
<pT>, ftpc
TH2F * m_petaF
eta, tpc,tpc+svt
TH1F * m_prim_zfFE
z-coord. of first hit on trk, ftpc
TH2F * m_pnt_timeT
z dist. of hits, svt
TH1F * m_pTFW
pT, ftpc east
TH2F * m_ist_nhit_tpc_mult
of hits vs ladder vs sensor
TH1F * m_global_hft_hit
of hits in ITS vs hits in ITS
TH2F * m_pfpoint_lengthTS
num fit points vs length, tpc
TH1F * m_primary_hft_hit
HFT hits per global track.
TH1F * m_glb_yfT
x-coord. of first hit on trk, ftpc west
TH2F * m_pnt_xyFE
barrel dist. of hits, svt
TH1F * m_prim_xfFW
x-coord. of first hit on trk, ftpc east
TH2F * m_chisq0_phiT
chisq0 vs zfirst, tpc+svt
TH2F * m_pmomF
momentum, tpc+svt
TH2F * m_glb_impactT
chi square [1], ftpc west
TH2F * m_chisq0_momTS
chisq0 vs momentum, tpc
TH1F * m_prim_xf0TS
x-coord. of first hit - at start of helix
TH2F * m_primtrk_xf_yfTE
pT versus eta, ftpcW
TH1F * m_prim_radfFE
radial (xy) coordinate of first hit, ftpc
TH1F * m_glb_radfT
transverse & longitudinal helix center of last hit - helix hit+svt
TH2F * m_pT_eta_recTS
pT versus eta, tpc
TH1F * m_primtrk_good_sm
tracks in table with iflag>0
TH1F * m_prim_radfFW
radial (xy) coordinate of first hit, ftpc east
TH1F * m_prim_xfFE
x-coord. of first hit on trk, ftpc
TH2F * m_pT_eta_recFE
pT versus eta, tpc+svt
TH1F * m_pTT
momentum, ftpc west
TH1F * m_prim_radfTS
radial (xy) coordinate of first hit, tpc
TH2F * m_prim_rzf0TS
transverse & longitudinal helix center of last hit - helix hit
TH2F * m_ppTF
pT, tpc,tpc+svt
TH1F * m_global_pxl_hit
of hits in inner vs outer PIXEL layer
TH1F * m_v_z
STAR reference.
TH2F * m_glb_impactTTS
signed impact parameter from primary vertex, tpc+svt
TH1F * m_ppointFW
number of points on the track - ftpc east
TH2F * m_tof_match_module
of matched hits vs tray
TH2F * m_npoint_lengthT
num fit points vs length, tpc,tpc+svt
TH2F * m_tof_hit_tray
of IST hits per primary track
TH1F * m_dcaToBeamZ2
z-dca to beam axis -210 to -105
TH2F * m_ppT_eta_recTS
pT versus eta, tpc
TH2F * m_chisq0_phiTS
chisq0 vs phi, tpc
TH1F * m_prim_z0T
azimuthal angle at start (deg), tpc+svt
TH1F * m_emc_point_sigeta[4]
Point Phi spectra.
TH1F * m_chisq1FE
chi square [1], ftpc
TH2F * m_ppsi_phiT
chisq0 vs zfirst, tpc+svt
TH1F * m_glb_thetaT
tan(dip) =pz/pt at start, tpc+svt
TH2F * m_prim_ratiomF
ratio of n fit pnts over max n pnts - tpc+svt
TH1F * m_pchisq0FW
chi square [0], ftpc east
TH1F * m_v0
z distribution of V0s relative to primVtx
TH1F * m_pmax_pointTS
number of max possible track points - tpc
TH1F * m_v_num_sm
number of vertices
TH1F * m_prim_yf0TS
y-coord. of first hit - at start of helix
TH2F * m_pchisq0_dipTS
chisq0 vs dip angle, tpc
TH2F * m_mom_trklengthTS
mom vs. trk length, tpc
TH1F * m_chisq0FE
chi square [0], ftpc
TH2F * m_ppT_eta_recFW
pT versus eta, ftpcE
TH2F * m_chisq0F
chi square [0], tpc,svt
TH2F * m_eta_trklengthFE
trk length vs. eta, tpc+svt
TH1F * m_dedx1T
dE/dx [0], tpc
TH2F * m_pxl_nhit_Pxl1_Pxl2
of hits in outer PIXEL layer vs TOF track multiplicity
TH2F * m_globtrk_xf_yfTE
pT versus eta, ftpcW
TH2F * m_pchisq1F
chi square [0], ftpc west
TH1F * m_emc_point_delphi[4]
Point DeltaEta spectra.
TH1F * m_ist_hit_ladder
of hits vs phi vs z position
TH1F * m_psiTS
psi reconstructed, tpc
TH1F * m_pnt_tpc
detector ID of the hit
TH1F * m_vtx_phi_dist
radius to primary vertex
TH2F * m_globtrk_padtimeFW
padlength vs timelength of hits on track, ftpcE
TH2F * m_pnt_rpTE
rphi dist. of hits, tpcW
TH1F * m_primtrk_tot_sm
tracks in table
TH1F * m_v_y
vertex coordinates in
TH2F * m_pfpoint_lengthTTS
num fit points vs length, tpc+svt
TH2F * m_petaTTS
eta, tpc+svt
TH2F * m_tof_vtx_z
of tof hits vs # of vpd hits
TH2F * m_fit_pointTTS
number of track points used for fitting - tpc+svt
TH2F * m_primtrk_meanptF
<pT>, tpc, tpc+svt
TH2F * m_chisq0_dipT
chisq0 vs eta, tpc+svt
TH1F * m_glb_xfFW
x-coord. of first hit on trk, ftpc east
TH2F * m_ptanl_zfTS
tanl(dip angle) vs zfirst-zvtx, tpc
TH2F * m_geant_reco_vtx_z_z
prim vtx z, diff geant - reco
TH1F * m_pchisq0TS
chi square [0], tpc
TH1F * m_pv_y
row1-vertex coordinates in
TH2F * m_chisq0_dipTS
chisq0 vs dip angle, tpc
TH1F * m_prim_zf0
y-coord. of first hit - at start of helix+svt
TH1F * m_pfit_pointTS
number of track points used for fitting - tpc
TH1F * m_prim_radfT
transverse & longitudinal helix center of last hit - helix hit+svt
TH1F * m_glb_zfFW
z-coord. of first hit on trk, ftpc east
TH2F * m_vtx_FtpcEastTpc_xy
number of kinks
TH1F * m_glb_zf0TS
z-coord. of first hit - at start of helix
TH1F * m_dcaToBeamZ3
z-dca to beam axis -105 to 0
TH2F * m_fpoint_lengthT
trk length vs. eta, ftpc west
TH1F * m_pv_pchi2
x versus y
TH1F * m_glb_phi0TS
azimuthal angle at start (deg), tpc
TH1F * m_pfit_pointT
number of max possible track points - ftpc west
TH1F * m_max_pointFW
number of max possible track points - ftpc east
TH2F * m_glb_radfF
impact parameter from primary vertex, ftpc
TH1F * m_prim_chargeFW
particle charge in units of |e| - ftpc east
TH1F * m_pointFW
number of points on the track - ftpc east
TH1F * m_pnt_zS
z dist. of hits, tpc
TH2F * m_pxl_nhit_Pxl1_tof_mult
of hits in outer PIXEL layer vs TPC track multiplicity
TH2F * m_chisq0_momT
tot num points vs length, ftpc west
TH1F * m_chisq1FW
chi square [1], ftpc east
TH1F * m_glb_ratiomT
ratio of n fit pnts over max n pnts - tpc,tpc+svt
TH2F * m_p_dedx_rec
psi vs phi, tpc+svt
TH1F * m_glb_padfT
z-coord at start (cm), tpc+svt
TH1F * m_ppTT
momentum, ftpc west
TH2F * m_tof_PID
vertex z from vpd vs verex z from TPC
TH2F * m_glb_simpactT
log impact parameter from primary vertex, tpc
TH2F * m_pchisq0_dipT
chisq0 vs eta, tpc+svt
TH1F * m_globtrk_fit_prob
TH2F * m_ppTTTS
pT, tpc+svt
TH2F * m_glb_padfTEW
padrow of first hit on trk, tpc
TH2F * m_glb_impactTS
chi square [0], tpc+svt
TH2F * m_prim_radfF
radial (xy) coordinate of first hit, tpc+svt
TH1F * m_pnt_barrelS
padrow dist. of hits, tpc
TH1F * m_ppsiFW
psi reconstructed, ftpc east
TH2F * m_ist_hit_ladder_sensor
of hits vs x ladder (integrated over ladder, per event)
TH1F * m_prim_yfFE
y-coord. of first hit on trk, ftpc
TH2F * m_pnt_planeF
padlength vs timelength of hits, ftpcW
TH1F * m_prim_curvTS
curvature (1/cm), tpc
TH1F * m_prim_zfTS
z-coord. of first hit on trk, tpc
TH1F * m_prim_zfFW
z-coord. of first hit on trk, ftpc east
TH2F * m_glb_chargeF
particle charge in units of |e| - tpc+svt
TH1F * m_ppointT
<eta>, ftpc
TH1F * m_glb_xfT
curvature (1/cm), tpc+svt
TH1F * m_prim_chargeFE
particle charge in units of |e| - ftpc
TH1F * m_glb_phifTS
phi dist. of first point on trk, tpc
TH2F * m_ppT_eta_recFE
pT versus eta, tpc+svt
TH2F * m_vtx_FtpcWestTpc_xy
FtpcEast prim vtx z - TPC prim vtx z.
TH2F * m_glb_rzf0TS
transverse & longitudinal helix center of last hit - helix hit
TH2F * m_glb_impactrTTS
impact parameter from primary vertex, tpc+svt
TH2F * m_glb_zfF
z-coord. of first hit on trk, tpc+svt
TH1F * m_pointFE
number of points on the track - ftpc
TH2F * m_eta_trklengthTS
trk length vs. eta, tpc
TH1F * m_globtrk_goodTTS
tracks in table with iflag>0,small range
TH1F * m_psiFE
psi reconstructed, ftpc
TH1F * m_glb_phi0T
radius at start (cm), tpc+svt
TH2F * m_eta_trklengthT
mom vs. trk length, tpc+svt
TH1F * m_glb_xfFE
x-coord. of first hit on trk, ftpc
TH2F * m_eta_trklengthFW
trk length vs. eta, ftpc east
TH2F * m_pnt_padtimeFW
padlength vs timelength of hits, ftpcE
TH1F * m_prim_ratiomT
ratio of n fit pnts over tot n pnts - ftpc west
TH1F * m_prim_xf0
z-coord. of first hit on trk, ftpc west
TH1F * m_pnt_ssd
SVT drift.
TH2F * m_pchisq0F
chi square [0], tpc,tpc+svt
TH2F * m_pointF
number of points on the track - tpc+svt
TH1F * m_pmax_pointFE
number of max possible track points - ftpc
TH2F * m_glb_f0
impact parameter from primary vertex, tpc
TH2F * m_globtrk_xf_yfFE
Y vs X of first hit on trk, tpc+svt.
TH1F * m_glb_xfTS
x-coord. of first hit on trk, tpc
TH1F * m_prim_zf0TS
z-coord. of first hit - at start of helix
TH1F * m_dedx0T
number of point to find dE/dx, tpc
TH1F * m_pmomTS
momentum, tpc
TH2F * m_primtrk_goodF
tracks in table with iflag>0, tpc,svt
TH2F * m_ptanl_zfT
Y vs X of first hit on trk, ftpc west.
TH1F * m_pnt_phiS
phi dist. of hits, tpc
TH1F * m_glb_radfFW
radial (xy) coordinate of first hit, ftpc east
TH1F * m_prim_xfTS
x-coord. of first hit on trk, tpc
TH1F * m_glb_yf0
x-coord. of first hit - at start of helix+svt
TH1F * m_glb_impactrT
signed impact parameter from primary vertex, tpc
TH1F * m_ndedxFE
dE/dx [0], ftpc
TH2F * m_prim_f0TS
overlayed hist of first point - helix point
TH1F * m_momT
eta, ftpc west
TH1F * m_glb_chargeTS
particle charge in units of |e| - tpc
TH1F * m_ev0_k0ma_hist
Lambda mass.
TH2F * m_glb_rzl0
transverse & longitudinal helix center of first hit - helix hit
TH2F * m_peta_trklengthTS
trk length vs. eta, tpc
TH2F * m_pxl_nhit_Pxl2_tof_mult
of hits in inner PIXEL layer vs TOF track multiplicity
TH1F * m_chisq0TS
length of track, tpc+svt
TH2F * m_primtrk_xf_yfFW
Y vs X of first hit on trk, ftpc east.
TH1F * m_globtrk_iflag
tracks in table with iflag>0, ftpc
TH1F * m_pmomT
eta, ftpc west
TH2F * m_zDcaZf
z-dca vs tanl
TH1F * m_petaFW
eta, ftpc east
TH1F * m_glb_radfTS
radial (xy) coordinate of first hit, tpc
TH1F * m_emc_points[4]
Point TrMom spectra.
TH2F * m_ppsiTTS
psi reconstructed, tpc+svt
TH1F * m_primtrk_goodTTS
tracks in table with iflag>0, small range
TH1F * m_vtx_z_dist
radial distribution of V0s relative to primVtx
TH1F * m_pmax_pointT
number of points on the track - ftpc west
TH2F * m_tof_match_tray
of hits vs module # (int over trays, east side use -32 - -1)
TH1F * m_tanlTS
tan(dip) =pz/pt at start, tpc
TH2F * m_chisq0_etaTS
chisq0 vs eta, tpc
TH2F * m_fpoint_lengthTS
num fit points vs length, tpc
TH2F * m_prim_yfF
y-coord. of first hit on trk, tpc+svt
TH2F * m_pnt_ftpc
number of hits tpc
TH1F * m_prim_thetaT
tan(dip) =pz/pt at start, tpc+svt
TH1F * m_max_pointFE
number of max possible track points - ftpc
TH2F * m_glb_f0TS
overlayed hist of first point - helix point
TH1F * m_glb_zf0
y-coord. of first hit - at start of helix+svt
TH2F * m_glb_impactrF
impact parameter from primary vertex, tpc,tpc+svt
TH1F * m_det_id
iflag value
TH1F * m_pnt_phiSSD
number of hits ssd
TH2F * m_pnt_xyFW
xy dist. of hits, ftpcE
TH1F * m_pv_z
STAR reference.
TH2F * m_prim_rzl0TS
transverse & longitudinal helix center of first hit - helix hit+svt
TH2F * m_dcaToBeamXY
detector id of track
TH1F * m_pv_vtxid
radius to vertex
TH1F * m_pointTS
number of points on the track - tpc
TH2F * m_psiTTS
psi reconstructed, tpc+svt
TH1F * m_glb_yfTS
y-coord. of first hit on trk, tpc
TH2F * m_prim_chargeF
particle charge in units of |e| - tpc+svt
TH2F * m_zDcaTanl
z-dca to beam axis 0 to 105
TH1F * m_etaFW
eta, ftpc east
TH2F * m_glb_phifT
radial (xy) coordinate of first hit, tpc+svt
TH2F * m_chisq0_zfT
chisq0 vs dip angle, tpc+svt
TH2F * m_pT_eta_recFW
pT versus eta, ftpcE
TH1F * m_pmomFE
momentum, ftpc
TH2F * m_prim_rzl0
transverse & longitudinal helix center of first hit - helix hit
TH1F * m_momTS
momentum, tpc
TH1F * m_pchisq1FE
chi square [1], ftpc
TH1F * m_glb_ratiomFE
ratio of n fit pnts over max n pnts - ftpc
TH1F * m_glb_curvTS
curvature (1/cm), tpc
TH1F * m_plengthFW
length of track, ftpc east
TH1F * m_momFW
momentum, ftpc east
TH1F * m_ev0_lama_hist
v0 vertices
TH2F * m_pchisq0_etaTS
chisq0 vs eta, tpc
TH2F * m_primtrk_xf_yfTS
Y vs X of first hit on trk, tpc west.
TH2F * m_pxl_hit_phi_z_Pxl1
HFT hits per primary track.
TH1F * m_plengthTS
length of track, tpc
TH2F * m_pxl_hit_sector_sensor_Pxl1
of hits vs x ladder (integrated over ladder, per event)
TH1F * m_lengthFW
length of track, ftpc east
TH1F * m_primary_pxl_hit
of PIXEL hits per global track
TH1F * m_glb_curvT
padrow of first hit on trk, tpc east and west
TH1F * m_pnt_tot
plane dist. of hits, ftpc
TH1F * m_prim_chargeT
ratio of n fit pnts over max n pnts - ftpc west
TH2F * m_globtrk_goodF
tracks in table with iflag>0, tpc, tpc+svt
TH1F * m_primtrk_good
tracks in table, small range
TH2F * m_pnpoint_lengthTS
tot num points vs length, tpc
TH2F * m_nhit_Pxl_Sst
of hits in PIXEL vs hits in IST
TH1F * m_glb_ratiomFW
ratio of n fit pnts over max n pnts - ftpc east
TH1F * m_prim_yfTS
y-coord. of first hit on trk, tpc
TH2F * m_glb_yfF
y-coord. of first hit on trk, tpc+svt
TH2F * m_pnpoint_lengthT
num fit points vs length, tpc
TH2F * m_svt_loc
diff of laser spots in svt
TH2F * m_pnpoint_lengthFE
tot num points vs length, tpc+svt
TH1F * m_pdet_id
iflag value
TH2F * m_glb_rzl0TS
transverse & longitudinal helix center of first hit - helix hit+svt
TH1F * m_ppointTS
number of points on the track - tpc
TH1F * m_glb_r0TS
radius at start (cm), tpc
TH1F * m_pchisq1FW
chi square [1], ftpc east
TH1F * m_primtrk_tot
dE/dx [0], ftpcW
TH2F * m_nhit_Ist_Sst
of hits in PIXEL vs hits in SST
TH1F * m_prim_r0TS
radius at start (cm), tpc
TH2F * m_pchisq0_etaT
chisq0 vs momentum, tpc+svt
TH1F * m_glb_yfFE
y-coord. of first hit on trk, ftpc
TH2F * m_ppsi_phiTS
psi vs phi, tpc
TH1F * m_emc_point_trmom[4]
Point DeltaPhi spectra.
TH1F * m_pmomFW
momentum, ftpc east
TH1F * m_pointT
z-dca vs phi0
TH1F * m_ptanlTS
tan(dip) =pz/pt at start, tpc
TH1F * m_glb_zfT
y-coord. of first hit on trk, ftpc west
TH2F * m_prim_zfF
z-coord. of first hit on trk, tpc+svt
TH2F * m_pnt_xyS
time dist. of hits, tpc
TH2F * m_glb_planefF
z-coord. of first hit on trk, ftpc west
TH2F * m_pxl_nhit_Pxl2_tpc_mult
of hits in inner PIXEL layer vs TPC track multiplicity
TH2F * m_lengthF
length of track, tpc
TH2F * m_ppsiF
psi reconstructed, tpc,tpc+svt
TH1F * m_glb_thetaTS
theta - tpc
TH2F * m_ppT_eta_recT
overlayed hist of first point - helix point
TH2F * m_chisq0_etaT
chisq0 vs momentum, tpc+svt
TH2F * m_etaF
eta, tpc,tpc+svt
TH2F * m_etaTTS
eta, tpc+svt
TH2F * m_chisq0TTS
chi square [0], tpc
TH1F * m_xi_ma_hist
number of xi vertices
TH1F * m_prim_yfFW
y-coord. of first hit on trk, ftpc east
TH1F * m_ppTFW
pT, ftpc east
TH2F * m_plengthF
length of track, tpc+svt
TH2F * m_pmax_pointF
number of max possible track points - tpc+svt
TH2F * m_tanl_zfTS
tanl(dip angle) vs zfirst-zvtx, tpc
TH2F * m_pchisq0_momT
tot num points vs length, ftpc west
TH2F * m_pmom_trklengthTS
mom vs. trk length, tpc
TH1F * m_ptanlT
psi reconstructed, ftpc west
TH2F * m_psi_phiT
chisq0 vs phi, tpc+svt
TH1F * m_emc_point_sigphi[4]
Point SigmaEta spectra.
TH2F * m_dedx0F
number of point to find dE/dx, ftpc
TH2F * m_pchisq0_zfT
chisq0 vs dip angle, tpc+svt
TH1F * m_lengthTS
phi dist. of first point on trk, tpc+svt
TH1F * m_prim_phi0TS
azimuthal angle at start (deg), tpc
TH1F * m_pnt_ftpcE
number of hits ftpc
TH1F * m_psiT
radial (xy) coordinate of first hit, ftpc west
TH1F * m_primtrk_iflag
fit points ratio prim/glob, all detectors
TH2F * m_max_pointF
number of max possible track points - tpc+svt
TH1F * m_prim_yf0
x-coord. of first hit - at start of helix+svt
TH1F * m_prim_curvT
z-coord at start (cm), tpc+svt
TH2F * m_glb_ratiomF
ratio of n fit pnts over max n pnts - tpc+svt
TH1F * m_ndedxT
number of tracks with dedx info
TH2F * m_primtrk_meanetaF
<eta>, tpc, tpc+svt
TH1F * m_globtrk_good_sm
tracks in table with iflag>0
TH2F * m_prim_rzf0
z-coord. of first hit - at start of helix+svt
TH2F * m_pnt_svtLaser
number of hits svt
TH2F * m_tof_vpd_hit
of matched hits vs module
TH1F * m_lengthT
pT, ftpc west
TH1F * m_dedxTTS
dE/dx [1], tpc
TH2F * m_fpoint_lengthTTS
num fit points vs length, tpc+svt
TH1F * m_glb_chargeT
ratio of n fit pnts over max n pnts - ftpc west
TH1F * m_prim_ratiomTS
ratio of n fit pnts over max n pnts - tpc
TH1F * m_plengthT
pT, ftpc west
TH1F * m_prim_ratiomFE
ratio of n fit pnts over max n pnts - ftpc
TH1F * m_dedx0FE
number of point to find dE/dx, ftpcE
TH2F * m_tanl_zfT
padlength vs timelength of hits on track, ftpcW
TH2F * m_ist_hit_phi_z
of PIXEL hits per primary track
TH2F * m_emc_nhit
prim vtx z, diff geant - reco vs reco z
TH1F * m_glb_zfFE
z-coord. of first hit on trk, ftpc
TH2F * m_bbc_adc[4]
Point Flag spectra.
TH1F * m_emc_point_deleta[4]
Point SigmaPhi spectra.
TH1F * m_pmax_pointFW
number of max possible track points - ftpc east
TH1F * m_glb_z0T
azimuthal angle at start (deg), tpc+svt
TH2F * m_glb_ratiomTTS
number of track pnts on svt - tpc+svt
TH2F * m_pxl_hit_sector_sensor_Pxl2
of hits vs sector vs sensor, inner PIXEL layer
TH1F * m_glb_zfTS
z-coord. of first hit on trk, tpc
TH1F * m_glb_r0T
particle charge in units of |e| - ftpc west
TH1F * m_z_hits
FtpcWest prim vtx z - TPC prim vtx z.
TH1F * m_primglob_fit
trks w/ iflag>0, ratio prim/glob, all detectors
TH2F * m_globtrk_xf_yfTW
Y vs X of first hit on trk, tpc east.
TH1F * m_prim_r0T
particle charge in units of |e| - ftpc west
TH1F * m_glb_chargeFE
particle charge in units of |e| - ftpc
TH2F * m_pnt_phiT
rphi dist. of hits, tpcE
TH1F * m_emc_point_phi[4]
Point Eta spectra.
TH1F * m_pnt_id
number of hits total, small range
TH1F * m_max_pointT
number of points on the track - ftpc west
TH1F * m_MtdNHits
TOF PID: 1/beta vs p.
TH1F * m_tanlT
psi reconstructed, ftpc west
TH2F * m_mom_trklengthT
tanl(dip angle) vs zfirst-zvtx, tpc+svt
TH2F * m_npoint_lengthFE
tot num points vs length, tpc+svt
TH2F * m_pnt_padrowT
phi dist. of hits, svt
TH1F * m_geant_reco_pvtx_y
prim vtx x, diff geant - reco
TH1F * m_pchisq0FE
chi square [0], ftpc
TH1F * m_momFE
momentum, ftpc
TH1F * m_fit_pointTS
number of track points used for fitting - tpc
TH2F * m_peta_trklengthFE
trk length vs. eta, tpc+svt
TH1F * m_v_vtxid
number of vertices,small range
TH1F * m_glb_radfFE
radial (xy) coordinate of first hit, ftpc
TH2F * m_npoint_lengthFW
tot num points vs length, ftpc east
TH1F * m_fit_pointT
number of max possible track points - ftpc west
TH1F * m_pnt_phiSST
number of hits sst
TH2F * m_globtrk_xf_yfTS
Y vs X of first hit on trk, tpc west.
TH2F * m_pTF
pT, tpc,tpc+svt
TH2F * m_prim_xfF
x-coord. of first hit on trk, tpc+svt
TH2F * m_ppointF
number of points on the track - tpc+svt
TH1F * m_ppsiFE
psi reconstructed, ftpc
TH2F * m_ist_nhit_tof_mult
of hits in IST vs TPC track multiplicity
TH1F * m_ppsiTS
psi reconstructed, tpc
TH1F * m_max_pointTS
number of max possible track points - tpc
TH1F * m_prim_chargeTS
particle charge in units of |e| - tpc
TH2F * m_zDcaPhi0
z-dca vs psi
TH2F * m_pchisq0_momTS
chisq0 vs momentum, tpc
TH1F * m_prim_thetaTS
theta - tpc
TH2F * m_psiF
psi reconstructed, tpc,tpc+svt
TH1F * m_etaT
theta - tpc+svt
TH1F * m_prim_yfT
x-coord. of first hit on trk, ftpc west
TH2F * m_glb_simpactTS
log impact parameter from primary vertex, tpc+svt
TH2F * m_pfpoint_lengthT
trk length vs. eta, ftpc west
TH1F * m_pnt_ftpcW
number of hits ftpcE
TH2F * m_psi_phiTS
psi vs phi, tpc
TH2F * m_zDcaPsi
z-dca vs z-first
TH2F * m_pnt_padtimeFE
xy dist. of hits, ftpcW
TH1F * m_glb_z0TS
z-coord at start (cm), tpc
TH1F * m_emc_point_eta[4]
Point Energy spectra.
TH2F * m_peta_trklengthFW
trk length vs. eta, ftpc east
TH1F * m_vtx_r_dist
azimuthal distribution of V0s relative to primVtx
TH1F * m_pv_r
row1-chisq per dof of vertex fit
TH1F * m_prim_z0TS
z-coord at start (cm), tpc
TH2F * m_pnpoint_lengthFW
tot num points vs length, ftpc east
TH2F * m_globtrk_xf_yfFW
Y vs X of first hit on trk, ftpc east.
TH1F * m_prim_ratiomFW
ratio of n fit pnts over max n pnts - ftpc east
TH2F * m_pchisq0TTS
chi square [0], tpc+svt
TH1F * m_chisq0FW
chi square [0], ftpc east
TH1F * m_primary_ist_hit
of IST hits per global track
TH1F * m_globtrk_good
tracks in table
TH1F * m_pxl_hit_ladder
of hits vs phi vs z position, outer PIXEL layer
TH1F * m_psiFW
psi reconstructed, ftpc east
TH2F * m_pxl_hit_phi_z_Pxl2
of hits vs phi vs z position, inner PIXEL layer
TH2F * m_ndedxF
<dE/dx>/(Bichsel <dE/dx>), tpc,tpc+svt
TH2F * m_pnt_rpTW
xy dist. of hits, svt
TH1F * m_emc_point_flag
Emc Point multiplicity.
TH1F * m_glb_ratiomTS
ratio of n fit pnts over max n pnts - tpc
TH2F * m_momF
momentum, tpc+svt
TH1F * m_geant_reco_pvtx_x
number of hits ftpcW
TH1F * m_glb_yf0TS
y-coord. of first hit - at start of helix
TH1F * m_glb_sptsTS
number of track points used for fitting - tpc,tpc+svt
TH2F * m_pxl_nhit_Pxl1_tpc_mult
of hits vs sector vs sensor, outer PIXEL layer
TH1F * m_prim_phi0T
radius at start (cm), tpc+svt
TH1F * m_prim_xfT
curvature (1/cm), tpc+svt
TH2F * m_chisq1F
chi square [0], ftpc west
TH2F * m_primtrk_xf_yfFE
Y vs X of first hit on trk, tpc+svt.
TH1F * m_v_r
chisq per dof of vertex fit