#!/bin/tcsh -f

goto get_vals
echo -n  "11111 or 12660  --> [11111] " 
set wave = ( $< ) 

if (("$wave" != "") && ("$wave" !=  "11111") )then
echo "using 12660 ev"
   set a = -0.07191
   set m = 0.02031
   goto  get_percent
endif 
echo "using 11111 ev"
set m               = 0.0205
set a               = -0.053  

get_vals:
set m               = 0.0202366    
set a               = -0.0522611  



get_percent:

set percent = $1
if ("$percent" != "") goto calculate
echo  "    based on full beam with slits wide open, what percent of the flux would you like:" 
echo -n "->"

set percent = ( $< ) 

calculate:
set hor_div = `echo $percent | awk -v m=$m -v a=$a '{printf"%.2f\n", $1 * m + a }'`
echo " Get approx. ${percent}%  of full beam (2.0 X .35 divergence) by doing this-->>    divergence.com $hor_div .35"

exit




set divergence = `divergence.com | awk '{printf "%.2f %.2f\n", $3, $4}'`
set m    = 0.0202366
set a = -0.0522611
set percent = echo $divergence[1] | awk -v m=$m -v a=$a '{printf"%.2f\n", ($1 -a) / m }'`
echo "a: $a ; m: $m divergence: $divergence[1] ; percent: $percent"

diode.com in
shutter.com open
set divergence = "2.00 0.35"
divergence.com  $divergence
sleep 1
set FLUX = ` flux.com | awk '/Flux/{print $3}'`

echo " flux with full beam divergence =  $divergence is $FLUX nAmp"

  divergence.com $divergence >>& /dev/null
   set divergence = `divergence.com | awk '{printf "%.2f %.2f\n", $3, $4}'`
   set foo = `flux.com | awk '{if ($1 == "Flux") print $3, $4}'`
   set flux = $foo[1]
   set percent = `echo $flux $FLUX | awk '{print $1/$2 * 100}'`
   echo "divergence: $divergence, flux: $foo,  full beam percentage: $percent"





 fit f(x) "div35percent.log" using 1:2 via m,a

set m               = 0.0202366    
set a               = -0.0522611  
set percent = 50

set hor_div = `echo $percent | awk -v m=$m -v a=$a '{printf"%.2f\n", $1 * m + a }'`
echo " for $percent percent attenuation cut and paste the following:"
echo "  divergence.com $hor_div .35"
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$




cat atten_12660.txt | awk '/diver/{print $NF, $2}' > ! div12660ev.log

gnuplot> f(x) = m * x + a
 fit f(x) "div12660ev.log" using 1:2 via m,a

 gnuplot> print a
-0.0719177754652525
gnuplot> print m
0.0203150026785394

plot f(x), "div12660ev.log"
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
cat atten_12660.txt
divergence: 2.00 0.35, flux: 224.75 uA,  full beam percentage: 99.6541
divergence: 1.75 0.35, flux: 200.44 uA,  full beam percentage: 88.8751
divergence: 1.50 0.35, flux: 174.48 uA,  full beam percentage: 77.3644
divergence: 1.25 0.35, flux: 149.16 uA,  full beam percentage: 66.1375
divergence: 1.00 0.35, flux: 123.30 uA,  full beam percentage: 54.6712
divergence: 0.80 0.35, flux: 101.29 uA,  full beam percentage: 44.912
divergence: 0.70 0.35, flux: 89.67 uA,  full beam percentage: 39.7597
divergence: 0.60 0.35, flux: 77.39 uA,  full beam percentage: 34.3147
divergence: 0.50 0.35, flux: 64.15 uA,  full beam percentage: 28.4441
divergence: 0.40 0.35, flux: 51.27 uA,  full beam percentage: 22.7331
divergence: 0.35 0.35, flux: 45.24 uA,  full beam percentage: 20.0594
divergence: 0.25 0.35, flux: 31.85 uA,  full beam percentage: 14.1223
divergence: 0.25 0.25, flux: 30.84 uA,  full beam percentage: 13.6745






$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$4
    set info = `tac /data/log/change.log | awk '/Iend/{print;exit}'`
    set cpna = 3000000
    set Iend  = `echo "$info" | awk '{print $23/1000}'`
    set flux = `echo $Iend $cpna | awk '{printf "%.4g", $1*$2*1000}'`

    set energy = "$info[6]"
    set pinhole = 100

    set doserate = `echo $flux $pinhole $energy | awk '{printf "%.4g", $1/(3.14*($2/2)**2)/(2000/(12398.42/$3))}'`  

    echo "Iend flux energy pinhole doserate  $Iend $flux $energy $pinhole $doserate"


                                                                                                                                                                      
../enumerate_fields.com $info
The number of arguments in argv is: 23
Argument 1  is Aug
Argument 2  is 8
Argument 3  is 15:33:49
Argument 4  is 2016
Argument 5  is energy:
Argument 6  is 12660.0
Argument 7  is Hdiv:
Argument 8  is 0.3
Argument 9  is Vdiv:
Argument 10  is 0.3
Argument 11  is Iring:
Argument 12  is 500.327
Argument 13  is Iin:
Argument 14  is 2.2792
Argument 15  is Iout:
Argument 16  is 0.8379
Argument 17  is Izero:
Argument 18  is 0.5270
Argument 19  is blBL-831
Argument 20  is Chi2
Argument 21  is 1.28043
Argument 22  is Iend:
Argument 23  is 38016.5