#!/bin/bash

gnuplot<<EOF
set terminal postscript enhanced color landscape
set title "H_2 + Cu(111) \n Ratios of rotationally inelastic scattering probabilities \n P({/Symbol n} = 1, J = 2) / P({/Symbol n} = 1, J = 0)" font "Helvetica,18"
set ylabel "P({/Symbol n} = 1, J = 2) / P({/Symbol n} = 1, J = 0)" font "Helvetica, 14"
set xlabel "translational energy [eV]" font "Helvetica, 14"
set output "test.ps"
set key left font "Helvetica, 14"
set tics font "Helvetica, 14"
set xrange[0:0.3]
#plot "experiment" u (\$1 / 100):2 w lp lw 2 title"exp." , "state-to-state_v1j0-0-2_WP_B" u 1:(\$3/\$2) w l lw 3 title "QD B86SRP68-DF2",  "" u (\$1 - 0.04):(\$3/\$2) w l lw 3 title"QD B86SRP68-DF2 - 0.04 eV)" , "SRP48" u (\$1 / 100):2 w l lw 3 title"QD SRP48", "SRP48" u ((\$1 / 100)-0.04):2 w l lw 3 title"QD SRP48 - 0.04 eV", "v1j0mj0-to-0_2_WP_B_PBEa57-DF2" u 1:(\$3/\$2) w l lw 3 title "QD PBE{/Symbol a}57-DF2"

#plot "experiment" u (\$1 / 100):2 w lp lw 2 title"exp." , "state-to-state_v1j0-0-2_WP_B" u 1:(\$3/\$2) w l lw 3 title "QD B86SRP68-DF2", "SRP48" u (\$1 / 100):2 w l lw 3 title"QD SRP48", "v1j0mj0-to-0_2_WP_B_PBEa57-DF2" u 1:(\$3/\$2) w l lw 3 title "QD PBE{/Symbol a}57-DF2"
#plot "experiment" u (\$1 / 100):2 w lp lw 2 title"exp." , "state-to-state_v1j0-0-2_WP_B" u (\$1 - 0.04):(\$3/\$2) w l lw 3 title"QD B86SRP68-DF2 - 0.04 eV)", "SRP48" u ((\$1 / 100)-0.04):2 w l lw 3 title"QD SRP48 - 0.04 eV" 

plot "experiment" u (\$1 / 100):2 w lp lw 2 title"exp." , "state-to-state_v1j0-0-2_WP_B" u (\$1 ):(\$3/\$2) w l lw 3 lt 2 title"QD B86SRP68-DF2 )", "state-to-state_v1j0-0-2_WP_C" u (\$1 ):(\$3/\$2) w l lw 3 lt 2 title"", "SRP48" u (\$1 / 100):2 w l lw 3 title"QD SRP48"


EOF
ps2pdf test.ps
rm test.ps
evince test.pdf