set title "Material Depolarization Correction for Gaussian/Exponential Model"
set size 1.618,1.
set key spacing 1.5
set ylabel "P_{/Symbol m}^A{/Symbol x}" 
set border 3
set ytics nomirror
set xtics nomirror
set x2tics (1050.,9000.)
set grid x2
set term postscript eps enhanced color "Times-Roman" 21
set out "gauexp.eps"
#set grid y lt 18
set yrange [0.98:1.04]
set ytics ( 0.98,0.99, 1.00, "data" 1.0036,"sim." 1.0106,1.02)
set xlabel 'muon decay time (ns)'
#set key 5200.,0.992
f(x) = 1.0048*exp( -0.0000015496*x )
g(x) = 1.0024*exp( -(x**2/2./(51111.**2)) )
fg(x) = 0.5*f(x) + 0.5*g(x)
f2(x) = 1.0106*exp( -0.0000015496*x )
g2(x) = 1.0106*exp( -(x**2/2./(51111.**2)) )
fg2(x) = 0.5*f2(x) + 0.5*g2(x)

plot [0.:10000.] f(x) title 'data, exponen. P_{/Symbol m}^A{/Symbol x}(0)=1.0048 <P_{/Symbol m}^A{/Symbol x}>=1.0001' lw 3,\
                fg(x) title 'data, exp+Gaus P_{/Symbol m}^A{/Symbol x}(0)=1.0036 <P_{/Symbol m}^A{/Symbol x}>=1.0001' lw 3,\
                 g(x) title 'data, Gaussian P_{/Symbol m}^A{/Symbol x}(0)=1.0024 <P_{/Symbol m}^A{/Symbol x}>=1.0001' lw 3,\
                f2(x) title 'sim., exponen. P_{/Symbol m}^A{/Symbol x}(0)=1.0106 <P_{/Symbol m}^A{/Symbol x}>=1.0059' lw 3,\
               fg2(x) title 'sim., exp+Gaus P_{/Symbol m}^A{/Symbol x}(0)=1.0106 <P_{/Symbol m}^A{/Symbol x}>=1.0071' lw 3,\
                g2(x) title 'sim., Gaussian P_{/Symbol m}^A{/Symbol x}(0)=1.0106 <P_{/Symbol m}^A{/Symbol x}>=1.0082' lt 7 lw 3