* macro to plot neutrino fluxes properly normalized on 
* same plot for several polarizations and 
* neutrino energies
* all are plotted on axis at 1 km 



exec flux#numu1 
wait
exec flux#numu0 
wait
exec flux#numum1
wait
exec flux#nue1
wait
exec flux#nue0
wait
exec flux#nuem1 
wait

exec flux#enumu1 
wait
exec flux#enumu0 
wait
exec flux#enumum1
wait
exec flux#enue1
wait
exec flux#enue0
wait
exec flux#enuem1 
wait
exec flux#numu0-wdbnd
wait
exec flux#nue0-wdbnd
wait
exec flux#efxnumu0
wait
exec flux#efxnue0 
wait

return

macro  numu1 

for/file 63 numu1.eps  
meta 63 -113  

exec pretty


title 'mu- decays polarization= +1'
null 0. 20. 0. 0.0007  
itx 4. 0.00066 '[n]?[m]! flux at 1 km on axis' 
itx 4. 0.00060 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/muon' 

func/plot flxnumu.f(x,1.,20.) 0. 20. s 

func/plot flxnumu.f(x,1.,7.) 0. 7. s 

func/plot flxnumu.f(x,1.,4.) 0. 4. s 

close 63 

return  

macro numu0


for/file 63 numu0.eps  
meta 63 -113  

exec pretty


title 'mu- decays polarization= 0'
null 0. 20. 0. 0.0007  
itx 4. 0.00066 '[n]?[m]! flux at 1 km on axis' 
itx 4. 0.00060 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/muon' 

func/plot flxnumu.f(x,0.,20.) 0. 20. s 

func/plot flxnumu.f(x,0.,7.) 0. 7. s 

func/plot flxnumu.f(x,0.,4.) 0. 4. s 
close 63 


return 


macro numu0-wdbnd


for/file 63 numu0-wdbnd.eps  
meta 63 -113  

exec pretty

hist/file 1 ../../beam/numuv14.hst 

hist/opera/add 40143 40143 9000 3. 3. 

title 'mu- decays polarization= 0'
opt logy  

null 0. 15. 0. 0.0007  

itx 2. 0.0004 '[n]?[m]! flux at 1 km on axis' 
itx 2. 0.0002 'nufact flux for E?[m]! = 20, 7, 4 GeV' 

itx 2. 0.0001 'muon/proton = 1/6' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/proton' 

func/plot flxnumu.f(x,0.,20.)/6. 0. 20. s 

func/plot flxnumu.f(x,0.,7.)/6. 0. 7. s 

func/plot flxnumu.f(x,0.,4.)/6. 0. 4. s 

set hcol 2 
hist/pl 40143 s
set txci 2 
itx 2. 2.e-5 '28 GeV BNL wide band beam'
arrow 3. 2. 1.7e-5  1.2e-5
set hcol 1 

close 63 


return 




macro nue0-wdbnd


for/file 63 nue0-wdbnd.eps  
meta 63 -113  

exec pretty

hist/file 1 ../../beam/numuv14.hst 



title 'mu- decays polarization= 0'
opt logy  

null 0. 15. 0. 0.0007  

itx 2. 0.0004 'Anti-[n]?[e]! flux at 1 km on axis' 
itx 2. 0.0002 'nufact flux for E?[m]! = 20, 7, 4 GeV' 

itx 2. 0.0001 'muon/proton = 1/6' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/proton' 

func/plot flxnue.f(x,0.,20.)/6. 0. 20. s 

func/plot flxnue.f(x,0.,7.)/6. 0. 7. s 

func/plot flxnue.f(x,0.,4.)/6. 0. 4. s 

set hcol 2 
hist/pl 40143 s
set txci 2 
itx 2. 2.e-5 '28 GeV BNL wide band beam'
arrow 3. 2. 1.7e-5  1.2e-5
set hcol 1 

close 63 


return 






macro numum1


for/file 63 numum1.eps  
meta 63 -113  

exec pretty


title 'mu- decays polarization= -1'
null 0. 20. 0. 0.0007 
itx 4. 0.00066 '[n]?[m]! flux at 1 km on axis' 
itx 4. 0.00060 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/muon' 

func/plot flxnumu.f(x,-1.,20.) 0. 20. s 

func/plot flxnumu.f(x,-1.,7.) 0. 7. s 

func/plot flxnumu.f(x,-1.,4.) 0. 4. s 

close 63

return 



macro nue1

for/file 63 nue1.eps  
meta 63 -113  

exec pretty

title 'mu- decays polarization= 1'
null 0. 20. 0. 0.0007 
itx 4. 0.00066 'Anti-[n]?e! flux at 1 km on axis' 
itx 4. 0.00060 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/muon' 

func/plot flxnue.f(x,1.,20.) 0. 20. s 

func/plot flxnue.f(x,1.,7.) 0. 7. s 

func/plot flxnue.f(x,1.,4.) 0. 4. s 

close 63 

return 



macro nue0


for/file 63 nue0.eps  
meta 63 -113  
exec pretty

title 'mu- decays polarization= 0'
null 0. 20. 0. 0.0007 
itx 4. 0.00066 'Anti-[n]?e! flux at 1 km on axis' 
itx 4. 0.00060 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/muon' 

func/plot flxnue.f(x,0.,20.) 0. 20. s 

func/plot flxnue.f(x,0.,7.) 0. 7. s 

func/plot flxnue.f(x,0.,4.) 0. 4. s 

close 63 

return 


macro nuem1


for/file 63 nuem1.eps  
meta 63 -113  

exec pretty

title 'mu- decays polarization= -1'
null 0. 20. 0. 0.0007 
itx 4. 0.00066 'Anti-[n]?e! flux at 1 km on axis' 
itx 4. 0.00060 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'nus/GeV/m^2!/muon' 

func/plot flxnue.f(x,-1.,20.) 0. 20. s 

func/plot flxnue.f(x,-1.,7.) 0. 7. s 

func/plot flxnue.f(x,-1.,4.) 0. 4. s 


close 63 


return 




macro  enumu1 

for/file 63 numu1.evts.eps  
meta 63 -113  

exec pretty


title '1e20 mu- decays polarization= +1'
null 0. 20. 0. 50. 
itx 4. 45. '[n]?[m]! events at 2900 km on axis' 
itx 4. 40. 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnumu.f(x,1.,20.) 0. 20. s 

func/plot evtsnumu.f(x,1.,7.) 0. 7. s 

func/plot evtsnumu.f(x,1.,4.) 0. 4. s 

close 63 

return  

macro enumu0


for/file 63 numu0.evts.eps  
meta 63 -113  

exec pretty


title '1e20 mu- decays polarization= 0'
null 0. 20. 0. 50.  
itx 4. 45. '[n]?[m]! evts at 2900 km on axis' 
itx 4. 40. 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnumu.f(x,0.,20.) 0. 20. s 

func/plot evtsnumu.f(x,0.,7.) 0. 7. s 

func/plot evtsnumu.f(x,0.,4.) 0. 4. s 
close 63 


return 

macro enumum1


for/file 63 numum1.evts.eps  
meta 63 -113  

exec pretty


title '1e20 mu- decays polarization= -1'
null 0. 20. 0. 50. 
itx 4. 45. '[n]?[m]! evts at 2900 km on axis' 
itx 4. 40. 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnumu.f(x,-1.,20.) 0. 20. s 

func/plot evtsnumu.f(x,-1.,7.) 0. 7. s 

func/plot evtsnumu.f(x,-1.,4.) 0. 4. s 

close 63

return 



macro enue1

for/file 63 nue1.evts.eps  
meta 63 -113  

exec pretty

title '1e20 mu- decays polarization= 1'
null 0. 20. 0. 50. 
itx 4. 45. 'Anti-[n]?e! evts at 2900 km on axis' 
itx 4. 40. 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnue.f(x,1.,20.) 0. 20. s 

func/plot evtsnue.f(x,1.,7.) 0. 7. s 

func/plot evtsnue.f(x,1.,4.) 0. 4. s 

close 63 

return 



macro enue0


for/file 63 nue0.evts.eps  
meta 63 -113  
exec pretty

title '1e20 mu- decays polarization= 0'
null 0. 20. 0. 50. 
itx 4. 45. 'Anti-[n]?e! evts at 2900 km on axis' 
itx 4. 40. 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnue.f(x,0.,20.) 0. 20. s 

func/plot evtsnue.f(x,0.,7.) 0. 7. s 

func/plot evtsnue.f(x,0.,4.) 0. 4. s 

close 63 

return 


macro enuem1


for/file 63 nuem1.evts.eps  
meta 63 -113  

exec pretty

title '1e20 mu- decays polarization= -1'
null 0. 20. 0. 50. 
itx 4. 45 'Anti-[n]?e! evts at 2900 km on axis' 
itx 4. 40 'E?[m]! = 20, 7, 4 GeV' 

atitle 'E?[n]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnue.f(x,-1.,20.) 0. 20. s 

func/plot evtsnue.f(x,-1.,7.) 0. 7. s 

func/plot evtsnue.f(x,-1.,4.) 0. 4. s 


close 63 


return 

macro efxnumu0

for/file 63 numu0.e4.7.evts.eps  
meta 63 -113  

exec pretty

title '1e20 mu- decays polarization= 0'
opt logy
null 0. 20. 0. 10. 

itx 4. 4. '[n]?[m]! evts at 2900 km on axis' 
itx 4. 2.3 'E?[n]! = 9.4, 4.7, 1.6, 0.94 GeV' 

atitle 'E?[m]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnumu.f(9.4,0.,x) 9.4 20. s 

func/plot evtsnumu.f(4.7,0.,x) 4.7 20. s 

func/plot evtsnumu.f(1.6,0.,x) 1.6 20. s 

func/plot evtsnumu.f(0.94,0.,x) 0.94 20. s 


close 63 


return 




macro efxnue0


for/file 63 nue0.e4.7.evts.eps  
meta 63 -113  

exec pretty

title '1e20 mu- decays polarization= 0'
opt logy
null 0. 20. 0. 10. 
itx 4. 5. 'Anti-[n]?e! evts at 2900 km on axis' 
itx 4. 3. 'E?[n]! = 9.4, 4.7, 1.6, 0.94 GeV' 

atitle 'E?[m]! GeV' 'evts/GeV/1kT/1e20mu' 

func/plot evtsnue.f(9.4,0.,x) 9.4 20. s 

func/plot evtsnue.f(4.7,0.,x) 4.7 20. s 

func/plot evtsnue.f(1.6,0.,x) 1.6 20. s 

func/plot evtsnue.f(0.94,0.,x) 0.94 20. s 


close 63 


return