X.2.1 Event Filtering Failure Codes
X.2.2 Pattern recognition Failure Codes
X.2.4 Kalman Filtering Failure Codes
X.3.8 Kalman Filter Structures
X. Appendecies
X.1 Namelist Variables
MOFIA> show name
NAMELIST DESCRIPTION:
BATCH BATCH LOG control
parameters
HCUTS DPLOT user cuts
DCSET Drift Chamber
SETtings
PCSET Proportional
Chamber SETtings
HIST HISTogramming
parameters
PHOTO PHOTO flags
DCCFLAGS Drift Chamber
Calibration FLAGS
SCSET SCintillator
SETtings
SCCUTS SCintillator CUTS
RFCUTS RF CUTS
DCCUTS Drift Chamber CUTS
PCCUTS Proportional
Chamber CUTS
KPUNIT MOFIA print units
SCFLAGS SCintillator
FLAGS
KFLAGS MOFIA execution
control flags
PRCNTL Print control flags
GLOBAL GLOBAL settings
QOD QOD Monitor
params
STRSET STR SETtings
KalmanCuts Kalman
Tracking Cuts
Alignment Alignment
parameters
FirstGuess First Guess
parameters
HelixFit HelixFit
parameters
TimeZero Time zero fit
settings
Efficiency Efficiency
settings
REZS REZolutionS
controls
MOFIA> show name DCSET
NameList DCSET: Drift Chamber SETtings
FirstPlaneDC = 1 (dflt =
1): First DC plane
LastPlaneDC = 44 (dflt =
44): Last DC plane
MOFIA> show name PCSET
NameList
PCSET: Proportional
Chamber SETtings
FirstPlanePC = 1 (dflt =
1): First PC plane
LastPlanePC = 12 (dflt =
12): Last PC plane
MOFIA> show name hist
NameList HIST: HISTogramming parameters
nEVTprocessed = -1 (dflt = -1,OFF): Write hists every
nEVTprocessed events
FillRawHist = T
(dflt = T) : Fill raw histograms
FillHist = T
(dflt = T) : Fill histograms after initial filtering
FillTrackHist = T
(dflt = T) : Fill tracking histograms
FillPatternHist = T
(dflt = T) : Fill pattern recognition histograms
PulserHistToggle = T
(dflt = F) : Fill histograms with
random pulser data(T) or normal triggers(F)
FillPhysicsHist = T
(dflt = T) : Fill physics histograms
FillFirstGuessHist = F
(dflt = T) : Fill helix first guess histograms
FillXtalkHist = T
(dflt = T) : Fill cross talk histograms
PlaneHists = T
(dflt = T) : Generate individual plane histograms
Globalmem_toggle = T
(dflt = T) : Turn Global Memory Section on (T) /off (F)
TDC_MIN = 0. (dflt = 0.) : Lower limit on TDC
spectra
TDC_MAX = 6000. (dflt = 6000.) : Upper limit on TDC
spectra
Raw_XMI_tdcslot = 0. (dflt = 0.) : Lower limit on RAW TDC
plots
Raw_XMA_tdcslot = 30000. (dflt = 30000.) : Upper limit on RAW
TDC plots
Raw_NX_tdcslot = 3 (dflt = 3) : Number of RAW x Channels for
TDC plots
WEvent_per_plane = F
(dflt = F) : Turns on/off #of wires hit/event/plane hists
HMult_times = F
(dflt = F) : Turns on/off Time difference between multiple hits on a wire hists
MOFIA> show name
dccflags
NameList DCCFLAGS: Drift Chamber Calibration FLAGS
FindPlanePos = F (dflt =
F) : Find plane positions
FindWirePos = F (dflt =
F) : Find wire positions
FindPlaneRot = F (dflt =
F) : Find plane rotations
FindTDC0 = F (dflt =
F) : Find TDC time zero
FindResolution = F (dflt = F) : Find drift chamber resolution
MOFIA> show name dccuts
NameList DCCUTS: Drift Chamber CUTS
DC_MAXTDC_CUT = 20000.00 (dflt = -1, OFF): Max TDC
channel cut
DC_MINTDC_CUT = *******
(dflt = -1, OFF): Min TDC channel cut
DC_MAXWTDC_CUT = 150.00 (dflt = -1, OFF): Max TDC width cut
DC_MINWTDC_CUT = -1.00 (dflt = -1, OFF): Min TDC width cut
DC_MAX_HITS_IN_PLANE = -1
(dflt = -1, OFF): Max hit wires in plane cut
DC_MIN_PLANES = -1 (dflt = -1, OFF): Min planes cut
DC_XTALK_WCUT = -1.00 (dflt = -1, OFF): Min cross talk TDC
width cut
DC_NOISE_WCUT = -1.00 (dflt = -1, OFF): Min noise TDC width
cut
MOFIA> show name pccuts
NameList PCCUTS: Proportional Chamber CUTS
PC_MAXTDC_CUT = 6000.00 (dflt = -1, OFF): Max TDC channel
cut
PC_MINTDC_CUT = 0.00 (dflt = -1, OFF): Min TDC channel cut
PC_MAXWTDC_CUT = 150.00 (dflt = -1, OFF): Max TDC width cut
PC_MINWTDC_CUT = -1.00 (dflt = -1, OFF): Min TDC width cut
PC_MAX_HITS_IN_PLANE = -1
(dflt = -1, OFF): Max hit wires in plane cut
PC_MIN_PLANES = -1 (dflt = -1, OFF): Min planes cut
PC_XTALK_WCUT = -1.00 (dflt = -1, OFF): Min cross talk TDC
width cut
PC_NOISE_WCUT = -1.00 (dflt = -1, OFF): Min noise TDC width
cut
MOFIA> show name global
BField (def = 0.)
0.000000
UnpackMC = F
MOFIA> show name
alignment
NameList Alignment: Alignment parameters
nEventsPlane (def = 10000) =
10000 Total number of events per iteration for calculating plane
positions
nEventsWire(def = 100000) =
100000 Total number of events per iteration for calculating wire
positions
AlignAngleU
(def = 0.0) = 0.000 Angle of the U
planes transverse-alignment line with respect to the beam
AlignAngleV (def = 0.0) =
0.000 Angle of the V planes transverse-alignment line with respect to
the beam
FixPlanes (def = FALSE) =
F Choose a line for transverse alignment defined by any 2 U and 2 V
planes
FixedPlane1 (def = 1) =
1 Plane number for one of the four fixed planes used to define the
alignment line
FixedPlane1 (def = 1) =
2 Plane number for one of the four fixed planes used to define the
alignment line
FixedPlane1 (def = 1) =
7 Plane number for one of the four fixed planes used to define the
alignment line
FixedPlane1 (def = 1) =
8 Plane number for one of the four fixed planes used to define the
alignment line
MOFIA> show name efficiency
Namelist Efficiency: Efficiency calculation parameters
RadiusCutDense
= 5.000 (dflt = 15 cm): DC
dense stack radius cut
RadiusCutSparse
= 15.000 (dflt = 15 cm): DC
sparse stack radius cut
CellCut = 2 (dflt = 2): Max number of adjacent cells
to investigate on each side of the cell expected to have a hit.
CellCut = 10 (dflt = 10):
Minimum number of hit planes
FindEff = F (dflt =
.FALSE.): Calculate chamber efficiency
MOFIA> show name timezero
NameList TimeZero: T0 fit settings
T0_TDC_MIN (def = 3000.) =
3000.0 Min range for TDC spectra rising time
T0_TDC_MAX (def = 3080.) =
3080.0 Max range for TDC spectra rising time
TDCnsPerBin (def = 0.5) =
0.5 Number of bins per channel
FitT0 (def = FALSE) = F Accumulate histograms for fitting T0
spectra
TimeBackwards (def = TRUE) = T Time increases backwards
TriggerTimeDC = 0.0
TriggerTimePC = 0.0
TriggerTimeSC = 0.0
TCAP Cut Low Time = -1.0
TCAP Cut High Time = -1.0
MOFIA>
show
name rezs
Namelist REZS: Resolutions
nEventsMax = 5000 (dflt = 5000): # Events per Resolution
Iteration
MinHistFitEntries = 1000 (dflt = 1000): # histogram entries
necessary to do fit
nResidualBins = 300 (dflt = 300): # Bins in residual
histograms.
MinResidualValue = 0.30 (dflt = -0.3): Min histogram channel
(cm).
MaxResidualValue = 0.00 (dflt = 0.3): Max histogram channel
(cm).
PLEASE do NOT mess with the histogram definitions without THINKING
about the consequences. See resolution_mod
for details
calcResidualMethod = *** (dflt = 3): Residual calculation method.
1=> Gaussian Fit to residual histograms.
2=> SquaredSum calculation.
3=> FWHM calculation.
peakTOL =
MOFIA> show name sccuts
NameList SCCUTS: SCintillator CUTS
S1_MAX_NHITS = -1 (dflt = -1, OFF): Scint 1 max number of
hits cut
S2_MAX_NHITS = -1 (dflt = -1, OFF): Scint 2 max number of
hits cut
S3_MAX_NHITS = -1 (dflt = -1, OFF): Scint 3 max number of
hits cut
S1_WIDTH_CUT = -1.00
(dflt = -1, OFF): Scint 1 peak width cut
S2_WIDTH_CUT = -1.00
(dflt = -1, OFF): Scint 2 peak width cut
S3_WIDTH_CUT = -1.00
(dflt = -1, OFF): Scint 2 peak width cut
S1_MAX_TDC = *******
(dflt = -1, OFF): Scint 3 peak width cut
S2_MAX_TDC = *******
(dflt = -1, OFF): Scint 1 max TDC channel cut
S3_MAX_TDC = 2500.00
(dflt = -1, OFF): Scint 2 max TDC channel cut
S1_MIN_TDC = -1.00 (dflt = -1, OFF): Scint 1 min TDC
channel cut
S2_MIN_TDC = -1.00 (dflt = -1, OFF): Scint 2 min TDC
channel cut
S3_MIN_TDC = -1.00 (dflt = -1, OFF): Scint 2 min TDC
channel cut
S1_MAX_ADC = 5000.00
(dflt = -1, OFF): Scint 1 max ADC channel cut
S2_MAX_ADC = 5000.00
(dflt = -1, OFF): Scint 2 max ADC channel cut
S1_MIN_ADC = -1.00 (dflt = -1, OFF): Scint 1 min ADC
channel cut
S2_MIN_ADC = -1.00 (dflt = -1, OFF): Scint 2 min ADC
channel cut
MOFIA> show name rfcuts
NameList RFCUTS: RF CUTS
RF_MIN_TDC = -1.00 (dflt
= -1, OFF) : RF min TDC cut
RF_MAX_TDC = -1.00 (dflt
= -1, OFF) : RF max TDC cut
MOFIA> show name strset
use_cos_increments (def= T) T str_angle_inc (def= 0.05) 5.000000E-02
str_upper_angle_limit (def= 85.
degrees) 85.0000
MOFIA> show name FirstGuess
NameList FirstGuess: First Guess parameters
enableFirstGuess = T (logical: T or F)
enableFirstGuessNtuple = F (logical: T or F)
NameList FirstGuess: Time separation between windows
winPCthreshold = 500.000 (ns)
NameList FirstGuess: DC window start
winDCstart = -100.000 (ns)
NameList FirstGuess: DC window end
winDCend = 1000.000 (ns)
MOFIA> show name HelixFit
NameList HelixFit: HelixFit parameters
enableHelixFit= T
enableHelixNtuple= F
HelixFitVerbose= 0
HelixFitDPDS= 0.000000 energy loss, (MeV/cm)
HelixFitStartFG, StartMC= T F
HelixFitUpstream, HelixFitDownstream= T T
HelixFitCosTmin, HelixFitCosTmax= 0.0100
1.1000
HelixFitWireRes, HelixFitDriftRes,
HelixFitTimeRes= 0.1600 0.0300
50.0
HelixFitCutWC, CutDrift, CutTime, CutTref,
CutFit= -1.0000 -1.0000 -1.0000 0.0500
0.0200
HelixFitMaxIter= 20
MOFIA> Show name KalmanCuts
NameList KalmanCuts: Kalman tracking cuts
ChiDiffCluster (def = 1.E-02) = 0.10E-01 Chi2 convergence level
for cluster iteration
ChiDiffTime (def = 1.E-04) = 0.10E-03 Chi2 convergence level for
timing iteration
Chi2sCutCluster (def = 1.E05) = 0.10E+06 Chi2 cut for cluster iteration
Chi2sCutTime (def = 1.E05) = 0.10E+06 Chi2 cut for timing
iteration
MaxIterateCluster
(def = 50) = 50 Maximum number of
iterations for cluster fit
MaxIterateTime (def = 50) =
50 Maximum number of iterations for timing fit
EnableKalman (def = TRUE) =
T Execute Kalman filterng code
SwitchLR = (def = TRUE) =
T Attempt reducing the value of Chi2 by switching left and right
NoiseExcludeMax (def = 2) =
MOFIA> show name photo
NameList PHOTO: PHOTO FLAGS
IPIC = 1 (dflt = 1) : Chooses priority PHOTO view
HARD = F (dflt = F) :
Set TRUE for automatic PHOTO hardcopy
EDGR = F (dflt = F) :
Set TRUE to invoke EDGR picture editor
PRINTER = 2 (dflt = 2) : Printer is one of Printronix,
HPLaser HPThinkJet, La100, HPPaintJet
DURATION = 0.0 (dflt=0.0) : duration of automatic
photos
DETAIL = 0 (dflt = 0) : number of times to
automatically add detail to photo
SHOW_COORD = F (dflt = F) : Set TRUE to display counter
coordinates
DRAWFOILS = T (dflt = T) :
Set TRUE to display UTC superlayer boundaries
TrackedOnly = F (dflt = F)
: IF TRUE skip events with no tracks
MOFIA> show name qod
n_QOD_buffers (def = 4) 4
qod bufferlength (def =10000) 10000
short term warning message prob
1.000000E-03 long term warning message prob 1.000000E-04 Ratio of hot wire counts to
base wire counts 5.00000 baseline
filename (default )
Gives DC wire occupancies, INCLUDING all multiple hits.Set the variable
to the plane number you'd like to see.
0 (Default = 0, no hist)
MOFIA> show name
scflags
NameList SCFLAGS: SCintillator FLAGS
S1_signal = T (dflt = TRUE) : Require a signal in S1
S2_signal = T (dflt = TRUE) : Require a signal in S2
S3_signal = T (dflt = TRUE) : Require a signal in S3
X.2 Failure Codes
X.2.1 Event Filtering Failure Codes
|
Failure Code |
Description |
Associated Namelist |
|
1 |
successful
unpacking |
|
|
2 |
at least
one DC hit present |
|
|
3 |
there is
less than S1_MAX_NHITS in scint 1 |
SCCUTS |
|
4 |
there is
less than S2_MAX_NHITS in scint 2 |
SCCUTS |
|
5 |
require a
signal in scint 1 |
SCCUTS |
|
6 |
require a
signal in scint 2 |
SCCUTS |
|
7 |
|
|
|
8 |
|
|
|
9 |
time from
scint 1 is close to its peak location |
SCCUTS |
|
10 |
time from
scint 2 is close to its peak location |
SCCUTS |
|
11 |
number of
hits in DC plane is less than DC_MAX_HITS_IN_PLANE |
DCCUTS |
|
12 |
require
the numbr of DC hits in each plane to be less than DC_MIN_PLANES |
DCCUTS |
|
13 |
require
the total number of DC hits to be less than DCmaxAllowedHits |
DCCUTS |
|
14 |
number of
hits in PC plane is less than PC_MAX_HITS_IN_PLANE |
PCCUTS |
|
15 |
require
the numbr of PC hits in each plane to be less than PC_MIN_PLANES |
PCCUTS |
|
16 |
require
the RF time to be within the range [RF_MIN_TDC,RF_MAX_TDC] |
RFCUTS |
|
17 |
require the
total number of PC hits to be less than PCmaxAllowedHits |
PCCUTS |
|
18 |
there is
less than S3_MAX_NHITS in scint 3 |
SCCUTS |
|
19 |
require a
signal in scint 3 |
SCCUTS |
|
20 |
time from
scint 3 is close to its peak location |
SCCUTS |
Table 7 Event filtering failure
codes.
X.2.2 Pattern recognition Failure Codes
|
Failure Code |
Description |
Namelist |
|
1 |
track > MaxTmpTracks in InsertCluster |
|
|
2 |
iPair = PrevPair(iFGTrack) in InsertCluster |
|
|
3 |
FGCl(iFGTrack) % nCl >= MaxPairs in InsertCluster |
|
|
4 |
Cluster does not fit circle in InsertCluster |
|
|
5 |
Less than three clusters in ResolveCircle |
|
|
6 |
ChiSquare unacceptably large in ResolveCircle |
|
|
7 |
nRows in matrices A and B are unequal in MatSolv
(Matrix_mod) |
|
|
8 |
Unique solution does not exist for matrix equation
in MatSolv (Matrix_mod) |
|
|
9 |
0 row in matrix in Triangulate (Matrix_mod) |
|
|
10 |
Failed to find pivot in Triangulate (Matrix_mod) |
|
|
11 |
ChiSquare unacceptably large in ResolvePhiLambda |
|
|
12 |
ChiSquare unacceptably large in ResolvePhiLambda |
|
Table 8 Pattern recognition failure codes.
X.2.3 c2 Fit Failure codes
X.2.4 Kalman Filtering Failure Codes
|
Failure Code |
Description |
Namelist |
|
1 |
Not enough hits for tracking |
|
|
2 |
Filter failed in cluster fit |
|
|
3 |
Smoother failed in cluster fit |
|
|
4 |
Bad track Chi2 – cluster iteration |
|
|
5 |
Failed computing hit position |
|
|
6 |
Filter failed in timing fit |
|
|
7 |
Smoother failed in timing fit |
|
|
8 |
Bad tracking Chi2 - timing iteration |
|
|
9 |
Failed computing physics paramters |
|
Table 9 Kalman filter failure codes.
X.3 Data Structures
X.3.1 Geometry Structures
X.3.2 TDC Structures
X.3.3 Calibrations Structures
X.3.4 Windowing Structures
X.3.5 Clustering Structures
X.3.6 First Guess Structures
X.3.7c2 Helix Fit Structures
X.3.8 Kalman Filter Structures
X.3.9 MC Banks Structures
|
|
X.4 Flowcharts
X.4.1 Initialization Branch
X.4.2 Analysis Branch
