Tutorial: MOFIA kcm files

K.Olchanski, 17-Jan-2002

What's a KCM file?

A KCM file is a list commands telling MOFIA to use certain configuration files (geometry, fastbus maps, etc), to set the analysis parameters ("namelist variables"), open input and output data files, and what to do with the data in general.

Corresponding to the 4 types of data (magnet on or off, and real or simulated data), there are 4 families of kcm files: depending on real vs simulated data, we use different configuration files, and depending on magnet on vs off we enable different reconstruction code.

Basic rules: one command per line, blank lines are okey, anything that starts with a bang (!) is a comment. To execute a kcm file from the MOFIA> use the @filename command. Use the HELP command to obtain a complete list of MOFIA commands.

To have MOFIA always execute a kcm file on startup, set the environment variable MOFIA_INIT to the name of that kcm file: setenv MOFIA_INIT filename

Example kcm file: real data, magnet on

This file was used during the November run to run "mofia" and "photo" on the data files. It can be found at /home/twistonl/online_anal/raw.kcm.

namelist hist NEVTPROCESSED = 2000      ! save histograms every 2000 events
set dc/cuts on                          ! turn on DC cuts (RTFM!)
name dccuts DC_MINTDC_CUT=-10000.0
name dccuts DC_MAXTDC_CUT=20000.
name dccuts dc_xtalk_wcut = 50.0
name pccuts pc_xtalk_wcut = 50.0
name firstguess winpcthreshold = 500.0  ! set the windowing
name firstguess windcstart = -100.0     ! parameters (RTFM!)
name firstguess windcend = 1000.0
set photo on                            ! enable the graphics
                                        ! if running the "photo"
                                        ! executable.
name global     BField = 2.0            ! set the magnetic field
name global     unpackmc = F            ! this is NOT simulated data
name firstguess enableWindowing  = T    ! enable windowing
name firstguess enableClustering = T    ! enable clustering
name KalmanCuts enableKalman     = F    ! disable the straight-track
                                        ! reconstruction
name firstguess enableFirstGuess = T    ! enable first-guess
name firstguess enableFirstGuessNtuple = F ! do not make f-g ntuples
name helixfit   enableHelixFit    = T   ! enable the helix fitter
name helixfit   enableHelixNtuple = F   ! do not make helix ntuples

                                        ! these commands are commented out:
!name timezero   triggertime = 10181.0    ! set the global time offset
!mtin "/twist/data_onl/2001/run02673.ybs" ! open an input file
!ana                                      ! analyze all events
!save                                     ! save the histograms
                                          ! (remember that they are also
                                          ! saved after every 2000 events!)

Example kcm file: simulated data, magnet on

This is the helix_mc.kcm file from mofia/source/user CVS repository, recommended for analyzing the geant files.

! Use the same CFM files as were used by GEANT (see the geant e614.com file)
use dt geo   "/home/e614/e614soft/caldb_ascii/dt_geo.00021"
use fbc1 map "/home/e614/e614soft/caldb_ascii/fbc1_map.00007"
use fbc2 map "/home/e614/e614soft/caldb_ascii/fbc2_map.00007"
use dc str   "/home/e614/e614soft/caldb_ascii/dc_drift_times.00001"
use dc t0  "/home/e614/e614soft/caldb_ascii/dc_t0.mc"
use pc t0  "/home/e614/e614soft/caldb_ascii/pc_t0.mc"
use sc t0  "/home/e614/e614soft/caldb_ascii/sc_t0.mc"
use pc adc "/home/e614/e614soft/caldb_ascii/pc_adc.mc"
use sc adc "/home/e614/e614soft/caldb_ascii/sc_adc.mc"
use dc res  "/home/e614/e614soft/caldb_ascii/dc_res.mc"

set dc/cuts on                          ! enable DC cuts (RTFM!)
namelist hist NEVTPROCESSED = 2000      ! save histograms every 2000 events
name global     BField = 2.2            ! set the magnetic field
                                        ! (should match the geant FFCARD)
name timezero triggertime = 0.00001     ! set the global time offset
                                        ! (required for geant data!)
set photo on                            ! enable the "photo" graphics
name global     unpackmc = T            ! tell MOFIA: this is simulated data
name KalmanCuts enableKalman     = F    ! disable magnet-off code

! Setup first guess
name firstguess enableFirstGuess   = T  ! enable the code
name firstguess enableFirstGuessNtuple = F ! disable ntuples
name firstguess enableFirstGuessJim= F  ! select which f-g code to run
name firstguess enableFirstGuessArt= F
name firstguess enableFirstGuessKO = T
name firstguess FirstGuessVerbose  = 0  ! do not produce debug output

! Setup thehelix fitter
name helixfit   enableHelixFit    = T   ! enable the code
name helixfit   enableHelixNtuple = T   ! enable ntuple output
name helixfit   HelixFitDPDS      = 0.001 ! set the energy loss (should
                                          ! match the geant LOSS ffcard)
name helixfit   HelixFitVerbose   = 0   ! disable the debug output

                                    ! open the input file
mtin "/twist/data4/olchansk/twist/mc/mc-clean.ffcards.0001/run1.dat"
ana                                 ! analyze all events
save                                ! save the histograms
exit                                ! stop mofia
Konstantin Olchanski
Last modified: Thu Jun 20 14:43:37 PDT 2002