set terminal post eps enhanced solid color "Times-Roman" 20
set nokey

set style line 1 lt 1 lw 4 pt 5  ps 1.5  
set style line 2 lt 2 lw 4 pt 7  ps 1.5 
set style line 3 lt 3 pt 9  ps 1.2
set style line 4 lt 4 pt 13 ps 1.2

# set grid 5 5 

rhoaup = 0.01487
rhobup = -0.00927
rhoadn = 0.3459
rhobdn = -0.02506

delaup = 0.09159
delbup = -0.06224
deladn = -0.03697
delbdn = -0.03401

xiaup = -0.1053
xibup = 0.0782
xiadn = 0.00284
xibdn = -0.05566

set output "aup_dn_all.eps"

set size 2,3
set origin 0,0
set multiplot

set size 2, 1
set origin 0, 2
set ylabel "/Symbol{r} x 10^4"
set xlabel "Data Set Number"
set xrange [66:95]
set yrange [-10:10]

plot "ecal_par_bias.txt" us 1:(rhoaup*$2 + rhobup*$4 + rhoadn*$6 + rhobdn*$8)*10000:\
     sqrt((rhoaup*$3)^2 + (rhobup*$5)^2 + (rhoadn*$7)^2 + (rhobdn*$9)^2)*10000 \
     linestyle 1 w errorbars 


set size 2, 1
set origin 0, 1
set ylabel "/Symbol{d} x 10^4"
set xlabel "Data Set Number"
set xrange [66:95]
set yrange [-10:10]

plot "ecal_par_bias.txt" us 1:(delaup*$2 + delbup*$4 + deladn*$6 + delbdn*$8)*10000:\
     sqrt((delaup*$3)^2 + (delbup*$5)^2 + (deladn*$7)^2 + (delbdn*$9)^2)*10000\
     linestyle 1