set term pdf fontscale 0.5 size 6in,4in

file1='tbt_9.1584e-06_bl_bz_0._br_bz_2._0_-7000000._eta_0._radial_frac_2.'
outfile = 'spin_vs_fields.pdf'
set output outfile

set timestamp 
set label 1 '$d9/g-2/examples/example10/edm/plotting_scripts/spin_vs_fields.gnu' at graph 0.6 ,-0.12 font 'Verdana,6' noenhance
!pwd > directory
!sed -n -e 's/^/name="/' -e 's/$/"/p' directory > name
load 'name'
print 'directory='.name
set label 2 at graph 0.0,1.02 name font 'Verdana,10' noenhance
set title ' '
set xrange [0:10]
set xlabel 'time [us]'
set ylabel 's_{long}'
set y2label 's_{vert}'
set yrange [-1:1]
set ytics -1,.2
set  y2range [-1.:1.]
set y2tics -1.,0.2
set label 3 at graph 0.05,0.95 '{/Symbol h}=10^{-6}'
set label 5 at graph 0.05,0.9 'd = 4.67 X 10^{-20} e-cm'
factory2=1000.
print log10(factory2)
#set label 1 at graph 1, 0 sprintf("× 10^{%d}",log10(factorX)) offset 0.5,0.3
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./10000.)) offset 0,0.7
plot 'tbt_0.0_bl_bz_0._br_bz_0._0_-7000000._eta_1.e-6' u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'

factory2=10.
set label 3 at graph 0.05,0.95 '{/Symbol h}=1.158 X 10^{-4}'
set label 5 at graph 0.05,0.9 'd = 5.4 X 10^{-18} e-cm'
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.99999999)) offset 0,0.7
plot 'tbt_0.0_bl_bz_0._br_bz_0._0_-7000000._eta_1.15848e-4_radial_frac_1.'  u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'

set xrange [0:10]
factory2=100.
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set label 3 at graph 0.05,0.95 '{/Symbol q}_{pitch}=1 mrad'
unset label 5
plot 'tbt_0.001_bl_bz_0._br_bz_0._0_-7000000._eta_0._radial_frac_1.' u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'
factory2=10000.
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set label 6 at graph 0.05,0.85 'Filtered'
plot 'tbt_0.001_bl_bz_0._br_bz_0._0_-7000000._eta_0._radial_frac_1._Filt' u 2:($6*factory2) skip 1 w l axes x1y2 t 's_{vert}', '' u 2:5 skip 1 w l dt 2 t 's_{long}''

factory2=10000.
unset label 6
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set label 3 at graph 0.05,0.95 'B_{long}=10^{-6} B_y'
plot 'tbt_0.0_bl_bz_1.e-6_br_bz_0._0_-7000000._eta_0.' u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'

factory2=100000.
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set label 3 at graph 0.05,0.95 'B_{radial}=10^{-6} B_y'
set label 5 at graph 0.05,0.9 'B_{radial} uniform around ring'
plot 'tbt_0.0_bl_bz_0._br_bz_1.e-6_0_-7000000._eta_0._radial_frac_1.' u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'
plot 'tbt_0.0_bl_bz_0._br_bz_1.e-6_0_-7000000._eta_0._radial_frac_1._Filt' u 2:($6*factory2) w l axes x1y2 t 's_{vert}', '' u 2:5 w l dt 2 t 's_{long}'

#set xrange [0:5]
#plot 'tbt_0.0_bl_bz_0._br_bz_1.e-6_0_-7000000._eta_0._radial_frac_1.' u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'

set label 5 at graph 0.05,0.9 'B_{radial} limited to 72 deg arc'
factory2=100000.
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set xrange [0:20]
plot 'tbt_0.0_bl_bz_0._br_bz_1.e-6_0_-7000000._eta_0._radial_frac_0.2' u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'
set xrange [0:5]
plot 'tbt_0.0_bl_bz_0._br_bz_1.e-6_0_-7000000._eta_0._radial_frac_0.2' u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'

factory2=10000.
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set label 3 at graph 0.05,0.95 'B_{radial}=4.3 X 10^{-5} B_y'
unset label 5
set xrange [0:20]
plot 'tbt_0.0_bl_bz_0._br_bz_43.e-6_0_-7000000._eta_0._radial_frac_1.' u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'
reset
set xlabel 'time [us]'
set xrange [0:20]
set ylabel 'vertical displacement [mm]'
set yrange [-20:20]
set label 3 at graph 0.05,0.95 'B_{radial}=4.3 X 10^{-5} B_y'
plot 'tbt_0.0_bl_bz_0._br_bz_43.e-6_0_-7000000._eta_0._radial_frac_1.' u ($14*1.e6):($6*1000) w l not




set xlabel 'time [us]'
set ylabel 's_{long}'
set y2label 's_{vert}'
set yrange [-1:1]
set ytics -1,.2
set  y2range [-1.:1.]
set y2tics -1.,0.2

factory2=10000.
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set label 3 at graph 0.05,0.95 'B_{radial}=10^{-5} B_y'
set label 5 at graph 0.05,0.9 'B_{radial} limited to 72 deg arc'
plot 'tbt_0.0_bl_bz_0._br_bz_1.e-5_0_-7000000._eta_0._radial_frac_0.2' u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'

reset
set xlabel 'time [us]'
set xrange [0:10]
set ylabel 'vertical displacement [mm]'
set autoscale y
set label 3 at graph 0.05,0.95 '{/Symbol q}_{pitch}=1 mrad'
plot 'tbt_0.001_bl_bz_0._br_bz_0._0_-7000000._eta_0.' u ($14*1.e6):($6*1000) w l not



set xlabel 'time [us]'
set ylabel 's_{long}'
set y2label 's_{vert}'
set yrange [-1:1]
set ytics -1,.2
set  y2range [-1.:1.]
set y2tics -1.,0.2

factory2=10000.
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.999999)) offset 0,0.7
set label 3 at graph 0.05,0.95 'B_{radial}=10^{-5} B_y'
set label 5 at graph 0.05,0.9 'B_{radial} limited to 72 deg arc'
plot 'tbt_0.0_bl_bz_0._br_bz_1.e-5_0_-7000000._eta_0._radial_frac_0.2' u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'

factory2=1000.
set label 3 at graph 0.05,0.95 '{/Symbol h}=1.502 X 10^{-6}'
set label 5 at graph 0.05,0.9 'd = 7 X 10^{-20} e-cm'
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.99999999)) offset 0,0.7
plot 'tbt_0._bl_bz_0._br_bz_0._0_-7000000._eta_1.502e-6_radial_frac_1.' u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'

factory2=10000.
set label 3 at graph 0.05,0.95 'B_{radial}=22 X 10^{-6} B_{vert}'
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.99999999)) offset 0,0.7
set label 5 at graph 0.05,0.9 'Oscillating about closed orbit'
plot 'tbt_0.0_bl_bz_0._br_bz_22.e-6_0_-7000000._eta_0._radial_frac_1._voff_0.'  u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'
plot 'tbt_0.0_bl_bz_0._br_bz_22.e-6_0_-7000000._eta_0._radial_frac_1._voff_0._Filt'  u 2:($6*factory2) w l axes x1y2 t 's_{vert}', '' u 2:5 w l dt 2 t 's_{long}'

set label 5 at graph 0.05,0.9 'Vertical closed orbit +1.37 mm'
plot 'tbt_0._bl_bz_0._br_bz_22.e-6_0_-7000000._eta_0._radial_frac_1._voff_1.367'  u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'

factory2=10000.
set label 3 at graph 0.05,0.95 'B_{radial}= measured'
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.99999999)) offset 0,0.7
set label 5 at graph 0.05,0.9 'On closed orbit'
plot 'tbt_9.1584E-06_bl_bz_0._br_bz_1._0_-7000000._eta_0._radial_frac_1.'  u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'
plot 'tbt_9.1584E-06_bl_bz_0._br_bz_1._0_-7000000._eta_0._radial_frac_1._Filt' u 2:($6*factory2) w l axes x1y2 t 's_{vert}', '' u 2:5 w l dt 2 t 's_{long}'


set label 3 at graph 0.05,0.95 'B_{radial}= B_{measured} - <B_{measured}>'
set label 5 at graph 0.05,0.9 'On closed orbit -0.259 mm'
plot file1  u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'
plot file1.'_Filt' u 2:($6*factory2) w l axes x1y2 t 's_{vert}', '' u 2:5 w l dt 2 t 's_{long}'



#plot 'tbt_9.1584E-06_bl_bz_0._br_bz_1._0_-7000000._eta_0._radial_frac_1.'   u ($14*1.e6):($2 == 17 ? $11*factory2:1/0) w lp t '0 deg',\
# ''  u ($14*1.e6):($2 == 1 ? $11*factory2:1/0) w lp t '22.5 deg', ''  u ($14*1.e6):($2 == 3 ? $11*factory2:1/0) w lp t '45 deg',\
# ''  u ($14*1.e6):($2 == 5 ? $11*factory2:1/0) w lp t '90 deg' ,\
# ''  u ($14*1.e6):($2 == 7 ? $11*factory2:1/0) w lp t '135 deg', ''  u ($14*1.e6):($2 == 9 ? $11*factory2:1/0) w lp t '180 deg', ''  u ($14*1.e6):($2 == 11 ? $11*factory2:1/0) w lp t '225 deg'

#plot 'tbt_9.1584E-06_bl_bz_0._br_bz_1._0_-7000000._eta_0._radial_frac_1.'   u ($14*1.e6):($2 == 17 ? $11*factory2:1/0) w lp t '0 deg', '' u ($14*1.e6):($2 == 17 ? $12:1/0) w lp t '0 deg',\
# ''  u ($14*1.e6):($2 == 1 ? $11*factory2:1/0) w lp t '22.5 deg', ''  u ($14*1.e6):($2 == 1 ? $12:1/0) w lp t '22.5 deg', ''  u ($14*1.e6):($2 == 3 ? $11*factory2:1/0) w lp t '45 deg',\
# ''  u ($14*1.e6):($2 == 3 ? $12:1/0) w lp t '45 deg'

reset
set label 3 at graph 0.05,0.95 'B_{radial}= measured'
set xtics 0,0.015
set xrange [0:0.15]
set yrange [3:5]
set ytics 0,0.2
set xlabel 'time [us]'
set ylabel 'vertical displacement [mm]'
set label 5 at graph 0.05,0.9 'Closed orbit'
plot 'tbt_9.1584E-06_bl_bz_0._br_bz_1._0_-7000000._eta_0._radial_frac_1.'  u ($14*1.e6):($6*1000) w l not
set label 3 at graph 0.05,0.95 'B_{radial}= B_{measured} - <B_{measured}>'
set yrange [-1:1]
set ytics -1,0.2

plot file1  u ($14*1.e6):($6*1000) w l not

reset
set label 3 at graph 0.05,0.95 'B_{radial}= measured'
set xrange [0:10]
set yrange [3:5]
set ytics 0,0.2
set xlabel 'time [us]'
set ylabel 'vertical displacement [mm]'
set label 5 at graph 0.05,0.9 'Closed orbit'
plot 'tbt_9.1584E-06_bl_bz_0._br_bz_1._0_-7000000._eta_0._radial_frac_1.'  u ($14*1.e6):($6*1000) w l not

set label 3 at graph 0.05,0.95 'B_{radial}= B_{measured} - <B_{measured}>'
set yrange [-1:1]
set ytics -1, 0.2

plot file1  u ($14*1.e6):($6*1000) w l not

reset
set xrange [0:360]
set xtics 0,60
set title 'B_{radial}'
set xlabel 'Azimuth [deg]'
set ylabel 'B_{rad} [ppm]'
plot 'fort.99' u ($1/4.4686E+01*360):($2*1.e6) i 0:0 w p pt 6 ps .5 not

print outfile