# C-shell script to run STR job for B=10
#
source /usr/local/triumf/setup.com
cd ~chensm/utc/tdc
setup kofia 2.1
date
setenv KOFIALOG /home/chensm/utc/tdc/findstr1998.klog
setenv MTIN /home/chensm/utc/tdc/run1998.skim1
unsetenv TRAP_FPE
#
# The parameters for FUNC are
#  1: Save output histograms
#  2: Clear the histograms
#  3: Calculate a new STR
#  4: Write the STR's to a file
#  5: Read histograms from a file. This also sets the max times for which
#      the STR's are valid.
#  6: Modify the HBOOK fitting options
#  7: Modify the window for accumulating avg. residuals
#  8: Modify the STR fitting mode
#      0: One STR per layer, same for left and right side of cell
#      1: STR left and STR right are computed separately
#      2: Same as 0, except mean residual per time bin is fit for
#          rather than taken from a profile histogram
#      3: Same as 1, except mean residual per time bin is fit
#  9: Set the mode for calculating residuals
#      1: Full track fit. Include 'residual' hit in fit.
#      2: Full track fit. Exclude 'residual' hit from fit.
#      3: Local track fit. Include 'residual' hit.
#      4: Local track fit. Exclude 'residual' hit.
# 10: Set the number of time bins, the max TDC time and the max
#      residual for the time vs residual histograms
# 11: Set the fraction of the residual to be applied as an STR correction
# 12: Set the minimum number of entries required in a time bin in
#      order to compute an avg. residual
# 13: Set the max. drift distance allowed for a hit to be used in a
#      fit.  (The units are in half-cell widths.) Also set the
#      max. drift time (for each layer) allowed in a fit.
#      CURRENTLY NOT IMPLEMENTED FOR RESMODE=1.
# 14: Set the min. number of layers required in the original track fit
#      and the min. number of layers required in subsequent track refits
# 15: If .eq. 1 then for each layer included in the fit, use only the
#      hit with the smallest drift time.
# 16: If .eq.1 then compute a STR in angular bins.
#     Also set the minimum number of entries required in a time/angle bin.
# 17: Set the cell fraction above which to compute an STR correction.
#     Only meaningful when computing the STR in angle bins.
# 18: If .eq. 1 then force the STR to be monotonically increasing.
# 19: If .eq. 1 then constrain the STR's at t=0.
#
findstr.exe <<+
use UT PED '/home/chensm/utc/tdc/t0fit1998.txt'
@findstr_b.kcm
namelist utxytrk
utcxyres=2.0
&
func 6
Q
func 7
3.0
func 8
3
func 9
2
func 10
50
250
3.0
func 11
0.5
func 12
50.
func 13
2.0
500.,500.,500.,500.,500.,500.,500.,500.,500.,500.,500.,500.
func 14
10
8
func 15
0
func 16
1
10
func 17
0.0
func 18
0
func 19
1
files 100
an
func 3
func 4
/home/chensm/utc/tdc/findstr1998.txt
func 1
/home/chensm/utc/tdc/findstr1998.paw
exit
+
date
#mt -t /dev/mx6v rewind