The OMNTUPLE Ntuple
Version 9

This page describes the OMntuple that Evgeny Popkov uses for his analysis. It has been made specially for the Squark and Gluino search. Adam Lyon, Joey Thompson and Gervasio Gomez are the authors of the  ntuple.
Version 8 has many new fields in the muon section and two new blocks. Many fields in the EM section have been renamed to make more sense.
Version 9 has several new fields in FILT, JET and DTRK sections and new QSIM block.

Fields in the ntuple

The omntuple is a large column wise ntuple with many blocks. Blocks can be enabled or disabled when the ntuple is produced, so an ntuple you are examining may not have all of the blocks listed below.
 
 
EVENT (glob bank info) 
TRIG (Triggers passed) 
FILT (Filters passed) 
AID (Hot cell [AIDA]) 
VERTEX (Vertices) 
LVL0 (Level zero and luminosity) 
EM (EM objects) 
EMT (Tracks associated with PELCs) 
EMV (Vertex corrected EM objects) 
MUON (Muons) 
JET (Jets) 
JNEP (Jets without electron energy) 
UJET (Uncorrected jets) 
L1L2JET (Level 1, Large Tile and Level 2 jets**) 
MET (Missing ET
HCIJ (Hot cells in jets**) 
TAU (Taus) 
DTRK CDC track information 
ISAJET (Isajet info) 
PJET (Pjets) 
OMET (Missing ET from other vertices) 
SCANMET (Missing ET at various vertices along z
MTC (MTC) 
ASTUB (A-PHI prototype counter data) 
QSIM (Quick Simulation Information)
NOTE: Here are some codes you'll see when describing the ntuple fields:

EVENT block:

Basic event information
run             I run number (from RUNNO)
event           I event number (from EVENT)
reco_ver        I reco version (from HEAD/RECO/HSTR)
cal_sumet     ! R total scalar Et in calorimeter (from GLOB)
cal_sume      ! R Total energy in Calorimeter (from GLOB)
ecs_sume        R Total energy in south EC (")
ccup_sume       R Total energy in upper CC
ccdn_sume       R Total energy in lower CC
ecn_sume        R Total energy in north EC
ec_sumet        R Total scalar Et in EC
cc_sumet        R Total scalar Et in CC
mr_sume         R Total energy along Main Ring (from GLOB)
mr_sumet        R Total scalar Et along Main Ring (from GLOB)
hot_sume        R Total energy of HOT cells  (from GLOB)
hot_vsumet      R Total vectorial sum of Et of HOT cells (from GLOB)
ev29            R Time since MR event 29 (seconds) (from GLOB)
l1met           R Level 1 missing Et (from ESUM)
l1metphi        R Phi of level 1 missing Et (from ESUM)
l2met           R Level 2 missing Et (from ESUM)
l2metphi        R Phi of level 2 missing Et (from ESUM)
l2vtxz          R Z position of vertex from level 2 (from ESUM)
recovtxz      ! R Z position of vertex from reco (from ESUM)
mi_tool         I Results from the multiple interaction tool
mi_tool1        I Results from the NEW multiple_interaction_tool_run1
nege_cal        R Negative Energy in the calorimeter (GLOB)
neget_cal       R Negative Et in the calorimeter
mrbs_loss       L True if event was taken during MRBS_LOSS time
micro_blank     L True if event was taken during MICRO_BLANK time
mask_badrun   * I The bad run mask.

TRIG block:

This block contains info about which level 1 triggers were passed. The variables are all logical and are true if that trigger was passed. The trigger is searched for in the TSUM bank by its name.
trmet           L Passed missing_et trigger 
trj3mis         L Passed jet_3_miss trigger
trj1mis         L Passed jet_1_miss trigger
trem1hi         L Passed em 1 high
trem2med        L Passed em 2 med
trem1mon        L Passed em monitor trigger
trjetmin        L Passed jet min trigger
trjt35          L Passed jet 35 trigger (single interaction required!)
trjt60          L Passed jet 60 trigger
trjetmax        L Passed jet max trigger
trljt25         L Passed jet 25 trigger
trljt15         L Passed jet 15 trigger
trljt10nol0     L Passed jet 10 no level zero required trigger
trljt15nol0     L Passed jet 15 no level zero required trigger
trjet2med     * L Passed jet 2 med

FILT block:

This block contians info about which level 2 filters were passed. The variables are all logical and are true if that filter was passed. The filter is searched for in the TSUM bank by its name.
flmet           L Passed missing_et filter
flmetmax        L Passed missing_et_max filter
flj3mis         L Passed jet_3_miss filter
flj2mis         L Passed jet_2_miss filter
flj2mismon  *** L Passed jet_2_miss_mon filter
flj1veto        L Passed jet_1_veto filter
fljmele         L Passed jet_miss_ele
flj1mis         L Passed jet_1_miss filter
fljmis_e        
fle1_eitkcc_ms
fle1_eitkcc_esc
fle1_esc
fle1_e_ms_m
fle1_e_m        L Passed em1_ele_mon filter
fle2_ei_ms
fle2_ei2_hi
fle2_ei_esc     L Passed em2_eis_esc filter
fle2_ei_i
fle2_ei_e
fle2_ei_ei
fle2_eitkcc_ei
fle1m_etk1_mon
fle1m_e1_mon    L Passed ele_1_mon filter
fle1m_e_ms_mon
flminbias       L Passed minbias filter
fljetmin        L Passed jet_min filter
fljet50         L Passed jet_50
fljet85         L Passed jet_85
fljetmax        L Passed jet_max
fljet140        L Passed jet_140
fljet35         L Passed jet_35
fljet60         L Passed jet_60
fljet30         L Passed jet_30
fljet20nol0     L Passed jet_20_nol0
fljet30nol0     L Passed jet_30_nol0
fljet15nol0     L Passed jet_15_nol0
fljet3mon     * L Passed jet 3 mon
fljet4mon     * L Passed jet 4 mon

AID block:

This block contains information returned by the AIDA hot cell killer.

Each array has naid elements, one for each hot cell reported by AIDA. The order of the elements are that returned by AIDA. These values are all from the CAID bank.

naid                    I Number of hot cells flagged
aid_ieta(naid)          R Detector eta of hot cell
aid_iphi(naid)          F Detector phi of hot cell
aid_ilyr(naid)          F Layer of hot cell
aid_scratio(naid)       F Scalar ratio (0-255)
aid_e(naid)             F Energy of hot cell

VERTEX block:

This block contains information about the vertices found by reco. The arrays have nvtx elements. The first element should be the primary and can be checked by vtx_main. These values are all obtained from the VERT bank. Microdsts only contain information about the z position.
nvtx                 I Number of primary vertices (FROM VERH)
vtx_x(nvtx)        ! R X position (ALL FROM VERT)
vtx_y(nvtx)        ! R Y position
vtx_z(nvtx)          R Z position
vtx_dx(nvtx)       ! R Uncert in X
vtx_dy(nvtx)       ! R           Y
vtx_dz(nvtx)         R           Z
vtx_wt(nvtx)       ! R Weight (% of tracks from this vertex)
vtx_ntrks(nvtx)    ! I Number of tracks contributing to this vertex
vtx_main(nvtx)     ! L Main primary vertex? [1]
vtx_bycdc(nvtx)    ! L Vertex found by CDC method 1?
vtx_byfdc(nvtx)    ! L Vertex found by FDC?
vtx_bycdc2(nvtx)   ! L Vertex found by CDC method 2?
vtx_bycdc3(nvtx)   ! L Vertex found by CDC method 3?

Notes:
[1] vtx_main is true if this vertex is flagged by reco as the primary
(should always be the first entry in arrays).

LVL0 Block:

This block contains level zero and instantaneous luminosity information.
l0_fastzvtx     R Vertex position from fast z
l0_slowzvtx     R Vertex position from slow z
l0_slowzqual    R Slow z quality factor
l0_goodfz       L Flag if fast z is good value
l0_miflag(4)    L Logical array with multiple interaction flags [1]
l0_goodsz       L Flag if slow z is good value
l0_phalo        L Flag if event during P halo
l0_pbhalo       L Flag if event during Pbar halo
l0_slowzint     L Flag for slow z interaction
lum_instlum(2)  R Instantaneous Luminosity [2]
lum_age(2)      R Time between last luminosity reading and event time (hrs)
lum_err         I Error code returned by GETLUM_CORR [3]
lum_vaxtime     I VAXtime used to look up luminosity
lum_beginrun  * R Not filled
lum_endrun    * R Not filled
lum_avelum    * R Average luminosity for that run

Notes:
[1] l0_miflag always has four elements 
    miflag(1) = multiple interaction flag 0
    miflag(2) = muitiple interaction flag 1
    ...
    miflag(4) = muitiple interaction flag 3

[2] lum_instlum and lum_age always have two elements.  The first is
the luminosity reading from before the event occured.  The second is
the luminosity reading after the event occured.  The luminosity is
obtained from GETLUM_CORR which contains D0LUMC to D0LUMSD
corrrections.

[3] See d0$physics_util$general:getlum_corr.for
(uff $d0physics_util/general/getlum_corr on unix) for values of lum_err.

EM block:

EM object information (Electrons and photons). EM objects are sorted by pT (highest first).
nemobj                  I Number of em objs in ntuple
nemtot                  I Number of em objs in ntuple   
em_id(nemobj)           I 11=PPHO bank, 12=PELC bank

em_etl1(nemobj)         R\ L1 ET     !from uDST using UDST_VALUE
em_etl2(nemobj)         R \Raw L2 ET !or from GET_L1L2_EMCLUS
em_etl2corr(nemobj)     R /Corr L2 ET (L2 cuts on this)
em_l2mask(nemobj)       I/ Bitwise pass/fail mask in L2 [1]

em_qual(nemobj)         I +1 pass RCP mask/-1 fail RCP mask
em_status(nemobj)       I bitwise status returned from CLEANEM
em_jetoverlap(nemobj)   I ID# (in NTUPLE) of overlapping jet
em_ex(nemobj)           R q(lclus+3)  (err=-999)
em_ey(nemobj)           R q(lclus+4)  (err=-999)
em_ez(nemobj)           R q(lclus+5)  (err=-999)
em_etot(nemobj)         R q(lclus+6)  CAFIX corrected energy
em_et(nemobj)           R q(lclus+7)  (err=-999)
em_theta(nemobj)        R q(lclus+8)  (err=-999)
em_eta(nemobj)          R q(lclus+9)  (err=-999)
em_phi(nemobj)          R q(lclus+10) (err=-999)        
em_chisq(nemobj)        R cquan(4)  H-Matrix chi-square
em_deteta(nemobj)       R from TEMP_FIX_THETA, based on PELC 
                               track if available (err=-999)
em_ieta(nemobj)         I cquan(5)  ieta of hottest cell in em3
em_emfrac(nemobj)       R cquan(9)  
em_iso1(nemobj)         R cquan(13) frac clus energy in cone .4
em_iso2(nemobj)         R cquan(14) frac clus energy in cone .7
em_etiso1(nemobj)       R cquan(13) frac clus et in cone .4
em_etiso2(nemobj)       R cquan(13) frac clus et in cone .7
em_phiso(nemobj)      * R Isolation as per photon
em_match(nemobj)        R tquan(12) PELC trk match signif
em_dcrk(nemobj)         R cquan(20) angular dist to phi crack
em_ncells(nemobj)       R cquan(21) number of cells in clust
em_cdcmip(nemobj)       R tquan(13) Frac MIP in CDC track
em_fdcmip(nemobj)       R tquan(14) Frac MIP in FDC track
em_vtxmip(nemobj)       R tquan(15) Frac MIP in VTX track
                                trd info
em_trd_MEAN(nemobj)     R tquan(19)
em_trd_ACCP(nemobj)     R tquan(22)
em_trd_EFFC(nemobj)     R tquan(23)
                                info from Clean_Photon
em_flgcd(nemobj)        R tquan(24)
em_rhvtxw(nemobj)       R tquan(25)
em_rhcdcw(nemobj)       R tquan(26)
em_nhvtxxy(nemobj)      R tquan(27)
em_nhvtx3d(nemobj)      R tquan(28)
em_nhcdcxy(nemobj)      R tquan(29)
em_nhcdc3d(nemobj)      R tquan(30)
em_nhcdczs(nemobj)      R tquan(31)
em_ccflag(nemobj)       L passed ieta and crack cuts in CC
em_ecflag(nemobj)       L passed ieta cuts in EC

nemobj2                 I Identical to nemobj1 [2]
em_ntrks(nemobj2)       I  q(lclus+21) Num of tracks in road
em_em1(nemobj2)         R\  From CEMENR.FOR
em_em2(nemobj2)         R \   fractional energies in each layer of
em_em3(nemobj2)         R  |  cluster shower
em_em4(nemobj2)         R /
em_fh1(nemobj2)         R/
                        electron likelihoods and errors
em_5varlkl(nemobj2)     R<-|TRD+chisq+emf+trksig+dedx 
em_5varlklerr(nemobj2)  R  |
em_4varlkl(nemobj2)     R  |chisq+emf+trksig+dedx
em_4varlklerr(nemobj2)  R  |
em_3varlkl(nemobj2)     R  |chisq+emf+trksig
em_3varlklerr(nemobj2)  R  |
em_2varlkl(nemobj2)     R  |trd+dedx
em_2varlklerr(nemobj2)  R  |
em_2ntvarlkl(nemobj2)   R  |chisq+emf
em_2ntvarlklerr(nemobj2)R<-|

Notes:
[1] 1st bit - passes ELE shape cuts,
    2nd bit - passes ELE shape cuts and track match,
    3rd bit - passes ELE shape cuts and 0.4 isolation
    4th bit - passes ELE shape cuts and 0.6 isolation
    5th bit - version of L2EM bank is >= 3
[2] This block is actually split among to two blocks, thus the
number of EM objects must be repeated.

EMT block:

PELC tracks for EM objects. Entry i corresponds to the same ith entry in the EM block.
nemtrk                  I  Number of EM objects (= nemobj)
em_dtexist(nemtrk)      L  associated CDC track?
em_vtexist(nemtrk)      L  associated VTX track?
em_ftexist(nemtrk)      L  associated FDC track?
em_dtphi(nemtrk)        R  q(ldtrk+6)  CDC phi
em_dtx0(nemtrk)         R  q(ldtrk+7)      x c.o.g. of track
em_dty0(nemtrk)         R  q(ldtrk+8)      y c.o.g. of track
em_dttheta(nemtrk)      R  q(ldtrk+9)      theta
em_dtz0atr0(nemtrk)     R  q(ldtrk+11)     z0 at R0 point
em_dtchixy(nemtrk)      R  q(ldtrk+12)     chisq of xy fit
em_dtchirz(nemtrk)      R  q(ldtrk+13)     chisq of rz fit
em_dtndof(nemtrk)       I  iq(ldtrk+14)    total num of deg of free
em_dtphierr(nemtrk)     R  q(ldtrk+16)     error on phi
em_dtthetaerr(nemtrk)   R  q(ldtrk+18)     error on theta
em_ftstatus(nemtrk)     I  iq(lfdct+1) FDC bitwise track status
em_ftx0(nemtrk)         R  q(lfdct+4)      x at z=z0
em_fty0(nemtrk)         R  q(lfdct+5)      y at z=z0
em_ftphi(nemtrk)        R  q(lfdct+6)      phi
em_ftchi(nemtrk)        R  q(lfdct+19)     chisq
em_fttheta(nemtrk)      R  q(lfdct+22)     theta
em_ftphierr(nemtrk)     R  q(lfdct+23)     error on phi
em_ftthetaerr(nemtrk)   R  q(lfdct+24)     error on theta
em_ftnptsfit(nemtrk)    R  iq(lfdct+25)    num pts used in track fit
em_vtphi(nemtrk)        R  q(lvtxt+6)  VTX phi
em_vtx0(nemtrk)         R  q(lvtxt+7)      x c.o.g. 
em_vty0(nemtrk)         R  q(lvtxt+8)      y c.o.g.
em_vttheta(nemtrk)      R  q(lvtxt+9)      theta
em_vtvgthcov(nemtrk)    R  q(lvtxt+10)     
em_vtzg(nemtrk)         R  q(lvtxt+11)
em_vtdzdr(nemtrk)       R  q(lvtxt+14)     dzdr of road center
em_vtzvtx(nemtrk)       R  q(lvtxt+15)     zvtx used in road
em_vtphierr(nemtrk)     R  q(lvtxt+16)     error of phi
em_vtxyerr(nemtrk)      R  q(lvtxt+17)     error of cog xy proj
em_vtthetaerr(nemtrk)   R  q(lvtxt+18)     error of theta
em_vtrzerr(nemtrk)      R  q(lvtxt+19)     error of cog rz proj

EMV block:

Correction information for EM objects.
nemcor                  I  Number of EM objects (= nemobj)

            *** Cafixed and uncorrected comparison ***
em_vcdpx(nemcor)        R  q(lvcor+3)  ExCafixed - ExUncorrected
em_vcdpy(nemcor)        R  q(lvcor+4)  Ey
em_vcdpz(nemcor)        R  q(lvcor+5)  Ez
em_vcdp(nemcor)         R  q(lvcor+6)  E
em_vcdet(nemcor)        R  q(lvcor+11) Et

            *** Revertexing with PELC track ***
em_pvtx_id(nemcor)      I  revertexing ID [1]
em_pvtx_tktyp(nemcor)   I  6=VTX,7=CDC,8=FDC
em_pvtx_tkmch(nemcor)   R  tquan(12)  tk mch sig 
em_pvtx_dedx(nemcor)    R  q(ltrak+20) dEdx of track
em_pvtx_theta(nemcor)   R  <-| 
em_pvtx_deta(nemcor)    R    | Detector eta
em_pvtx_et(nemcor)      R    | New Cluster quantities from
em_pvtx_vert(nemcor)    R    | repositioned vertex (em_pvtx_vert)
em_cm3posx(nemcor)      R    | 
em_cm3posy(nemcor)      R    |
em_cm3posz(nemcor)      R  <-|

           *** Revertexing with a track with the best track match
               significance to the EM cluster (PELCs and PPHOs)   ***
em_cvtx_id(nemcor)      I  revertexing ID [2]
em_cvtx_tktyp(nemcor)   I  6=VTX,7=CDC,8=FDC
em_cvtx_tkmch(nemcor)   R  from ELE_TRACKS  tk mch sig 
em_cvtx_dedx(nemcor)    R  from ELE_TRACKS
em_cvtx_theta(nemcor)   R  <-| 
em_cvtx_deta(nemcor)    R    | New cluster quantities from 
em_cvtx_et(nemcor)      R    | repositioned vertex (em_cvtx_vert)
em_cvtx_vert(nemcor)    R    |
em_cvtx_x_pv(nemcor)    R    |
em_cvtx_y_pv(nemcor)    R    |
em_cvtx_z_pv(nemcor)    R  <-|

Notes:
[1] Names of many entries changed compared to OMNtuple version 7
[2] New vertex is assigned a unique, arbitrary id number.
    Is 0 if error or no track (PPHO).
[3] New vertex is assigned a unique, arbitrary id number.
    Is 0 if error or no track found.
    If the new vertex is the same as the one discovered with 
    the PELC track, em_cvtx_id(i) = -em_pvtx_id(i)

MUON block:

Muon quantities (each array has nmuo entries). Entries are from the PMUO and MUOT banks and are sorted by muon pT (highest first).
nmuo                    Number of pmuo banks
muo_chg(numo)         I Charge of muon (-14 = -mu, +14 = +mu)
muo_fit(nmuo)         I Fit status (pmuo+4)
muo_cdtr(nmuo)        I Number of CD tracks in cone
muo_quad(nmuo)        I Quadrant
muo_ifw4(nmuo)        I Quality flag (IFW4)
muo_p(nmuo)           R Total momentum of muon
muo_pt(nmuo)          R Transverse momentum
muo_theta(nmuo)       R Theta direction
muo_eta(nmuo)         R Eta
muo_phi(nmuo)         R Phi
muo_chi(nmuo)         R Chi squared / degrees of freedom
muo_flt0(nmuo)        R Floating T0 offset (ns)
muo_isoc(nmuo)        R Isolation parameter (cells hit only)
muo_isocn(nmuo)       R Isolation parameter (cells hit + neighbors)
muo_isoc2n(nmuo)      R Isolation parameter (cells hit + 2 neighbors)
muo_iso4(nmuo)        R Isolation parameter (cone size 0.4)
muo_elecal(nmuo)      R E loss expected in calorimeter
muo_ecalc1(nmuo)      R E observed in calorimeter (cells hit + one neighbor)
muo_ecalc2(nmuo)      R E observed in cal (cells hit + two neighbors)
muo_ecal4(nmuo)       R E observed in cal (cone size 0.4)
muo_dphi(nmuo)        R Change in phi (degrees)
muo_dth(nmuo)         R Change in theta (degrees)
muo_impv(nmuo)        R Impact parameter from vertex
muo_impf(nmuo)        R Impact parameter after global fit
muo_impbm(nmuo)       R Impact parameter bend view, muon only
muo_impnbm(nmuo)      R Impact parameter non-bend view, muon only
muo_impbc(nmuo)       R Impact parameter bend view, CD
muo_impnbc(nmuo)      R Impact parameter, non-bend view, CD
muo_hont(nmuo)        I Hits on track, A,B,C
muo_hontf(nmuo)       I Hits on track fit, A,B,C
muo_vert(nmuo)        I Vertex used
muo_bdl(nmuo)         R B . dl
muo_l1conf(nmuo)    ! I Level 1 confirmation
muo_l15conf(nmuo)   ! I Level 1.5 confirmation
muo_l2pt(nmuo)      ! R Level 2 Pt
muo_l2phi(nmuo)     ! R Level 2 Phi
muo_alayerx(nmuo)   ! R x position of a-layer hit
muo_alayery(nmuo)   ! R y position of a-layer hit
muo_alayerz(nmuo)   ! R z position of a-layer hit
muo_ifw1(nmuo)     *! I Flag word one
muo_hfrac(nmuo)    *  R MTC Hfrac value
muo_frac(nmuo)     *  R MTC Frac value
muo_efrac(nmuo)    *  R MTC efrac_h(1) value
muo_drjet(nmuo)    *  R Distance in eta x phi of closest jet
muo_cljet(nmuo)    *  I Jet number (from jet block) of the closest jet

JET block:

Jets information. Only one jet algorithm/cone size is supported. Usual cone size is DR = 0.5. The jets are sorted by descending ET.
njets                 I Number of jet banks
jet_emoverlap(njets)  I Entry number in EM block of EM object overlapping jet
jet_jnepn(njets)      I Entry number in JNEP block of JNEP associated with jet
jet_hasjnep(njets)    L True if this jet has a jnep bank associated with it
jet_px(njets)         R Px of jet (from JETS)
jet_py(njets)         R Py
jet_pz(njets)         R Pz
jet_e(njets)          R Energy of jet
jet_et(njets)         R Et of jet
jet_eta(njets)        R Eta of jet
jet_ieta(njets)       I detector eta of jet
jet_phi(njets)        R phi of jet
jet_sigex2(njets)   ! R SIGMA(Ex)**2
jet_sigey2(njets)   ! R
jet_sigez2(njets)   ! R
jet_etaw(njets)       R RMS eta width
jet_phiw(njets)       R RMS phi width
jet_emfrac(njets)     R EM fraction 
jet_msflag(njets)     I Flag for merging and splitting
jet_ncath(njets)      I Number of cells above threshold
jet_ficdmg(njets)     R Fraction of ICD/MG Et(EtICD+EtMG/TOTAL_ET)
jet_fchet(njets)      R Fraction of CH Et (Et CH/TOTAL_ET)
jet_rhtnh(njets)      R Ratio of hottest to next hottest cell
jet_ntow(njets)       I Number of TOWERS comprising 90% of jet Et
jet_cone(njets)       R ConeSize (-1.0 for Nearest Neighbor)
jet_metphi(njets)     R Angle between missing Et (PNUT(4)) and jet 
jet_ecorwd(njets)   ! I Energy correction status word
jet_vtxwd(njets)      I Vertex matching word - 4 bytes
jet_ecor(njets)     ! R Energy Correction Factor
jet_etcor(njets)    ! R Transverse Energy Correction Factor
jet_etacor(njets)   ! R Eta correction made to jet
jet_emfcor(njets)   ! R Emf correction made to jet
jet_phicor(njets)   ! R Phi correction made to jet
jet_vtxok(njets) ***! R Error code returned by vertex_confirm_jet [1] 
jet_zvtx(njets)  ***! R primary vertex position found by the jet axis pointing
jet_etax(njets)  ***! R eta of jet axis in the eta-z plane
Notes:
[1] jet_vtxok codes are as follows:
    0 = jet axis pointing is ok 
    1 = jet 10<Et<=25, otherwise ok
    2 = there are less then 4 layers with at least 2 GeV deposited energy in jet 
    3 = jet is not central or jet Et < 10 GeV or nvtx=0

JNEP block:

Information on jets with EM object cells excised. To find the entry in the jnep block corresponding to a jet, do the following:
C*****Example of accessing JNEP information for jet 3.      
      IF ( jet_hasjnep(3) ) THEN   ! Does jet have a jnep entry?        
        w = jet_jnepn(3) ! w holds the entry number in JNEP for this jet
        CALL HFILL(100, jnep_et(w), 0., 1.) ! Plot ET without electron
      ENDIF
There is one entry for each JNEP bank
njnep                 I Number of jnep entries
jnep_px(njnep)        R Px of non/pelc/ppho part of jet (FROM JNEP)
jnep_py(njnep)        R
jnep_pz(njnep)        R
jnet_e(njnep)         R
jnet_et(njnep)        R
jnep_eta(njnep)       R Similar to above but from (JNEP)
jnep_phi(njnep)       R 
jnep_emfrac(njnep)    R

UJET block:

Uncorrected jets (same algorithm and cone size as used for the JETS block).
nujets                  I Number of jets (same as njets)
ujet_px(nujets)         R Px of uncorrected jet
ujet_py(nujets)         R Py
ujet_pz(nujets)         R Pz
ujet_e(nujets)          R Energy of uncorrected jet
ujet_et(nujets)         R Et of uncorrected jet
ujet_eta(nujets)        R Eta of uncorrected jet
ujet_ieta(nujets)       I Detector eta of uncorrected jet
ujet_phi(nujets)        R Phi of uncorrected jet
ujet_emfrac(nujets)     R EM fraction of uncorrected jet
ujet_ficdmg(nujets)     R Fraction if ICD+MG Et
ujet_cone(nujets)       R Cone size for uncorrected jet
ujet_metphi(nujets)     R Angle between missing Et (always PNUT 2) and 
                             uncorrected jet

L1L2JET block:

Information about level 1, large tile and level 2 jets. If a micro-dst, information comes from JUTL bank, otherwise from ESUM.
nl1jets              * I Number of level 1 jets
l1jet_et(nl1jets)    * R Et of level 1 jet
l1jet_eta(nl1jets)   * R Eta of level 1 jet
l1jet_phi(nl1jets)   * R Phi of level 1 jet [1]
nltjets              * I Number of large tile jets
ltjet_et(nltjets)    * R Et of large tile jet
ltjet_eta(nltjets)   * R Eta
ltjet_phi(nltjets)   * R Phi [1]
nl2jets              * I Number of level 2 jets (total over all cones)
l2jet_cone(nl2jets)  * R Cone size of level 2 jet
l2jet_et(nl2jets)    * R Et of level 2 jet
l2jet_eta(nl2jets)   * R Eta
l2jet_phi(nl2jets)   * R Phi [1]

Notes:
[1] = Phi valuses of trigger/filter jets may be corrupted. 
      A -1 is stored for the phi value if that is the case.

MET block:

Missing ET information. There is one entry for each PNUT bank. For example, met_et(2) corresponds to PNUT(2); met_et(4) corresponds to PNUT(4), if it exists.
nmet               I Number of PNUT banks (should be 5 for CAFIXed)
met_ex(nmet)       R Px of missing energy (From PNUT)
met_ey(nmet)       R
met_ez(nmet)       R
met_e(nmet)        R Total missing energy
met_et(nmet)       R Missing transverse energy
met_eta(nmet)      R Eta of missing Et
met_phi(nmet)      R Phi of missing Et
met_siget(nmet)  ! R Sigma(MEt)**2
met_etscal(nmet)   R Scalar Et

HCIJ block:

Information about hot cells in jets. Any hot cell within R=0.5 of a reconstructed jet is declared to be "within a jet".
nhcij                   * I Number of hot cells found within jets
hcij_met                * R Missing Et with hot cells in jets added back in
hcij_metx               * R x component of missing Et above
hcij_mety               * R y component of missing Et above
hcij_metphi             * R Phi of missing Et above
hcij_njets              * I Number of jets with hot cells found within
hcij_tojet(hcij_njets)  * I Jet in jets block that has a hot cell found within
hcij_n(hcij_njets)      * I Number of hot cells found within jet
hcij_dele(hcij_njets)   * R Total energy of hot cells found within this jet
hcij_delet(hcij_njets)  * R Total Et of hot cells found within this jet
hcij_et(hcij_njets)     * R New jet Et with hot cells added back in [1]
hcij_e(hcij_njets)      * R New jet E "" [1]
hcij_emfrac(hcij_njets) * R New EM fraction for jet with hot cells added back [1]
hcij_fchet(hcij_njets)  * R New CH fraction for jet with hot cells added back [1]

Notes:
[1] Cafix corrections are applied, so all quantities are corrected (but correction
    is based on original jet Et).

TAU block:

Tau information the tau group uses.
ntau               I Number of ptau banks
tau_px(ntau)       R Px of 0.7 cone tau jet (from JETS bank, *not* PTAU)
tau_py(ntau)       R Py
tau_pz(ntau)       R Pz
tau_et(ntau)       R Et of 0.7 cone tau jet (from JETS bank, *not* PTAU)
tau_eta(ntau)      R Eta of 0.7 cone tau jet (from JETS bank, *not* PTAU)
tau_phi(ntau)      R Phi of 0.7 cone tau jet (from JETS bank, *not* PTAU)
tau_emfrac(ntau)   R EM fraction of 0.7 cone tau jet (from JETS bank, *not* PTAU)
tau_ethot1(ntau)   R Et of hottest cell in tau jet (PTAU bank)
tau_ethot2(ntau)   R Et of 2nd hottest cell in tau jet (PTAU bank)
tau_etawid(ntau)   R Eta width of 0.7 cone tau jet (from JETS bank, *not* PTAU)
tau_phiwid(ntau)   R Phi width of 0.7 cone tau jet (from JETS bnak, *not* PTAU)
tau_r(ntau)        R Width of tau jet (R = sqrt(deta**2 + dphi**2)) (from PTAU)

DTRK block:

Information on CDC tracks. Note that microdsts only have tracks in roads! DSTs and STAs contain all CDC tracks.
ndtrk                   I Number of CDC tracks (150 max)
dtrk_ok(ndtrk)          I Track quality code [1]
dtrk_x0(ndtrk)          R x position of track's center of gravity
dtrk_y0(ndtrk)          R y position of track's center of gravity
dtrk_r0(ndtrk)          R r (sqrt(x^2+y^2))  position of track's center of gravity
dtrk_z0(ndtrk)          R z position of track's center of gravity
dtrk_theta(ndtrk)       R Theta direction of track
dtrk_phi(ndtrk)         R Phi direction of track
dtrk_dedx(ndtrk)        R dE/dx of track
dtrk_zcross(ndtrk)      R z position where track crosses R=0
dtrk_xyimp(ndtrk)       R Track impact parameter in x-y plane
dtrk_vert(ndtrk)        I Entry in VERTEX block of vertex within 10 cm of track
dtrk_jet(ndtrk)         I Entry in JET block of jet associated with this track [2]
dtrk_jetdr(ndtrk)       R Distance in eta-phi space of track to jet axis
dtrk_em(ndtrk)          I Entry in EM block of EM object associated with this track
dtrk_emdr(ndtrk)        R Distance in eta-phi space of track to EM cluster center
dtrk_trkmat(ndtrk)      R track match?
dtrk_nhxy(ndtrk)    *** I Number of track hits in X-Y plane
dtrk_nhrz(ndtrk)    *** I Number of track hits in R-Z plane 
dtrk_chixy(ndtrk)   *** R chi2/(Nhxy-2) 
dtrk_chirz(ndtrk)   *** R chi2/(Nhrz-2)
dtrk_eta(ndtrk)     *** R track eta
dtrk_deta(ndtrk)    *** R track detector eta  

Notes:
[1] dtrk_ok codes are as follows:
    0 = Track is ok
    1 = No R-Z information
    2 = Too few hits
    3 = Bad chi square
    4 = XY impact parameter to far from r=0

[2] Looks within a cone DR from the jet axis for tracks.  The
    distance is set in the omntuple.rcp file.  Similar for dtrk_em.
    If no jet or EM object can be associated with the track, the 
    respective field (dtrk_jet or dtrk_em) is 0.

ISAJET block:

Some ISAJET information
nisv1              I Number of isv1 banks (up to 4)
isv1_z(nisv1)      R Z vertex given from isv1 bank (ISV1)
ipart(2)           I ISAJET particle number of ISAJ particle (primary interaction) [1]
iscolor(2)         L True if ISAJ particle is SUSY and carries color (squark or gluino)
imass(2)           R Mass of primary ISAJ particle
nislep             I Number of isajet leptons (from ISAL)
islep_id(nislep)   I ID of isajet lepton
islep_pt(nislep)   R Pt of isajet lepton
islep_phi(nislep)  R Phi of isajet lepton
lslep_eta(nislep)  R Eta of isajet lepton
isamet             R Missing Et from isajet particles [1]
isametphi          R Phi of isajet missing Et

Notes:
[1] Always two entries, one for each produced primary isajet particle.

[2] The isajet missing ET is calculated from all particles that interact in
    the calorimeter in the ISP1 banks (muons, neutrinos, sneutrinos and LSPs
    are excluded).

PJET block:

PJet information. The first PJET algorithm is used.
npjet             I Number of pjets
pjet_px(npjet)    R Px of Pjet
pjet_py(npjet)    R Py
pjet_pz(npjet)    R Pz
pjet_e(npjet)     R E
pjet_et(npjet)    R ET
pjet_eta(npjet)   R eta
pjet_phi(npjet)   R phi

OMET block:

Missing ET calculated from other vertices. CAEP or CEAQ banks are required, so nothing filled if run on microdsts. The order of the entries are as follows: The vertex number corresponds to the entry in the VERTEX block.
nomet                   ! I Number of missing_et's calculated
omet_vtxn(nomet)        ! I Number of vertex this MEt is calc for (entry in VERTEX block)
omet_met(nomet)         ! R Magintude of missing transverse energy
omet_metphi(nomet)      ! R Phi direction of MEt
omet_mez(nomet)         ! R z component of missing energy (watch out!)
omet_scet(nomet)        ! R Scalar Et for this vertex

SCANMET block:

Missing ET is calculated at various "fake" vertices along the z axis, starting at -70 cm to +70 cm in steps of 10 cm. This block also requires CAEP or CAEQ banks.
scanmet(15)     ! R Missing Et at vertex at z = -80+10*i cm, i=1...15

MTC block:

Muon MIP tracks found in the calorimeter by MTC
mtc_ntrks               I Number of tracks found
mtc_ierror(mtc_ntrks)   I Error report from MTC [1]
mtc_eta(mtc_ntrks)      R Eta of track
mtc_phi(mtc_ntrks)      R Phi of track
mtc_e(mtc_ntrks)        R Energy of track
mtc_fr(mtc_ntrks)       R Fraction of layers utilized with nonzero energy
mtc_hdfr(mtc_ntrks)     R Fraction of hadronic layers
mtc_tres(mtc_ntrks)     R Track residual

Notes:
[1] 0 = ok, -1 = failed, +1 = too many tracks

ASTUB block:

Prototype A-phi scintillator counter information.
as_npmt              I Number of counters examined (same for every event)
as_nhit              I Number of counters actually hit in the event
as_addr(as_npmt)     I Address of hit counter (0 if counter not hit) [1]
as_time(as_npmt)     I Time of hit (-1000 if counter not hit)
       
Notes:
[1] as_addr(i) will be the same for each event if the counter was hit, and 0
    if not hit.

 QSIM block:
Quick SIMulation information.
 ISAJET information:
 iparton1                I 1st primary particle id   
 iparton2                I 2nd primary particle id [1] 
 nele_gen                I number of electrons   
 ele_pt_gen(nele_gen)    R Et of electron (E*sin(Thet))
 ele_eta_gen(nele_gen)   R eta of electron 
 ele_phi_gen(nele_gen)   R phi of electron  
 nmu_gen                 I number of muons 
 muo_pt_gen(nmu_gen)     R Et of muon (E*sin(Thet))
 muo_eta_gen(nmu_gen)    R eta of muon
 muo_phi_gen(nmu_gen)    R phi of muon
 PJET information (d0pjet):
 njet_gen                I number of jets 
 jet_pt_gen(njet_gen)    R Pt of jet (sqrt(Px**2+Py**2))
 jet_eta_gen(njet_gen)   R eta of jet 
 jet_phi_gen(njet_gen)   R phi of jet
 jet_et_gen(njet_gen)    R Et of jet (E*sin(Thet))
 Calculated by QSIM using ISAJET and PJET information:
 met_gen                 R missing Et using Pt=sqrt(sum(Px)^2+sum(Py)^2) of 
                           isv1 particles excluding neutrinos and LSP's
 met_phi_gen             R phi of missing Et of jets
 scalar_et_gen           R scalar sum of Pt=sqrt(Px^2+Py^2) of
                           isv1 particles excluding neutrinos and LSP's
 met_mu_gen              R is not implemented
 met_phi_mu_gen          R is not implemented 
 scalar_et_mu_gen        R is not implemented
 ISAJET and PJET information smeared by QSIM:
 nele_sm                 I number of electrons (the same as unsmeared one)
 ele_pt_sm(nele_sm)      R Et of electron smeared by d0smear 
 ele_eta_sm(nele_sm)     R eta of electron (the same as unsmeared one)
 ele_phi_sm(nele_sm)     R phi of electron (the same as unsmeared one)
 ele_deta_sm(nele_sm)    R detector eta of electron (the same as unsmeared one)
 ele_dr_jet(nele_sm)     R delta R between electron and the closest jet
 iele_jet(nele_sm)       R closest jet number 
 nmu_sm                  I number of muons (the same as unsmeared one)
 muo_pt_sm(nmu_sm)       R Et of muon smeared by d0smear
 muo_eta_sm(nmu_sm)      R eta of electron (the same as unsmeared one)
 muo_phi_sm(nmu_sm)      R phi of electron (the same as unsmeared one) 
 njet_sm                 I number of jets (the same as unsmeared one)
 jet_pt_sm(njet_sm)      R Et of jet smeared by d0smear
 jet_eta_sm(njet_sm)     R eta of jet (the same as unsmeared one)
 jet_phi_sm(njet_sm)     R phi of jet (the same as unsmeared one)
 jet_deta_sm(njet_sm)    R detector eta of jet (the same as unsmeared one)
 jet_dr_ele(njet_sm)     R delta R between jet and the closest electron
 ijet_ele(njet_sm)       R delta R between jet and the closest electron
 met_sm                  R smeared missing Et 
 met_phi_sm              R phi of smeared missing Et
 scalar_et_sm            R scalar sum of jet_pt_sm (and ele_pt_sm if it is close to jet)
 met_mu_sm               R is not implemented
 met_phi_mu_sm           R is not implemented
 scalar_et_mu_sm         R is not implemented
 zvert_gen               R z-position of ISAJET isv1 primary vertex
 zvert_sm                R the same as zvert_gen
 Are not used:

 nele_par                   
 ele_par1(nele_par)
 ele_par2(nele_par)
 ele_par3(nele_par)
 ele_par4(nele_par)
 ele_par5(nele_par)
 nmu_par
 muo_par1(nmu_par)
 muo_par2(nmu_par)
 muo_par3(nmu_par)        
 muo_par4(nmu_par)
 muo_par5(nmu_par)
 njet_par
 jet_par1(njet_par)
 jet_par2(njet_par)
 jet_par3(njet_par)
 jet_par4(njet_par)
 jet_par5(njet_par)
 met_par1
 met_par2
 met_par3
 met_par4
 met_par5
 Notes:
[1] Always two primary isajet particle (iparton1 >= iparton2)