(Patrice Verdier) |
The code SQUARKANA is used for the squark analysis
(stop and sbottom) : stop -> c neutralino1
and stop -> c gluino.
The code is written in the SKELANA2.0 framework.
REMCLU is not used for the moment, but the code is ready
to use it.
The jet algorithm used is DURHAM in the YKERN package : The events are
forced in 2 and 4 jets.
SKELANA flags :
IFLTRA = 1 ! Track information COMMON Filled (1 value = default)
IFLODR = 1 ! Ordered track first charged then neutral (1= default)
IFLVEC = 2 ! Fill all tracks except the charged outgoing ones and PA's related to REMCLU clusters
IFLSTR = 1 ! Selection applied and bad tracks removed form VECSUB
IFLCUT = 3 ! Recommended Track selection tuning
IFLRVR = 000! Recover charged tracks as neutral.
IFLSIM = 1 ! Simulation info filled
IFLBSP = 2 ! Beamspot info red from beamspot file
IFLBTG = 0 ! Btag filled from XshortDST
IFLPVT = 1 ! Btagging primary vertex
IFLVDR = 1 ! VD only refitting
IFLFCT = 1 ! FCA/FCB refitting with PV constraint
IFLRNQ = 1 ! Run quality selection
IFLBHP = 1 ! Skip the BAD 1997 HPC
IFLUTE = 1 ! Unassociated TE bank filled
IFLVDH = 1 ! Vertex detector hits filled
IFLMUO = 1 ! Muon id filled from Xshort DST
IFLECL = 0 ! No remclu
IFLELE = 1 ! Elephant elec id
IFLEMC = 1 ! Electromagnetic calorimetry
IFLPHO = 1 ! Photon id
IFLPHC = 1 ! Photon Conversion
IFLSTC = 1 ! Stic info
IFLHAC = 1 ! Hadron calo
IFLHAD = 1 ! Hadron id
IFLRV0 = 1 ! V0 reconstruction
IFLJET = 0 ! No jet reconstruction
Preselection :
(1) at least 6 charged particles (2) transverse energy >= 4 GeV (3) at least 1 particle with Pt >= 1.5 GeV/c
Ntuple description :
****************************************************************** * Ntuple ID = 100 Entries = 113 Squark analysis ****************************************************************** * Var numb * Type * Block * Name * ****************************************************************** * 1 * I*4 * EVENT * NORUN * # run * 2 * I*4 * EVENT * NOEVE * # event * 3 * I*4 * EVENT * NOFIL * # fill * 4 * R*4 * EVENT * ECENT * # sqrt(s) * 5 * I*4 * EVENT * MSQU * Generated squark Mass * 6 * I*4 * EVENT * MLSP * Generated LSP Mass ****************************************************************** * 1 * R*4 * SNNSO * SBOGDM(3) * NN output sbot. Dm > 30 * 2 * R*4 * SNNSO * SBOMDM(3) * NN output sbot. Dm < 30 * 3 * R*4 * SNNSO * STOGDM(3) * NN output stop. Dm > 30 * 4 * R*4 * SNNSO * STOMDM(3) * NN output stop. Dm < 30 * 1 * I*4 * ENERG * NCHAR * # charged * 2 * I*4 * ENERG * NNEUT * # neutral * 3 * I*4 * ENERG * NSCHA * sum(charge) * 4 * I*4 * ENERG * NTOTA * # particles * 5 * R*4 * ENERG * EVISI * Evis. * 6 * R*4 * ENERG * ECHAP * Echar. * 7 * R*4 * ENERG * DCHAR * error(Echar) * 8 * R*4 * ENERG * ENEUT * Eneut. * 9 * R*4 * ENERG * DNEUT * error(Eneut) * 10 * R*4 * ENERG * ETSCA * E transverse * 11 * R*4 * ENERG * ETCHA * E transverse charged * 12 * R*4 * ENERG * PTMIS * Pt missing * 13 * R*4 * ENERG * FMVIS * Mvis. * 14 * R*4 * ENERG * FMINV * Minv. * 15 * R*4 * ENERG * EELEC * E(electron+gamma) * 16 * R*4 * ENERG * ECO15 * E cone 15 around z axis * 17 * R*4 * ENERG * ECO20 * E cone 20 around z axis * 18 * R*4 * ENERG * ECO30 * E cone 30 around z axis * 19 * R*4 * ENERG * ECO45 * E cone 45 around z axis * 20 * R*4 * ENERG * ENE15 * Eneut. cone 15 around z axis * 21 * R*4 * ENERG * ENE20 * Eneut. cone 20 around z axis * 22 * R*4 * ENERG * ENE30 * Eneut. cone 30 around z axis * 23 * R*4 * ENERG * ENE45 * Eneut. cone 45 around z axis * 24 * R*4 * ENERG * PTMAX * Max Pt of a part. * 25 * R*4 * ENERG * PTMCH * Max Pt of a charged part. * 26 * R*4 * ENERG * PTMNE * Max Pt of a neutral part. * 27 * R*4 * ENERG * PMISS(6) * Missing P info * 28 * R*4 * ENERG * SPRIM * sqrt(s') * 29 * R*4 * ENERG * SPR3C * sqrt(s') 3C FIT ****************************************************************** * 1 * I*4 * DETEC * NVD * # part. VD only * 2 * I*4 * DETEC * NVI * # part. VD-ID only * 3 * I*4 * DETEC * NVIA * # part. VD-ID only with Ecalo>0 * 4 * I*4 * DETEC * NVIT * # part. VD-ID-TPC only * 5 * I*4 * DETEC * NVITO * # part. VD-ID-TPC-OD only * 6 * R*4 * DETEC * EVD * E of VD only * 7 * R*4 * DETEC * EVI * E of VD-ID only * 8 * R*4 * DETEC * DVI * error(E of VD-ID only) * 9 * R*4 * DETEC * EVIA * E of VD-ID only with Ecalo>0 * 10 * R*4 * DETEC * DVIA * error(E of VD-ID only with Ecalo>0) * 11 * R*4 * DETEC * EVIT * E of VD-ID-TPC only * 12 * R*4 * DETEC * DVIT * error(E of VD-ID-TPC only) * 13 * R*4 * DETEC * EVITO * E of VD-ID-TPC-OD only * 14 * R*4 * DETEC * DVITO * error(E of VD-ID-TPC-OD only) ****************************************************************** * 1 * I*4 * CALOM * MHPC * # of part. Ehpc>0. * 2 * I*4 * CALOM * NNHPC * # of neut. part. Ehpc>0. * 3 * I*4 * CALOM * MEMF * # of part. Eemf>0. * 4 * I*4 * CALOM * NNEMF * # of neut part. Eemf>0. * 5 * I*4 * CALOM * MHAC * # of part. Ehac>0. * 6 * I*4 * CALOM * NNHAC * # of neut part. Ehac>0. * 7 * I*4 * CALOM * MSTC * # of part. Estic>0. * 8 * I*4 * CALOM * NNSTC * # of neut. part. Estic>0. * 9 * R*4 * CALOM * EHPC * E in HPC * 10 * R*4 * CALOM * EEMF * E in EMF * 11 * R*4 * CALOM * EHAC * E in HAC * 12 * R*4 * CALOM * ESTC * E in STIC * 13 * R*4 * CALOM * ENHPC * neut. E in HPC * 14 * R*4 * CALOM * ENEMF * neut. E in EMF * 15 * R*4 * CALOM * ENHAC * neut. E in HAC * 16 * R*4 * CALOM * ENSTC * neut. E in STIC ****************************************************************** * 1 * R*4 * BTAGG * PROBAP * btag pos. * 2 * R*4 * BTAGG * PROBAN * btag neg. * 3 * R*4 * BTAGG * PROBAS * btag tot. * 4 * R*4 * BTAGG * PROBZ0 * combined Z0 * 5 * R*4 * BTAGG * PROBHZ * combined HZ * 6 * R*4 * BTAGG * PROBHA * combined HA ****************************************************************** * 1 * R*4 * JETII * Y2 * Ycut 1-2 jet * 2 * R*4 * JETII * Y3 * Ycut 2-3 jet * 3 * R*4 * JETII * P2J(2,6) * 2 Jets 4-vecotr * 4 * R*4 * JETII * PLEAD2(2,8)* leading part. of jets * 5 * R*4 * JETII * ETRA2(2) * E transverse * 6 * R*4 * JETII * ECHA2(2) * E charged * 7 * R*4 * JETII * ETCH2(2) * E transverse charged * 8 * R*4 * JETII * EELE2(2) * E electromagnetic * 9 * R*4 * JETII * EHAD2(2) * E hadronic * 10 * R*4 * JETII * T2J12 * theta jet 1-2 * 11 * I*4 * JETII * NPAR2(2) * # part. * 12 * I*4 * JETII * NCHA2(2) * # char. part. * 13 * R*4 * JETII * PRZ02 * batg Comb. Z0 with jet * 14 * R*4 * JETII * PRHZ2 * btag Comb. HZ with jet * 15 * R*4 * JETII * PRHA2 * btag Comb. HA with jet * 16 * R*4 * JETII * PRJP2(2) * batg. pos ****************************************************************** * 1 * R*4 * JETIV * Y4 * * 2 * R*4 * JETIV * Y5 * * 3 * R*4 * JETIV * P4J(4,6) * * 4 * R*4 * JETIV * PLEAD4(4,8)* * 5 * R*4 * JETIV * ETRA4(4) * * 6 * R*4 * JETIV * ECHA4(4) * * 7 * R*4 * JETIV * ETCH4(4) * * 8 * R*4 * JETIV * EELE4(4) * * 9 * R*4 * JETIV * EHAD4(4) * * 10 * R*4 * JETIV * T4J12 * * 11 * R*4 * JETIV * T4J13 * idem 4 jets * 12 * R*4 * JETIV * T4J14 * * 13 * R*4 * JETIV * T4J23 * * 14 * R*4 * JETIV * T4J24 * * 15 * R*4 * JETIV * T4J34 * * 16 * I*4 * JETIV * NPAR4(4) * * 17 * I*4 * JETIV * NCHA4(4) * * 18 * R*4 * JETIV * PRZ04 * * 19 * R*4 * JETIV * PRHZ4 * * 20 * R*4 * JETIV * PRHA4 * * 21 * R*4 * JETIV * PRJP4(4) * ****************************************************************** * 1 * R*4 * TOPOL * COPLA * theta jet 1-2 in trans. plane * 2 * R*4 * TOPOL * ACOPL * acoplanarity (PUTHRU) * 3 * R*4 * TOPOL * SPHER * sphericity * 4 * R*4 * TOPOL * APLAN * aplanarity * 5 * R*4 * TOPOL * OBLAT * oblatness * 6 * R*4 * TOPOL * ACOLI * acolinearity * 7 * R*4 * TOPOL * THRUS(4) * thrust coordinate * 8 * R*4 * TOPOL * H10 * FOX WOLFRAM 1 * 9 * R*4 * TOPOL * H20 * " 2 * 10 * R*4 * TOPOL * H30 * " 3 * 11 * R*4 * TOPOL * H40 * " 4 ****************************************************************** * 1 * I*4 * PARTI * NELEC * # of electrons * 2 * I*4 * PARTI * NMUON * # of muons * 3 * R*4 * PARTI * PELE(9,NELEC)* electrons infos * 4 * R*4 * PARTI * PMUO(9,NMUON)* muons infoss * 5 * I*4 * PARTI * NGAT * # of gamma * 6 * I*4 * PARTI * NGAM * # of gamma <=2 * 7 * R*4 * PARTI * PGAM(9,NGAM)* info of the 2 most E. gamma * 8 * R*4 * PARTI * PCH1(6) * most E. ch. part. 4-vector * 9 * R*4 * PARTI * PCH2(6) * 2nd most E. ch. part. 4-vector * 10 * R*4 * PARTI * PNE1(6) * most E. ne. part. 4-vector * 11 * R*4 * PARTI * PNE2(6) * 2nd most E. ne. part. 4-vector ****************************************************************** * 1 * R*4 * HERMI * A40JJ * * 2 * R*4 * HERMI * A40JP * * 3 * R*4 * HERMI * A40PP * * 4 * R*4 * HERMI * A40PJ * * 5 * R*4 * HERMI * A90JJ * * 6 * R*4 * HERMI * A90JP * * 7 * R*4 * HERMI * A90PP * TAGGERS INFO * 8 * R*4 * HERMI * A90PJ * * 9 * R*4 * HERMI * APHJJ * * 10 * R*4 * HERMI * APHJP * * 11 * R*4 * HERMI * APHPP * * 12 * R*4 * HERMI * APHPJ * * 13 * I*4 * HERMI * NT40 * * 14 * I*4 * HERMI * NT90 * * 15 * I*4 * HERMI * NTPH * ****************************************************************** * 1 * I*4 * GAMMA * NTPC * # of TPC tracks * 2 * I*4 * GAMMA * NTPC4 * # of TPC tracks 4 pads * 3 * I*4 * GAMMA * NTPC10 * # of TPC tracks 10 pads * 4 * I*4 * GAMMA * NTPC16 * # of TPC tracks 16 pads * 5 * R*4 * GAMMA * PTPC * P of TPC tracks * 6 * R*4 * GAMMA * PTPC4 * P of TPC tracks 4 pads * 7 * R*4 * GAMMA * PTPC10 * P of TPC tracks 10 pads * 8 * R*4 * GAMMA * PTPC16 * P of TPC tracks 16 pads * 9 * R*4 * GAMMA * PVDID * P of VD-ID only * 10 * R*4 * GAMMA * EHEM1 * E hemisphere 1 * 11 * R*4 * GAMMA * EHADR1 * E hadronic hem. 1 * 12 * R*4 * GAMMA * EELEC1 * E electromagnetic hem. 1 * 13 * R*4 * GAMMA * EHEM2 * E hemisphere 2 * 14 * R*4 * GAMMA * EHADR2 * E hadronic hem. 2 * 15 * R*4 * GAMMA * EELEC2 * E electromagnetic hem. 2 * 16 * R*4 * GAMMA * EFOR10 * E cone 10 around z axis * 17 * R*4 * GAMMA * PCHT * P charged particles * 18 * I*4 * GAMMA * ISTAG * dafne tag fo single tagged gg * 19 * I*4 * GAMMA * IUTAG * dafne tag fo untagged gg ****************************************************************** * 1 * I*4 * HEAVY * NHEAV * # of heavy tracks * 2 * R*4 * HEAVY * PHEAV(11,NHEAV)* info of heavy tracks ******************************************************************