# C-shell script to run Kofia
#
source /usr/local/triumf/setup.com
cd ~chensm/utc/tdc
setenv scratch /usr/people/chensm/utc/tdc
setup kofia 2.1
date
setenv KOFIALOG $scratch/basic_1998.klog
#setenv MTIN $scratch/km21_utcadc100.dat

setenv MTIN /dev/mx1d
#
# Parameters are:
#  PARAM(65): If non-zero, then diagnostic printing is turned on
#  PARAM(66): Nblank_max for clustering; see utc_z_cuts.cmn
#  PARAM(67): Low edge of the finely binned TDC histograms
#  PARAM(68): Low edge of the finely binned anode ADC histograms
#  PARAM(69): High edge of the finely binned anode ADC histograms
#  PARAM(70): Low edge of the finely binned cathode ADC histos
#  PARAM(71): High edge
#  PARAM(72): Min. pulse height for ADC profile plots
#  PARAM(73): Max. pulse height for ADC profile plots
#  PARAM(74): If 1 then turn on TDC pedestal subtraction
#  PARAM(75): If 1 then turn on ADC ped. subtraction
#  PARAM(76): If 1 then just look at the first TDC hit on each
#              wire/strip
#  PARAM(77): If .gt.0. then this becomes the number of sigma(noise)
#              cut that is applied before an entry is made in an ADC
#              histogram
#  PARAM(78): Min time for TDC profile plots
#  PARAM(79): Max time for TDC profile plots
#  PARAM(80): If 1 then order the anode channels according to the TDC
#              numbering.
#  PARAM(81): Min TDC required for making an entry in an anode ADC plot
#  PARAM(82): Max TDC
#  PARAM(83): Same as PARAM(81) except for strips
#  PARAM(84): Same as PARAM(82) except for strips
#  PARAM(85): N sigma cut for determining noisy anode TDC channels
#  PARAM(86): Same as PARAM(85) except for strips
#  PARAM(87): Same as PARAM(85) for ADC's
#  PARAM(88): Same as PARAM(86) for ADC's
#  If the thresholds param(85) through param(88) are set to 0 sigma
#  then we will flag any channel with occupancy less than 50% of the
#  average or more than twice the average.
#  PARAM(89): 0 if we want to look at all events
#             1 if we want to run RD_TRK (ADC version) first and look
#                only at events containing a track with energy
#                between param(90) and param(91)
#             2 same as 1 except TD version of RD_TRK is used
#  PARAM(90): Min. RD_TRK energy
#  PARAM(91): Max. RD_TRK energy
#
#  use param(68)=300, param(69)=800, param(70)=300, param(71)=800
#   if not ADC ped subtracted
#
# FUNC values:
#   1: Analyze histograms and print out list of dead/noisy channels
#   2: Output the histograms to a file
#   3: Reset the histograms
#
nice +4 basic.exe <<+
namelist hcuts
param(65)=0
param(66)=2
param(67)=0.
param(68)=0.
param(69)=4000.
param(70)=0.
param(71)=4000.
param(72)=0.
param(73)=2000.
param(74)=0.
param(75)=1.
param(76)=1.
param(77)=3.
param(78)=0.
param(79)=500.
param(80)=0.
param(81)=10.
param(82)=60.
param(83)=10.
param(84)=60.
param(85)=5.
param(86)=5.
param(87)=5.
param(88)=5.
param(89)=0.
param(90)=88.
param(91)=120.
&
namelist utxypat
utangwind=t
utanglo=-15.
utanghi=15.
&
use UA GAI "/usr/people/chensm/utc/tdc/ua_gai.00001"
use UC GAI "/usr/people/chensm/utc/tdc/uc_gai.00001"
files 100
an 
show errors
func 2
/usr/people/chensm/utc/tdc/basic1998_3.paw
func 1
exit
+
date
mt -t /dev/mx1d rewind