Final AG Ntuple Description Andrew Godley, 21st June 2007 NOTES: Charged track variables. NCand = NCht = Number of charge tracks, in, near or pointing too (hangers) the primary vertex. Prefered variable to use. NPrim = Same as NCand except no hangers that point to primary vertex. NTkFit = Same as NPrim, but tracks must be "good quality" tnchgd = Charged primary tracks - different, and unknown, definition of primary. Help with Variable Names: S,g = Variable is Simulated/generated R = Reconstructed l,lep = charged track, formerly implied lepton candidate Padova versus Harvard ECAL: One can select between Padova and Harvard ECAL clusterisation routines using the calisto datacard. If you choose Harvard, TimeCluster, Chi2Xlep and Chi2Ylep will not be filled, but padova ecal energy will. If you choose Padova, no Harvard variables will be filled, and standard energy eg (ECALE) will be Padova. Also, Calisto correctly removes out of time ECAL information. ****************************************************************** * Ntuple ID = 500 Entries = 44035 Numu->Nue, Selected events ****************************************************************** * Var numb * Type * Packing * Range * Block * Name * ****************************************************************** * 1 * I*4 * * * HEADER * Run Run Number * 2 * I*4 * * * HEADER * Event Event Number, within a run * 3 * I*4 * * * HEADER * EvType OLD ************************************************************************** MCGENE block is only filled for Monte Carlo events * 1 * R*4 * * * MCGENE * XVS * 2 * R*4 * * * MCGENE * YVS * 3 * R*4 * * * MCGENE * ZVS X, Y and Z position of primary vertex * 4 * R*4 * * * MCGENE * PxLS * 5 * R*4 * * * MCGENE * PyLS * 6 * R*4 * * * MCGENE * PzLS X, Y and Z momentum of tagging lepton * 7 * R*4 * * * MCGENE * PxHPChgS * 8 * R*4 * * * MCGENE * PyHPChgS * 9 * R*4 * * * MCGENE * PzHPChgS Summed X, Y and Z momentum of postive hadrons * 10 * R*4 * * * MCGENE * PxHNChgS * 11 * R*4 * * * MCGENE * PyHNChgS * 12 * R*4 * * * MCGENE * PzHNChgS Summed X, Y and Z momentum of negative hadrons * 13 * R*4 * * * MCGENE * PxHNeuS * 14 * R*4 * * * MCGENE * PyHNeuS * 15 * R*4 * * * MCGENE * PzHNeuS Summed X, Y and Z momentum of neutral hadrons * 16 * R*4 * * * MCGENE * ENuS Visible energy * 17 * R*4 * * * MCGENE * Q2S Q^2 momentum transfer * 18 * R*4 * * * MCGENE * W2S W^2 invariant mass of hadronic system * 19 * R*4 * * * MCGENE * NuS Hadronic Energy * 20 * R*4 * * * MCGENE * XBjS * 21 * R*4 * * * MCGENE * YBjS X and Y Bjorken * 22 * R*4 * * * MCGENE * PxFermi * 23 * R*4 * * * MCGENE * PyFermi * 24 * R*4 * * * MCGENE * PzFermi * 25 * R*4 * * * MCGENE * EFermi Fermi momentum 4 vector * 26 * I*4 * * * MCGENE * ParentID Neutrino Parent * 27 * I*4 * * * MCGENE * LeadS Flag of true leading lepton ************************************************************************** GLOBREC block contains general information about an event * 1 * I*4 * * * GLOBREC * ISele OLD * 2 * I*4 * * * GLOBREC * tnchgd Total number of primary charged particles (from DST) * 3 * I*4 * * * GLOBREC * tnneut Total number of primary neutrals * 4 * I*4 * * * GLOBREC * NUnH Number of unused driftchamber hits * 5 * R*4 * * * GLOBREC * FracUnH Fraction of unused driftchamber hits * 6 * I*4 * * * GLOBREC * NTkFit Number of tracks fitted to the primary vertex * 7 * I*4 * * * GLOBREC * NPrim Number of Charged tracks in/near the primary vertex * 8 * I*4 * * * GLOBREC * NCht Number of Charged tracks * 9 * I*4 * * * GLOBREC * nvxneu Number of vertices not close to primary with no incoming tracks * 10 * I*4 * * * GLOBREC * PzBack Number of backwards going tracks * 11 * R*4 * * * GLOBREC * VtxChg Summed charge at the primary vertex * 12 * R*4 * * * GLOBREC * VtxChgNH Summed charge at the primary vertex excluding hangers * 13 * R*4 * * * GLOBREC * XVR * 14 * R*4 * * * GLOBREC * YVR * 15 * R*4 * * * GLOBREC * ZVR X, Y and Z position on vertex * 16 * I*4 * * * GLOBREC * unusedver Unused hits in the vicinity of the primary vertex * 17 * R*4 * * * GLOBREC * ECalOut Total energy that is flagged as out of time by ECAL * 18 * R*4 * * * GLOBREC * HCalNeu Neutral HCAL energy ************************************************************************** Information for verticies with one or two charged tracks * 19 * I*4 * * [0,25] * GLOBREC * nvzero Number of secondary verticies * 20 * R*4 * * * GLOBREC * v0posx(nvzero) * 21 * R*4 * * * GLOBREC * v0posy(nvzero) * 22 * R*4 * * * GLOBREC * v0posz(nvzero) Position of vertex * 20 * R*4 * * * GLOBREC * v0gposx(nvzero) * 21 * R*4 * * * GLOBREC * v0gposy(nvzero) * 22 * R*4 * * * GLOBREC * v0gposz(nvzero) Generated position of vertex * 23 * I*4 * * * GLOBREC * v0type(nvzero) Geant ID of particle that decayed/interacted at vertex * 24 * I*4 * * * GLOBREC * v0hits(nvzero) Unused hits around vertex * 25 * R*4 * * * GLOBREC * v0chi2(nvzero) Chi2 of vertex * 26 * I*4 * * * GLOBREC * decay(nvzero) 10 one positive track -10 one negative track 0 a positive and a negative track 1 two postive tracks -1 two negative tracks ************************************************************************** Information about positive track from above vertex * 27 * I*4 * * * GLOBREC * pv0typ(nvzero) Geant ID of positive track * 28 * R*4 * * * GLOBREC * pv0px(nvzero) * 29 * R*4 * * * GLOBREC * pv0py(nvzero) * 30 * R*4 * * * GLOBREC * pv0pz(nvzero) Momentum * 28 * R*4 * * * GLOBREC * pv0gpx(nvzero) * 29 * R*4 * * * GLOBREC * pv0gpy(nvzero) * 30 * R*4 * * * GLOBREC * pv0gpz(nvzero) Generated Momentum * 31 * R*4 * * * GLOBREC * pv0epxpx(nvzero) * 32 * R*4 * * * GLOBREC * pv0epypy(nvzero) * 33 * R*4 * * * GLOBREC * pv0epzpz(nvzero) * 34 * R*4 * * * GLOBREC * pv0epxpy(nvzero) * 35 * R*4 * * * GLOBREC * pv0epypz(nvzero) * 36 * R*4 * * * GLOBREC * pv0epxpz(nvzero) Momentum error rror matrix * 37 * R*4 * * * GLOBREC * pv0startx(nvzero) * 38 * R*4 * * * GLOBREC * pv0starty(nvzero) * 39 * R*4 * * * GLOBREC * pv0startz(nvzero) First hit position * 40 * R*4 * * * GLOBREC * pv0endpx(nvzero) * 41 * R*4 * * * GLOBREC * pv0endpy(nvzero) * 42 * R*4 * * * GLOBREC * pv0endpz(nvzero) Momentum at track end * 43 * R*4 * * * GLOBREC * pv0endx(nvzero) * 44 * R*4 * * * GLOBREC * pv0endy(nvzero) * 45 * R*4 * * * GLOBREC * pv0endz(nvzero) Last hit of track * 46 * I*4 * * * GLOBREC * phits(nvzero) Number of hits in track * 47 * R*4 * * * GLOBREC * plength(nvzero) Length of track * 48 * R*4 * * * GLOBREC * ptrdpicon(nvzero) * 49 * R*4 * * * GLOBREC * ptrdelacc(nvzero) TRD info * 50 * R*4 * * * GLOBREC * pecale(nvzero) Associated ECAL energy Harvard * 51 * R*4 * * * GLOBREC * pecalpad(nvzero) Associated ECAL energy Padova * 52 * R*4 * * * GLOBREC * peradx(nvzero) * 53 * R*4 * * * GLOBREC * perady(nvzero) X,Y radius of ECAl hit Harvard * 54 * R*4 * * * GLOBREC * pprsx(nvzero) * 55 * R*4 * * * GLOBREC * pprsy(nvzero) PRS info * 56 * R*4 * * * GLOBREC * pchi2(nvzero) Track Chi2 * 57 * I*4 * * * GLOBREC * pphaseii(nvzero) Is track a Phase II muon * 58 * I*4 * * * GLOBREC * pflag(nvzero) Positron? * 59 * R*4 * * * GLOBREC * pbrem(nvzero) Associated brem energy ************************************************************************** Information for the negative track * 60 * I*4 * * * GLOBREC * nv0typ(nvzero) * 61 * R*4 * * * GLOBREC * nv0px(nvzero) * 62 * R*4 * * * GLOBREC * nv0py(nvzero) * 63 * R*4 * * * GLOBREC * nv0pz(nvzero) * 61 * R*4 * * * GLOBREC * nv0gpx(nvzero) * 62 * R*4 * * * GLOBREC * nv0gpy(nvzero) * 63 * R*4 * * * GLOBREC * nv0gpz(nvzero) * 64 * R*4 * * * GLOBREC * nv0epxpx(nvzero) * 65 * R*4 * * * GLOBREC * nv0epypy(nvzero) * 66 * R*4 * * * GLOBREC * nv0epzpz(nvzero) * 67 * R*4 * * * GLOBREC * nv0epxpy(nvzero) * 68 * R*4 * * * GLOBREC * nv0epypz(nvzero) * 69 * R*4 * * * GLOBREC * nv0epxpz(nvzero) * 70 * R*4 * * * GLOBREC * nv0startx(nvzero) * 71 * R*4 * * * GLOBREC * nv0starty(nvzero) * 72 * R*4 * * * GLOBREC * nv0startz(nvzero) * 73 * R*4 * * * GLOBREC * nv0endpx(nvzero) * 74 * R*4 * * * GLOBREC * nv0endpy(nvzero) * 75 * R*4 * * * GLOBREC * nv0endpz(nvzero) * 76 * R*4 * * * GLOBREC * nv0endx(nvzero) * 77 * R*4 * * * GLOBREC * nv0endy(nvzero) * 78 * R*4 * * * GLOBREC * nv0endz(nvzero) * 79 * I*4 * * * GLOBREC * nhits(nvzero) * 80 * R*4 * * * GLOBREC * nlength(nvzero) * 81 * R*4 * * * GLOBREC * ntrdpicon(nvzero) * 82 * R*4 * * * GLOBREC * ntrdelacc(nvzero) * 83 * R*4 * * * GLOBREC * necale(nvzero) * 84 * R*4 * * * GLOBREC * necalpad(nvzero) * 85 * R*4 * * * GLOBREC * neradx(nvzero) * 86 * R*4 * * * GLOBREC * nerady(nvzero) * 87 * R*4 * * * GLOBREC * nprsx(nvzero) * 88 * R*4 * * * GLOBREC * nprsy(nvzero) * 89 * R*4 * * * GLOBREC * nchi2(nvzero) * 90 * I*4 * * * GLOBREC * nphaseii(nvzero) * 91 * I*4 * * * GLOBREC * nflag(nvzero) * 92 * R*4 * * * GLOBREC * nbrem(nvzero) ************************************************************************** Information for secondary tracks coming from vertices with 3 or greater tracks. * 102 * I*4 * * [0,25] * GLOBREC * nsecond Number of secondaries from 3+ track verticies * 103 * I*4 * * * GLOBREC * secmcid(nsecond) Geant ID of track * 104 * I*4 * * * GLOBREC * secq(nsecond) Charge of track The rest of the variables are as described for the positive track above. * 105 * R*4 * * * GLOBREC * secpx(nsecond) * 106 * R*4 * * * GLOBREC * secpy(nsecond) * 107 * R*4 * * * GLOBREC * secpz(nsecond) * 108 * R*4 * * * GLOBREC * secgpx(nsecond) * 109 * R*4 * * * GLOBREC * secgpy(nsecond) * 110 * R*4 * * * GLOBREC * secgpz(nsecond) * 111 * R*4 * * * GLOBREC * secepxpx(nsecond) * 112 * R*4 * * * GLOBREC * secepypy(nsecond) * 113 * R*4 * * * GLOBREC * secepzpz(nsecond) * 114 * R*4 * * * GLOBREC * secepxpy(nsecond) * 115 * R*4 * * * GLOBREC * secepypz(nsecond) * 116 * R*4 * * * GLOBREC * secepxpz(nsecond) * 117 * R*4 * * * GLOBREC * secstartx(nsecond) * 118 * R*4 * * * GLOBREC * secstarty(nsecond) * 119 * R*4 * * * GLOBREC * secstartz(nsecond) * 120 * R*4 * * * GLOBREC * secendpx(nsecond) * 121 * R*4 * * * GLOBREC * secendpy(nsecond) * 122 * R*4 * * * GLOBREC * secendpz(nsecond) * 123 * R*4 * * * GLOBREC * secendx(nsecond) * 124 * R*4 * * * GLOBREC * secendy(nsecond) * 125 * R*4 * * * GLOBREC * secendz(nsecond) * 126 * I*4 * * * GLOBREC * sechits(nsecond) * 127 * R*4 * * * GLOBREC * seclength(nsecond) * 128 * R*4 * * * GLOBREC * sectrdpicon(nsecond) * 129 * R*4 * * * GLOBREC * sectrdelacc(nsecond) * 130 * R*4 * * * GLOBREC * sececale(nsecond) * 131 * R*4 * * * GLOBREC * sececalpad(nsecond) * 132 * R*4 * * * GLOBREC * secradx(nsecond) * 133 * R*4 * * * GLOBREC * secrady(nsecond) * 134 * R*4 * * * GLOBREC * secprsx(nsecond) * 135 * R*4 * * * GLOBREC * secprsy(nsecond) * 136 * R*4 * * * GLOBREC * secchi2(nsecond) * 137 * I*4 * * * GLOBREC * secphaseii(nsecond) ************************************************************************** Neutral cluster information * 93 * I*4 * * [0,25] * GLOBREC * nclu Number of Clusters * 94 * I*4 * * * GLOBREC * TimeCluster(nclu) Timing - Padova only * 95 * I*4 * * * GLOBREC * clustyp(nclu) Geant ID of cluster * 96 * R*4 * * * GLOBREC * ECluster(nclu) Energy * 97 * R*4 * * * GLOBREC * RadXCluster(nclu) * 98 * R*4 * * * GLOBREC * RadYCluster(nclu) X,Y radius Harvard * 99 * R*4 * * * GLOBREC * XCluster(nclu) * 100 * R*4 * * * GLOBREC * YCluster(nclu) Cluster position * 101 * I*4 * * * GLOBREC * ncells(nclu) Number of ECAL cells in cluster ************************************************************************** Primary track information * 1 * I*4 * * [0,25] * CCEVEREC * ncand Number of tracks * 2 * I*4 * * * CCEVEREC * LepType(ncand) Type of tagging lepton - OLD * 3 * I*4 * * * CCEVEREC * isprimary(ncand) OLD Is reconstructed primary lepton according AG criteria * 4 * I*4 * * * CCEVEREC * istruelep(ncand) MC is this the true leading lepton * 5 * I*4 * * * CCEVEREC * mctrackid(ncand) Geant ID of track * 6 * I*4 * * * CCEVEREC * parenttype(ncand) Geant ID of parent * 7 * I*4 * * * CCEVEREC * LeadLepR(ncand) OLD * 8 * I*4 * * * CCEVEREC * LepQ(ncand) Charge * 9 * I*4 * * * CCEVEREC * TotVeto(ncand) Hits in first plane of DC * 10 * I*4 * * * CCEVEREC * TotTube(ncand) Hits in tube back to first plane * 11 * I*4 * * * CCEVEREC * trackhits(ncand) Hits in track * 12 * R*4 * * * CCEVEREC * mctrackPx(ncand) * 13 * R*4 * * * CCEVEREC * mctrackPy(ncand) * 14 * R*4 * * * CCEVEREC * mctrackPz(ncand) MC momentum of track * 15 * R*4 * * * CCEVEREC * DeltaP(ncand) Error on Momentum * 16 * R*4 * * * CCEVEREC * epxpx(ncand) * 17 * R*4 * * * CCEVEREC * epypy(ncand) * 18 * R*4 * * * CCEVEREC * epzpz(ncand) * 19 * R*4 * * * CCEVEREC * epxpy(ncand) * 20 * R*4 * * * CCEVEREC * epypz(ncand) * 21 * R*4 * * * CCEVEREC * epxpz(ncand) Momentum error matrix * 16 * R*4 * * * CCEVEREC * rqtr(ncand) RQT kinematic variable * 17 * R*4 * * * CCEVEREC * XEndLep(ncand) * 18 * R*4 * * * CCEVEREC * YEndLep(ncand) * 19 * R*4 * * * CCEVEREC * ZEndLep(ncand) X, Y and Z position of last hit * 20 * R*4 * * * CCEVEREC * XStartLep(ncand) * 21 * R*4 * * * CCEVEREC * YStartLep(ncand) * 22 * R*4 * * * CCEVEREC * ZStartLep(ncand) X, Y and Z position of first hit * 23 * R*4 * * * CCEVEREC * PXEndlep(ncand) * 24 * R*4 * * * CCEVEREC * PYEndLep(ncand) * 25 * R*4 * * * CCEVEREC * PZEndLep(ncand) Momentum at last hit * 26 * R*4 * * * CCEVEREC * tracklength(ncand) Length of track * 27 * R*4 * * * CCEVEREC * PxLR(ncand) * 28 * R*4 * * * CCEVEREC * PyLR(ncand) * 29 * R*4 * * * CCEVEREC * PzLR(ncand) Momentum * 30 * R*4 * * * CCEVEREC * PxHNeuR * 31 * R*4 * * * CCEVEREC * PyHNeuR * 32 * R*4 * * * CCEVEREC * PzHNeuR Summed neutral hadron momentum * 33 * R*4 * * * CCEVEREC * Evis(ncand) Visible energy * 34 * R*4 * * * CCEVEREC * Q2R(ncand) * 35 * R*4 * * * CCEVEREC * W2R(ncand) * 36 * R*4 * * * CCEVEREC * NuR(ncand) * 37 * R*4 * * * CCEVEREC * XBjR(ncand) * 38 * R*4 * * * CCEVEREC * YBjR(ncand) Kinematic variables * 39 * R*4 * * * CCEVEREC * PtMis(ncand) Magnitude of missing transverse momentum * 40 * R*4 * * * CCEVEREC * PtLepMis(ncand) P_T of lepton with respect to missing P_T * 41 * R*4 * * * CCEVEREC * NuHat(ncand) (P_T lepton - P_T missing)/P_T Hadron * 42 * R*4 * * * CCEVEREC * HCalE(ncand) Associated HCAL energy * 43 * R*4 * * * CCEVEREC * DeltaZ(ncand) Z distance of start of lepton track to primary vertex * 44 * R*4 * * * CCEVEREC * PiConLep(ncand) * 45 * R*4 * * * CCEVEREC * ElAccLep(ncand) * 46 * R*4 * * * CCEVEREC * praccppi(ncand) * 47 * R*4 * * * CCEVEREC * piconppi(ncand) * 48 * R*4 * * * CCEVEREC * praccpk(ncand) * 49 * R*4 * * * CCEVEREC * kconppi(ncand) TRD info * 50 * R*4 * * * CCEVEREC * EPrsXLep(ncand) * 51 * R*4 * * * CCEVEREC * EPrsYLep(ncand) PRS info * 52 * R*4 * * * CCEVEREC * EBremLep(ncand) Associated Brem strip energy - Harvard * 53 * R*4 * * * CCEVEREC * Chi2XLep(ncand) * 54 * R*4 * * * CCEVEREC * Chi2YLep(ncand) Chi2 of ECAL X and Y shower profile - Padova only * 55 * R*4 * * * CCEVEREC * EOPLep(ncand) E-P / E+P calculated with BREM and Pion first hit momentum * 56 * R*4 * * * CCEVEREC * EOPLepSan(ncand) E-P / E+P calculated with ECAL and end momentu * 57 * R*4 * * * CCEVEREC * ECalE(ncand) Associated ECAL energy Harvard * 58 * R*4 * * * CCEVEREC * ECalPadova(ncand) Associated ECAL energy Padova * 59 * R*4 * * * CCEVEREC * ECalRadX(ncand) * 60 * R*4 * * * CCEVEREC * ECalRadY(ncand) X and Y radius of ECAL energy cluster * 61 * R*4 * * * CCEVEREC * ECalX(ncand) * 62 * R*4 * * * CCEVEREC * ECalY(ncand) X and Y extrapolated to ECAL * 63 * R*4 * * * CCEVEREC * PhiLepHad(ncand) Angle between lepton and hadron shower * 64 * R*4 * * * CCEVEREC * PhiMisHad(ncand) Angle between missing P_T and hadron shower * 65 * R*4 * * * CCEVEREC * PhiMinCh(ncand) Minimum phi (transverse angle) between lepton and any charged track * 66 * R*4 * * * CCEVEREC * ThMinCh(ncand) Minimum theta (longitudinal angle) between lepton and any charged track * 67 * R*4 * * * CCEVEREC * QlepLep(ncand) P_T of lepton with respect to the hadron shower * 68 * R*4 * * * CCEVEREC * QtLep(ncand) P_T of lepton with respect to the event * 69 * R*4 * * * CCEVEREC * probmu(ncand) Probability to be a muon from the muon chambers * 70 * R*4 * * * CCEVEREC * probmucal(ncand) Probability to be a muon from the calorimiters * 71 * R*4 * * * CCEVEREC * mutaglike(ncand) Muon likelihood * 72 * R*4 * * * CCEVEREC * likenccc(ncand) Event CC/NC likelihood * 73 * I*4 * * * CCEVEREC * nx1(ncand) * 74 * I*4 * * * CCEVEREC * ny1(ncand) * 75 * I*4 * * * CCEVEREC * nx2(ncand) * 76 * I*4 * * * CCEVEREC * ny2(ncand) Hits in muon chamber stations 1 and 2, x and y projections * 77 * I*4 * * * CCEVEREC * isstub(ncand) Is a muon track stub * 78 * I*4 * * * CCEVEREC * nveto(ncand) Hits in veto postion of muon chambers * 79 * R*4 * * * CCEVEREC * chi2mu1(ncand) Chi2 of match to muon chamber station 1 * 80 * R*4 * * * CCEVEREC * chi2mu2(ncand) Chi2 of match to muon chamber station 2 * 81 * I*4 * * * CCEVEREC * NUpHanS Number of hangers starting before the primary vertex * 82 * I*4 * * * CCEVEREC * NUpHanE Number of hangers ending before the primary vertex * 83 * I*4 * * * CCEVEREC * NUpHan2 Number of hangers starting and ending before the primary vertex * 84 * I*4 * * * CCEVEREC * NUpBadS * 85 * I*4 * * * CCEVEREC * NUpBadE * 86 * I*4 * * * CCEVEREC * NUpBad2 Number of bad tracks starting and ending before the primary vertex * 87 * I*4 * * * CCEVEREC * MuPV Number of mu+ in the event * 88 * I*4 * * * CCEVEREC * MuPPrimV Number of mu+ from primary * 89 * I*4 * * * CCEVEREC * MuMV Number of mu- in the event * 90 * I*4 * * * CCEVEREC * MuMPrimV Number of mu- from primary ****************************************************************** * Block * Entries * Unpacked * Packed * Packing Factor * ****************************************************************** * HEADER * 44035 * 12 * 12 * 1.000 * * MCGENE * 44035 * 108 * 108 * 1.000 * * GLOBREC * 44035 * 8180 * Var. * Variable * * CCEVEREC * 44035 * 7656 * Var. * Variable * * Total * --- * 15956 * Var. * Variable * ****************************************************************** * Blocks = 4 Variables = 221 Max. Columns = 3989 * ******************************************************************