close 1
opt stat 
opt fit
nrun=13
v/cr run([nrun]) i 2385 2386 2387 2388 2389 2390 2391 2392 2393 2394 2395 2396 2397
datafile=/w/work506/e94104/espace/harp/bpm_mar03_raw_new.dat
hbookfile=/w/work506/e94104/espace/hbook/e94104_harp_mar03_new.hbook
psfile=/w/work506/e94104/espace/harp/bpm_mar03_raw_new.ps
shell rm -f [psfile]
*shell rm -f  [datafile]
hist/file 1 [hbookfile]
for/file 99 [psfile]
meta 99 -111
opt stat
opt fit
opt nbox
set hcol 2 
do i=1,[nrun]
zone 2 4
h/del 999
v/del par       
	irun=$eval(run([i]))
        title runnumber"J#[irun]
	mess [irun]
	nt/plot 1.bpmaxraw n_number=[irun] ! -999  
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
	opt nstat
        h/fit 999 g  ! 3 par 
	axr=$eval(par(2))
	daxr=$eval(par(3))
	
        
	h/del 999 
	v/del par
        opt stat 
	nt/plot 1.bpmayraw n_number=[irun] ! -999 
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
	opt nstat
        h/fit 999 g  ! 3 par 
	ayr=$eval(par(2))
	dayr=$eval(par(3))
        
        opt stat
	h/del 999 
	v/del par 
	nt/plot 1.bpmbxraw n_number=[irun] ! -999
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
	opt nstat
        h/fit 999 g  ! 3 par 
	bxr=$eval(par(2))
	dbxr=$eval(par(3))

        opt stat
	h/del 999 
	v/del par 
	nt/plot 1.bpmbyraw n_number=[irun] ! -999
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
        opt nstat
	h/fit 999 g  ! 3 par 
	byr=$eval(par(2))
	dbyr=$eval(par(3))

        opt stat
	h/del 999 
	v/del par 
	nt/plot 1.beam_x n_number=[irun] ! -999 
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
        opt nstat
	h/fit 999 g  ! 3 par 
	x=$eval(par(2))
	dx=$eval(par(3))

        opt stat
	h/del 999 
	v/del par 
	nt/plot 1.beam_y n_number=[irun] ! -999 
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
        opt nstat
	h/fit 999 g  ! 3 par 
	y=$eval(par(2))
	dy=$eval(par(3))

        opt stat
	h/del 999 
	v/del par 
	nt/plot 1.beam_ph n_number=[irun] ! -999 
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
        opt nstat
	h/fit 999 g  ! 3 par 
	ph=$eval(par(2))
	dph=$eval(par(3))

        opt stat
	h/del 999 
	v/del par 
	nt/plot 1.beam_th n_number=[irun] ! -999 
	v/cr par(3) R $hinfo(999,'max') $hinfo(999,'mean') $hinfo(999,'rms')  
        opt nstat
	h/fit 999 g  ! 3 par 
	th=$eval(par(2))
	dth=$eval(par(3))
    

	shell echo [irun] [axr]	[daxr]	[ayr]	[dayr]	[bxr]	[dbxr]	[byr]	[dbyr]	[x]	[dx]	[y] [dy] [ph] [dph] [th] [dth] >> [datafile]
wait
enddo        
        close 1
        close 99