Mofia Ntuples in TWIST

Last updated 03.05.09 by Robert MacDonald.

The most current description of the ntuples available for use in Mofia can be found in the comment blocks in the file helixntuple_mod.f90 in Mofia's source/user directory. NTuple type is selected by setting a Mofia namelist variable HelixNtupleType in the HelixFit namelist. For example, the command:

name HelixFit HelixNtupleType = 8
selects Ntuple type 8.

As of this writing, Ntuple type 7 requires some "fancy" indexing, since it contains information about all tracks in all windows. Blair Jamieson explains in a Forum posting:

I also modified ntyple type 7 to contain information about all windows
and detector stream fits.  This could be used if we want to compare
fits for multiple beam positrons or decay positrons in an event.  Each
muon track has one track (upstream), each decay positron has two
tracks (upstream, and downstream), and each beam positron has three
tracks (upstream, bothstream...probably bogus, and downstream).  I
broke up the indicies into fit types for this ntuple type into 3
indicies: muidx, dkidx, and eidx.  For example if there was two
muons(fits 1and 3), two decay positrons(fits 2,3, and,5,6), and a beam
positron (fits 7,8,9), then the index arrays would be:

muidx ( 1, 3, 0, 0, 0, 0, 0, 0, 0)
dkidx ( 2, 3, 5, 6, 0, 0, 0, 0, 0)
eidx  ( 7, 8, 9, 0, 0, 0, 0, 0, 0)

ROOT users (and potential ROOT users!) will also be interested in Konstantin's ROOT Tree functionality. For more information, see Mofia's source/user/ko_tdctree.f90 file.


Return to the TWIST Mofia home page.