Tags:
create new tag
view all tags
Introduction

The "mn" version of ntuple is designed for keeping as a general ntuple for AMS analysis, therefore it is keeping as much information as possible. The anverage size for a typical run would be about 300MiB, and total size for 5 year data would be 40TiB.

In practise we only select a subset of the events in ISS data. To be clear, the file name of the output also shows the basic selection used. The naming convection for makeNtuple is: mn-{datatype}-{dataid}.eminiXX.root, in which XX is 2 hex digits combined from the seleection bits. The bits on selection is: (from lower bit to higher)

  • bit 0(0x1): fullspan only
  • bit 1(0x2): above 30GeV
  • bit 2(0x4): positive prescale
  • bit 3(0x8): shower event only
  • bit 4(0x10): magical
  • bit 5(0x20): nuclei (q_inner>1.5)

Therefore, emini08 would keep only events with shower, emini04 would apply prescale to positive events, and emini03 would be events with fullspan track and above 30GV.

Beside that, events without any shower nor any reconstructed particle are also skipped.

In the ntuple, basicly only 1 reconstructed object in each detector is kept. When multiple showers exist, the one with highest energy is selected. Then tracker track, TRD track and TOF beta objects are selected by matching them with the leading shower. If there is no shower in the event, the objects assosiated with particle #0 are selected.

Now The data could be found in directory

/afs/cern.ch/work/t/tangzhch/ams/linkmn/ISS.B950/pass6/

which containts links to files on NAS.

For electron/positron analysis, a preselected vertion of mn ntuple was produced and saved on NAS. See PositronAnalysis for details.

Variable desciption

Here list the short description for variables in the current common ntuple, produced with program "makeNtuple"

General

  • Int_t info_evvalid
    • There are valid values in the event info values if valid==0.
  • UInt_t info_run
    • The run number
  • UInt_t info_event
    • The event eumber
  • UInt_t info_utime
    • The unix time recorded in JMDC.
  • Float_t info_lf
    • The livetime of the event.
  • Float_t info_maxcf[4][2]
    • The maximum cutoff calculated with Stormer model. The first index for 25, 30, 35, 40 degree FOV, and then index of 0, 1 for negative and positive events.
  • Float_t info_maxcfIGRF[4][2]
    • The maximum cutoff calculated with IGRF model. The index follow the same rule with info_maxcf.
  • Float_t info_cfIGRF_ev_cfi
    • The
  • Float_t info_cfIGRF_ev_cfineg
    • The
  • Float_t info_cfStormer
    • The Stormer cutoff of the given particle direction.
  • Float_t info_cf
    • The other Stormer cutoff stored in the event.
  • Int_t info_btres
    • The back tracing result of IGRF model: the higher bits(info_btres>>8) is back tracing status (1: under cutoff, 2: above cutoff, 3: trapped, other: error); the lower 8bits(info_bres&256) is the back tracking result showing which position and time infomation are used in the calculation. For detail see AMSEventR::DoBacktracing();
  • Float_t info_theta
    • The
  • Float_t info_phi
    • The
  • Float_t info_glong
    • The
  • Float_t info_glat
    • The
  • Float_t info_zenith
    • The
  • Float_t info_ThetaS
    • The
  • Float_t info_PhiS
    • The
  • Float_t info_ThetaM
    • The
  • Float_t info_PhiM
    • The
  • Bool_t info_inTrd
    • The
  • Bool_t info_inSAR
    • Whether the particle pass through solar array on space station.
  • Int_t info_nAnti
    • The
  • Int_t info_ipart
    • The id of selected particle object.
  • Float_t mcinfo_p
    • The generated momentum of MC particle. For data it's 0.
  • Float_t mcinfo_theta
    • The theta of MC generated particle.
  • Float_t mcinfo_phi
    • The phi of MC generated particle.
  • Float_t mcinfo_x
    • The
  • Float_t mcinfo_y
    • The
  • Int_t mcinfo_inUpperTOF
    • The
  • Int_t mcinfo_inLowerTOF
    • The
  • Int_t mcinfo_inEcalExit
    • The
  • Int_t mcinfo_inEcalEntry
    • The
  • Int_t mcinfo_isPassTOF[4]
    • The
  • Float_t mcinfo_ecalentry[2]
    • The
  • Float_t mcinfo_ecalexit[2]
    • The
  • Float_t mcmisalign[2][2]
    • The simulated external tracker misalignment on layer1(0) and layer9(1): the x(0) and y(1) value.
  • Int_t lvl1_valid
    • The level1 trigger variables is valid (==0).
  • Int_t lvl1_JMembPatt
  • Int_t lvl1_PhysBPatt
  • Int_t lvl1_JMembPattMC
    • The
  • Int_t lvl1_PhysBPattMC
    • The
  • Int_t lvl1_rret
    • The
  • Float_t lvl1_ftrate
    • The

ECAL

  • Int_t ecal_ish
    • The selected shower id.
  • Int_t ecal_nsh
    • The number of showers in the event.
  • Float_t ecal_theta
    • The
  • Float_t ecal_phi
    • The
  • Int_t ecal_nhits
    • The
  • Int_t ecal_stat
    • The shower status bit 25, or -1 if no shower.
  • Float_t ecal_ene
    • The reconstructed energy <nop>EnergyE.
  • Float_t ecal_ene14
    • The new version of <nop>EnergyE.
  • Float_t ecal_edep
    • The sum of hit energy in the shower, in <nop>GeV
  • Float_t ecal_enc
    • The reconstructed energy <nop>EnergyC.
  • Float_t ecal_enp
    • The reconstructed energy <nop>EnergyP, which is the corrected <nop>EnergyC.
  • Float_t ecal_cog[3]
    • The simple center of gravity for the shower.
  • Float_t ecal_entry[3]
    • The entry position of the shower, Z at the center of the first ECAL layer.
  • Float_t ecal_entry_old[3]
    • The old version based on simple COG.
  • Float_t ecal_exit[3]
    • The exit position of the shower, Z at the center of the last ECAL layer.
  • Float_t ecal_ehit[18][72]
    • The deposited energy for each cell.
  • Float_t ecal_adcl[18][72]
    • The ADC value from low gain channel for each cell.
  • Float_t ecal_adch[18][72]
    • The ADC value from high gain channel for each cell.
  • Float_t ecal_orp
    • The orphan energy in ECAL.
  • Float_t ecal_s3
    • The product of Energy3C [0] [1] and [2] (3cm, 5cm, 7cm).
  • Float_t ecal_e3c0
    • The energy deposition within 3cm of shower axis (Energy3C [0])
  • Float_t ecal_e3c1
    • The energy deposition within 5cm of shower axis (Energy3C [1])
  • Float_t ecal_rl
    • The
  • Float_t ecal_sl
    • The
  • Float_t ecal_e1st
    • The
  • Float_t ecal_bdtp
    • The most recent Pisa BDT value (currently version 7)
  • Float_t ecal_bdtp5
    • The version 5 Pisa BDT value.
  • Float_t ecal_bdtchi2
    • The
  • Float_t ecal_bdtp5s
    • The
  • Float_t ecal_bdtchi2s
    • The
  • Float_t ecal_bdtl
    • The
  • Float_t ecal_bdtl3
    • The
  • Float_t ecal_bdtl3eff
    • The
  • Float_t ecal_nchi2
    • The

Tracker

  • Int_t trk_itrk
    • The
  • Int_t trk_ntrk
    • The number of tracker track in the event.
  • Int_t trk_nhitx
    • The number of 2D(xy)-hits in the track.
  • Int_t trk_nhity
    • The number of hits in the track.
  • Int_t trk_pat
    • The tracker hits pattern in the track. Bit 0~8: layer1 to layer9
  • Int_t trk_patx
    • The 2D hit pattern in the track. Bit 0~8: layer1 to layer9
  • Float_t trk_theta
    • The
  • Float_t trk_phi
    • The
  • Float_t trk_bcorr
    • The
  • Float_t trk_b_corr[6]
    • The
  • Int_t trk_b_corrflag[6]
    • The
  • Float_t trk_rig[6]
    • The reconstructed regidity for different tracker hit pattern: 0: inner, 1: maxspan, 2:inner+L1, 3:inner+L9, 4:inner upper half, 5: inner lower half. The rigidities are obtained from fitid returned from TrTrackR::iTrTrackPar() with fit pattern 3, 0, 5, 6, 1, 2.
  • Float_t trk_rigerr[6]
    • The error of rigidity fitting. Index convention is the same as trk_rig[]
  • Float_t trk_rig_al[6]
    • The rigidity calculated with Alcaraz method.
  • Float_t trk_chi2x[6]
    • The tracker fitting chi2 on x projection
  • Float_t trk_chi2y[6]
    • The tracker fitting chi2 on y projection
  • Float_t trk_chi2y_op[6]
    • The chi2 given opposite rigidity sign in fitting.
  • Float_t trk_res1[6][2]
    • The tracker fitting residual on layer 1 with different fitting pattern, on x and y.
  • Float_t trk_res9[6][2]
    • The tracker fitting residual on layer 1 with different fitting pattern, on x and y.
  • Float_t trk_qin
    • The tracker charge calculated from inner tracker.
  • Float_t trk_qinrms
    • The RMS of inner tracker charge.
  • Float_t trk_qinup
    • The tracker charge using layer 2 to layer 4
  • Float_t trk_qinlow
    • The tracker charge using layer 5 to layer 8
  • Float_t trk_qlay[9]
    • The tracker charge on each layer.
  • Float_t trk_ecalentry[3]
    • The coordinate of tracker extrapolate to ECAL entry (Z=-143.2cm)
  • Float_t trk_ecalexit[3]
    • The coordinate of tracker extrapolate to ECAL exit (Z=-158.92cm)
  • Float_t trk_L1theta
    • The
  • Float_t trk_L1phi
    • The
  • Float_t trk_L9theta
    • The
  • Float_t trk_L9phi
    • The
  • Float_t trk_mindisty
    • The
  • Int_t trk_nnoise
    • The
  • Float_t trk_laypos[9][2]
    • The
  • Float_t trk_aligndiff[2][3]
    • The (external) alignment difference on layer 1 and 9, with x, y, z value in centimeter.
  • Float_t trk_aligndiffRTI[2][3]
    • The same as above, but this is the average value in one second (from RTI)
  • Float_t trk_edeplay_offtrack[9]
    • The
  • Int_t trk_nhitlay_offtrack[9]
    • The
  • Float_t trk_edeplay_neighbor[9]
    • The
  • Int_t trk_nhitlay_neighbor[9]
    • The
  • Float_t trk_mindisty_lay[9]
    • The
  • Float_t trk_mindisty_sort[9]
    • The
  • Int_t trk_nsumhit_radius[5]
    • The
  • Float_t trk_esumhit_radius[5]
    • The
  • Float_t trk_lkhd
    • The
  • Float_t trk_ccbdt
    • The
  • Float_t trk_layresx[9][2]
    • The
  • Float_t trk_layresy[9][2]
    • The

TOF

  • Int_t tof_ibeta
    • The
  • Int_t tof_ibetah
    • The
  • Float_t tof_beta
    • The
  • Float_t tof_q
    • The
  • Float_t tof_qlow
    • The
  • Float_t tof_qup
    • The
  • Int_t tof_match
    • The
  • Int_t tof_patth
  • Int_t tof_nlay
    • The number of layers used in tof_qh reconstruction.
  • Int_t tof_nclus
    • The
  • Int_t tof_nclusH
    • The
  • Float_t tof_betah
    • The
  • Float_t tof_qh
    • The TOF charge
  • Float_t tof_qhrms
    • The
  • Float_t tof_qhlay[4]
    • The
  • Float_t tof_chi2c
    • The TOF betaH coordinate fitting Chi2.
  • Float_t tof_chi2t
    • The TOF BetaH time fitting Chi2.
  • Float_t tof_edeplay[4]
    • The
  • Int_t tof_nclslay_offtrack[4]
    • The
  • Float_t tof_edeplay_offtrack[4]
    • The
  • Float_t tof_edeplay_2nd[4]
    • The
  • Int_t tof_nclslay_2nd[4]
    • The

TRD

  • Int_t trd_itrd
  • Int_t trd_ihtrd
  • Int_t trd_type
    • The
  • Int_t trd_ntrdhits
    • The number of total TRD hits
  • Int_t trd_nseg
    • The number of TRD segments
  • Float_t trd_Q
    • The charge estimated from TRD track
  • Int_t trd_nhitsR
    • The number of recorded trd hits.
  • Float_t trd_hitamp[300] //[td_nhitsR]
    • The amplitude of recorded (max 300) trd hits.
  • Short_t trd_tubeid[300] //[td_nhitsR]
    • The tube id of recorded (max 300) hits.
  • Int_t trd_nhtrd
    • The (historical) combined number of TRD tracks: trd_nhtrd/100: number of TrdTrack. trd_nhtrd%100: number of TrdHTrack.
  • Int_t trd_nhits[3]
    • The number of hits used in likelihood calculation for inner(0), maxspan(1) and trd refit(2)
  • Int_t trd_offhits[3]
    • Then umber of off-track(nearby) hits used in likelihood for inner(0), maxspan(1) and trd refit(2)
  • Int_t trd_fitcode[3]
    • The (tracker) fit code used in likelihood calcluation for inner(0) and maxspan(1) .
  • Float_t trd_kraw[3][3]
    • The TRD likelihood from inner(0), maxspan(1) and trd refit(2) track position (1st dimension) and for electron(0), proton(1) and helium(2) particle (2nd dimension), calculated from rigidity.
  • Float_t trd_klkhd[3][3]
    • The TRD likelihood ratio from inner(0), maxspan(1) and trd refit(2) track postion (1st dimension) and for e/p(0), e/He(1), p/He(2) ratios (2nd dimension), calculated from rigidity.
  • Float_t trd_klkhd_e[3][3]
    • The likelihood ratio calculated from ECAL energy estimation.
  • Float_t trd_theta
    • The theta of TRD track.
  • Float_t trd_phi
    • The
  • Float_t trd_ecalentry[3]
    • The
  • Float_t trd_ecalexit[3]
    • The
  • Float_t trd_ipchi2
    • The
  • Float_t trd_trds[3]
    • The trds likelihood estimation (Aachen version) for e/p, e/He, p/He ratios.
  • Float_t trd_asyd
    • The
  • Float_t trd_asye
    • The

RICH

  • Int_t rich_inv
    • The
  • Int_t rich_jjok
    • The
  • Short_t rich_aero
    • The
  • Short_t rich_good
    • The
  • Short_t rich_clean
    • The
  • Float_t rich_beta
    • The
  • Float_t rich_errb
    • The
  • Float_t rich_prob
    • The
  • Float_t rich_udist
    • The
  • Short_t rich_used
    • The
  • Short_t rich_usedm
    • The
  • Short_t rich_pmts
    • The
  • Float_t rich_npes
    • The
  • Float_t rich_npese
    • The
  • Float_t rich_rinz
    • The
  • Int_t rich_status
    • The
  • Float_t rich_betap
    • The
  • Float_t rich_rinzp
    • The
  • Float_t rich_lkhd
    • The
  • Float_t rich_empnt[3]
    • The

Comments:


Edit | Attach | Watch | Print version | History: r15 < r14 < r13 < r12 < r11 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r15 - 2017-04-13 - TangZhicheng
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2018 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback