Data ntuple ID: 202; MC ntuple ID: 203
Blocks in data ntuples:
EVENT block
TRACK block
PIDTRK block
NDSTRK block
NVTX2 block
GAMMA block
VTXPID block
BEAMLIST block
NPART block
PILEUP block
GAMPID block
SPECASSC block
MC ntuples have the blocks above, plus:
GBEAM block
GPART block
EVENT block: deals with information at the
event-level (other blocks at track-level), and beam info
| RUN | run number |
| ISP | spill number |
| IEV | event number |
| EV_STATUS(3) | event status values (not used) |
| LOBITS | trigger bits - before prescaling |
| HIBITS | trigger bits - after prescaling |
| SPTIME | time in the spill (ms) |
| S1TIME | beam counter time start (in 50ps bins) |
| SF_NPL1 | scint fiber # hits on plane 1 |
| SF_NPL2 | scint fiber # hits on plane 2 |
| SFMASK1(3) | scint fiber plane 1 hit mask |
| SFMASK2(3) | scint fiber plane 2 hit mask |
| BESTATUS | |
| BCKVADC(3) | beam ckov adc's 1-3 |
| BCKVTDC(3) | beam ckov tdc's |
| A14STAT | A14 reconstruction status |
| A56STAT | A56 reconstruction status |
| PTOTBEAM | total beam momentum |
| XBEAM(3) | |
| DXDZBEAM | |
| DYDZBEAM | |
| V0NTRKS | |
| V0CHISQ | chi2 for primary vertex fitting |
| V0STATUS | |
| V0X | x-coord of primary vertex |
| V0Y | y-coord of primary vertex |
| V0Z | z-coord of primary vertex (targte centered at 88cm) |
| V0DX(3) | |
| GPEVTYPE | event type (0=inelastic, 1=beam-like, 2=QE, 3=junk) |
| GPEVSTAT | |
| NSECOND | # of secondaries |
| NDELTAE | # of delta electrons (removed from secondaries) |
| NPOS | number of positive tracks |
| NNEG | negative tracks |
| NNEUT | neutral tracks |
| SUMPX | sum of px components of all secondaries |
| SUMPY | py |
| SUMPZ | pz |
| SUMPT | sum of transverse momentum components of all secondaries |
| SUMPTOT | total momentum sum |
| DSUMPX | |
| DSUMPY | |
| DSUMPZ | |
| DSUMPT | |
| DSUMPTOT | |
| XMAX | |
| XMAXIDX | |
| NDCPLHIT(7,5) |
TRACK block: some TPC information
| NTRACK | number of tracks in TRACK block |
| SPECID(NTRACK) | track ID |
| STATUS(NTRACK) | bit mask for track status |
| NHITTPC(NTRACK) | number of raw TPC hits |
| NUNSHTPC(NTRACK) | |
| XBSPECX(NTRACK) | |
| XBSPECY(NTRACK) | |
| XBSPECZ(NTRACK) | |
| XESPECX(NTRACK) | |
| XESPECY(NTRACK) | |
| XESPECZ(NTRACK) | |
| S1ST(NTRACK) | |
| SBFREE(NTRACK) | |
| PSPECX(NTRACK) | px from spectrometer |
| PSPECY(NTRACK) | py |
| PSPECZ(NTRACK) | pz |
| DPOVERPX(NTRACK) | |
| DPOVERPY(NTRACK) | |
| DPOVERPZ(NTRACK) | |
| X0SPEC(NTRACK) | x-comp of primary vertex from spectrometer |
| Y0SPEC(NTRACK) | y-comp |
| Z0SPEC(NTRACK) | z-comp |
| HELICITY(NTRACK) | electric charge (helicity=+1 is negative track...) |
| LASTCHI2(NTRACK) | |
| TPCCHI2X(NTRACK) | |
| TPCCHI2Y(NTRACK) | |
| LASTMETH(NTRACK) | |
| DEDX(NTRACK) | dE/dx (arbitrary units) |
| NUSED(NTRACK) | # of TPC hits used for dE/dx calculation |
| TRDEIDEL(NTRACK) | dE/dx deviation (# of sigmas), assuming pcle is an electron |
| TRDEIDPI(NTRACK) | assuming pcle is pion |
| TRDEIDK(NTRACK) | assuming pcle is kaon |
| TRDEIDP(NTRACK) | assuming pcle is proton |
| TRDEIDD(NTRACK) | assuming pcle is deuteron |
| TRDEIDT(NTRACK) | assuming pcle is triton |
| TRDEIDH3(NTRACK) | assuming pcle is He3 nucleus |
| TRDEIDH4(NTRACK) | assuming pcle is He4 nucleus |
PIDTRK block: some information useful for
particle ID from all subdetectors. Other blocks have more specific info
| NPID | # of pcles in PIDTRK block |
| PIDSPID(NPID) | track ID |
| PIDTRST(NPID) | bit mask for track status |
| PIDCHG(NPID) | charge |
| PTOT(NPID) | total momentum |
| DPTOT(NPID) | |
| PIDNDCHT(NPID) | |
| PIDDEDX(NPID) | dE/dx |
| PIDNUSED(NPID) | # of hits in the TPC |
| PIDSLAT(NPID) | TOF slat number (between 1 and 32) |
| INVBETA(NPID) | 1/beta from TOF (beta=v/c) |
| CKMIRR0(NPID) | # of Ckov mirrors fired |
| CKNPE(NPID) | # of Ckov PEs |
| DEIDE(NPID) | dE/dx deviation (# of sigmas), assuming pcle is an electron |
| DEIDPI(NPID) | assuming pcle is a pion |
| DEIDK(NPID) | assuming pcle is a kaon |
| DEIDP(NPID) | assuming pcle is a proton |
| DEIDD(NPID) | assuming pcle is a deuteron |
| DEIDT(NPID) | assuming pcle is a triton |
| DEIDHE3(NPID) | assuming pcle is a He3 nucleus |
| DEIDHE4(NPID) | assuming pcle is a He4 nucleus |
| TFIDE(NPID) | TOF deviation (# of sigmas), assuming pcle is an electron |
| TFIDPI(NPID) | assuming pcle is a pion |
| TFIDK(NPID) | assuming pcle is a kaon |
| TFIDP(NPID) | assuming pcle is a proton |
| TFIDD(NPID) | assuming pcle is a deuteron |
| TFIDT(NPID) | assuming pcle is a triton |
| TFIDHE3(NPID) | assuming pcle is a He3 nucleus |
| TFIDHE4(NPID) | assuming pcle is a He4 nucleus |
| CKIDE(NPID) | Ckov # of PEs (# of sigmas), assuming pcle is an electron |
| CKIDPI(NPID) | assuming pcle is a pion |
| CKIDKA(NPID) | assuming pcle is a kaon |
| CKIDP(NPID) | assuming pcle is a proton |
| GPIDE(NPID) | |
| GPIDPI(NPID) | |
| GPIDK(NPID) | |
| GPIDP(NPID) | |
| GPIDD(NPID) | |
| GPIDT(NPID) | |
| PROBE(NPID) | global probability that pcle is an electron (from all subdets) |
| PROBPI(NPID) | prob to be a pion |
| PROBK(NPID) | prob to be a kaon |
| PROBP(NPID) | prob to be a proton |
| PROBD(NPID) | prob to be a deuteron |
| PROBT(NPID) | prob to be a triton |
| PROBHE3(NPID) | prob to be a He3 nucleus |
| PROBHE4(NPID) | prob to be a He4 nucleus |
| E910PID(NPID) | E910 pcle numbering scheme (-1=ambiguous, 0=e, 1=mu, 2=pi,
3=k, 4=p, 5=d,...) |
| BESTPID(NPID) | |
| GPIDSTAT(NPID) | bit mask for particle ID (particle ID mask, detector mask) |
NDSTRK block: block with info on downstream
tracks (drift chambers, Cherenkov, TOF)
| NDSTRK | # of pcles in NDSTRK block |
| DSSPECID(NDSTRK) | track ID |
| DSTRSTAT(NDSTRK) | bit mask for spectrack status |
| DSHELIC(NDSTRK) | electric charge (helicity=+1 is negative) |
| DSPTOT(NDSTRK) | total momentum |
| DSDEDXEL(NDSTRK) | dE/dx deviation (# of sigmas) assuming pcle is an electron |
| DSDEDXPI(NDSTRK) | assuming pcle is a pion |
| DSDEDXK(NDSTRK) | assuming pcle is a kaon |
| DSDEDXPR(NDSTRK) | assuming pcle is a proton |
| DSTFSTAT(NDSTRK) | TOF status (-1=initial value, 0=rejected track, 1=OK, 2double track hit) |
| DSTFSLAT(NDSTRK) | TOF slat # |
| PATHL(NDSTRK) | track pathlength to TOF (about 8m from target) |
| TOFCORR(NDSTRK) | corrected TOF |
| DE(NDSTRK) | measures TOF pulse height, for slewing corrections |
| DXTOF(NDSTRK) | x-position in TOF slat (useful for TOF association cuts) |
| DYTOF(NDSTRK) | y-position in TOF slat |
| TOFPID(9,NDSTRK) | TOF deviation (# of sigmas) for 9 different pcle assumptions |
| DSCKSTAT(NDSTRK) | status bits for Ckov detector |
| NMIRROR(NDSTRK) | # of Ckov mirrors hit |
| DSMIRR0(NDSTRK) | |
| DSMIRR1(NDSTRK) | |
| DSMIRR2(NDSTRK) | |
| DSMIRR3(NDSTRK) | |
| DSMNPE0(NDSTRK) | |
| DSMNPE1(NDSTRK) | # of PEs from mirror 0 |
| DSMNPE2(NDSTRK) | from mirror 1 |
| DSMNPE3(NDSTRK) | from mirror 2 |
| DCASQ0(NDSTRK) | from mirror 3 |
| TOTNPE(NDSTRK) | total # of Ckov PEs |
| CKPID(9,NDSTRK) | Ckov deviation (# of sigmas) for 9 different pcle assumptions |
| DSV0CHI2(NDSTRK) | chi2 for primary vertex fit |
| DSUSCHI2(NDSTRK) | upstream chi2 |
| DSDSCHI2(NDSTRK) | downstream chi2 |
| NDCUSHIT(NDSTRK) | # of planes hit in upstream drift chambers (3 US chambers) |
| NDCDSHIT(NDSTRK) | # of planes hit in downstream drift chambers (2 DS chambers) |
NVTX2 block: block to reconstruct 2-body decays
of neutral partciles (K_shorts, lambdas, phis). TPC information in here.
| NVTX2 | # of tracks in NVTX2 block |
| V2STATUS(NVTX2) | bit mask for track status |
| TRACKID1(NVTX2) | track ID for one daughter pcle |
| TRACKID2(NVTX2) | track ID for other daughter pcle |
| HELICITY1(NVTX2) | charge of daughter 1 |
| HELICITY2(NVTX2) | charge of daughter 2 |
| V2P1SPEC(3,NVTX2) | |
| V2P2SPEC(3,NVTX2) | |
| V2S1(NVTX2) | |
| V2S1(NVTX2) | |
| V2XSPEC(NVTX2) | |
| V2YSPEC(NVTX2) | |
| V2ZSPEC(NVTX2) | |
| V2PXSPEC(NVTX2) | |
| V2PYSPEC(NVTX2) | |
| V2PZSPEC(NVTX2) | |
| VTDIRX(NVTX2) | |
| VTDIRY(NVTX2) | |
| VTDIRZ(NVTX2) | |
| DIST2TGT(NVTX2) | distance between pcle decay and target (cm) |
| V2CHISQ(NVTX2) | |
| PPER(NVTX2) | pt of charged daughter tracks wrt parent direction |
| PPAR1(NVTX2) | |
| PPAR2(NVTX2) | |
| ALFA(NVTX2) | measures aymmetry of decay |
GAMMA block: deals with gammas...
| NGAMMA | |
| POSSPID(NGAMMA) | |
| NEGSPID(NGAMMA) | |
| POSTRIDX(NGAMMA) | |
| NEGTRIDX(NGAMMA) | |
| POSORDER(NGAMMA) | |
| NEGORDER(NGAMMA) | |
| GAMMAVX(NGAMMA) | |
| GAMMAVY(NGAMMA) | |
| GAMMAVZ(NGAMMA) | |
| GAMMAPX(NGAMMA) | |
| GAMMAPY(NGAMMA) | |
| GAMMAPZ(NGAMMA) | |
| GAMMAP(NGAMMA) | |
| G1ELECID(NGAMMA) | |
| G2ELECID(NGAMMA) | |
| ANGLE(NGAMMA) | |
| QT(NGAMMA) | |
| DCA(NGAMMA) | |
| GAMMAST(NGAMMA) |
VTXPID block: similar to VTX2, but with more
information
| NVTXPID | number of pcles in VTXPID block (most of the times, 0) |
| VIDTYPE(NVTXPID) | pcle type (0 or 2=lambda, 1=k0) |
| VGLBVID(NVTXPID) | |
| VTX2ID(NVTXPID) | |
| VIDSTAT(NVTXPID) | |
| VIDTRID1(NVTXPID) | |
| VIDTRID2(NVTXPID) | |
| VIDNHIT1(3,NVTXPID) | |
| VIDNHIT2(3,NVTXPID) | |
| VIDCHG1(NVTXPID) | |
| VIDCHG2(NVTXPID) | |
| INVMASS(NVTXPID) | invariant mass |
| PREINVM(NVTXPID) | |
| VIDALFA(NVTXPID) | measures decay asymmetry |
| VIDPPER(NVTXPID) | pt of charged dughter tracks wrt parent's direction |
| VIDD2TGT(NVTXPID) | |
| ALIGNANG(NVTXPID) | |
| OPENANG(NVTXPID) | |
| VIDVX(NVTXPID) | |
| VIDVY(NVTXPID) | |
| VIDVZ(NVTXPID) | |
| VIDPX(NVTXPID) | |
| VIDPY(NVTXPID) | |
| VIDPZ(NVTXPID) | |
| VIDP(NVTXPID) | |
| VIDP1(3,NVTXPID) | |
| VIDP2(3,NVTXPID) | |
| KCM(NVTXPID) | |
| COSTHCM(NVTXPID) | |
| DIRCOSCM(3,NVTXPID) | |
| VIDPID1(3,NVTXPID) | |
| VIDPID2(3,NVTXPID) | |
| VIDXB1(3,NVTXPID) | |
| VIDXE1(3,NVTXPID) | |
| VIDXB2(3,NVTXPID) | |
| VIDXE2(3,NVTXPID) | |
| STPATH1(NVTXPID) | |
| STPATH2(NVTXPID) | |
| VIDCHI2(NVTXPID) | |
| TRK1CHI2(NVTXPID) | |
| TRK2CHI2(NVTXPID) | |
| MASSPENA(NVTXPID) | |
| ANGPENA(NVTXPID) | |
| VIDNDF(NVTXPID) | |
| TRK1NDF(NVTXPID) | |
| TRKNDF(NVTXPID) | |
| VIDCONF(NVTXPID) |
| NBEAMLIST | |
| BMLSTAT(NMBEAMLIST) | |
| TDCY(NMBEAMLIST) | |
| MATCHST(NMBEAMLIST) | |
| MATCHID(NMBEAMLIST) |
NPART block: block with very processed info on
momenta, positions, and best guess for particle ID. Final summary of
event.
| NPART | # of pcles in NPART block; best estimate for pcle multiplicity |
| PARTTYPE(NPART) | particle type (0=good, other=beam-like pcles, or short tracks) |
| PARTSRC(NPART) | track ID |
| PSRCIDX(NPART) | different track ID |
| PARTCHG(NPART) | charge |
| PDGID(NPART) | GEANT pcle ID |
| PE910ID(NPART) | E910 pcle ID |
| PARTPROB(NPART) | probability that track is correctly identified |
| PARTPX(NPART) | px |
| PARTPY(NPART) | py |
| PARTPZ(NPART) | pz |
| PARTX(NPART) | x-coordinate |
| PARTY(NPART) | y-coordinate |
| PARTZ(NPART) | z-coordinate |
| SFIBX(NPART) | scintillating fiber x (trigger detector) |
| SFIBY(NPART) | scintillating fiber y |
| SFIBDXDZ(NPART) |
| NPILEUP | |
| PUX0(NPILEUP) | |
| PUY0(NPILEUP) | |
| PUPTOT(NPILEUP) | |
| PUDXDZ(NPILEUP) | |
| PUDYDZ(NPILEUP) | |
| PUDEIDP(NPILEUP) | |
| PUDEIDPI(NPILEUP) | |
| BMLSTIDX(NPILEUP) |
| NGAMPID | |
| GGLBVID(NGAMPID) | |
| GAMMAIDX(NGAMPID) | |
| GIDSTAT(NGAMPID) | |
| GIDTRIDP(NGAMPID) | |
| GIDTRIDN(NGAMPID) | |
| GIDNHITP(3,NGAMPID) | |
| GIDNHITN(3,NGAMPID) | |
| GIDANG(NGAMPID) | |
| GIDPRANG(NGAMPID) | |
| GIDQT(NGAMPID) | |
| GIDPRQT(NGAMPID) | |
| GIDALFA(NGAMPID) | |
| GIDDIST(NGAMPID) | |
| GIDVX(NGAMPID) | |
| GIDVY(NGAMPID) | |
| GIDVZ(NGAMPID) | |
| GIDPX(NGAMPID) | |
| GIDPY(NGAMPID) | |
| GIDPZ(NGAMPID) | |
| GIDPTOT(NGAMPID) | |
| GIDPP(3,NGAMPID) | |
| GIDPN(3,NGAMPID) | |
| GIDXB1(3,NGAMPID) | |
| GIDXE1(3,NGAMPID) | |
| GIDXB2(3,NGAMPID) | |
| GIDXE2(3,NGAMPID) | |
| GIDSTPP(NGAMPID) | |
| GIDSTPN(NGAMPID) | |
| GIDCHI2(NGAMPID) | |
| GIDCHI2P(NGAMPID) | |
| GIDCHI2N(NGAMPID) | |
| QTPENA(NGAMPID) | |
| GIDNDF(NGAMPID) | |
| GIDNDFP(NGAMPID) | |
| GIDNDFN(NGAMPID) | |
| GIDCONF(NGAMPID) |
SPECASSC block: IDs to match tracks in
different blocks (spectrometer association)
| NSPECTRK | |
| TRACKIDX(NSPECTRK) | |
| DSTRKIDX(NSPECTRK) | |
| PIDTRIDX(NSPECTRK) | |
| PARTIDX(NSPECTRK) |
| GBEVMODE | |
| GBNCOL | |
| GNCOLEL | |
| GBIMP | |
| GBMOM |
GPART block (MC only): this has the generated
("true") pcle info
| NGPART | # of pcles in GPART block |
| GPDGID(NGPART) | PDG ID |
| GEANTID(NGPART) | GEANT ID |
| GPNUMBER(NGPART) | |
| PARENTID(NGPART) | parent track ID |
| GTRACKID(NGPART) | track ID |
| GPMASS(NGPART) | mass |
| GPCHARGE(NGPART) | charge |
| GPX(NGPART) | x-coord. |
| GPY(NGPART) | y-coord. |
| GPZ(NGPART) | z-coord. |
| GPPX(NGPART) | generatedpx |
| GPPY(NGPART) | generated py |
| GPPZ(NGPART) | generated pz |
| GPX_SPEC(NGPART) | |
| GPY_SPEC(NGPART) | |
| GPZ_SPEC(NGPART) | |
| GPPX_SPEC(NGPART) | |
| GPPY_SPEC(NGPART) | |
| GPPZ_SPEC(NGPART) |
GPIDSTAT bit mask
definition
0-8 bits for particle type mask
9-11 bits for detector mask ( dE/dx, TOF, Ckov in
order)
Examples:
IAND(NINT(GPIDSTAT(IPID)/2.**9), 4) .NE. 0 , where 1<IPID<NPID, means that Ckov PEs
callculation was successful
IAND (NINT(GPIDSTAT(IPID)/2.**9),7 ) .EQ. 0 , where 1<IPID<NPID, means that at least one
among TPC, TOF, or Ckov
subdetectors did not have useful
information for this track
PIDTRST, DSTRSTAT bit
mask definitions
The meaning of each bit is defined as follows:
| TK_SUCCESSED | 0x00000000 | good track |
| TK_FAILED | 0x80000000 | failed in track fitting |
| TK_BM_PILEUP | 0x00000001 | beam pileup track |
| TK_FROM_V0 | 0x00000002 | track starts from prim vertex |
| TK_DECAY | 0x00000004 | track starts from secondary decay |
| TK_HAS_VTX | 0x00000008 | track has been used for vertex reconstruction |
| TK_DEDX | 0x00000010 | dE/dX calculation successful |
| TK_LINE | 0x00000020 | staright line track |
| TK_NOFIELD | 0x00000040 | zero field |
| TK_FITCORR | 0x00000080 | done circle fitting correction |
| TK_DEVMEAN | 0x00000100 | has large xdev or ydev mean |
| TK_MERGED | 0x00000200 | merged track |
| TK_ORG_INTPC | 0x00000400 | track originated from inside the TPC |
| TK_TER_INTPC | 0x00000800 | track terminated inside the TPC |
| TK_NEAR_V0 | 0x00001000 | track is close to primary vertex, but wasn't used in the vertex fit |
| TK_BADV0FIT | 0x00002000 | vertex fit has large chisquare |
| TK_DCHIT_MIN | 0x00004000 | has min number of dchits from elnino |
| TK_RESERVE | 0x00008000 | reserved bit for obsolete track flag |
| TK_SPLIT | 0x00010000 | track has been produced by splitting |
| TK_DISTCORR | 0x00020000 | distortion correction recalculated |
| TK_BEAMGAS | 0x00040000 | from beam gas interaction |
| TK_TRACKFIT | 0x00080000 | track has been successfully fitted with variable curvature fitting method |
| TK_DELTA_E | 0x00100000 | track is delta electron |
| TK_QUASI_EL | 0x00200000 | track is from quasi elastic collision |
| TK_BEAM | 0x00400000 | track is uninteracted beam pcle |
| TK_GHOST | 0x00800000 | determined to be duplicate of another track |
Examples:
IAND( DSTRSTAT(IDSTRK), 2 ) .NE. 0 , where 1<IDSTRK<NDSTRK, means that the track is
coming from primary vertex
IAND( ISHFT(PIDTRST(IPID),-1), 1) .EQ. 0 ,
where 1<IPID<NPID, means that the track is not coming from
primary vertex
IAND( ISHFT(PIDTRST(IPID),-4), 1 ) .NE. 0
, where 1<IPID<NPID, means that dE/dx calculation was successful
The last 8 bits are reserved for error states (when TK_FAILED is
true):
| TK_ERR_MASK | 0x7F000000 | error mask |
| TK_ERR_DISASTER | 0x01000000 | reconstruction failed completely |
| TK_ERR_CHISQ | 0x02000000 | failed chisq cut |
| TK_ERR_RIGID | 0x03000000 | failed rigidity cut |
| TK_ERR_NHITS | 0x04000000 | failed nhits cut |
| TK_ERR_NSHARED | 0x05000000 | failed minimum unshared hits cut |
| TK_ERR_HITSDEV | 0x06000000 | bad xdev or ydev distribution |
| TK_ERR_MOMERR | 0x07000000 | momentum error too large |
| TK_ERR_HIGHMOM | 0x08000000 | momentum magnitude too large |
DSCKSTAT bit mask
definition
Status bits for Ckov detector. The meaning of each bit is defined as
follows:
| CKPID_BAD | 0x10000000 | bad pid |
| CKPID_OK | 0x00000001 | good |
| CKPID_NODATA | 0x00000002 | no data in ckovhits_cal |
| CKPID_NAN | 0x00000004 | error in npe calculation |
| CKPID_CENT_BAD | 0x00000008 | bad center mirror |
| CKPID_FIRED | 0x00000010 | ckov fired |
| CKPID_MIR_SHARED | 0x00000020 | light is shared by mirrors |
| CKPID_TRK_SHARED | 0x00000040 | track shared by mirrors |
| CKPID_SMALL | 0x00000080 | small region |