set size 1.4,1.0
set out 'mlrs_limits_pmuxi.eps'
set term postscript enhanced eps color "Times-Roman" 16
set grid y2 lt 18
set grid x2 lt 18
set key spacing 1.5
set title 'Manifest LRS Model Limits from {/Symbol D}P_{/Symbol m}{/Symbol x}'
set ylabel '{/Symbol z}'
set xlabel 'm_2 (GeV/c^2)'

set samples 500000
set border 3
set xtics nomirror
set ytics nomirror
mw = 80.390
dpmuxi1 = 0.008
fp1(x) = 0.5*(-(mw/x)**2+sqrt( dpmuxi1-3.*(mw/x)**4 ) ) 
fm1(x) = 0.5*(-(mw/x)**2-sqrt( dpmuxi1-3.*(mw/x)**4 ) ) 
dpmuxi2 = 0.004
fp2(x) = 0.5*(-(mw/x)**2+sqrt( dpmuxi2-3.*(mw/x)**4 ) ) 
fm2(x) = 0.5*(-(mw/x)**2-sqrt( dpmuxi2-3.*(mw/x)**4 ) ) 
dpmuxi3 = 0.001
fp3(x) = 0.5*(-(mw/x)**2+sqrt( dpmuxi3-3.*(mw/x)**4 ) ) 
fm3(x) = 0.5*(-(mw/x)**2-sqrt( dpmuxi3-3.*(mw/x)**4 ) )
dpmuxi4 = 0.0003
fp4(x) = 0.5*(-(mw/x)**2+sqrt( dpmuxi4-3.*(mw/x)**4 ) ) 
fm4(x) = 0.5*(-(mw/x)**2-sqrt( dpmuxi4-3.*(mw/x)**4 ) )
set yrange [-0.06:0.1]
set xtics 300,100.
set ytics -0.1,0.02
set y2tics (-0.052,-0.037,-0.018,-0.010,0.006,0.013,0.029,0.042)
set x2tics (352.,419.,594.,805.)
set label 'Allowed' at 1000.,-0.005
plot [300.:1200.] fp1(x) lt 2 title 'Pre TWIST {/Symbol D}P_{/Symbol m}{/Symbol x} = 0.0080',fm1(x) lt 2 notitle, \
                fp2(x) lt 1 title 'Current {/Symbol D}P_{/Symbol m}{/Symbol x} = 0.0040',fm2(x) lt 1 notitle, \
                fp3(x) lt 3 title 'Next {/Symbol D}P_{/Symbol m}{/Symbol x} = 0.0010',fm3(x) lt 3 notitle, \
                fp4(x) lt 4 title 'TWIST goal {/Symbol D}P_{/Symbol m}{/Symbol x} = 0.0003',fm4(x) lt 4 notitle

set out 'mlrs_limits_pmuxi.ps'
set term postscript enhanced color "Times-Roman" 16
replot