macro load_vecs dist=2540 nomatter=5e+09 matter=4e+09
   ve/del *
   ve/read cpm,eem,aam cpnue-prob-matter-[dist]/cpnue-[nomatter].vec
   ve/read cpn,een,aan cpnue-prob-nomatter-[dist]/cpnue-[matter].vec
   ve/read en,an,bn,cn prob/nomatter_numu_[dist]_0_0.8_1.0_0.06_6.0_2.5.vec
   ve/read em,am,bm,cm prob/matter_numu_[dist]_0_0.8_1.0_0.06_6.0_2.5.vec
 return
 
 macro plot_nomatter_vs_energy dist=2540
   sigma gev=en/1e9
   title 'No matter, CP=0, '//[dist]//' km'
   ve/pl an%gev
   atitle 'E?[n]! (GeV)' 'Prob([n]?[m]! "5# [n]?e!)' ' ' 321
 return
 macro plot_matter_vs_energy dist=2540
   sigma gev=en/1e9
   title 'Matter Effects, CP=0, '//[dist]//' km'
   ve/pl am%gev
   atitle 'E?[n]! (GeV)' 'Prob([n]?[m]! "5# [n]?e!)' ' ' 321
 return
 macro plot_nomatter_vs_cp energy dist=2540
   title 'No matter, E?[n]! = '//[energy]//' eV, '//[dist]//' km'
   ve/pl aan%cpn
   atitle 'CP phase (degree)' 'Prob([n]?[m]! "5# [n]?e!)' ' ' 321
 return
 macro plot_matter_vs_cp energy dist=2540 
   title 'Matter Effects, E?[n]! = '//[energy]//' eV, '//[dist]//' km'
   ve/pl aam%cpm
   atitle 'CP phase (degree)' 'Prob([n]?[m]! "5# [n]?e!)' ' ' 321
 return
 
 macro plot_dist dist=2540 nomatter=5e+09 matter=4e+09
   exec load_vecs [dist] [nomatter] [matter]
   exec cpprob#plot_nomatter_vs_energy  [dist]
   exec text#footer
   exec cpprob#plot_matter_vs_energy  [dist]
   exec text#footer
   exec cpprob#plot_nomatter_vs_cp [nomatter] [dist]
   exec text#footer
   exec cpprob#plot_matter_vs_cp  [matter] [dist]
   exec text#footer
 return
 
 macro plot_all
   exec plot_dist 2540    5e+09    4e+09
   exec plot_dist  735  1.8e+09  1.8e+09
 return
 
 
 macro print_all 
   exec print cpprob.ps -114 cpprob#plot_all
 return