opt nstat
opt nfit 
set ltyp 1
set pmci
set mtyp

ve/cr error_adc(30) R 30*15.

ve/cr res_5_om(36) I 36*0
ve/cr res_5_a(36) R 36*0.
ve/cr res_5_b(36) R 36*0.
ve/cr res_5_c(36) R 36*0.
ve/cr res_5_d(36) R 36*0.
ve/cr res_5_e(36) R 36*0.
ve/cr res_5_ll(36) R 36*0.
ve/cr res_5_ul(36) R 36*0.

ve/cr res_6_om(36) I 36*0
ve/cr res_6_a(36) R 36*0.
ve/cr res_6_b(36) R 36*0.
ve/cr res_6_c(36) R 36*0.
ve/cr res_6_d(36) R 36*0.
ve/cr res_6_e(36) R 36*0.
ve/cr res_6_ll(36) R 36*0.
ve/cr res_6_ul(36) R 36*0.

ve/cr res_7_om(36) I 36*0
ve/cr res_7_a(36) R 36*0.
ve/cr res_7_b(36) R 36*0.
ve/cr res_7_c(36) R 36*0.
ve/cr res_7_d(36) R 36*0.
ve/cr res_7_e(36) R 36*0.
ve/cr res_7_ll(36) R 36*0.
ve/cr res_7_ul(36) R 36*0.

ve/cr res_8_om(36) I 36*0
ve/cr res_8_a(36) R 36*0.
ve/cr res_8_b(36) R 36*0.
ve/cr res_8_c(36) R 36*0.
ve/cr res_8_d(36) R 36*0.
ve/cr res_8_e(36) R 36*0.
ve/cr res_8_ll(36) R 36*0.
ve/cr res_8_ul(36) R 36*0.

ve/cr res_9_om(36) I 36*0
ve/cr res_9_a(36) R 36*0.
ve/cr res_9_b(36) R 36*0.
ve/cr res_9_c(36) R 36*0.
ve/cr res_9_d(36) R 36*0.
ve/cr res_9_e(36) R 36*0.
ve/cr res_9_ll(36) R 36*0.
ve/cr res_9_ul(36) R 36*0.

ve/cr res_10_om(36) I 36*0
ve/cr res_10_a(36) R 36*0.
ve/cr res_10_b(36) R 36*0.
ve/cr res_10_c(36) R 36*0.
ve/cr res_10_d(36) R 36*0.
ve/cr res_10_e(36) R 36*0.
ve/cr res_10_ll(36) R 36*0.
ve/cr res_10_ul(36) R 36*0.

* The below hybrid designations come from a list of hybrid OM's.
* Ex: for string 5, the first OM position is a hybrid (1), the
* second position is standard (0), etc.  By the end of string 5,
* the last 25 OM's are standards (0).   NOTE: in the finished
* .result file, hybrid is designated as 0.0, while standards 
* are designated as 1.0.
ve/cr hybrids_5(36) I 1 0 1 0 1 0 1 0 1 0 1 25*0
ve/cr hybrids_6(36) I 1 0 1 0 1 0 1 0 1 0 1 25*0
ve/cr hybrids_7(36) I 31*0 1 0 0 0 1
ve/cr hybrids_8(36) I 1 0 1 0 1 0 1 0 1 0 1 25*0
ve/cr hybrids_9(36) I 36*0
ve/cr hybrids_10(36) I 36*0

* OPEN THE HBOOK FILE
*hi/file 1 test4.small.out.partcalib.string10.hbook
ntuple/chain -blah
ntuple/chain 'blah' /home/aheffwong/crosstalk/hbook/ab_2002_170_5781_011.data.mu.min_bias.string9.hbook
ntuple/chain 'blah' /home/aheffwong/crosstalk/hbook/ab_2002_170_5781_012.data.mu.min_bias.string9.hbook
ntuple/chain 'blah' /home/aheffwong/crosstalk/hbook/ab_2002_170_5781_087.data.mu.min_bias.string9.hbook
ntuple/chain 'blah' /home/aheffwong/crosstalk/hbook/ab_2002_170_5781_088.data.mu.min_bias.string9.hbook
ntuple/chain 'blah' /home/aheffwong/crosstalk/hbook/ab_2002_277_5895_005.data.mu.min_bias.string9.hbook
ntuple/chain 'blah' /home/aheffwong/crosstalk/hbook/ab_2002_277_5895_105.data.mu.min_bias.string9.hbook


* WE ITERATE OVER STRING AND OM NUMBERS
for string in  '9'
*'5' '6' '7' '8' '9'  '10'
for pos in _
'01' '02' '03' '04' '05' '06' '07' '08' '09' '10' _
'11' '12' '13' '14' '15' '16' '17' '18' '19' '20' _
'21' '22' '23' '24' '25' '26' '27' '28' '29' '30' _
'31' '32' '33' '34' '35' '36'


* DETERMINE BASIC INFO ABOUT THE OM
is_hybrid=$EVAL(hybrids_[string]([pos]))
om=$EVAL(([string]-5)*36+86+[pos])
*om=[pos]
title 'ADC-TOT FIT string'//$EVAL([string])//' OM'//$EVAL([om])//' POS'//$EVAL([pos])
mess POS [pos] OM [om] STRING [string]man.string[string].adc_tot_fit.result 

* WE WANT TO FIND THE PEAK OF TDC TO CLEAN THE DATA
1d 30 '' 132 0. 33000.
nt/proj 30 //blah/100.time ch=[pos].and.0<time<33000
mean_tmp = $hinfo(30,'mean')
1d 31 '' 200 $EVAL([mean_tmp]-2000.) $EVAL([mean_tmp]+2000.)
nt/proj 31 //blah/100.time ch=[pos].and.$EVAL([mean_tmp]-2000.)<time<$EVAL([mean_tmp]+2000.) 
tdc_mean= $hinfo(31,'mean')

nt/pl //blah/100.time ch=[pos].and.$EVAL([tdc_mean]-200.)<time<$EVAL([tdc_mean]+250.)
mess 'Now we have gotten here.'


* NOW WE OBTAIN THE EDGE BY HAND
nt/pl //blah/100.amp%tot ch=[pos].and.$EVAL([tdc_mean]-200.)<time<$EVAL([tdc_mean]+250.).and.0<tot<1000.and.0<amp<2000

mess '...and there.'

set pmci 2
vlocate vtot vadc *S

* THE SIZE OF THE VECTORS IS:
adc_dim=$vdim(vadc)

mess 'About to fit the vector...'
* Now we fit the the vector
IF ([is_hybrid].eq.0) then
   ve/cr par(4) R 0.5 0.02 15 180
   ve/fit vtot vadc error_adc adc_tot_fit.f 'WQN0' 4 par
*   ve/fit vtot vadc error_adc /home/aheffwong/crosstalk/adc_tot_fit.f 'WQN0' 4 par

ELSEIF ([is_hybrid].eq.1) then
   ve/cr par(4) R 4.0 0.005 0. 0.
   ve/fit vtot vadc error_adc adc_tot_fit_hyb.f 'WQN0' 2 par(1:2) 
*   ve/fit vtot vadc error_adc /home/aheffwong/crosstalk/adc_tot_fit_hyb.f 'WQN0' 2 par(1:2)

ENDIF
mess 'Fitting is done.'

set mtyp 
set pmci
nt/pl //blah/100.amp%tot ch=[pos].and.$EVAL([tdc_mean]-200.)<time<$EVAL([tdc_mean]+250.).and.0<tot<1000.and.0<amp<2000

mess 'Done making another plot'

set pmci 2
set mtyp 8
pmarker [adc_dim] vtot vadc
set mtyp
set pmci

set ltyp
IF ([is_hybrid].eq.0) then
   set plci 2
   fun/pl adc_tot_fit.f 0 2000 s
   fun/pl adc_tot_fit.f 0 2000 s
   fun/pl adc_tot_fit.f 0 2000 s
   fun/pl adc_tot_fit.f 0 2000 s
   fun/pl adc_tot_fit.f 0 2000 s
   set plci 3
   fun/pl adc_tot_fit_s.f 0 2000 s
   fun/pl adc_tot_fit_s.f 0 2000 s
   fun/pl adc_tot_fit_s.f 0 2000 s
   fun/pl adc_tot_fit_s.f 0 2000 s
   fun/pl adc_tot_fit_s.f 0 2000 s
   set plci
   set ltyp 2
   line 125 0 125 2000
ELSEIF ([is_hybrid].eq.1) then
   set plci 2
   fun/pl adc_tot_fit_hyb.f 0 2000 s
   fun/pl adc_tot_fit_hyb.f 0 2000 s
   fun/pl adc_tot_fit_hyb.f 0 2000 s
   fun/pl adc_tot_fit_hyb.f 0 2000 s
   fun/pl adc_tot_fit_hyb.f 0 2000 s
   set plci 3
   fun/pl adc_tot_fit_h_s.f 0 2000 s
   fun/pl adc_tot_fit_h_s.f 0 2000 s
   fun/pl adc_tot_fit_h_s.f 0 2000 s
   fun/pl adc_tot_fit_h_s.f 0 2000 s
   fun/pl adc_tot_fit_h_s.f 0 2000 s
   set plci
   set ltyp 2
   line 200 0 200 2000  
ENDIF
mess 'All those fun/pl commands done'

*pic/print string_[string]/man.adc_tot_fit.string[string]_om[om].gif
pic/print /home/aheffwong/crosstalk/paw_results/man.adc_tot_fit.string[string]_om[om].gif


set plci

* Fill out the results

ve/input res_[string]_om($EVAL($UNQUOTE([pos]))) $EVAL($UNQUOTE([om]))
ve/input res_[string]_a($EVAL($UNQUOTE([pos]))) $EVAL(par(1))
ve/input res_[string]_b($EVAL($UNQUOTE([pos]))) $EVAL(par(2))
ve/input res_[string]_d($EVAL($UNQUOTE([pos]))) $EVAL(par(3))
ve/input res_[string]_e($EVAL($UNQUOTE([pos]))) $EVAL(par(4))
IF ([is_hybrid].eq.0) then
   ve/input res_[string]_c($EVAL($UNQUOTE([pos]))) 1.
ELSEIF  ([is_hybrid].eq.1) then
   ve/input res_[string]_c($EVAL($UNQUOTE([pos]))) 0.
ENDIF
ve/input res_[string]_ll($EVAL($UNQUOTE([pos]))) vtot(1)
ve/input res_[string]_ul($EVAL($UNQUOTE([pos]))) vtot([adc_dim])
ve/print par

* PRINT OUT THE RESULTS

list= vtot,vadc
*fname=string[string].om[om].adc_tot_edge
fname=/home/aheffwong/crosstalk/paw_results/string[string].om[om].adc_tot_edge
form='2(1x,F8.3)'
ve/write [list] [fname] [form]

ve/del vtot
ve/del vadc
ve/del par
close 1

endfor
endfor

for string in  '9'
*'6' '7' '8' '9'  '10'
vaulist = res_[string]_om,res_[string]_a,res_[string]_b,res_[string]_c,res_[string]_d,res_[string]_e,res_[string]_ll,res_[string]_ul

form = '(1x,I3),2(1x,F9.6),(1x,F3.1),2(1x,F9.3),2(1x,F5.1)'
ve/write [vaulist] ' ' [form]

fname=/home/aheffwong/crosstalk/paw_results/man.string[string].adc_tot_fit.result
ve/write [vaulist] [fname] [form]

endfor