The combined ntuple is a set of blocks. All blocks can be set on or off via datacards. A group of block (tracking information,  jet information,....) can be set off with the datacards BASIC (TRAC.BASIC, JFNT.BASIC,....).  The maximum number of objects in each block can also be set via datacards. It is up to the user to chose the blocks he wants for a specific analysis (but the maximum number of columns should not exceed 50000 (Max. Columns, bottom right when typing nt/pri 3333) .

It is up to the user to make sure the information he requires has been reconstructed by having the appropriate reconstruction packages called (documentation on atrecon is here).

The following information is available: KINE, Inner Detector Tracks, Conversions, General Calorimeter, Missing Energy,  Jets, Tile Calorimeter Cells, Calorimeter projected matrix, Electromagnetic clusters, Muons, MUon IDentification, ATLFAST , User

Run, event info

always present.

Variables:

  • RUN           run number
  • EVENT     event number

    • KINE/GENZ information.

    The GENZ bank contain all the tracks generated by the generator (e.g PYTHIA). All the final GENZ tracks are copied in the KINE bank. In addition, the KINE bank contain all the particles created by GEANT (interaction or decay) (depending on how the generator was run strange particle maybe decayed by the generator or by GEANT (more correct)). Only a selection of these particles are saved in the KINE/GENZ ntuple block. A final GENZ particle being duplicated in a KINE particle, they are saved as a single particle.

    Driving data cards:

  • KINE.BASIC    1 (default) fill KINE+GENZ blocks (if 0 do not)
  • KINE.NPARMAX  600 (default,<10000) maximum number of particles (if 0 do not fill particle block)
  • KINE.NVERMAX  100 (default,<10000) maximum number of particles (if 0 do not fill vertex block)
  • KINE.PILEUP 0(default) no-pile-up, 1 + SCT+PIXEL pile-up (~24 Min. Bias), 2 + TRT pile-up (~32 Min. Bias), 3 full pile-up (~48 Min. Bias) (only KINE pile-up is saved, not GENZ)
  • KINE.GENZ  0 nothing  1 (default) only hard scatter and track mother of saved vertices , 2 plus hard scatter daughters 3 full GENZ bank
  • KINE.BPHYS 0 (default) no 1 all B,D hadrons and tau are saved 2 plus all daughters down to primary KINE tracks
  • MUPTMIN  (3. default) minimum PT to save a KINE muon (if negative do not save muons), with |eta|<2.75
  • ELPTMIN (0.5 default) minimum PT to save a KINE electron (if negative do not save electrons) , with |eta|<2.55.
  • GAPTMIN (10. default) minimum PT to save a gamma (if negative do not save gamma ), with | eta|<2.55.
  • JETS  0 (default) no 1 save all KINE tracks contributing to jets i.e all final KINE tracks (no further daughters) with JFNT.ETAMIN<eta<JFNT.ETAMAX which is -5<eta<5 by default
  • TRACK 0  no 1 (default) save all potentially vsible charged KINE tracks: PT>0.4 GeV and |ETA|<2.55 and inside a reconstruction road with pT>pTmin-0.1, or referred to (KINEREF) by any reconstructed charged track  2 do not check inside reconstruction road
  • MOTHER 0 no 1 (default)  save mother of  KINE particle selected above up to primary KINE track 2 idem but go up to parent string (excluded)
  • DAUGHTER 0 (default) 1 save daughter of selected KINE particle down to the last 2 save also daughter of the mothers (selected with MOTHER flag)
  • VERT 0 save only primary vertex 1 (default) save all vertices referred to by one of the particle selected above, as origin or end vertex
  • DIRVERT 0 (default) 1 copy vertex info in vertex block in new variables in kine/genz trackblock. This allows easier interactive use of vertex info but does not add information (see below) and takes more room if lots of tracks have to be saved

    •  

     

    By default are saved visible electron, photon, muon and charged tracks, with their mother up to the primary.  If you want to save all the KINE
    bank (ony particle with eta>5 will be missing) , then use JETS=1 MOTHER=1.
     If no GENZ bank is present or if it contain only 1 track, KINE track number 1 is saved anyway, as it is probably a single track event.
    Genz particle are saved with increasing GENZ number, KINE only particles are saved at the end with increasing KINE number.  Final GENZ/ primary KINE particles duplicates are interleaved with GENZ only particles. The same is true for vertices.
     
     

    Useful formulae: PX=PTGEN*COS(PHIGEN),
                                     PY=PTGEN*SIN(PHIGEN),
                                     PZ=PTGEN*SINH(ETAGEN),
                                     P=PTGEN*COSH(ETAGEN),
                                     COTAN=SINH(ETAGEN),
                                     A0=RVGEN*SIN(PHIVGEN-PHIGEN)

     

    Variables:

    The following block is created as soon as NPARMAX>0
  •  NPAR   number of true particles  saved
  • IFGENZ index of first GENZ particle (1 if any GENZ particle 0 if not )
  • ILGENZ index of last GENZ particle (-1 if no GENZ particle) DO I=IFGENZ,ILGENZ  will loop on all GENZ particle
  • IFKINE index of first KINE particle (0 if no KINE particle)
  • ILKINE index of last KINE particle (NPAR if some KINE particle, -1 if no KINE particle ) WARNING all tracks between IFKINE and ILKINE are not necessary KINE particles (a check on GENREF should be done (see below))
  •  TYPE(NPAR) 10000 *PDG code of the particle +GEANT code of particle (lots of unstable particle have 0 GEANT code, very few particles (nuclei) have 0 PDG code). Useful formula: GEANT=ABS(MOD(TYPE,10000)), PDG=INT(TYPE/10000)
  •  PTGEN(NPAR)  Pt
  •  PHIGEN(NPAR) phi
  •  ETAGEN(NPAR) eta
  •  MGEN(NPAR)  mass
  •  CHARGE(NPAR)  charge
  • GENSTAT(NPAR) for a GENZ only particle GENSTAT=GENZ  status code (IPSTAT=3 history particle, 2 for decayed particle,....).

    •    For a KINE particle or a stable GENZ particle (MOD(GENREF,100000)>0) GENSTAT=(10000000*BCOD+1000000*DCOD+100000*TAUCOD+10000*SECOCOD+GENZ status code)*SGN, where
             BCOD is the PDG code of the second quark in the B hadron (plus 5 if it is a baryon ) (1 B0 2 B+ 3 Bs 4 Bc 5 bottomonium >=6 b-baryon), which is the ancestor of the particle
             DCOD is the same for D hadron (1 D+ 2 D0 3 Ds 4 charmonia >=6 c-baryon),
             TAUCOD is 1 if one ancestor to the particle is  a tau lepton,
             SECOCOD=1 for bremstrahlung or photon conversion,
                                      2 K0 decay,
                                      3 /\ decay,
                                      4 other strange baryon decay,
                                      5 K+ decay/interaction energy conserved,
                                      6 pi+ decay/interaction  energy conserved,
                                      7 any interaction (energy not conserved)
                                      8 mu "decay" copy in spectrometer,
                                      9 anything else
                                   (if several cascade interaction/decay only the last one is used),
             SGN is the sign of the PDG code of the heaviest particle which was used in the calculation of genstat (e.g gives the flavor of the B hadron if BCOD>0)
  •  GENREF(NPAR) 100000*(original GENZ number)+original KINE number. Useful formula GENZ=INT(GENREF/100000)., KINE=MOD(GENREF,100000),.  Unstable particles generated by the generator (e.g PYTHIA)  have GENZ>0 and KINE=0. Stable particles generated by the generator have GENZ>0 and KINE>0, since they are duplicated in GENZ and KINE banks, except if they were filtered out at DICE level (usually if eta>|2.7|) in which case they have KINE=0. Particles generated by GEANT (interaction with matter or decay of long-lived particle) have GENZ=0 and KINE>0.
  • KMOTHNT(NPAR) index of mother in the Ntuple. If the mother was not saved, KMOTHNT=-GENREF of the mother.
  • KDAUNT(NPAR) 100000*(index of first daughter in Ntuple)+(index of last daughter in ntuple) WARNING as in the original GENZ and KINE banks all tracks between first and last daughters are not necessarily daughters. This must be checked. A loop on daughters of track ITK should read

    •         IF (KDAUNT(ITK).GT.0) THEN
                DO I=INT(KDAUNT/100000), MOD(KDAUNT,100000)
                      IF (KMOTHNT(I).EQ.ITK) THEN
                         XXX
                     ENDIF
                ENDDO
            ENDIF
           If no daughters are saved in the ntuple, then KDAUNT is -GENREF of the first daughter.
  • KVNT(NPAR) 100000*(index of  origin vertex in ntuple)+(index of end vertex in ntuple)

    •  

     

    The following variables are saved only if DIRVERT=1. Note that if the vertex block is saved, this is only a duplication of the information contained in the vertex block (see below), so that RVGEN(ITK)=RV(INT(KVNT(ITK)/100000). This is done only because nt/plo RV(INT(KVNT/100000)) is not possible in paw.

  •  RVGEN(NPAR)  radius of origin vertex
  •  PHIVGEN(NPAR) Phi of origin vertex
  •  ZVGEN(NPAR)    Z of origin vertex
  •  BUNCHNUM(NPAR)   Bunch number (main event 0, pile-up -6 to -1, in average 8 event per bunch number). Note that  since August 1998 (new calorimeter pile-up scheme) , some of the pile-up events have 0 time of flight in KINE bank. To avoid confusion these are labelled -4 in the ntuple.It is also the case for low-luminosity pile-up.

    •  

     

     The following block is created as soon as NVERMAX>0

  • NVER number of true vertices saved
  • IFVGENZ index of first GENZ vertex  1 if any GENZ vertex 0 if not (a GENZ (KINE) vertex is a vertex which mother is a GENZ  (KINE) track.)
  • ILVGENZ index of last GENZ vertex  -1 if no GENZ vertex
  • IFVKINE index of first KINE vertex 0 if no KINE vertex
  • ILVKINE index of last KINE vertex  -1 if no KINE vertex
  •  RV(NVER)  radius of vertex
  •  PHIV(NVER) Phi of vertex
  •  ZV(NVER)    Z of vertex
  •  BUNCHV(NVER)   Bunch number (main event 0, pile-up -6 to -1, in average 8 event per bunch number)Note that  since August 1998 (new calorimeter pile-up scheme) , some of the pile-up events have 0 time of flight in KINE bank. To avoid confusion these are labelled -4 in the ntuple..It is also the case for low-luminosity pile-up.
  • KORINTV (NVER): index of the mother in track ntuple block. If mother is not in ntuple block -GENREF of the mother.

    •  

     

    Inner Detector Tracks

    Tracks and seeds from iPatrec, Pixlrec, xKalman and xHourec are available. An interface to CTVMFT (CDF vertexing package adapted to Atlas by Francesco Tartarelli) has been done by Paolo Iengo, as documented here by Simon Robins.

    Driving Data cards

  • TRAC.BASIC  1 (default)  if information to be saved
  • TRAC.NTRACMAX 300 (default) maximum number of track to be saved (track may be counted several times in case of multiple fit, multiple pattern recognition used)
  • TRAC.ROAD 1 (default) fill seed=road information
  • TRAC.NSEEDMAX 10 (default) maximum number of road saved
  • TRAC.IPAT 1 (default) save IPATREC tracks
  • TRAC.PIXL 1 (default) save PIXLREC tracks
  • TRAC.XKAL 1 (default) save XKALMAN tracks
  • TRAC.XHOU 1 (0 default) save XHOUREC tracks
  • TRAC.FITCOV 2 (default) save full covariance matrix at vertex  1 only diagonal terms 0 nothing
  • TRAC.FITHIT 1(default) save detailed fit information (in particular hit information)
  • TRAC.FITEND 1 (0 default) save fit parameters at TRT exit
  • TRAC.FITCVE 2 (0 default) save full covariance matrices at TRT exit 1 only diagonal terms 0 nothing
  • TRAC.FITVER 1 (0 default) save vertex fit  (note that PTINV and PHI with vertex constraint are saved anyway in VPTINVVERT and VPHIVERT)
  • TRAC.FITSEED 1 (0 default) save seed fit

    • To keep minimal tracking info  use TRAC.FITEND= 0 FITHIT= 0 ROAD= 0 FITCOV= 0 FITSEED= 0 ROAD=0

    Variables

     
  •  NSEEDS   number of seeds saved
  •  NTRKSEED(NSEEDS) number of track in seed
  •  ROADTYPE(NSEEDS) type of seed
  •  ISEEDref(NSEEDS)=Ipro*1000000+ISeed where IPRO=1 for ipatrec track, 2 pixlrec, 3 xkalman, 4 xhourec, ISeed is the Seed number,
  •  IndexSEED(NSEEDS) index of object used as a seed (e.g number of e.m. cluster)
  •   PHISEED(NSEEDS) phi of seed
  •  ETASEED(NSEEDS) eta of seed
  •  WPHISEED(NSEEDS) phi width of seed
  •  WETASEED(NSEEDS) eta width of seed

    •  
    Useful formulae: PX=COS(PHIVERT)/ABS(PTINVVERT)
                                     PY=SIN(PHIVERT)/ABS(PTINVVERT),
                                     PZ=COTTHVERT/ABS(PTINVVERT),
                                     PT=1/ABS(PTINVVERT),
                                     ETA=-LOG(TAN(ATAN2(1.,COTTHVERT)/2.)), CHARGE=SIGN(1.,PTINVVERT)

     
  •  NHEL    number of tracks ( possibly repeated with different  pattern recognition and/or fits)
  •  ITRKRef(NHEL) =IPro*1000000+ISeed*10000+ITRK where

    •               IPRO=1 for Ipatrec track, 2 pixlrec, 3 xkalman, 4 xhourec
                  ISeed is the Seed number,
                  ITRK the track number for this seed.
                  ITRKRef is an unambiguous track number to be used for reference.
  • ItrkSeedNT index of seed in SEED ntuple block, if any
  •  BVSFIT(NHEL)=1 if basic fit, 2 if vertex fit, 3 if seed fit
  •  CHI2(NHEL)  chi2 of fit
  •  A0VERT(NHEL)  impact parameter
  •  ZVERT(NHEL)  Z at closest approach
  •  PHIVERT(NHEL) Phi at closet approach
  •  COTTHVERT(NHEL) cottan(theta)
  •  PTINVVERT(NHEL) signed inverse pt
  • KINEREF(NHEL)  original KINE track matching this track (correspond to KINEpar in KINE ntuple block)
  • KINENT(NHEL) index of KINE track in KINE ntuple block
  • TKQUAL(NHEL) Track quality flag (HIT+100*TR) HIT=1 if passes b-tag cuts 2 plus NSTRAW>=20, TR=1 if passes low threshold TR electron id cuts 2 if passes high threshold TR cuts

    •  

     

    The following variables are filled only if TRAC.FITHIT=1 (default)

  • TKFITFLAG(NHEL) Track fit flag (XKALMAN: -2 electron fit successful, -1 normal fit)
  •   NSIHITS(NHEL) number of sct+pixel hit
  •  NSIHOLES(NHEL) number of sct+pixel holes
  •  PATTERN(NHEL) bit pattern of sct and pixel hits associated
  •  NSTRAWHITS(NHEL) number of straws associated
  •  NSTRAWHOLE(NHEL) number of straw holes
  •  NSTRAWTIME(NHEL) number of straw drift time hit
  •  NTRHits(NHEL) number of Transition Radiation hits
  •  VPTINVVERT(NHEL) PT at vertex
  •  VPHI(NHEL)      phi at vertex
  •  RSTART(NHEL)   radius of innermost hit
    • Following number of hits are stored as follows: 100000*N b+1000*N SCT+100*N PIXL+N TRT (not fully implemented in all packages)
  •  UNIQUEHITS(NHEL)  number of UNIQUE hits
  •  SPOILTHITS(NHEL)  number of SPOILT hits
  •  WRONGHITS(NHEL)  number of WRONG hits
  •  SHAREDHITS(NHEL)  number of hits SHARED  with another reco track

    •  

     

    Track covariance matrix  on (A0, Z,, PHI, COTTH, PTINV) fully filled if  TRAC.FITCOV=2 (default) (only diagonal if  TRAC.FITCOV=1)

  • NHELCOV  = NHEL (repeated for technical reasons)
  •  COVVERT11(NHELCOV)
  •  COVVERT21(NHELCOV)
  •  COVVERT22(NHELCOV)
  •  COVVERT31(NHELCOV)
  •  COVVERT32(NHELCOV)
  •  COVVERT33(NHELCOV)
  •  COVVERT41(NHELCOV)
  •  COVVERT42(NHELCOV)
  •  COVVERT43(NHELCOV)
  •  COVVERT51(NHELCOV)
  •  COVVERT52(NHELCOV)
  •  COVVERT53(NHELCOV)
  •  COVVERT54(NHELCOV)
  •  COVVERT55(NHELCOV)

    •  

     

    Tracks parameters at ID exit  only filled if  TRAC.FITEND=1

  • NHelEnd  = NHEL
  • REnd(NHelEnd)      radius at exit of Inner detector
  • PhiPEnd(NHelEnd)   phi
  • ZEnd(NHelEnd)      Z
  • DPhiEnd(NHelEnd)   dphi
  • CotThEnd(NHelEnd)  cottan(theta)
  • PTInvEnd(NHelEnd)  1/PT
  • BarEnd(NHelEnd)    1 if leave ID in the barrel 2 if leave ID in the endcap

    •  

     

    Covariance on  tracks parameters at ID exit  only filled if  TRAC.FITEND=1 and TRAC.FITCVE=2 (diagonal only if  TRAC.FITCVE=1)

  • NHELCVE=NHEL
  • COVEND11(NHELCVE)
  •  COVEND21(NHELCVE)
  •  COVEND22(NHELCVE)
  •  COVEND31(NHELCVE)
  •  COVEND32(NHELCVE)
  •  COVEND33(NHELCVE)
  •  COVEND41(NHELCVE)
  •  COVEND42(NHELCVE)
  •  COVEND43(NHELCVE)
  •  COVEND51(NHELCVE)
  •  COVEND52(NHELCVE)
  •  COVEND53(NHELCVE)
  •  COVEND54(NHELCVE)
  •  COVEND55(NHELCVE)

    •  

    Conversions

    Conversions are reconstructed using module XCONVER, which combines IPATREC, XKALMAN and XHOUREC tracks
     

    Driving data cards

  • CVNT.Basic  1 (0 default)  XCONVER information
  • CVNT.NCONVMAX  (10 default ,<1000) maximum number of converisons saved
  • CVNT.FITCOV 0 (default) no covariance matrix 1 diagonal elements 2 full covariance matrix

    • Variables

  • NCONV   Number of conversions saved
  • CVITRP(NCONV)   Track number (a la ITRKREF) of the positive track
  • CVITRM(NCONV)  Track number (a la ITRKREF) of the negative track
  • CVRTRP(NCONV)  Radius of innermost hit of the positive track
  • CVRTRM(NCONV)  Radius of innermost hit of the negative track
  • CVPTSUM(NCONV) refitted Pt of the conversion
  • CVCHI2(NCONV)  Chi2 of the conversion
  • CVR(NCONV)   Radius of the conversion point
  • CVZ(NCONV)   Z of the conversion point
  • CVPHI(NCONV) Phi of the conversion at the conversion point
  • CVCOT(NCONV) Cotth(theta) of the conversion at the conversion point
  • CVPTINVP(NCONV)  Inverse pt of positive track at conversion point
  • CVPTINVM(NCONV)  Inverse pt of negative track at conversion point

    •  

     

    Conversion covariance matrix, involving parameters R (1)  to PTINVM (6).  Not written out if CVNT.FITCOV=0, only diagonal terms if CVNT.FITCOV=1, complete if CVNT.FITCOV=2

  • NConvCov  == NCONV
  • Cvcov11(NConvCov)
  • Cvcov21(NConvCov)
  • Cvcov22(NConvCov)
  • Cvcov31(NConvCov)
  • Cvcov32(NConvCov)
  • Cvcov33 (NConvCov)
  • Cvcov41(NConvCov)
  • Cvcov42(NConvCov)
  • Cvcov43(NConvCov)
  • Cvcov44(NConvCov)
  • Cvcov51(NConvCov)
  • Cvcov52(NConvCov)
  • Cvcov53(NConvCov)
  • Cvcov54(NConvCov)
  • Cvcov55(NConvCov)
  • Cvcov61(NConvCov)
  • Cvcov62(NConvCov)
  • Cvcov63(NConvCov)
  • Cvcov64(NConvCov)
  • Cvcov65(NConvCov)
  • Cvcov66(NConvCov)

    •  

     
     
     

    Total energy per calorimeter region

     

    Driving data cards

  • CALO.Basic  1 (default)  Calorimeter information

    • Variables

    (these variables are currently filled from DETM bank, not from RECB)
  • Nh_Calo Number of hits in full calo
  • Eh_Calo  Total energy in full calo
  • Nh_EM Number of hits in  em calo
  • Eh_EM  Total energy in em calo
  • Nh_HAD Number of hits in had calo
  • Eh_HAD Total energy in had calo
  • Nh_Dead Number of hits in dead material
  • Eh_Dead Total energy in dead material

    •  

     

    In the following XXX stands for COPS (Presampler), ACCB (Em barrel), ENDE (EM endcap), TILE (Tile barrel calo), HEND (Hadronic endcap), FWDC( Forward calorimeter), CBEA (Barrel cryostat), CEEA (Endcap Cryostat), CRAC( Barrel Crack), CREC(Endcap Crack), CEPS (Endcap Presampler)
     

  • NregXXX  Number of region in calorimeter XXX
  • E_XXX(NregXXX) Energy in region of calorimeter XXX
  • N_XXX(NregXXX) Number of cell in calorimeter XXX

    •  

    Missing energy

    Driving data cards

  • MISS.BASIC 1 (default) Missing energy info
  • MISS.INFOCAL 1 (default) missing Et per calo (COPS, ACCB, ENDE, TILE, HEND, FWDC, see above)
  • MISS.INFOETA 1 (default) missing Et  per eta region  (0<eta<1.5, 1.5<eta<3, eta>3)

    • Variables

    Always there:
  • ETMISS  Missing energy
  • EXMISS X component of missing energy
  • EYMISS Y component of missing energy

    •  
    IF MISS.INFOCAL=1
  • ETMISS _CAL(6) Missing energy in given calo
  • EXMISS _CAL(6) X component of missing energy in given calo
  • EYMISS _CAL(6) Y component of missing energy in given calo

    •  
    IF MISS.INFOETA=1
  • ETMISS _ETA(3) Missing energy in given eta region (0<eta<1.5, 1.5<eta<3, eta>3)
  • EXMISS _ETA(3) X component of missing energy in given eta region
  • EYMISS _ETA(3) Y component of missing energy in given eta region

    • Hadronic jets

    Driving data cards

  • JFNT.Basic 1 (default) if jet information
  • JFNT.Njet   20 (default) Maximum number of jets
  • JFNT.InfoCal 1 (default) if detailed jet calo info
  • JFNT.InfoReg 1 (default) if detailed calo region info
  • JFNT.InfoTau  1 ( default) 1 if tau info to be filled
  • JFNT.InfoKin  1 (default) if KINE jets info
  • JFNT.Neta     100 (default)     define a 2D...
  • JFNT.EtaMin -5 (default)    ...histogram used...
  • JFNT.EtaMax 5 (default)    ...to calculate true...
  • JFNT.Nphi      64 (default)     ...jet energy and direction

    • Variables

    Jets are ordered by decreasing ET.
  • NUMJ number of jet
  • IJETREF(j) : original jet bank number (to be used for reference)
  • EntJ(j)         : ET of  the j'th jet
  • EtaJ(j)         : pseudorapidity of the j'th jet
  • PhiJ(j)         : azimuthal angle of the j'th jet

    •  

     

    From JETS/CISO bank,  only if JFNT.Infocal=1
     

  • EneJ_CAL(11,numj) : Total energy in i'th calorimeter of j'th jet

    •  

     

                   from JETS bank , only if JFNT.InfoReg=1
     

  • EneJ_XXX(i,j)  : Energy in i'th region of Calo XXX  (see Total energy section) of j'th jet

    •  

     

                   finally from Kine information,  only if JFNT.InfoKin=1  (KINE energy summed in cone around reconstructed jet direction . THIS IS NOT THE PARTON ENERGY!!!)
     

  • EntJ_Kin(j)     : Kine. ET of  the j'th jet
  • EtaJ_Kin(j)     : Kine. pseudorapidity of the j'th jet
  • PhiJ_Kin(j)     : Kine. azimuthal angle of the j'th jet
    •  
    Tau identification (only if JFNT.InfoTau=1)
  • NUMD number of jet (=NUMJ)
  • itau(j) 1 if jet is a tau using calo criteria and reco tracks criteria
  • itaucal(j) 1 if jet is a tau using calo criteria only
  • itauk(j) 1 if jet is a tau using calo criteria and kine tracks criteriarecal (j) em radius
  • tra2k(j) number of associated tracks pt>2 geV from KINE bank
  • ntra2r(j) number of associated tracks pt>2 geV reconstructed
  • recal(j) electromagnetic radius
  • econe(j) fraction of transverse energy in cone .1-.2

    • Tile calorimeter cells

    Driving data cards

  • CENT.Basic   1 (0 default) Calorimeter cell info to be filled
  • CENT.NCEL 512 (default, <=4096) Maximum number of cells to be filled
  • CENT.ETTH 0 (default) Threshold in Et applied

    • Variables

     
    From bank  HCAL/TICE
  • NCELT Number of cells
  • ENTT(NCELT) ET in the cell
  • EPDT(NCELT)  condensed info abaout eta, phi, deta: (Sampling Nb + 100 * width in eta+ 10000 phi module number (1->64)+ 1000000 abs(eta) of cell )* sign (eta)

    •  

     

    Useful formulae to decode EPDT:
    * sampling number : 1 to 3 for barrel; 4 to 6 for ext.barrel
       samp = abs(mod(mod(epdt,100000.),10.))
    * eta bin size (10 for deta=0.1, 20 for deta=0.2)
       etasize = mod((abs(epdt)-samp),100.)
    * phi bin number 1 to 64
       phibin = mod((abs(epdt)-etasize-samp),10000.)/100.
     *value of center of eta bin
       etabin = sign(1.,epdt)*(abs(epdt)-phibin*100.-etasize-samp)/1000000.
     

    Calorimeter projected matrix

    From bank HCAL/ETPN. Transverse Energy per bin eta*phi=0.1*0.1 in all calorimeters with E.M weight

    Driving data cards

  • ETNT.Basic   1 (0 default) Projected matrix to be filled
  • ETNT.NCEL 512 (default, <=6400) Maximum number of cells to be filled
  • ETNT.ETTH 0 (default) Threshold in Et applied at ntuple filling time

    • Variables

     
  • NCELPM Number of cells above threshold
  • ENTPM(NCELPM) ET in the cell
  • EPPM(NCELPM)  Phi/eta map encoded: eta=(int(eppm)/10)/100.,phi=mod(abs(eppm),1.)*10
  • NCELPM_KIN Number of cells above threshold (particle level) (fill only if jets info was required)
  • ENTPM_KIN(NCELPM) ET in the cell: sum of pt of particles. magnetic field applied as in ATLFAST
  • EPPM_KIN(NCELPM)  Phi/eta map encoded: eta=(int(eppm_kin)/10)/100., phi=mod(abs(eppm),1.)*10

    •  

    Electromagnetic clusters

    Driving data cards

  • EMCA.Basic   1 (default) E.m cluster info to be filled (eta*phi=3x5 )
  • EMCA.Wind      10 (default)   Switch for additional cluster sizes.Wind=xxxxxx with x=0 or 1 to switch on 7x7, 5x7, 5x5, 3x7, 3x5, 3x3.  (switch for 3x5 is not used because this block is always filled)
  • EMCA.ID        1(default) if fill identification block

    • Variables

     Clusters are ordered by decreasing ET
    From bank CTYP(2)/EMCL (3*5 cluster)
  • NEMCL  Number  of cluster
  • IEMCLREF  Number of the cluster bank (to be used for reference)
  • PHI2     Phi in second sample
  • ETA1   Eta in first sample
  • ETA2  Eta in second sample
  • ZREC  Z vertex
  • EZREC  Error on z vertex
  • DEPTH Shower depth
  • ETA  Eta  pointing
  • E0_35 Energy in presampler with 3*5 cluster
  • E1_35  Energy in first sample with 3*5 cluster
  • E2_35 Energy in second sample with 3*5 cluster
  • E3_35Energy in third sample with 3*5 cluster
  • Es_35Energy in scintillator with 3*5 cluster
    •  
    If required in EMCA.Wind (000001)
  • E0_33 Energy in presampler with 3*3 cluster
  • E1_33  Energy in first sample with 3*3 cluster
  • E2_33 Energy in second sample with 3*3 cluster
  • E3_33Energy in third sample with 3*3 cluster
  • Es_33Energy in scintillator with 3*3 cluster
    •  
    If required in EMCA.Wind (000100)
  • E0_37 Energy in presampler with 3*7 cluster
  • E1_37  Energy in first sample with 3*7 cluster
  • E2_37 Energy in second sample with 3*7 cluster
  • E3_37Energy in third sample with 3*7 cluster
  • Es_37Energy in scintillator with 3*7 cluster
    •  
    If required in EMCA.Wind (001000)
  • E0_55 Energy in presampler with 5*5 cluster
  • E1_55  Energy in first sample with 5*5 cluster
  • E2_55 Energy in second sample with 5*5 cluster
  • E3_55Energy in third sample with 5*5 cluster
  • Es_55Energy in scintillator with 5*5 cluster
    •  
    If required in EMCA.Wind (010000)
  • E0_57 Energy in presampler with 5*7 cluster
  • E1_57  Energy in first sample with 5*7 cluster
  • E2_57 Energy in second sample with 5*7 cluster
  • E3_57Energy in third sample with 5*7 cluster
  • Es_57Energy in scintillator with 5*7 cluster
    •  
    If required in EMCA.Wind (100000)
  • E0_77 Energy in presampler with 7*7 cluster
  • E1_77  Energy in first sample with 7*7 cluster
  • E2_77 Energy in second sample with 7*7 cluster
  • E3_77Energy in third sample with 7*7 cluster
  • Es_77Energy in scintillator with 7*7cluster
    •  
     
    From bank  EMCL/EMID  (Particle Identification) (if emca.ID=1)
  • NEMID  Number of clusters==NEMCL
  • Etha1     ET in first sampling of Had Calo (0.2x0.2)
  • F1        E1/Etot
  • F3        E3/Etot
  • E233      E in 3x3 in sampling 2 - not calibrated for containment
  • E237      E in 3x7 in sampling 2 - not calibrated for containment
  • E277      E in 7x7 in sampling 2 - not calibrated for containment
  • Weta1     sigma(3 strips) in sampling 1 - corrected
  • Weta2     sigma(3x5) in sampling 2 - corrected
  • E2Ts1     ET of 2nd max (3 strips)
  • E2Tsts1   ET on central strip of 2nd max
  • Widths1   sigma(3 strips) - not corrected
  • Barys1    position in strip
  • Wtots1    sigma(40 strips)
  • Emins1    minimum E between 2 maxima
  • Fracs1    (E7-E3)/E3
  • IsEM      flag = 0 if passes EM cuts

    •                  = 1 if |eta| > 2.47
                     = 2 if fails cut on Etha1
                     = 3 if fails cut on sampling 2
                     = 4 if fails cut on sampling 1

    Electron photon identification (from bank EGAM/PHOT and EGAM/ELEC (ON TEST!)

    Driving data cards

  • EGNT.Basic   1 (0 default) fill e/gamma block
  • EGNT.NEGMAX (20 default, <100) maximum number of electronand photon

    • Variables

    Electron (egCH=+/-1) and photons (egCH=0) are put in the same block
  • NEGAM Number  of electron/photon
  • egICLUS(NEGAM) cluster number (a la IEMCLREF)
  • egISEED(NEGAM) seed number for best xkalman (if positive) or pixlrec (if negative) track (if electron) or photon (if conversion)
  • egITRC(NEGAM) track number for best xkalman (if positive) or pixlrec (if negative) track (if electron) or photon (if conversion)
  • egENE(NEGAM)  e/gamma energy
  • egETA(NEGAM) e/gamma eta
  • egPHI(NEGAM) e/gamma phi
  • egZV(NEGAM) Z vertex (estimated from calo)
  • egEZV(NEGAM) error on z vertex estimated from calo
  • egCH(NEGAM) charge (+/-1 for electrons, 0 for photons)

    • Muon in spectrometer

    Driving data cards

  • MUON.Basic   1 (default) Muon spectrometer info to be filled
  • MUON.NTRACKRR     Maximum number of high PT muons (block off if <=0)
  • MUON.NTRFDRR Maximum number of muons fit with straight segments  (block off if <=0)
  • MUON.NSGINSR Maximum Number of Inner station segment (block off if <=0)

    • Variables

    BLOCK MUONFULL(reconstructed tracks at Muon Spectrometer entrance)
==============
NTRACKR[0,100] nber of tracks
CHIR_MU(NTRACKR) chi2
NDOF_MU(NTRACKR) n.o.d.f
IMTH_MU(NTRACKR) obsolete
THPT_MU(NTRACKR) Theta of the crossing point
PHPT_MU(NTRACKR) Phi   of the crossing point
THVE_MU(NTRACKR) Theta of the unit vector
PHVE_MU(NTRACKR) Phi   of the unit vector
PTIR_MU(NTRACKR) Inverse momentum (not transverse momentum)
ER11_MU(NTRACKR) elements of cov.mat.
ER21_MU(NTRACKR)
ER31_MU(NTRACKR)
ER41_MU(NTRACKR)
ER51_MU(NTRACKR)
ER22_MU(NTRACKR)
ER32_MU(NTRACKR)
ER42_MU(NTRACKR)
ER52_MU(NTRACKR)
ER33_MU(NTRACKR)
ER43_MU(NTRACKR)
ER53_MU(NTRACKR)
ER44_MU(NTRACKR)
ER54_MU(NTRACKR)
ER55_MU(NTRACKR)
XCPT_MU(NTRACKR) the x coordinate of the crossing point
YCPT_MU(NTRACKR) the y coordinate of the crossing point
ZCPT_MU(NTRACKR) the z coordinate of the crossing point
XCVE_MU(NTRACKR) the x coordinate of the unit vector
YCVE_MU(NTRACKR) the y coordinate of the unit vector
ZCVE_MU(NTRACKR) the z coordinate of the unit vector
PMON_MU(NTRACKR) momentum
CHAR_MU(NTRACKR) charge
NFI_MU             Total nb. of fits performed
NWI_MU(10,NTRACKR) Statistics of hit origin
MWX_MU(10,NTRACKR) Statistics of hit origin
NWX_MU(10,NTRACKR) Statistics of hit origin
ITI_MU(12,NTRACKR) Codes of stations used in the track
ITX_MU(12,NTRACKR) Codes of stations used in the track
PES_MU(NTRACKR)   Estimated momentum when starting the pattern
RMI_MU(NTRACKR)   minimum distance to the orign (strai. line extr.)
ZMI_MU(NTRACKR)   Z at minimum distance to the orign (strai. line extr.)
CHI_MU(NTRACKR)   chi2
VMU_MU(8,NTRACKR) x, y, z, vx, vy, vz, Preconstructed,charge
CXX_MU(NTRACKR)   chi2
DXX_MU(NTRACKR)   chi2
VXX_MU(8,NTRACKR) x, y, z, vx, vy, vz, Preconstructed, charge

Details on MUONFULL:
  Reconstructed tracks are given at the entrance of MUON system,
 i.e at cylinder of radius 425cm and length 2 times 682cm centered
 on origin.
  One gives the crossing point, i.e the point where
 the track crosses the cylinder, the unit vector tangent to the 
 track and the inverse of the momentum at this point as well as 
 the elements of the covariance matrice on these parameters
    BLOCK MUONTREK(reconstructed tracks at Calorimeter entrance)
==============
NTRAEK[0,100] nber of tracks
THPT_MUK(NTRAEK) Theta of the crossing point
PHPT_MUK(NTRAEK) Phi   of the crossing point
THVE_MUK(NTRAEK) Theta of the unit vector
PHVE_MUK(NTRAEK) Phi   of the unit vector
PTIR_MUK(NTRAEK) Inverse momentum (not transverse momentum)
ER11_MUK(NTRAEK) elements of cov.mat.
ER21_MUK(NTRAEK)
ER31_MUK(NTRAEK)
ER41_MUK(NTRAEK)
ER51_MUK(NTRAEK)
ER22_MUK(NTRAEK)
ER32_MUK(NTRAEK)
ER42_MUK(NTRAEK)
ER52_MUK(NTRAEK)
ER33_MUK(NTRAEK)
ER43_MUK(NTRAEK)
ER53_MUK(NTRAEK)
ER44_MUK(NTRAEK)
ER54_MUK(NTRAEK)
ER55_MUK(NTRAEK)
XCPT_MUK(NTRAEK) the x coordinate of the crossing point
YCPT_MUK(NTRAEK) the y coordinate of the crossing point
ZCPT_MUK(NTRAEK) the z coordinate of the crossing point
XCVE_MUK(NTRAEK) the x coordinate of the unit vector
YCVE_MUK(NTRAEK) the y coordinate of the unit vector
ZCVE_MUK(NTRAEK) the z coordinate of the unit vector
PMON_MUK(NTRAEK) momentum

Details on MUONTREK:
  Reconstructed tracks are given at the entrance of MUON system,
 i.e at cylinder of radius 105.cm and length 2 times 320cm centered
 on origin.
  One gives the crossing point, i.e the point where
 the track crosses the cylinder, the unit vector tangent to the 
 track and the inverse of the momentum at this point as well as 
 the elements of the covariance matrice on these parameters
    BLOCK MUONTRVT(reconstructed tracks at vertex)
==============
NTRAVT[0,100] nber of tracks
A0_MUV(NTRAVT)   Impact parameter of the track at D.C.A.  
Z0_MUV(NTRAVT)   Z   of the track at D.C.A.
PHI_MUV(NTRAVT)  Phi of the track at D.C.A
COTH_MUV(NTRAVT) Cotg of the theta of the track at D.C.A.
PTIR_MUV(NTRAVT) Inverse of transverse momentum
ER11_MUV(NTRAVT) elements of cov.mat.
ER21_MUV(NTRAVT)
ER31_MUV(NTRAVT)
ER41_MUV(NTRAVT)
ER51_MUV(NTRAVT)
ER22_MUV(NTRAVT)
ER32_MUV(NTRAVT)
ER42_MUV(NTRAVT)
ER52_MUV(NTRAVT)
ER33_MUV(NTRAVT)
ER43_MUV(NTRAVT)
ER53_MUV(NTRAVT)
ER44_MUV(NTRAVT)
ER54_MUV(NTRAVT)
ER55_MUV(NTRAVT)
XCPT_MUV(NTRAVT) x  of the track at D.C.A. 
YCPT_MUV(NTRAVT) y  of the track at D.C.A. 
ZCPT_MUV(NTRAVT) z  of the track at D.C.A. 
XCVE_MUV(NTRAVT) ux of the track at D.C.A. 
YCVE_MUV(NTRAVT) uy of the track at D.C.A.
ZCVE_MUV(NTRAVT) uz of the track at D.C.A.
PMON_MUV(NTRAVT) momentum

Details on MUONTRVT:
  Reconstructed tracks are given at the the D.C.A. point
    BLOCK MUONSEGF (results of the fit using straight track segments)
==============
NTRFDR[0,100] nber of tracks
CHIRF_MU(NTRFDR) chi2
THPTF_MU(NTRFDR) Theta of the crossing point
PHPTF_MU(NTRFDR) Phi of the crossing point
THVEF_MU(NTRFDR) Theta of unit vector along the track at crossing point
PHVEF_MU(NTRFDR) Phi of unit vector along the track at crossing point
PTIRF_MU(NTRFDR) Inverse momentum (not transverse momentum)
XCPTF_MU(NTRFDR) the x coordinate of the crossing point
YCPTF_MU(NTRFDR) the y coordinate of the crossing point
ZCPTF_MU(NTRFDR) the z coordinate of the crossing point
XCVEF_MU(NTRFDR) the x coordinate of the unit vector
YCVEF_MU(NTRFDR) the y coordinate of the unit vector
ZCVEF_MU(NTRFDR) the x coordinate of the unit vector
PMONF_MU(NTRFDR) momentum
CHARF_MU(NTRFDR) charge

Details on MUONSEGF:
  Reconstructed tracks are given at the entrance of MUON system,
 i.e at cylinder of radius 425cm and length 2 times 682cm centered
 on origin.
  One gives the crossing point, i.e the point where the track crosses
 the cylinder, the unit vector tangent to the track
 and the inverse of
 the momentum at this point as well as the elements
 of the covariance
 matrice on these parameters
    BLOCK MUONINST (track segments in Inner station of Muon system)
==============
NSGINS[0,100] nber of track segments
MLUI_MU(NSGINS) nber of multilayers used to build the segment
XPTI_MU(NSGINS) the x coordinate of the point defining the segment
YPTI_MU(NSGINS) the y coordinate of the point defining the segment
ZPTI_MU(NSGINS) the z coordinate of the point defining the segment
XVEI_MU(NSGINS) the x coordinate of the vector defining the segment
YVEI_MU(NSGINS) the y coordinate of the vector defining the segment
ZVEI_MU(NSGINS) the z coordinate of the vector defining the segment
QFAI_MU(NSGINS) quality factor

Details on MUONINST:
  The track segments in inner stations of Muon system:
 are defined by a point and a vector 
 are qualified by a quality factor and the number of multilayer used
 The precision of the above informations of the track segments
 is very poor in orthoradial direction

    Muon Identification (MUID)

    Datacards and variable list is available here.
     

    ATLFAST ntuple blocks

    The standard ATLFAST ntuple is provided (current dev version). Blocks can be switched on and off using datacards. ATLFAST driving parameters are provided through the usual  atlfast.dat file, which should be present in the working directory. This allows an easy comparison between fast and full simulation. Specific ATLFAST documentation is here.
     
     

    Driving data cards

  • ATLF.Basic   0 ( default) ATLFAST ntuple blocks  are filled
  • ATLF.INFO  1 (default) general information
  • ATLF.MISS 1 (default) missing energy information
  • ATLF.NLEPM 12 (default) maximum number of isolated leptons(if <=0 no block)
  • ATLF.NPHOM 12 (default) maximum number of isolated photons (if <=0 no block)
  • ATLF.NJETAM 20 (default) maximum number of  jet (if <=0 no block)
  • ATLF.NONMUXM 12 (default) maximum number of  non-isolated muon (if <=0 no block)
  • ATLF.NTRAM 300 (default) maximum number of  tracks (if <=0 no block)
  • ATLF.NPARTM 10 (default) maximum number of particle in history block (if <=0 no block)
  • ATLF.NBPHYSM 40 (default) maximum number of B0,/\b,D0,Ds,J/Y,K0s,/\,e,mu (if <=0 no block)
  • ATLF.TRIG 1 (default) trigger information

    • Variables

    Block ATLFINFO
  • ISUB Process simulated
  • NEL number of isolated electron
  • NMU number of isolated muon
  • NMUX number of non isolated muons
  • NPH number of isolated photons
  • JETB number of b jet
  • JETC number of c jet
  • JETL number of light jet
  • CIRCJ jets circularity
  • CIRCE event circularity
  • THRUST event thrust
  • OBLAT event oblateness
    •  
    Block PMISSING
  • PMISS(2) missing px, py
    •  
  • Block PLEPTONS
  • NLEP number of leptons
  • KFLEP(NLEP) lepton PDG number
  • KTRGLEP(NLEP) 1 if trigger
  • PXLEP(NLEP) lepton 4 momentum
  • PYLEP(NLEP)
  • PZLEP(NLEP)
  • EELEP(NLEP)
    •  
    Block PPHOTONS
  • NPHO number of isolated photons
  • KFPHO(NPHO) photon PDG number
  • PXPHO(NPHO) photon 4 momentum
  • PYPHO(NPHO)
  • PZPHO(NPHO)
  • EEPHO(NPHO)
    •  
    Block PPJETS
  • NJETA number of jets
  • KFJET(NJETA)  jet PDG number
  • PXJET(NJETA) jet 4 momentum
  • PYJET(NJETA)
  • PZJET(NJETA)
  • EEJET(NJETA)
    •  
    Block PMUXS
  • NMUXS number of non-isolated muons
  • KFMUXS(NMUXS) muon PDG number
  • KTRGMUX(NMUXS) 0 no trigger 1 low threshold 2 high threshold
  • PXMUX(NMUXS) muon 4 momentum
  • PYMUX(NMUXS)
  • PZMUX(NMUXS)
  • EEMUX(NMUXS)
    •  
    Block PTRACKS
  • NTRA number of tracks
  • KPTRA(NTRA) GENZ track number
  • KFTRA(NTRA) PDG track number
  • KPM1TRA(NTRA) track 1st ancestorGENZ number
  • KFM1TRA(NTRA) track 1st ancestor PDG number
  • KPM2TRA(NTRA) track 2nd ancestorGENZ number
  • KFM2TRA(NTRA) track 2ndancestor PDG number
  • KPM3TRA(NTRA) track 3rd ancestorGENZ number
  • KFM3TRA(NTRA) track 3rd ancestor PDG number
  • KPM4TRA(NTRA) track 4th ancestorGENZ number
  • KFM4TRA(NTRA) track 4th ancestor PDG number
  • KPM5TRA(NTRA) track 5th ancestorGENZ number
  • KFM5TRA(NTRA) track 5th ancestor PDG number
  • KPM6TRA(NTRA) track 6th ancestorGENZ number
  • KFM6TRA(NTRA) track 6th ancestor PDG number
  • D0TRACRU(NTRA) true D0
  • Z0TRACRU(NTRA) true Z0
  • PHITRACRU(NTRA) true phi
  • COTTRACRU(NTRA) true cottan theta
  • PTINVTRACRU(NTRA) true inverse pT
  • D0TRAC(NTRA) simulated D0
  • Z0TRAC(NTRA) simulated Z0
  • PHITRAC(NTRA) simulated phi
  • COTTRAC(NTRA) simulatedcottan theta
  • PTINVTRAC(NTRA) simulated inverse pT
  • CORR11(NTRA)  simulated covariance matrix
  • CORR21(NTRA)  simulated covariance matrix
  • CORR31(NTRA)  simulated covariance matrix
  • CORR41(NTRA)  simulated covariance matrix
  • CORR51(NTRA)  simulated covariance matrix
  • CORR22(NTRA)  simulated covariance matrix
  • CORR32(NTRA)  simulated covariance matrix
  • CORR42(NTRA)  simulated covariance matrix
  • CORR52(NTRA)  simulated covariance matrix
  • CORR33(NTRA)  simulated covariance matrix
  • CORR43(NTRA)  simulated covariance matrix
  • CORR53(NTRA)  simulated covariance matrix
  • CORR44(NTRA)  simulated covariance matrix
  • CORR54(NTRA)  simulated covariance matrix
  • CORR55(NTRA)  simulated covariance matrix
  • EFFTRA(NTRA) track efficiency
  • ISTATRA(NTRA) track status (1 if correctly reconstructed)
    •  
    Block PHISTORY : particles from the original hard process (kept for historical reasons, please use KINEGENZ block in CBNT)
  • NPART number of particle
  • KFPAR(NPART) particle PDG number
  • PXPAR(NPART) particle 4 momentum
  • PYPAR(NPART)
  • PZPAR(NPART)
  • EEPAR(NPART)
    •  
    Block BPHYSICS (kept for historical reasons, please use KINEGENZ block in CBNT with BPHYS=1 option)
  • NBPHYS number of B0,/\b,D0,Ds,K0s,/\,J/Y,e,mu in GENZ bank
  • KPBPHYS(NBPHYS) GENZ number
  • KFBPHYS(NBPHYS) PDG number
  • PXBPHYS(NBPHYS) 4-momentum
  • PYBPHYS(NBPHYS) 4-momentum
  • PZBPHYS(NBPHYS) 4-momentum
  • EEBPHYS(NBPHYS) 4-momentum
  • VXBPHYS(NBPHYS)  origin vertex
  • VYBPHYS(NBPHYS)  origin vertex
  • VZBPHYS(NBPHYS)  origin vertex
    •  
    Block PTRIGGER
  • TGALL 1  if any trigger
  • TGEM1 1 if one isolated electron
  • TGPH1 1 if one isolated photon
  • TGEM2 1 if two isolated electron/photon
  • TGMU1 1 if one isolated muon
  • TGMU2 1 if two isolated muons
  • TGEMU 1 if electron/muon pair
  • TGJT1 1 if one jet
  • TGJT2 1 if two jets
  • TGJT4 1 if four jets

    • User block

    A simple user routine (cbntuser) is provided just as an example. It is up to the user to extend this routine (see instructions)  to provide whatever additional information he needs (specific GENZ information for example). More driving datacards can also be defined.

    Driving data cards

  • USER.Basic   1 (0 default) user block to be filled (should be left)
  • USER.NOBJECTMAX      just an example

    • Variables

  • NOBJECT number of object (just an example)
  • UOBJECT(NOBJECT) properties of object (just an example)

    • Comments/suggestions/assistance on this page or the CBNT itself  David.Rousseau@cern.ch