COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cpbarpTotXs.f
Go to the documentation of this file.
1 ! pbarp total xsection, pbap ela
2  subroutine cpbarptotxs(p, xs)
3  implicit none
4 #include "Zmass.h"
5  real*8 p ! input. momentum of pbar. in GeV
6  real*8 xs ! output. total pbarp cross section in mb
7  integer np, i, m
8  real*8 error
9  parameter(np=34, m=5)
10  real*8 px(np), mb(np)
11  real*8 Z, Y1, Y2, B, s, rts, s0, eta1, eta2
12  parameter(z= 35.45, y1=42.53, y2=33.34, b=0.308,
13  * s0=5.38**2, eta1=0.458, eta2=0.545)
14 
15  data ( px(i), i= 1 , np )/
16  1 9.7289e-02, 0.1485 , 0.1925 , 0.2330 ,
17  2 0.2821 , 0.3558 , 0.4426 , 0.5736 ,
18  3 0.6850 , 0.8180 , 1.004 , 1.319 ,
19  4 1.686 , 2.098 , 2.911 , 4.207 ,
20  5 5.918 , 9.286 , 13.79 , 21.64 ,
21  6 32.15 , 49.08 , 70.95 , 109.8 ,
22  7 144.3 , 192.2 , 274.0 , 418.3 ,
23  8 702.7 , 1133. , 1877. , 2986. ,
24  9 5015. , 1.5568e+04
25  * /
26  data ( mb(i), i= 1 , np )/
27  1 445.7 , 380.5 , 321.4 , 283.2 ,
28  2 249.5 , 214.1 , 179.9 , 153.6 ,
29  3 134.6 , 123.0 , 113.7 , 103.9 ,
30  4 95.54 , 86.43 , 75.75 , 66.40 ,
31  5 59.76 , 54.64 , 51.30 , 47.16 ,
32  6 45.22 , 43.82 , 42.92 , 42.03 ,
33  7 41.82 , 41.39 , 41.84 , 42.07 ,
34  8 42.31 , 43.22 , 43.70 , 45.12 ,
35  9 46.59 , 50.73
36  * /
37  if(p .gt. 15000.) then
38 ! s =(m+E)^2 - p^2 = m^2 + p^2+m^2+2mE-p^2
39 ! = 2m(m+E)
40  s = (masp + sqrt(p**2+masp**2))*masp*2.
41  rts = sqrt(s)
42  xs =z + b*log(s/s0)**2 + y1*(1./s)**eta1 +y2*(1./s)**eta2
43  elseif( p .gt. 0.1) then
44  call kpolintplogxyfe(px, 1, mb, 1, np, m, 3, p, xs, error)
45  else
46  xs =mb(1)
47  endif
48  end
49  subroutine cpbarpelaxs(p, xs)
50 ! pbarp elastic cross section in mb
51  implicit none
52  real*8 p ! input. momentum of pbar in GeV
53  real*8 xs ! output pbarp elastic xs. mb.
54 
55  integer np, m, i
56  parameter(np=35, m=5)
57  real*8 px(np), mb(np)
58  real*8 error
59  real*8 xssave/-1.d0/
60 
61 
62  data ( px(i), i= 1 , np )/
63  1 9.8633e-02, 0.1485 , 0.2299 , 0.3062 ,
64  2 0.4078 , 0.5358 , 0.6850 , 0.8758 ,
65  3 1.166 , 1.575 , 1.986 , 2.505 ,
66  4 3.572 , 4.957 , 6.784 , 9.286 ,
67  5 11.87 , 15.18 , 21.64 , 32.60 ,
68  6 56.27 , 83.59 , 127.6 , 289.4 ,
69  7 486.1 , 936.0 , 1551. , 3241. ,
70  8 6680. , 1.4542e+04, 3.3434e+04, 6.4372e+04,
71  9 1.5845e+05, 6.2031e+05, 1.7030e+06
72  * /
73  data ( mb(i), i= 1 , np )/
74  1 106.0 , 96.90 , 86.76 , 78.90 ,
75  2 68.79 , 60.62 , 53.69 , 47.81 ,
76  3 43.25 , 37.51 , 32.19 , 26.90 ,
77  4 20.89 , 16.91 , 13.91 , 11.81 ,
78  5 10.80 , 9.772 , 9.078 , 8.568 ,
79  6 7.960 , 7.513 , 7.358 , 7.170 ,
80  7 7.135 , 7.252 , 7.488 , 8.239 ,
81  8 9.064 , 10.13 , 11.21 , 12.07 ,
82  9 13.43 , 15.42 , 17.42
83  * /
84  save
85  if( p .gt. 2000.) then
86 ! assume prop. to total
87  if(xssave .lt. 0.) then
88  call cpbarptotxs(px(27), xssave)
89  endif
90  call cpbarptotxs(p, xs)
91  xs = xs * mb(27)/xssave
92  elseif(p .gt. 0.1) then
93  call kpolintplogxyfe(px, 1, mb, 1, np, m, 3, p, xs, error)
94  else
95 ! get value at ~0.1
96  xs = 106.
97  endif
98  end
99  subroutine cpbarpinelaxs(p, xs)
100  implicit none
101  real(8),intent(in)::p
102  real(8),intent(out)::xs
103 
104  real(8)::txs, exs
105  call cpbarptotxs(p, txs)
106  call cpbarpelaxs(p, exs)
107  xs =max( txs - exs, 0.d0)
108  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
cpbarptotxs
subroutine cpbarptotxs(p, xs)
Definition: cpbarpTotXs.f:9
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
cpbarpelaxs
subroutine cpbarpelaxs(p, xs)
Definition: cpbarpTotXs.f:51
masp
masp
Definition: Zmass.h:5
cpbarpinelaxs
subroutine cpbarpinelaxs(p, xs)
Definition: cpbarpTotXs.f:116