set terminal postscript eps color enhanced font "Helvetica,14pt" size 9cm,6cm 

#set terminal pngcairo enhanced font "Helvetica,42pt" size 4000,3000

set output "comparison.eps"
#set output "comparison.png"

# common settings
set logscale xy
set format xy "10^{%L}"
set tics scale 2
set tics nomirror
set xr [0.00000001:500] #nanometres
set yr [100:1E26] #counts/sec
set border linewidth 1

set logscale x2 10
set format x2 "10^{%L}"
set x2tics nomirror
set x2r [0.0001:0.000000000000002]
set x2label "Coherence time (s)"

set xlabel "Spectral width {/Symbol Dl} (nm)"
set ylabel "Photoevents per second per mode" offset 0.8,0


#this work
set label "This work" at 0.0000005,1E18

#previous work - sub-threshold laser
set label "STL^{12}" at 2E-5,8E14

#thermal sources
set label "MS^{13}" at 4E-8,3E7
set label "^1^4Sunlight" at 0.00002,800000
set label "Hg^1^5" at 0.0015,9000
set label "RGG^1^6" at 0.02,8000
set label "EDFA^1^7" at 1,1E16
set label "SLD^1^8" at 8,3E19

#spdc sources
set label "MRR^{19}" at 4E-4,4E8
set label "^2^0FWM" at 0.09,600000
set label "PPLN^2^1" at 8,60000000
set label "BBO^2^2" at 15,1200000
set label "PPKTP^2^3" at 3,1400


M(x)=4E21*x
set label "4{/Symbol \264}10^2^1 photoevents/sec/nm" at 0.00005,5E18 rotate by 15

N(x)=40000000.*x
set label "4{/Symbol \264}10^7 photoevents/sec/nm" at 0.005,5000000. rotate by 15

#H(x)=1000000000000.
#set arrow from 0.07, graph 0 to 0.07, graph 1 nohead lw 2 lc "black" dt 2 


set label "Thermal" at 0.0000005,2E23
set label "SPDC" at 0.0000005,5E21

plot 'thermal.txt' using 3:2 notitle "Thermal" with points ps 1.5 lw 2 pt 7 lc rgb "red", 'spdc.txt' using 3:2 notitle "SPDC" with points ps 1.5 lw 2 pt 5 lc rgb "blue", M(x) lw 2 lc "black" notitle, N(x) lw 2 lc "black" notitle

#, H(x) lw 2 lc "black" dt 2 notitle

reset