macro fit_12c_shape field=para
vec/del *
nseg=1
* 43774 12C 12C para
nrun=43774
* 43775 12C+he para
 if ([field] .eq. 'para') then
nrun=43775
vec/create vthlo(4) r 11.25 12.25 12.8 13.5
vec/create vthhi(4) r 12.25 12.8  13.5 15.0
vec/create vth1lo(2) r 11.25 12.8 
vec/create vth1hi(2) r 12.8  15.0
endif
* 43229 12C+he perp
 if ([field] .eq. 'perp') then
nrun=43229
 vec/create vthlo(4) r 11.5 12.5  13.1 13.8
vec/create vthhi(4) r  12.5 13.1  13.8 16.0
vec/create vth1lo(2) r 11.5 13.1
vec/create vth1hi(2) r  13.1 16.0
endif
vec/create par(4) r
chain -h[field]
fname=/work/hallc/e93026/rss/pass3/ntup/hms[nrun].[nseg].hbook
while ($FEXIST([fname])) do
message add [fname] to chain
chain h[field] [fname]
nseg=[nseg]+1
fname=ntup/hms[nrun].[nseg].hbook
endwhile
zone 1 1
ttit= W  nrun=[nrun]
nbin=80
wlo=.7
whi=1.20
1dhisto 1011 [ttit] [nbin]  [wlo] [whi] 
ntu/proj 1011 //h[field]/9010.w 0.8<hsshtrk/hsp<1.4
*set/norm 1011 1.0
sum=$hinfo(1011,'sum')
hist/plot 1011 e
*vec/create x([nbin]) r
*vec/create y([nbin]) r
*vec/create yerr([nbin]) r
*hist/get_vect/contents 1011 y
*hist/get_vect/abscissa 1011 x
*hist/get_vect/errors 1011 yerr
*vec/create par(2) r 
*vec/fit x y yerr p1 SQ 2 par
hist/fit 1011 p3 SQ 2 par
vec/write par ' ' 4e15.5
sigma par=par/[sum]
vec/write par ' ' 4e15.5
wait
*
zone 2 2
l=1
do i=1,4
thlo=vthlo([i])
*i=[l]+1
*l=[l]+2
thhi=vthhi([i])
ntu/cut $1 [thlo]<hstheta<=[thhi]
message  [thlo]<hstheta<=[thhi]
1dhisto 1012 ' '  [nbin]  [wlo] [whi]
ntu/proj 1012 //h[field]/9010.w $1.and.0.8<hsshtrk/hsp<1.4
hist/plot 1012 e
vec/del par[i](4)
vec/create par[i](4) r
hist/fit 1012 p3 SQ 4 par[i]
vec/write par ' ' 4e15.5
sigma par[i]=par[i]/[sum]
vec/write par[i] ' ' 4e15.5
enddo
zone 2 2
wait
l=1
do i=1,2
thlo=vth1lo([i])
*i=[l]+1
*l=[l]+2
thhi=vth1hi([i])
ntu/cut $1 [thlo]<hstheta<=[thhi]
message  [thlo]<hstheta<=[thhi]
1dhisto 1012 ' '  [nbin]  [wlo] [whi]
ntu/proj 1012 //h[field]/9010.w $1.and.0.8<hsshtrk/hsp<1.4
hist/plot 1012 e
vec/del par[i](4)
vec/create par[i](4) r
hist/fit 1012 p3 SQ 4 par[i]
vec/write par ' ' 4e15.5
sigma par[i]=par[i]/[sum]
vec/write par[i] ' ' 4e15.5
enddo
chain -h[field]