COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Functions/Subroutines
cnpTotXs.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine cnptotxs (p, xs)
 
subroutine cnpelaxs (p, xs)
 
subroutine cnpinelaxs (p, xs)
 

Function/Subroutine Documentation

◆ cnpelaxs()

subroutine cnpelaxs ( real*8  p,
real*8  xs 
)

Definition at line 70 of file cnpTotXs.f.

References cnptotxs(), cppelaxs(), cpptotxs(), modpdgxs::csofstu(), e, false, kpolintpfe(), kpolintplogxyfe(), and parameter().

Referenced by cnpinelaxs(), and cnptotxs().

70  use modpdgxs, only : csofstu
71 ! np elastic cross section in mb
72  implicit none
73 #include "Zmass.h"
74  real*8 p ! input. momentum of n in GeV
75  real*8 xs ! output np elastic xs. mb.
76 
77  integer np, m, i
78  parameter(np=37, m=5)
79  real*8 px(np), mb(np)
80  real*8 error
81  real*8 xssave/-1./
82  real(8):: s
83  real(8),parameter:: pnorm=50.0
84  real(8),save:: xsn, norm
85  real(8):: xspel, xspt, xspit, ppp
86  real(8),save::xsnorm
87 
88  logical,save:: first = .true.
89 
90 
91  data ( px(i), i= 1 , np )/
92  1 0.9860e-02, 0.1271e-01, 0.1942e-01, 0.3230e-01,
93  2 0.4935e-01, 0.7025e-01, 0.1000 , 0.1187 ,
94  3 0.1294 , 0.1472 , 0.1692 , 0.1988 ,
95  4 0.2386 , 0.2833 , 0.3257 , 0.3785 ,
96  5 0.4803 , 0.5145 , 0.5818 , 0.6548 ,
97  6 0.7628 , 0.8543 , 0.9906 , 1.132 ,
98  7 1.243 , 1.413 , 1.575 , 1.712 ,
99  8 2.180 , 2.960 , 3.788 , 5.068 ,
100  9 8.427 , 16.91 , 37.07 , 66.05 ,
101  * 100.0
102  * /
103  data ( mb(i), i= 1 , np )/
104  1 8452. , 7556. , 6039. , 4437. ,
105  2 3306. , 2362. , 1549. , 1169. ,
106  3 1022. , 813.3 , 621.7 , 438.5 ,
107  4 301.1 , 198.5 , 136.3 , 100.1 ,
108  5 60.07 , 50.88 , 42.83 , 40.02 ,
109  6 34.96 , 32.67 , 31.88 , 30.92 ,
110  7 31.31 , 29.80 , 26.85 , 24.78 ,
111  8 21.65 , 18.23 , 15.92 , 12.92 ,
112  9 11.43 , 10.17 , 9.629 , 9.575 ,
113  * 9.578
114  * /
115  save
116  if( p .gt. pnorm) then
117  if( first ) then
118  ppp = pnorm
119  call cppelaxs(ppp, xspel)
120  call cpptotxs(ppp, xspt)
121  call cnptotxs(ppp, xspit)
122  call kpolintpfe(px, 1, mb, 1, np, m,
123  * ppp, xsnorm, error)
124  xs = xspit*xspel/xspt
125  norm = xs - xsnorm
126  first=.false.
127  endif
128  ppp = p
129  call cppelaxs(ppp, xspel)
130  call cpptotxs(ppp, xspt)
131  call cnptotxs(ppp, xspit)
132  xs = xspit*xspel/xspt
133  xs = xs - norm
134  elseif(p .gt. 0.01) then
135  call kpolintplogxyfe(px, 1, mb, 1, np, m, 3, p, xs, error)
136  else
137 ! get value at 0.01
138  xs = mb(1)
139  endif
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
norm
common ZdedxAir norm
Definition: ZdedxAir.h:2
e
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
cpptotxs
subroutine cpptotxs(p, xs)
Definition: cppTotXs.f:3
nrl_atmos::first
logical, save first
Definition: cNRLAtmos.f:8
i
nodes i
Definition: cstdatmos0.fNew2.h:93
cnptotxs
subroutine cnptotxs(p, xs)
Definition: cnpTotXs.f:3
true
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
Definition: cblkElemag.h:7
kpolintplogxyfe
subroutine kpolintplogxyfe(xa, xstep, ya, ystep, nt, m, logxy, x, y, error)
Definition: kpolintp.f:22
mb
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in mb
Definition: ZbpCnst.h:3
xs
! constants thru Cosmos real ! if multiplied to deg radian Torad ! light velocity m sec ! infinty ! kg m2 *Togpcm2 g cm2 ! g cm2 *Tokgpm2 kg m2 ! cm *Tom m ! m *Tocm cm ! g cm3 *Tokgpm3 kg m3 ! kg m3 *Togpcm3 g cm3 ! sec *Tonsec nsec ! Tesla m ! Avogadro *A2deninv ! mfp *n * xs
Definition: Zglobalc.h:18
p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
cppelaxs
subroutine cppelaxs(p, xs)
Definition: cppTotXs.f:68
modpdgxs::csofstu
real(8) function, public csofstu(p, Mp, Mt)
Definition: cpdgXs.f:21
m
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
kpolintpfe
subroutine kpolintpfe(xa, xstep, ya, ystep, nt, m, x, y, error)
Definition: kpolintp.f:134
false
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon ! knockon is considered Obsolete *PhotoProd false
Definition: cblkElemag.h:7
modpdgxs
Definition: cpdgXs.f:1
Here is the call graph for this function:
Here is the caller graph for this function:

◆ cnpinelaxs()

subroutine cnpinelaxs ( real(8), intent(in)  p,
real(8), intent(out)  xs 
)

Definition at line 142 of file cnpTotXs.f.

References cnpelaxs(), cnptotxs(), and d0.

Referenced by cinelx0().

142  implicit none
143  real(8),intent(in)::p
144  real(8),intent(out)::xs
145 
146  real(8)::txs, exs
147  call cnptotxs(p, txs)
148  call cnpelaxs(p, exs)
149  xs =max( txs - exs, 0.d0)
cnptotxs
subroutine cnptotxs(p, xs)
Definition: cnpTotXs.f:3
xs
! constants thru Cosmos real ! if multiplied to deg radian Torad ! light velocity m sec ! infinty ! kg m2 *Togpcm2 g cm2 ! g cm2 *Tokgpm2 kg m2 ! cm *Tom m ! m *Tocm cm ! g cm3 *Tokgpm3 kg m3 ! kg m3 *Togpcm3 g cm3 ! sec *Tonsec nsec ! Tesla m ! Avogadro *A2deninv ! mfp *n * xs
Definition: Zglobalc.h:18
p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
d0
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
cnpelaxs
subroutine cnpelaxs(p, xs)
Definition: cnpTotXs.f:70
Here is the call graph for this function:
Here is the caller graph for this function:

◆ cnptotxs()

subroutine cnptotxs ( real*8  p,
real*8  xs 
)

Definition at line 3 of file cnpTotXs.f.

References a, b, c, cnpelaxs(), modpdgxs::csigmat(), modpdgxs::csofstu(), d, e, kpolintpfe(), masn, masp, and parameter().

Referenced by cnpelaxs(), cnpinelaxs(), and ctotx0().

3  use modpdgxs, only : csigmat, csofstu
4  implicit none
5 #include "Zmass.h"
6 #include "Zevhnp.h"
7  real*8 p ! input. momentum of n. in GeV
8  real*8 xs ! output. total np cross section in mb
9  integer np, i, m
10  real*8 error, s
11  parameter(np=54, m=5)
12  real*8 px(np), mb(np)
13  data ( px(i), i= 1 , np )/
14  1 0.1001 , 0.1126 , 0.1224 , 0.1337 ,
15  2 0.1489 , 0.1716 , 0.1912 , 0.2119 ,
16  3 0.2395 , 0.2590 , 0.2856 , 0.3135 ,
17  4 0.3407 , 0.3739 , 0.4207 , 0.4595 ,
18  5 0.4996 , 0.5512 , 0.6082 , 0.6646 ,
19  6 0.7227 , 0.7781 , 0.8378 , 0.9246 ,
20  7 1.051 , 1.225 , 1.400 , 1.584 ,
21  8 1.827 , 2.161 , 2.506 , 2.949 ,
22  9 3.470 , 4.063 , 4.950 , 5.767 ,
23  a 6.854 , 8.432 , 10.07 , 12.09 ,
24  b 15.02 , 18.75 , 23.07 , 28.66 ,
25  c 36.68 , 47.41 , 61.27 , 78.80 ,
26  d 99.85 , 130.3 , 174.4 , 236.8 ,
27  e 318.3 , 395.5
28  * /
29  data ( mb(i), i= 1 , np )/
30  1 1465. , 1264. , 1097. , 946.4 ,
31  2 782.0 , 619.0 , 489.9 , 395.0 ,
32  3 301.3 , 242.9 , 197.0 , 152.1 ,
33  4 122.7 , 100.7 , 78.76 , 65.88 ,
34  5 56.14 , 48.73 , 44.17 , 42.57 ,
35  6 41.28 , 38.34 , 35.40 , 34.12 ,
36  7 34.12 , 36.07 , 38.13 , 40.30 ,
37  8 41.06 , 42.61 , 43.41 , 43.42 ,
38  9 43.42 , 42.64 , 42.12 , 42.13 ,
39  a 41.11 , 40.86 , 40.62 , 39.64 ,
40  b 39.16 , 39.90 , 39.66 , 39.67 ,
41  c 39.43 , 39.44 , 38.73 , 39.46 ,
42  d 39.23 , 39.48 , 40.47 , 39.99 ,
43  e 40.00 , 41.00
44  * /
45 
46  if(p .gt. 50. ) then
47  if( totxsopt == 1 ) then ! pdg parameterization (lowest value)
48  xs = csigmat('n','p', p)
49  else
50  s = csofstu(p, masn, masp)
51  if( totxsopt == 2 ) then ! midlle of 1 and 3
52  xs = 42.6*(s)**(-0.46) - 33.4*(s)**(-0.545) + 35.5 +
53  * 0.295*log(s/29.1)**2
54  elseif( totxsopt == 3) then ! TOTEM (COMPETE) fitting
55  xs = 42.6*(s)**(-0.46) - 33.4*(s)**(-0.545) + 35.5 +
56  * 0.307*log(s/29.1)**2
57  else
58  write(0,*) ' TotXSopt=', totxsopt, ' invalid'
59  stop
60  endif
61  endif
62  elseif( p .gt. 0.4) then
63  call kpolintpfe(px, 1, mb, 1, np, m, p, xs, error)
64  else
65  call cnpelaxs(p, xs)
66  endif
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
i
nodes i
Definition: cstdatmos0.fNew2.h:93
masn
masn
Definition: Zmass.h:5
modpdgxs::csigmat
real(8) function, public csigmat(a, b, p)
Definition: cpdgXs.f:31
mb
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in mb
Definition: ZbpCnst.h:3
xs
! constants thru Cosmos real ! if multiplied to deg radian Torad ! light velocity m sec ! infinty ! kg m2 *Togpcm2 g cm2 ! g cm2 *Tokgpm2 kg m2 ! cm *Tom m ! m *Tocm cm ! g cm3 *Tokgpm3 kg m3 ! kg m3 *Togpcm3 g cm3 ! sec *Tonsec nsec ! Tesla m ! Avogadro *A2deninv ! mfp *n * xs
Definition: Zglobalc.h:18
p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
cnpelaxs
subroutine cnpelaxs(p, xs)
Definition: cnpTotXs.f:70
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
modpdgxs::csofstu
real(8) function, public csofstu(p, Mp, Mt)
Definition: cpdgXs.f:21
m
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
kpolintpfe
subroutine kpolintpfe(xa, xstep, ya, ystep, nt, m, x, y, error)
Definition: kpolintp.f:134
nrl_atmos::a
real(4), save a
Definition: cNRLAtmos.f:20
masp
masp
Definition: Zmass.h:5
nrl_atmos::b
real(4), save b
Definition: cNRLAtmos.f:21
modpdgxs
Definition: cpdgXs.f:1
c
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130
Here is the call graph for this function:
Here is the caller graph for this function: