cinelx.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cinelx (pj, A, Z, xs)
subroutine	cinelx0 (pj, A, Z, xs)
real(8) function	cinelcosbypdg (A)

Function/Subroutine Documentation

cinelcosbypdg()

real(8) function cinelcosbypdg

(

real(8), intent(in)

A

)

Definition at line 231 of file cinelx.f.

!             sigmaInela(pA)cosmos/sigmaInela(pA)PDG
!        @ 200 GeV.  SigmaInela by cosmos is the one
!             obtained by cinelx0  before  
!             cxp2xAXsec  included cxAbyxpXsec.
!       
      real(8),intent(in):: a
      real(8):: ratio

      real(8):: temp

      if( a > 48. ) then
         if(a >150.) then
            ratio = 1.0
         else
            ratio = (1.006+0.0268*a/100.)*0.957
         endif
      elseif(a > 20.) then
         ratio = 0.975
      elseif( abs(a - 16.) < 0.1 ) then
         ratio = 0.9427
      elseif( abs( a- 7.0) < 0.1 ) then
         ratio = 0.902
      elseif( a > 7. ) then 
         temp = log(a)
         ratio = (-0.100*temp + 0.57311)*temp +0.1578
      elseif( abs( a- 4.) < 0.1  ) then
         ratio = 0.975
      elseif( abs( a- 2.) < 0.1 ) then
         ratio = 1.137
      elseif (abs(a-1.0) < 0.1 ) then
         ratio =1.0138
      else
         ratio = -0.025*(a-4) + 0.975
      endif

cinelx()

subroutine cinelx	(	type(ptcl)	pj,
		real(8), intent(in)	A,
		real(8), intent(in)	Z,
		real(8), intent(out)	xs
	)

Definition at line 4 of file cinelx.f.

References cinelx0(), cmkptc(), kgnuc, knuc, masn, and masp.

Referenced by caaxsec2(), caveintenuc(), cfakegh(), csampcollina(), csampcolln(), csmpcolina2(), cxsdpmjet3(), cxsgheisha(), cxsjam(), cxsother(), cxsphits(), and cxsqgsjet2().

      implicit none
#include "Zptcl.h"
#include "Zcode.h"
#include "Zmass.h"
#include "Zevhnp.h"
      type(ptcl)::pj   !input projectile particle
      real(8),intent(in)::  a     ! input. effective target mass number
      real(8),intent(in)::  z     ! input. effective target charge number
      real(8),intent(out):: xs    ! output. in mb. inelastic xs.

      real(8):: gamma
      type(ptcl):: target
      integer:: tcharge
      real(8):: pja, pjz


      if( a == 1.0 .and. pj%code == kgnuc  ) then ! target is p/n, pj is heavy ion
               ! we use proj. rest system ;  if not, inside cinelx0 
               ! cAAXsec2 is called and it calles cinelx and recursive
               ! call happens.
         tcharge = z
         call cmkptc(knuc, -1, tcharge, target)  ! p or n
         gamma = pj%fm%p(4)/pj%mass       
         pja = pj%subcode
         pjz = pj%charge   
         if( tcharge == 0 ) then
            target%fm%p(4) = gamma*masn
         else
            target%fm%p(4) = gamma*masp
         endif
         call cinelx0(target, pja, pjz, xs)   ! target is now proj. pjA is target

      else
         call cinelx0(pj, a, z,  xs)      
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinelx0()

subroutine cinelx0	(	type(ptcl)	pj,
		real(8), intent(in)	A,
		real(8), intent(in)	Z,
		real(8), intent(out)	xs
	)

Definition at line 42 of file cinelx.f.

References antip, caaxsec2(), cgpxsec(), ckmninelaxs(), ckmpinelaxs(), ckpninelaxs(), ckppinelaxs(), cnbarpinelaxs(), cnpinelaxs(), cpbarpinelaxs(), cpimpinelaxs(), cpippinelaxs(), cppinelaxs(), cxp2xaxsec(), d0, kbomega, kddb, kdmes, kelec, kgnuc, kgzai, kkaon, klambda, klambdac, kmuon, kneue, kneumu, kneutau, knnb, knuc, kphoton, kpion, ksigma, modcosxsec::nisdiff, and norm.

Referenced by cinelx().

      use modcosxsec  ! refer only nIsDiff
      implicit none
#include "Zptcl.h"
#include "Zcode.h"
#include "Zevhnp.h"
      type(ptcl)::pj   !input projectile particle
      real(8),intent(in)::  a          ! input. effective target mass number
      real(8),intent(in)::  z     ! input. effective target charge number
      real(8),intent(out):: xs          ! output. in mb. total xsection. 


      real(8) p
      real(8),external:: cpdgsigmainepa
      real(8),save:: aprev=-100., norm
      real(8):: xs200, xs200a


      real(8)::shp, shn

      if( pj%fm%p(4) .le. pj%mass) then
         p =0.
         if(pj%code .eq. knuc .and. pj%subcode .eq. antip) then
            xs = largexs
         elseif(pj%code .eq. kelec .and. pj%charge .eq. 1) then
            xs = largexs
         else
            xs  = smallxs
         endif
      else   
         p = sqrt(pj%fm%p(4)**2 - pj%mass**2)
         if(pj%code .eq. knuc) then
            if(pj%charge .eq. 1) then
               if( a == 1.0 ) then
                  if( z == 1.0 ) then
                     call cppinelaxs(p, xs)
                  elseif( z == 0.) then
                     call cnpinelaxs(p, xs)
                  else
                     write(0,*) 'A,Z=',a,z, 'for p pj'
                     write(0,*) 'in cinelx; invalid'
                     stop
                  endif
               elseif(a > 1.) then
                  call cppinelaxs(p, shp)
                  if( nisdiff ) then 
                     call cnpinelaxs(p, shn)
                     shp = (shp + shn)/2.0
                  endif
                  call cxp2xaxsec(a, shp, xs)
               else
                  write(0,*) ' A,Z for cinelx invalid'
                  write(0,*) ' pj is p '
                  stop
               endif
            elseif(pj%charge .eq. -1) then
               call cpbarpinelaxs(p, shp)
               if( a .gt. 1.) then
                  if(nisdiff) then
                     call cnbarpinelaxs(p, shn)
                     shp = (shp + shn)/2.
                  endif
                  call cxp2xaxsec(a, shp, xs)
               else
                  xs = shp
               endif
            elseif(pj%subcode .eq. antip) then
!              anti-neutron; assume the same one as pbar
               call cnbarpinelaxs(p, shp)
               if( a .gt. 1.) then
                  call cxp2xaxsec(a, shp, xs)
               else
                  xs = shp
               endif
            else
!               neutron
               call cnpinelaxs(p, shp)
               if(a .gt. 1.) then
                  call cxp2xaxsec(a, shp, xs)
               else
                  xs = shp
               endif
            endif
         elseif(pj%code == kpion ) then
            if(pj%charge .eq. 1) then
               call cpippinelaxs(p, shp)
            elseif(pj%charge .eq. -1) then
               call cpimpinelaxs(p, shp) 
            else 
                  ! pi0 
               call cpimpinelaxs(p, shp)
            endif
            if(a .ne. 1.0) then
               call cxp2xaxsec(a, shp, xs)
            else
               xs = shp
            endif
         elseif(pj%code .eq. kkaon) then
            if(pj%charge .eq. 1) then
               if(a == 1.0) then
                  if( z==1.0) then
                     call ckppinelaxs(p, xs)
                  else
                     call ckpninelaxs(p, xs)
                  endif
               elseif(a .gt. 1.) then
                  call ckppinelaxs(p, shp)
                  if( nisdiff ) then
                     call ckpninelaxs(p, shn)
                     shp = (shp + shn)/2.0
                  endif
                  call cxp2xaxsec(a, shp, xs)               
               else
                  write(0,*) 'A,Z=',a,z, ' invalid for cinel%Kaon'
                  stop
               endif
            elseif(pj%charge .eq. -1) then
               call ckmpinelaxs(p, shp)
               if(a .gt. 1.) then
                  if( nisdiff ) then
                     call ckmninelaxs(p, shn)
                     shp = (shp + shn)/2.
                  endif
                  call cxp2xaxsec(a, shp, xs)               
               else
                  xs = shp
               endif
            else
!             k0; don't worry so much 
               call ckmpinelaxs(p, shp)
!!               call ckmnInelaXs(p, shn)
!!               shp = (shp+ shn)/2.
               if(a .gt. 1.) then
                  call cxp2xaxsec(a, shp, xs)               
               else
                  xs = shp
               endif
            endif
         elseif(pj%code .eq. kgnuc ) then
!                heavy ion 
            call caaxsec2(pj, a, z,  xs)
         elseif(pj%code .eq. kdmes) then
            call ckppinelaxs(p, shp)
            call cxp2xaxsec(a, shp, xs)               
         elseif(pj%code .eq. kgzai .or. pj%code .eq. ksigma .or.
     *           pj%code .eq. kbomega .or. pj%code .eq. klambda .or.
     *           pj%code .eq. klambdac) then
!          don't worry, not used almost at all; use proton
!            shp = ctotpp1(p)
            call cppinelaxs(p, shp)
            call cxp2xaxsec(a, shp, xs)
         elseif(pj%code == kphoton ) then
            call cgpxsec(a, pj%fm%p(4), xs)
         elseif(pj%code .eq. kneumu) then
            xs= smallxs
         elseif(pj%code .eq. kneue) then
            xs= smallxs
         elseif(pj%code == kneutau) then
            xs= smallxs
         elseif(pj%code == knnb .or. pj%code == kddb ) then
            xs = smallxs
         elseif(pj%code == kmuon) then
            write(0,*) 'cinelx: code=',pj%code, ' should not come'
            stop
         else
!     eta, ds, Xic,  etc
           xs = smallxs 
!!          use pion
   !         call cpippInelaXs(p, shp)
   !         call cxp2xAXsec(A, shp, xs)               
         endif
      endif
!!!       normalize to pdg at 200 GeV. norm const is
!!!   applied to other collisions
      if( xs /= smallxs ) then
         if(a > 1. .and. a /=  aprev ) then
            call cppinelaxs(200.d0, xs200)
            call cxp2xaxsec(a, xs200, xs200a)
!!         write(0,*) ' xs200A=',xs200A, ' pgd=',  cPDGsigmaInepA(A)
            norm =  cpdgsigmainepa(a) /xs200a
!!         write(0,*) ' norm=', norm
            aprev = a
         elseif( a == 1.)  then
            norm = 1.
            aprev = a
         endif
         xs = xs*norm
      endif

Here is the call graph for this function:

Here is the caller graph for this function: