macro beam_comp nrun=1 fit=0
filecase keep
exec utilGenColormap.kumac
close 1
zone 1 1
gsry_adc_off = 4100
gsrx_adc_off = 3700
txt = run = [nrun]
title_gl [txt]
hist/file 1 /work/hallc/e93026/rss/firstpass/ntup_dpara/hms[nrun].1.hbook
if ( [fit] .eq. 1 ) then
ntu/gcut $1 9010.ntbpmx%calsrx
opt stat
zone 1 1
2dhist 1000 ' ' 100 -4000 4000 100 -20 20
ntu/proj 1000 9010.ntbpmx%gsry_adc-[gsry_adc_off] $1
hist/plot 1000 colz
message CLick at lower y and upper y limit
vlocate x y +ls
bpmstart=y(1)
bpmend=y(2)
bpmnum=20.
vec/create bpmc([bpmnum]) r
vec/create adcc([bpmnum]) r
bpmstep=([bpmend]-[bpmstart])/[bpmnum] 
do i=1,[bpmnum]
  bpmlo = [bpmstart] + ([i]-1)*[bpmstep]
  bpmhi = [bpmlo] + [bpmstep]
  message [bpmlo] [bpmhi] [bpmstep]
  1dhist 100 ' ' 100 -20 20
  ntu/proj 100 9010.ntbpmx $1.and.[bpmlo]<ntbpmx<[bpmhi]
  hrin 100
  mean=$hinfo(100,'mean')
  vec/input bpmc([i]) [mean] 
  1dhist 100 ' ' 100 -4000 4000
  ntu/proj 100 9010.gsry_adc-[gsry_adc_off] $1.and.[bpmlo]<ntbpmx<[bpmhi]
  hrin 100
  mean=$hinfo(100,'mean')
  vec/input adcc([i]) [mean] 
enddo
hist/plot 1000 colz
vec/create parx(2) r
set fcol 2.
set fwid 2.
vec/fit adcc bpmc ? p1 s 2 parx 
wait
endif
*
if ( [fit] .eq. 1 ) then
ntu/gcut $1 9010.ntbpmy%calsry
opt stat
zone 1 1
2dhist 1000 ' ' 100 -4000 4000 100 -20 20
ntu/proj 1000 9010.ntbpmy%gsrx_adc-[gsrx_adc_off] $1
hist/plot 1000 colz
message CLick at lower y and upper y limit
vlocate x y +ls
bpmstart=y(1)
bpmend=y(2)
bpmnum=20.
vec/create bpmc([bpmnum]) r
vec/create adcc([bpmnum]) r
bpmstep=([bpmend]-[bpmstart])/[bpmnum] 
do i=1,[bpmnum]
  bpmlo = [bpmstart] + ([i]-1)*[bpmstep]
  bpmhi = [bpmlo] + [bpmstep]
  message [bpmlo] [bpmhi] [bpmstep]
  1dhist 100 ' ' 100 -20 20
  ntu/proj 100 9010.ntbpmy $1.and.[bpmlo]<ntbpmy<[bpmhi]
  hrin 100
  mean=$hinfo(100,'mean')
  vec/input bpmc([i]) [mean] 
  1dhist 100 ' ' 100 -4000 4000
  ntu/proj 100 9010.gsrx_adc-[gsrx_adc_off] $1.and.[bpmlo]<ntbpmy<[bpmhi]
  hrin 100
  mean=$hinfo(100,'mean')
  vec/input adcc([i]) [mean] 
enddo
hist/plot 1000 colz
vec/create pary(2) r
vec/fit adcc bpmc ? p1 s 2 pary
vec/write parx,pary ' ' 2f15.5 
n_sr_adcx_zero=[gsry_adc_off]
n_sr_adcy_zero=[gsrx_adc_off]
n_sr_slopex=parx(2)
n_sr_slopey=pary(2)
n_sr_offsetx=parx(1)
n_sr_offsety=pary(1)
message n_sr_adcx_zero=[n_sr_adcx_zero] n_sr_slopex=[n_sr_slopex] n_sr_offsetx=[n_sr_offsetx]
message n_sr_adcy_zero=[n_sr_adcy_zero] n_sr_slopey=[n_sr_slopey] n_sr_offsety=[n_sr_offsety]
wait
endif
zone 2 2
* 
txt = run = [nrun]
title_gl [txt]
zone 2 2
1dhisto 1000 ' ntbpmx' 100 -2.0 2.0
ntu/proj 1000 9010.ntbpmx/10.
hist/plot 1000
set hcol 2
1dhisto 1001 ' calsrx' 100 -2.0 2.0
ntu/proj 1001 9010.calsrx 
hist/plot 1001 s
set hcol 1
1dhisto 1002 ' ntbpmx -calsrx' 100 -.5 .5
ntu/proj 1002 9010.ntbpmx/10.-calsrx
hist/plot 1002
2dhisto 1003 ' calsrx vs. ntbpmx-calsrx' 100 -.5 .5 100 -2. 2.
ntu/proj 1003 9010.calsrx%ntbpmx/10.-calsrx
hist/plot 1003 colz
wait
zone 2 2
1dhisto 1000 ' ntbpmy' 100 -2.0 2.0
ntu/proj 1000 9010.ntbpmy/10.
hist/plot 1000
set hcol 2
1dhisto 1001 ' calsry' 100 -2.0 2.0
ntu/proj 1001 9010.calsry 
hist/plot 1001 s
set hcol 1
1dhisto 1002 ' ntbpmy -calsry' 100 -5. 5.
ntu/proj 1002 9010.ntbpmy/10.-calsry
hist/plot 1002
2dhisto 1003 ' calsry vs. ntbpmy-calsry' 100 -.5 .5 100 -2. 2.
ntu/proj 1003 9010.calsry%ntbpmy/10.-calsry
hist/plot 1003 colz
wait