#!/bin/sh
~/bin/saverage.pl sig.inp.1?.1 > Sig.inp.average
W1=$(head -12 Sig.inp.average |tail -1 |strans '#1')
W2=$(head -13 Sig.inp.average |tail -1 |strans '#1')
W3=$(head -14 Sig.inp.average |tail -1 |strans '#1')
S1=$(head -12 Sig.inp.average |tail -1 |strans '#3*(-1)')
S2=$(head -13 Sig.inp.average |tail -1 |strans '#3*(-1)')
S3=$(head -14 Sig.inp.average |tail -1 |strans '#3*(-1)')
#echo $W1 $W2 $S1 $S2
S12=$(calc $S2-$S1)
#echo $S12
W12=$(calc $W2-$W1)
#echo $W12
S23=$(calc $S3-$S2)
W23=$(calc $W3-$W2)
SS=$(calc $S12/$W12 )
SS2=$(calc $S23/$W23 )
m=$(calc 1+$SS)
m2=$(calc 1+$SS2)
dm=$(calc $m-$m2 )
dt=$(calc $W1+$W2 )
dt2=$(calc $dt/2.0 )
dt3=$(calc $W23/$dt2 )
dm2=$(calc $dm/$dt3 )
mav2=$(calc $m+$dm2 )
SS0=$(calc $S1/$W1 )
m0=$(calc 1+$SS0)
mav=$(calc $m*2.0/3.0+$m0/3.0)
S0=$(calc $S1-$m*$W1)
GammaT=$(calc $S0/$m/$W1*3.1415926535)
beta=$(calc 3.1415926535/$W1)
mav3=$(calc $mav2/2.0+$m0/4.0+$m/4.0 )

echo $mav3 $mav2 $mav $m $m0 $S0 $beta $GammaT

W1=$(head -12 Sig.inp.average |tail -1 |strans '#1')
W2=$(head -13 Sig.inp.average |tail -1 |strans '#1')
W3=$(head -14 Sig.inp.average |tail -1 |strans '#1')
S1=$(head -12 Sig.inp.average |tail -1 |strans '#5*(-1)')
S2=$(head -13 Sig.inp.average |tail -1 |strans '#5*(-1)')
S3=$(head -14 Sig.inp.average |tail -1 |strans '#5*(-1)')
#echo $W1 $W2 $S1 $S2
S12=$(calc $S2-$S1)
#echo $S12
W12=$(calc $W2-$W1)
#echo $W12
S23=$(calc $S3-$S2)
W23=$(calc $W3-$W2)
SS=$(calc $S12/$W12 )
SS2=$(calc $S23/$W23 )
m=$(calc 1+$SS)
m2=$(calc 1+$SS2)
dm=$(calc $m-$m2 )
dt=$(calc $W1+$W2 )
dt2=$(calc $dt/2.0 )
dt3=$(calc $W23/$dt2 )
dm2=$(calc $dm/$dt3 )
mav2=$(calc $m+$dm2 )
SS0=$(calc $S1/$W1 )
m0=$(calc 1+$SS0)
mav=$(calc $m*2.0/3.0+$m0/3.0)
S0=$(calc $S1-$m*$W1)
GammaT=$(calc $S0/$m/$W1*3.1415926535)
beta=$(calc 3.1415926535/$W1)
mav3=$(calc $mav2/2.0+$m0/4.0+$m/4.0 )

echo $mav3 $mav2 $mav $m $m0 $S0 $beta $GammaT

W1=$(head -12 Sig.inp.average |tail -1 |strans '#1')
W2=$(head -13 Sig.inp.average |tail -1 |strans '#1')
W3=$(head -14 Sig.inp.average |tail -1 |strans '#1')
S1=$(head -12 Sig.inp.average |tail -1 |strans '#7*(-1)')
S2=$(head -13 Sig.inp.average |tail -1 |strans '#7*(-1)')
S3=$(head -14 Sig.inp.average |tail -1 |strans '#7*(-1)')
#echo $W1 $W2 $S1 $S2
S12=$(calc $S2-$S1)
#echo $S12
W12=$(calc $W2-$W1)
#echo $W12
S23=$(calc $S3-$S2)
W23=$(calc $W3-$W2)
SS=$(calc $S12/$W12 )
SS2=$(calc $S23/$W23 )
m=$(calc 1+$SS)
m2=$(calc 1+$SS2)
dm=$(calc $m-$m2 )
dt=$(calc $W1+$W2 )
dt2=$(calc $dt/2.0 )
dt3=$(calc $W23/$dt2 )
dm2=$(calc $dm/$dt3 )
mav2=$(calc $m+$dm2 )
SS0=$(calc $S1/$W1 )
m0=$(calc 1+$SS0)
mav=$(calc $m*2.0/3.0+$m0/3.0)
S0=$(calc $S1-$m*$W1)
GammaT=$(calc $S0/$m/$W1*3.1415926535)
beta=$(calc 3.1415926535/$W1)
mav3=$(calc $mav2/2.0+$m0/4.0+$m/4.0 )

echo $mav3 $mav2 $mav $m $m0 $S0 $beta $GammaT

W1=$(head -12 Sig.inp.average |tail -1 |strans '#1')
W2=$(head -13 Sig.inp.average |tail -1 |strans '#1')
W3=$(head -14 Sig.inp.average |tail -1 |strans '#1')
S1=$(head -12 Sig.inp.average |tail -1 |strans '#9*(-1)')
S2=$(head -13 Sig.inp.average |tail -1 |strans '#9*(-1)')
S3=$(head -14 Sig.inp.average |tail -1 |strans '#9*(-1)')
#echo $W1 $W2 $S1 $S2
S12=$(calc $S2-$S1)
#echo $S12
W12=$(calc $W2-$W1)
#echo $W12
S23=$(calc $S3-$S2)
W23=$(calc $W3-$W2)
SS=$(calc $S12/$W12 )
SS2=$(calc $S23/$W23 )
m=$(calc 1+$SS)
m2=$(calc 1+$SS2)
dm=$(calc $m-$m2 )
dt=$(calc $W1+$W2 )
dt2=$(calc $dt/2.0 )
dt3=$(calc $W23/$dt2 )
dm2=$(calc $dm/$dt3 )
mav2=$(calc $m+$dm2 )
SS0=$(calc $S1/$W1 )
m0=$(calc 1+$SS0)
mav=$(calc $m*2.0/3.0+$m0/3.0)
S0=$(calc $S1-$m*$W1)
GammaT=$(calc $S0/$m/$W1*3.1415926535)
beta=$(calc 3.1415926535/$W1)
mav3=$(calc $mav2/2.0+$m0/4.0+$m/4.0 )

echo $mav3 $mav2 $mav $m $m0 $S0 $beta $GammaT

W1=$(head -12 Sig.inp.average |tail -1 |strans '#1')
W2=$(head -13 Sig.inp.average |tail -1 |strans '#1')
W3=$(head -14 Sig.inp.average |tail -1 |strans '#1')
S1=$(head -12 Sig.inp.average |tail -1 |strans '#11*(-1)')
S2=$(head -13 Sig.inp.average |tail -1 |strans '#11*(-1)')
S3=$(head -14 Sig.inp.average |tail -1 |strans '#11*(-1)')
#echo $W1 $W2 $S1 $S2
S12=$(calc $S2-$S1)
#echo $S12
W12=$(calc $W2-$W1)
#echo $W12
S23=$(calc $S3-$S2)
W23=$(calc $W3-$W2)
SS=$(calc $S12/$W12 )
SS2=$(calc $S23/$W23 )
m=$(calc 1+$SS)
m2=$(calc 1+$SS2)
dm=$(calc $m-$m2 )
dt=$(calc $W1+$W2 )
dt2=$(calc $dt/2.0 )
dt3=$(calc $W23/$dt2 )
dm2=$(calc $dm/$dt3 )
mav2=$(calc $m+$dm2 )
SS0=$(calc $S1/$W1 )
m0=$(calc 1+$SS0)
mav=$(calc $m*2.0/3.0+$m0/3.0)
S0=$(calc $S1-$m*$W1)
GammaT=$(calc $S0/$m/$W1*3.1415926535)
beta=$(calc 3.1415926535/$W1)
mav3=$(calc $mav2/2.0+$m0/4.0+$m/4.0 )

echo $mav3 $mav2 $mav $m $m0 $S0 $beta $GammaT