set terminal postscript eps color enhanced solid font "Helvetica,22" size 15cm,17cm
set colors classic
set nobars

set output 'rates_efficiency.eps'
set nobars
set border 3
set multiplot

# first plot
set xtics nomirror
set ytics nomirror
set size 1,0.5
set origin 0,0.5
set xlabel "P_{776} (mW)"
set ylabel "pair rate ( r_{p} )"
set  key font ",15" at 16.5,-34000 vertical maxrows 1
set label 1 font "Helvetica,17" "P_{780} :" at first -2.5, first -33000
plot [0:16][:]\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):($2/2.5):(sqrt($2/2.5)) w yerrorbars title " 30 {/Symbol m}W" ps 1.5,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):($2/2.5):(sqrt($2/2.5)) w yerrorbars title "130 {/Symbol m}W" ps 1.5,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):($2/2.5):(sqrt($2/2.5)) w yerrorbars title "290 {/Symbol m}W" ps 1.5,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):($2/2.5):(sqrt($2/2.5)) w yerrorbars title "420 {/Symbol m}W" ps 1.5,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):($2/2.5):(sqrt($2/2.5)) w yerrorbars title "580 {/Symbol u}W" ps 1.5
#,  'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):($2/2.5)title "640 mW" ps 1.5 lc rgb "purple"

# second plot
unset label
set size 0.51,0.45
set origin 0.,0.05
set xtics nomirror
set ytics nomirror
set ytics 0.04
set xlabel "P_{776} (mW)"
set ylabel "signal heralding efficiency"
set key off

#set key bottom
plot [0:16][:]\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):($2/$3) title " 45 {/Symbol m}W" ps 1.5,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):($2/$3) title "170 {/Symbol m}W" ps 1.5,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):($2/$3) title "380 {/Symbol m}W" ps 1.5,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):($2/$3) title "580 {/Symbol m}W" ps 1.5,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):($2/$3) title "770 {/Symbol u}W" ps 1.5
#,  'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):($2/$3) title "830 mW" ps 1.5 lc rgb "purple"

# Third plot
unset label
set size 0.52,0.45
set origin 0.5,0.05
set xtics nomirror
set ytics nomirror
set format x "%2.0f"
set xtics 4
set ytics 0.02
set xlabel "P_{776} (mW)"
set ylabel "idler heralding efficiency"
set key off
unset title

plot [0:16][:]\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):($2/$4) title " 30 {/Symbol m}W" ps 1.5,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):($2/$4) title "130 {/Symbol m}W" ps 1.5,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):($2/$4) title "290 {/Symbol m}W" ps 1.5,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):($2/$4) title "420 {/Symbol m}W" ps 1.5,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):($2/$4) title "580 {/Symbol u}W" ps 1.5
#,  'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):($2/$4) title "640 mW" ps 1.5 lc rgb "purple"

unset multiplot