fortran/file 20 flux.ps
meta 20 -111



*label axis size
igset lasi 0.022

* marker scale factor
igset mscf 2.0

* Axis label and values size
set asiz 0.35
set vsiz 0.30
set tsiz 0.30
set ksiz .56

* Yaxis label position
set xlab 1.80
set xmgl 3.0
set xval 0.6

* number divisions
set ndvx 6.
set ndvy 6.

*line widths
set hwid 4
set bwid 4
*igset lwid 20.

*text thickness
set pass 9
set cshi 0.01
igset txfp 60

zone 1 1




*cdir //pawc

opt utit
opt nbox
opt tic

**************************************************

vector/create e(5) r 10. 20. 50. 100. 250.

* mu+ fluxes

vector/create nue_1(5)   r 0.0e+0 6.3e-6  1.8e+8  2.7e+9  8.9e+10
vector/create nue_0(5)   r 2.7e+9 1.1e+10 7.1e+10 2.8e+11 1.6e+12
vector/create nue_m1(5)  r 5.6e+9 2.3e+10 1.4e+11 5.6e+11 3.2e+12

vector/create anum_1(5)   r 3.8e+9 1.4e+10 9.4e+10 3.7e+11 2.2e+12
vector/create anum_0(5)   r 2.8e+9 1.1e+10 7.1e+10 2.8e+11 1.6e+12
vector/create anum_m1(5)  r 1.8e+9 7.5e+9  4.7e+10 1.8e+11 1.1e+12

* mu- fluxes

vector/create num_1(5)   r 1.8e+9 7.5e+9  4.7e+10 1.8e+11 1.1e+12
vector/create num_0(5)   r 2.6e+9 1.0e+10 7.1e+10 2.8e+11 1.6e+12
vector/create num_m1(5)  r 3.3e+9 1.4e+10 9.3e+10 3.7e+11 2.2e+12

vector/create anue_1(5)   r 5.2e+9 2.2e+10 1.4e+11 5.6e+11 3.2e+12
vector/create anue_0(5)   r 2.6e+9 1.1e+10 7.0e+10 2.8e+11 1.6e+12
vector/create anue_m1(5)  r 3.1e-6 9.5e-6  1.4e+8  2.4e+9  8.9e+10

**************************************************
* modify fluxes to take account of f_decay versus energy
* which was set = 0.3 for all energies in calculation
* (so multiply by 0.25/0.3 = 0.8333) 
* and correct for different flight distance (L = 8480 km in 
* calculation, but distance to Gran sasso = 9900 km) -> 
* divide by (8480/9900)**2 = 0.734)
* hence fcor = 0.8333*0.734 = 0.612

vector/create fcor(5) r 5*1.

application sigma
nue_1 = fcor * nue_1
nue_0 = fcor * nue_0
nue_m1 = fcor * nue_m1
anum_1 = fcor * anum_1
anum_0 = fcor * anum_0
anum_m1 = fcor * anum_m1
num_1 = fcor * num_1
num_0 = fcor * num_0
num_m1 = fcor * num_m1
anue_1 = fcor * anue_1
anue_0 = fcor * anue_0
anue_m1 = fcor * anue_m1
exit




**************************************************

zone 2 2
opt logy

set xwin 0.01
set ywin 0.1
set xlab 0.01
set vsiz 0.001
set xmgr 0.3
*set ymgl 0.01
set ygti 0.4
set ymgu 1.2
*set xsiz 20.
*set ysiz 20.

set gsiz 0.4
title 'Neutrino Fluxes at L = 10000 Km'

** num

graphics/hplot/null 0. 300. 1.e+8 2.e+13 
text 20. 5.0e+12 '[n]?[m]!' 0.7 0.


text -60. 1.e+13 '10^13!' 0.45 0.
text -60. 1.e+12 '10^12!' 0.45 0.
text -60. 1.e+11 '10^11!' 0.45 0.
text -60. 1.e+10 '10^10!' 0.45 0.
text -60. 1.e+9  '10^9!'  0.45 0.

text -75. 1.0e+8 'Flux  (m^-2! year^-1!)' 0.6 90. c

igset mtyp 20
igset pmci 4
graph 5 e num_0 pc

igset mtyp 22
igset pmci 2
graph 5 e num_m1 pc

igset mtyp 21
igset pmci 6
graph 5 e num_1 pc

** anti nue

graphics/hplot/null 0. 300. 1.e+8 2.e+13 
text 19. 9.0e+12 '-' 0.8 0.
text 20. 5.0e+12 '[n]?e!' 0.7 0.

igset mtyp 20
igset pmci 4
graph 5 e anue_0 pc

igset mtyp 22
igset pmci 2
graph 3 e(3:5) anue_m1(3:5) pc

igset mtyp 21
igset pmci 6
graph 5 e anue_1 pc

** nue

graphics/hplot/null 0. 300. 1.e+8 2.e+13 
text 20. 5.0e+12 '[n]?e!' 0.7 0.

text -60. 1.e+13 '10^13!' 0.45 0.
text -60. 1.e+12 '10^12!' 0.45 0.
text -60. 1.e+11 '10^11!' 0.45 0.
text -60. 1.e+10 '10^10!' 0.45 0.
text -60. 1.e+9  '10^9!'  0.45 0.
text 0. 4.0e+7 '0' 0.45 0. c
text 100. 4.0e+7 '100' 0.45 0. c
text 200. 4.0e+7 '200' 0.45 0. c

igset mtyp 20
igset pmci 4
graph 5 e nue_0 pc

igset mtyp 22
igset pmci 2
graph 5 e nue_m1 pc

igset mtyp 21
igset pmci 6
graph 3 e(3:5) nue_1(3:5) pc


** anti num

graphics/hplot/null 0. 300. 1.e+8 2.e+13 
text 19. 9.0e+12 '-' 0.8 0.
text 20. 5.0e+12 '[n]?[m]!' 0.7 0.

text 0. 4.0e+7 '0' 0.45 0. c
text 100. 4.0e+7 '100' 0.45 0. c
text 200. 4.0e+7 '200' 0.45 0. c

igset mtyp 20
igset pmci 4
graph 5 e anum_0 pc

igset mtyp 22
igset pmci 2
graph 5 e anum_m1 pc

igset mtyp 21
igset pmci 6
graph 5 e anum_1 pc

text 0. 9.0e+6 'Muon Beam Energy  (GeV)' 0.6 0. c

igset pmci 6
graphics/hplot/key 70. 1.5e+10 21 
text 90. 1.25e+10 'Polarization = +1' 0.4 0.
igset pmci 4
graphics/hplot/key 70. 0.375e+10 20 
text 90. 0.3125e+10 'Unpolarized' 0.4 0.
igset pmci 2
graphics/hplot/key 70. 0.09375e+10 22 
text 90. 0.078125e+10 'Polarization = -1' 0.4 0.

****************


***********************************************

close 20

MESSAGE 'Plots Done'