mc_exam ntuple description
LIST OF VARIABLES IMPLEMENTED IN THE NTUPLE: ============================================ updated 16 Mar 1999 Ch. Royon updated 23 Mar 1999 A. Zylberstejn updated 14 Jun 1999 A. Zylberstejn updated 17 Jan 2000 A. Zylbertejn updated 10 Feb 2000 O. Kuznetsov: add mass_gen in block VERGEN, update block VERTEX updated 16 Mar 2000 A. Zylbertejn: add pz and d(1/pt ) on block TRACK (for V-00-12-00 and above)
***************** * Ntuple ID = 2 * ***************** ************************** * Block * Data type ************************* * CAL Calorimeter digi * CFTC CFT clusters * CFTD CFT Digi * CFTH CFT hits * CPSC CPS clusters * CPSD CPS digi * CPSH CPS hits * ETMIS Missing Et * FPSC FPS clusters * FPSD FPS digi * JETS_C Jets cone on CAL towers * JETS_KI Jets kT on MC particles with preclustering * JETS_KT Jets kT on CAL towers with preclustering * JETS_NN CAL clusters Nearest Neighbour algorithm * JETS_SC Jets simple cone on CAL towers * KINE Informations on generated event (track vertex...) * MUOND Muon digi * MUONH Muon hits * SMTB SMT reco clusters barrel * SMTC SMT reco clusters disks * SMTD SMT digi * SMTH SMT hits * TRACK Information on reconstructed tracks * VERGEN Information on vertex (?) ****************************************************************** * Blocks = 24 Variables = 186 Max. Columns = 18497 * ************************************************************************ Block Variable Description ************************************************************************ CAL * ncal nb. of informations in calo digi CAL * e_cal(ncal) energy in cell CAL * ieta_cal(ncal) detector eta CAL * iphi_cal(ncal) detector phi CAL * lyr_cal(ncal) layer CAL * pt_cal(ncal) cell pT CFTC * ncftc nb. of informations for CFT/SFT clusters CFTC * lay_cftc(ncftc) layer number (0 to 15 ) CFTC * nst_cftc(ncftc) width of the cluster (in nb. of strips) CFTC * str_cftc(ncftc) center of the cluster (in nb. of strips) CFTD * ncftd nb. of informations in CFT/SFT digi CFTD * fib_cftd(ncftd) fiber nb. CFTD * lay_cftd(ncftd) layer number CFTH * ncfth nb. of informations in CFT/SFT layers CFTH * ed_cfth(ncfth) deposited energy (de/dx) CFTH * g3_cfth(ncfth) Geant code for particle having created this hit CFTH * itr_sfth(ncfth) track number " " CFTH * st_sfth(ncfth) Range in detector CFTH * xg_sfth(ncfth) X in general system CFTH * yg_sfth(ncfth) Y " " CFTH * zg_sfth(ncfth) Z " " CPSC * ncpsc nb. of informations for CPS clusters CPSC * e_cpsc(ncpsc) energy CPSC * x_cpsc(ncpsc) X in general system CPSC * y_cpsc(ncpsc) Y CPSC * z_cpsc(ncpsc) Z CPSD * ncpsd nb. of informations for CPS digi CPSD * e_cpsd(ncpsd) energy CPSD * lay_cpsd(ncpsd) layer CPSD * str_cpsd(ncpsd) strip number CPSH * ncpsh nb. of informations for CPS hits CPSH * ed_cpsh(ncpsh) CPSH * g3_cpsh(ncpsh) CPSH * itr_cpsh(ncpsh) see description for CFTH CPSH * st_cpsh(ncpsh) CPSH * xg_cpsh(ncpsh) CPSH * yg_cpsh(ncpsh) CPSH * zg_cpsh(ncpsh) ETMIS * netmis nb. of informations for reconstructed Etmiss (=1) ETMIS * et_etmis(netmis) Et ETMIS * ex_etmis(netmis) Ex ETMIS * ey_etmis(netmis) Ey FPSC * nfpsc nb. of informations for FPS clusters FPSC * e_fpsc(nfpsc) energy FPSC * x_fpsc(nfpsc) X in general system FPSC * y_fpsc(nfpsc) Y FPSC * z_fpsc(nfpsc) Z FPSD * nfpsd nb. of informations for FPS digi FPSD * e_fpsd(nfpsd) energy FPSD * lay_fpsd(nfpsd) layer FPSD * sla_fpsd(nfpsd) sub_layer FPSD * wed_fpsd(nfpsd) wedge FPSH * nfpsh FPSH * ed_fpsh(nfpsh) FPSH * g3_fpsh(nfpsh) FPSH * itr_fpsh(nfpsh) see description for CFTH FPSH * st_fpsh(nfpsh) FPSH * xg_fpsh(nfpsh) FPSH * yg_fpsh(nfpsh) FPSH * zg_fpsh(nfpsh)JETS_C * njets_c number of cone jets (on CAL towers) JETS_C * E_c(njets_c) energy JETS_C * emE_c(njets_c) EM energy JETS_C * eta_c(njets_c) eta JETS_C * nbTwr_c(njets_c) number of towers in that jet JETS_C * phi_c(njets_c) phi JETS_C * pt_c(njets_c) pT JETS_C * px_c(njets_c) px JETS_C * py_c(njets_c) py JETS_C * pz_c(njets_c) pz JETS_KI * njets_ki number of kT jets (on MC particles, with preclustering) JETS_KI * E_ki(njets_ki) JETS_KI * emE_ki(njets_ki) JETS_KI * eta_ki(njets_ki) JETS_KI * nbTwr_ki(njets_ki) JETS_KI * phi_ki(njets_ki) JETS_KI * pt_ki(njets_ki) see description for cone jets JETS_KI * px_ki(njets_ki) JETS_KI * py_ki(njets_ki) JETS_KI * pz_ki(njets_ki) JETS_KT * njets_kt number of kT jets (on towers, with preclustering) JETS_KT * E_kt(njets_kt) JETS_KT * emE_kt(njets_kt) JETS_KT * eta_kt(njets_kt) JETS_KT * nbTwr_kt(njets_kt) JETS_KT * phi_kt(njets_kt) JETS_KT * pt_kt(njets_kt) see description for cone jets JETS_KT * px_kt(njets_kt) JETS_KT * py_kt(njets_kt) JETS_KT * pz_kt(njets_kt) JETS_NN * njets_nn number of neareast neighbour CAL clusters JETS_NN * E_nn(njets_nn) JETS_NN * emE_nn(njets_nn) JETS_NN * eta_nn(njets_nn) JETS_NN * nbTwr_nn(njets_nn) JETS_NN * phi_nn(njets_nn) JETS_NN * pt_nn(njets_nn) see description for cone jets JETS_NN * px_nn(njets_nn) JETS_NN * py_nn(njets_nn) JETS_NN * pz_nn(njets_nn) JETS_SC * njets_sc number of simple cone jets (on CAL towers) JETS_SC * E_sc(njets_sc) JETS_SC * emE_sc(njets_sc) JETS_SC * eta_sc(njets_sc) JETS_SC * nbTwr_sc(njets_sc) JETS_SC * phi_sc(njets_sc) JETS_SC * pt_sc(njets_sc) see description for cone jets JETS_SC * px_sc(njets_sc) JETS_SC * py_sc(njets_sc) JETS_SC * pz_sc(njets_sc)KINE * mpart nb. of informations for generated particles KINE * cent_kin(mpart) =0 particle not in aceptance =1 part. in 8 layers of SFT not in SMT =2 part. in first 6 layers of SFT not in smt =4 part. in SMT only (forward ) =5 fully contained in SFT and SMT =6 in 6 SFT and SMT KINE * char_kin(mpart) = charge of the particle KINE * e_kin(mpart) energy KINE * eta_kin(mpart) eta KINE * id_kin(mpart) particle id KINE * idx_kin(mpart) part. index (within its chunk) for ref. from other blocks KINE * pid_kin(mpart) particle id of the parent KINE * phi_kin(mpart) phi KINE * prim_kin(mpart) =1 if primary vertex KINE * pt_kin(mpart) Pt KINE * stab_kin(mpart) =1 if stable particle KINE * th_kin(mpart) theta KINE * vtx_kin(mpart) vertex number where this particle was createdMUOND * nmuond nb. of informations for muon digi MUOND * eta_mud(nmuond) MUOND * nl_mud(nmuond) MUOND * nr_mud(nmuond) MUOND * nt_mud(nmuond) MUOND * phi_mud(nmuond) MUONH * nmuonh nb. of informations for muon hits MUONH * edep_muh(nmuonh) MUONH * g3_muh(nmuonh) MUONH * ind_muh(nmuonh) MUONH * itr_muh(nmuonh) MUONH * step_muh(nmuonh) MUONH * xg_muh(nmuonh) MUONH * yg_muh(nmuonh) MUONH * zg_muh(nmuonh)SMTB * nsmtb nb. of informations clusters SMT barrel SMTB * bar_smtb(nsmtb) barrel number SMTB * ex1_smtb(nsmtb) de/dx "n" side SMTB * ex2_smtb(nsmtb) de/dx "p" side SMTB * lay_smtb(nsmtb) layer nb. SMTB * phi_smtb(nsmtb) SMTB * st1_smtb(nsmtb) nb. of strips "n" side SMTB * st2_smtb(nsmtb) nb. of strips "p" side SMTB * xg_smtb(nsmtb) SMTB * yg_smtb(nsmtb) SMTB * zg_smtb(nsmtb) SMTC * nsmtc nb. of informations clusters SMT disks SMTC * dwe_smtc(nsmtc) wedge number SMTC * ex1_smtc(nsmtc) de/dx "n" side SMTC * ex2_smtc(nsmtc) de/dx "p" side SMTC * num_smtc(nsmtc) disk number SMTC * xg_smtc(nsmtc) SMTC * yg_smtc(nsmtc) SMTC * zg_smtc(nsmtc) SMTD * nsmtd nb. of informations for SMT digits SMTD * e_smtd(nsmtd) SMTD * iba_smtd(nsmtd) barrel number SMTD * ild_smtd(nsmtd) ladder number SMTD * ily_smtd(nsmtd) layer number SMTD * ist_smtd(nsmtd) strip number SMTD * iv_smtd(nsmtd) View (n or p side) SMTH * nsmth nb. of informations for SMT hits SMTH * ed_smth(nsmth) SMTH * g3_smth(nsmth) SMTH * ind_smth(nsmth) SMTH * itr_smth(nsmth) SMTH * st_smth(nsmth) see description for CFTH SMTH * xg_smth(nsmth) SMTH * yg_smth(nsmth) SMTH * zg_smth(nsmth)TRACK * ntra nb. of informations for reconstructed tracks TRACK * cft_tra(ntra) nb. of CFT clusters on track TRACK * cftb_tra(ntra) packed word giving hit layers in CFT TRACK * chi2_tra(ntra) Chi2 TRACK * dq_tra(ntra d(1/Pt) at DCA TRACK * eta_tra(ntra) eta at Distance of Closest Approach (DCA) TRACK * id_tra(ntra) track identifier (see : global tracking) = 101 CFT only (central and overlap) 111 central CFT only -no SMT 121 overlap CFT only- no SMT 201 central and overlap CFT tracks extended into the SMT. 301 SMT only 401 tracks in all eta SFT+SMT TRACK * phi8(ntra) phi at CFT layer X8 TRACK * phi_tra(ntra) phi ar DCA TRACK * pt_tra(ntra) Pt TRACK * pz_tra(ntra) Pz at DCA TRACK * smt_tra(ntra) nb. of SMT clusters on track TRACK * smtb_tra(ntra) packed word giving hit layers in SMT TRACK * theta8(ntra) theta at CFT layer X8 TRACK * x8(ntra) X at CFT layer X8 TRACK * x_tra(ntra) X at DCA TRACK * y8(ntra) Y at CFT layer X8 TRACK * y_tra(ntra) Y at DCA TRACK * z8(ntra) Z at CFT layer X8 TRACK * z_tra(ntra) Z at DCA the folllowing elements are written on the Ntuple if the option error_last="diag" has been selected in MC_exam.rcp Diagonal elements for the error matrix at the last plane TRACK * wfi0_l(ntra) TRACK * wimp_l(ntra) TRACK * wpti_l(ntra) TRACK * wtgl_l(ntra) TRACK * wz0_l(ntra) the following elements are written on the Ntuple if the option error_dca="diag" has been selected in MC_exam.rcp Diagonal elements for the error matrix at the Distance of Closest Approach (DCA) TRACK * wfi0_tra(ntra) TRACK * wimp_tra(ntra) TRACK * wpti_tra(ntra TRACK * wtgl_tra(ntra) TRACK * wz0_tra(ntra) The matrix element are for the following track parameters : r is fixed (except at the DCA) 0 - phi 1 - z (cm) 2 - alpha = phi_dir - phi; tan(alpha) = r*dphi/dr 3 - sin(lambda) = dz/ds; tan(lambda) = dz/dsT 4 - q/pT (1/GeV/c) (pT is component of p parallel to cylinder) in case option string error_dca ="none"/ is selected TRACK * dq_tra(ntra) // delta(charge/Pt) at DCAVERGEN * mverg nb. of informations for vertices at generation VERGEN * ifla_gen(mverg) Vertex number (from 1 to mverg-1) VERGEN * ityp_gen(mverg) =1 if primary vertex VERGEN * ntrk_gen(mverg) nb. of tracks above PT min (=100 MeV/c) at this vertex VERGEN * mass_gen(mverg) invariant mass of the "vertex" VERGEN * xver_gen(mverg) =0 if primary VERGEN * yver_gen(mverg) =0 if primary VERGEN * zver_gen(mverg)VERTEX * nvert Total nb. of the reconstructed (and already selected) vertices VERTEX * ityp_ver(nvert) =1 if primary =0 if secondary vertex VERTEX * chi2_ver(nvert) chi^2 of the vertex VERTEX * ntrk_ver(nvert) Total nb. of tracks at the vertex; "Ndof" for the vertex fit equals Ndof=2*ntrk_ver ==> for primary vtx Ndof=2*ntrk_ver-3 ===>for secondary vtx VERTEX * gtr_ver(nvert) packed word giving the numbers of tracks attached to the vertex VERTEX * x_ver(nvert) VERTEX * Errx_ver(nvert) VERTEX * y_ver(nvert) VERTEX * Erry_ver(nvert) VERTEX * z_ver(nvert) VERTEX * Errz_ver(nvert) ----------------------------------------------------------------------------