COSMOS v7.655  COSMOSv7655
(AirShowerMC)
ckppTotXs.f
Go to the documentation of this file.
1 ! K+ p total xsection , elastic xs
2  subroutine ckpptotxs(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=40, m=5)
10  real*8 px(np), mb(np)
11  real*8 Z, Y1, Y2, B, s, rts, s0, eta1, eta2
12 
13  parameter(z=17.91, y1=7.1, y2=13.45, b=0.308,
14  * s0=5.38**2, eta1=0.458, eta2=0.545)
15 
16  data ( px(i), i= 1 , np )/
17  1 0.1005 , 0.1706 , 0.2619 , 0.3586 ,
18  2 0.4511 , 0.5370 , 0.6587 , 0.7537 ,
19  3 0.8003 , 0.8846 , 0.9779 , 1.034 ,
20  4 1.065 , 1.131 , 1.166 , 1.263 ,
21  5 1.388 , 1.466 , 1.549 , 1.728 ,
22  6 1.938 , 2.173 , 2.377 , 2.653 ,
23  7 3.319 , 3.913 , 4.824 , 6.409 ,
24  8 8.305 , 10.19 , 12.50 , 15.33 ,
25  9 20.07 , 27.34 , 41.97 , 59.50 ,
26  a 85.63 , 126.3 , 203.9 , 302.5
27  * /
28 
29  data ( mb(i), i= 1 , np )/
30  1 10.98 , 11.55 , 11.91 , 12.33 ,
31  2 12.54 , 12.56 , 12.46 , 12.59 ,
32  3 13.13 , 14.05 , 15.20 , 16.17 ,
33  4 16.95 , 17.81 , 18.31 , 18.54 ,
34  5 18.33 , 18.02 , 17.81 , 17.89 ,
35  6 17.73 , 17.55 , 17.49 , 17.36 ,
36  7 17.28 , 17.20 , 17.25 , 17.25 ,
37  8 17.33 , 17.38 , 17.43 , 17.45 ,
38  9 17.50 , 17.71 , 17.97 , 18.30 ,
39  a 18.67 , 19.19 , 19.92 , 20.62
40  * /
41 
42  save
43  if(p .gt. 250.) then
44 ! s =(M+E)^2 - p^2 = M^2 + m^2 +2ME
45 ! =
46  s = masp**2 + maskc**2 + 2*masp*sqrt(p**2 + maskc**2)
47  rts = sqrt(s)
48  xs =z + b*log(s/s0)**2 + y1*(1./s)**eta1 +y2*(1./s)**eta2
49  elseif( p .gt. 0.6) then
50 ! take log x only
51  call kpolintplogxyfe(px, 1, mb, 1, np, m, 1, p, xs, error)
52  else
53  call ckppelaxs(p, xs)
54  endif
55  end
56 !
57  subroutine ckppelaxs(p, xs)
58 ! pi+ p elastic cross section in mb
59  implicit none
60  real*8 p ! input. momentum of n in GeV
61  real*8 xs ! output np elastic xs. mb.
62 
63  integer np, m, i
64  parameter(np=27, m=5)
65  real*8 px(np), mb(np)
66  real*8 error
67  real*8 xssave/-1./
68  data ( px(i), i= 1 , np )/
69  1 0.1005 , 0.1706 , 0.2619 , 0.3586 ,
70  2 0.4511 , 0.5370 , 0.6330 , 0.7725 ,
71  3 0.9106 , 1.011 , 1.203 , 1.410 ,
72  4 1.669 , 1.996 , 2.497 , 3.107 ,
73  5 4.066 , 5.427 , 7.428 , 9.819 ,
74  6 14.92 , 22.57 , 34.81 , 52.64 ,
75  7 83.26 , 133.7 , 250.5
76  * /
77  data ( mb(i), i= 1 , np )/
78  1 10.98 , 11.55 , 11.91 , 12.33 ,
79  2 12.54 , 12.56 , 12.46 , 12.17 ,
80  3 12.06 , 11.56 , 10.62 , 9.627 ,
81  4 8.319 , 6.958 , 5.701 , 4.811 ,
82  5 4.129 , 3.525 , 3.366 , 3.312 ,
83  6 3.387 , 2.992 , 2.492 , 2.358 ,
84  7 2.433 , 2.429 , 2.686
85  * /
86 
87 
88  save
89 !
90  if( p .gt. 150.) then
91 ! assume prop.to total
92  if( xssave .lt. 0.) then
93  call ckpptotxs(px(np), xssave)
94  endif
95  call ckpptotxs(p, xs)
96  xs = xs * mb(np)/xssave
97  elseif(p .gt. px(1)) then
98 ! take log x only
99  call kpolintplogxyfe(px, 1, mb, 1, np, m, 1, p, xs, error)
100  else
101 ! get value at 0.1
102  xs = mb(1)
103  endif
104  end
105  subroutine ckppinelaxs(p, xs)
106  implicit none
107  real(8),intent(in)::p
108  real(8),intent(out)::xs
109 
110  real(8)::txs, exs
111  call ckpptotxs(p, txs)
112  call ckppelaxs(p, exs)
113  xs =max( txs - exs, 0.d0)
114  end
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
ckppelaxs
subroutine ckppelaxs(p, xs)
Definition: ckppTotXs.f:51
ckpptotxs
subroutine ckpptotxs(p, xs)
Definition: ckppTotXs.f:3
kpolintplogxyfe
subroutine kpolintplogxyfe(xa, xstep, ya, ystep, nt, m, logxy, x, y, error)
Definition: kpolintp.f:22
maskc
maskc
Definition: Zmass.h:5
d0
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
a
nodes a
Definition: cstdatmos0.fNew2.h:93
masp
masp
Definition: Zmass.h:5
ckppinelaxs
subroutine ckppinelaxs(p, xs)
Definition: ckppTotXs.f:116