COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cpippTotXs.f
Go to the documentation of this file.
1 ! pi+p total xsection , elastic xs
2  subroutine cpipptotxs(p, xs)
3  implicit none
4 #include "Zmass.h"
5  real*8 p ! input. momentum of n. in GeV
6  real*8 xs ! output. total np cross section in mb
7  integer np, i, m
8  real*8 error
9  parameter(np=53, m=5)
10  real*8 px(np), mb(np)
11  real*8 Z, Y1, Y2, B, s, rts, s0, eta1, eta2
12 
13 
14  parameter(z=20.86, y1=19.24, y2=6.03, b=0.308,
15  * s0=5.38**2, eta1=0.458, eta2=0.545)
16  data ( px(i), i= 1 , np )/
17  1 0.1190 , 0.1275 , 0.1465 , 0.1649 ,
18  2 0.1811 , 0.2050 , 0.2197 , 0.2389 ,
19  3 0.2573 , 0.2703 , 0.2854 , 0.3058 ,
20  4 0.3180 , 0.3356 , 0.3595 , 0.3850 ,
21  5 0.4185 , 0.4594 , 0.5068 , 0.5645 ,
22  6 0.6198 , 0.6509 , 0.6737 , 0.7146 ,
23  7 0.7848 , 0.8327 , 0.8922 , 0.9944 ,
24  8 1.092 , 1.193 , 1.285 , 1.364 ,
25  9 1.439 , 1.542 , 1.636 , 1.718 ,
26  a 1.805 , 2.011 , 2.241 , 2.401 ,
27  b 2.573 , 2.825 , 3.242 , 4.066 ,
28  c 5.201 , 7.377 , 9.961 , 13.85 ,
29  d 19.08 , 27.20 , 39.73 , 59.20 ,
30  e 76.10
31  * /
32  data ( mb(i), i= 1 , np )/
33  1 9.845 , 13.81 , 19.99 , 30.48 ,
34  2 44.11 , 75.42 , 100.4 , 137.3 ,
35  3 174.4 , 193.5 , 200.7 , 190.5 ,
36  4 174.4 , 143.8 , 113.1 , 88.09 ,
37  5 62.45 , 45.45 , 33.24 , 24.57 ,
38  6 19.84 , 17.15 , 15.77 , 14.36 ,
39  7 15.21 , 17.41 , 20.14 , 23.30 ,
40  8 26.13 , 30.54 , 34.97 , 38.40 ,
41  9 40.67 , 38.20 , 35.14 , 32.84 ,
42  a 30.84 , 29.27 , 29.89 , 30.83 ,
43  b 30.99 , 30.19 , 29.10 , 27.61 ,
44  c 26.89 , 25.38 , 24.59 , 24.32 ,
45  d 23.81 , 23.43 , 23.17 , 23.40 ,
46  e 23.39
47  * /
48 
49 
50  save
51  if(p .gt. 30.) then
52 ! s =(M+E)^2 - p^2 = M^2 + m^2 +2ME
53 ! =
54  s = masp**2 + maspic**2 + 2*masp*sqrt(p**2 + maspic**2)
55  rts = sqrt(s)
56  xs =z + b*log(s/s0)**2 + y1*(1./s)**eta1 -y2*(1./s)**eta2
57  elseif( p .gt. 0.6) then
58  call kpolintplogxyfe(px, 1, mb, 1, np, m, 3, p, xs, error)
59  else
60  call cpippelaxs(p, xs)
61  endif
62  end
63 !
64  subroutine cpippelaxs(p, xs)
65 ! pi+ p elastic cross section in mb
66  implicit none
67  real*8 p ! input. momentum of n in GeV
68  real*8 xs ! output np elastic xs. mb.
69 
70  integer np, m, i
71  parameter(np=54, m=5)
72  real*8 px(np), mb(np)
73  real*8 error
74  real*8 xssave/-1./
75 
76  data ( px(i), i= 1 , np )/
77  1 0.1190 , 0.1275 , 0.1465 , 0.1649 ,
78  2 0.1811 , 0.2050 , 0.2197 , 0.2389 ,
79  3 0.2573 , 0.2703 , 0.2854 , 0.3058 ,
80  4 0.3180 , 0.3356 , 0.3595 , 0.3850 ,
81  5 0.4185 , 0.4594 , 0.5068 , 0.5645 ,
82  6 0.5900 , 0.6137 , 0.6414 , 0.6638 ,
83  7 0.6903 , 0.7395 , 0.7805 , 0.8039 ,
84  8 0.8444 , 0.8914 , 0.9273 , 0.9888 ,
85  9 1.075 , 1.146 , 1.241 , 1.349 ,
86  a 1.438 , 1.556 , 1.659 , 1.931 ,
87  b 2.283 , 2.821 , 3.716 , 5.167 ,
88  c 7.511 , 11.03 , 16.03 , 19.8157 ,
89  d 29.373 , 42.482 , 59.9484 , 84.1809,
90  e 134.34 , 240.072
91  * /
92 
93 
94  data ( mb(i), i= 1 , np )/
95  1 9.845 , 13.81 , 19.99 , 30.48 ,
96  2 44.11 , 75.42 , 100.4 , 137.3 ,
97  3 174.4 , 193.5 , 200.7 , 190.5 ,
98  4 174.4 , 143.8 , 113.1 , 88.09 ,
99  5 62.45 , 45.45 , 33.24 , 24.57 ,
100  6 21.01 , 18.74 , 16.11 , 14.29 ,
101  7 12.74 , 10.67 , 9.417 , 8.846 ,
102  8 8.396 , 8.752 , 9.612 , 11.30 ,
103  9 12.73 , 14.35 , 16.01 , 17.31 ,
104  a 17.40 , 16.34 , 13.83 , 10.99 ,
105  b 8.740 , 7.710 , 6.627 , 5.695 ,
106  c 5.075 , 4.740 , 4.359 , 3.94511 ,
107  d 3.46975 , 3.24744 , 3.20098 , 3.27183 ,
108  e 3.294 , 3.30042
109 
110  * /
111 
112  save
113 !
114  if( p .gt. 240.) then
115 ! assume prop.to total
116  if( xssave .lt. 0.) then
117  call cpipptotxs(px(np), xssave)
118  endif
119  call cpipptotxs(p, xs)
120  xs = xs * mb(np)/xssave
121  elseif(p .gt. 0.1) then
122  call kpolintplogxyfe(px, 1, mb, 1, np, m, 3, p, xs, error)
123  else
124 ! get value at 0.1
125  xs = mb(1)
126  endif
127  end
128  subroutine cpippinelaxs(p, xs)
129  implicit none
130  real(8),intent(in)::p
131  real(8),intent(out)::xs
132 
133  real(8)::txs, exs
134  call cpipptotxs(p, txs)
135  call cpippelaxs(p, exs)
136  xs =max( txs - exs, 0.d0)
137  end
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
e
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
kpolintplogxyfe
subroutine kpolintplogxyfe(xa, xstep, ya, ystep, nt, m, logxy, x, y, error)
Definition: kpolintp.f:22
d0
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
cpipptotxs
subroutine cpipptotxs(p, xs)
Definition: cpippTotXs.f:3
a
nodes a
Definition: cstdatmos0.fNew2.h:93
b
dE dx *! Nuc Int sampling table b
Definition: cblkMuInt.h:130
cpippelaxs
subroutine cpippelaxs(p, xs)
Definition: cpippTotXs.f:61
masp
masp
Definition: Zmass.h:5
maspic
maspic
Definition: Zmass.h:5
cpippinelaxs
subroutine cpippinelaxs(p, xs)
Definition: cpippTotXs.f:180
c
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130