set term png enhanced transparent font "Helvetica" 15
set encoding utf8
set auto
set border 15 lw 0.5
set grid

###################################################################
#	Scalar Amplitudes
###################################################################

#set xlabel 't (GeV^2)'
#set ylabel '(micro b/GeV^2)'
set output 'plotAp.png'
plot 	'output0.txt' using 1:3 w l ls 2 lw 3 lc rgb "blue" title "Re A^{(+)}", \
	'output0.txt' using 1:4 w l ls 2 lw 3 lc rgb "red" title "Im A^{(+)}"

#set xlabel 't (GeV^2)'
#set ylabel 'Dsig/Dt (micro b/GeV^2)'
set output 'plotBp.png'
plot 	'output0.txt' using 1:5 w l ls 2 lw 3 lc rgb "blue" title "Re B^{(+)}", \
	'output0.txt' using 1:8 w l ls 2 lw 3 lc rgb "red" title "Im B^{(+)}"

#set xlabel 't (GeV^2)'
#set ylabel 'Dsig/Dt (micro b/GeV^2)'
set output 'plotCp.png'
plot 	'output0.txt' using 1:7 w l ls 2 lw 3 lc rgb "blue" title "Re C^{(+)}", \
	'output0.txt' using 1:8 w l ls 2 lw 3 lc rgb "red" title "Im C^{(+)}"


#set xlabel 't (GeV^2)'
#set ylabel '(micro b/GeV^2)'
set output 'plotAm.png'
plot 	'output1.txt' using 1:3 w l ls 2 lw 3 lc rgb "blue" title "Re A^{(-)}", \
	'output1.txt' using 1:4 w l ls 2 lw 3 lc rgb "red" title "Im A^{(-)}"

#set xlabel 't (GeV^2)'
#set ylabel 'Dsig/Dt (micro b/GeV^2)'
set output 'plotBm.png'
plot 	'output1.txt' using 1:5 w l ls 2 lw 3 lc rgb "blue" title "Re B^{(-)}", \
	'output1.txt' using 1:8 w l ls 2 lw 3 lc rgb "red" title "Im B^{(-)}"

#set xlabel 't (GeV^2)'
#set ylabel 'Dsig/Dt (micro b/GeV^2)'
set output 'plotCm.png'
plot 	'output1.txt' using 1:7 w l ls 2 lw 3 lc rgb "blue" title "Re C^{(-)}", \
	'output1.txt' using 1:8 w l ls 2 lw 3 lc rgb "red" title "Im C^{(-)}"



###################################################################
#	Observables
###################################################################

set xlabel 'plab (GeV)'
set ylabel 'milli b'
set output 'plotSigTot.png'
plot 	'SigTot.txt' using 2:4 w l ls 2 lw 3 lc rgb "blue" title "\sigma(\pi^- p)", \
	'SigTot.txt' using 2:5 w l ls 2 lw 3 lc rgb "red" title "\sigma(\pi^+ p)"


set log y
set xlabel 't (GeV^2)'
set ylabel 'micro b.GeV^{-2}'
set output 'plotDSigDt.png'
plot 	'Observables0.txt' using 2:3 w l ls 2 lw 3 lc rgb "blue" title "Dsigma/Dt(\pi^- p -> pi^- p)", \
	'Observables1.txt' using 2:3 w l ls 2 lw 3 lc rgb "red" title "Dsigma/Dt(\pi^+ p -> pi^+ p)", \
	'Observables2.txt' using 2:3 w l ls 2 lw 3 lc rgb "green" title "Dsigma/Dt(\pi^-p —> pi^0 n)"

unset log y
set yrange [-1:1]
set xlabel 't (GeV^2)'
set output 'plotPol.png'
plot 	'Observables0.txt' using 2:4 w l ls 2 lw 3 lc rgb "blue" title "P(\pi^- p -> pi^- p)", \
	'Observables1.txt' using 2:4 w l ls 2 lw 3 lc rgb "red" title "P(\pi^+ p -> pi^+ p)", \
	'Observables2.txt' using 2:4 w l ls 2 lw 3 lc rgb "green" title "P(\pi^-p —> pi^0 n)"