51 if( pj.fm.
p(4) .le. pj.
mass)
then 60 p = sqrt(pj.fm.
p(4)**2 - pj.
mass**2)
81 elseif(pj.
charge .eq. -1)
then 116 elseif(pj.
charge .eq. -1)
then 135 shp = (shp + shn)/2.0
140 elseif(pj.
charge .eq. -1)
then 170 write(0,*)
'Sorry: ctotx is not usable for heavy ions' 171 write(0,*)
' only inelastic xs is used for heavy ions' 172 write(0,*)
' so that you may use cinela or cAAXsec2' max ptcl codes in the kgzai
subroutine cpbarptotxs(p, xs)
subroutine cpptotxs(p, xs)
max ptcl codes in the kgnuc
subroutine cnptotxs(p, xs)
max ptcl codes in the klambdac
subroutine ckpptotxs(p, xs)
max ptcl codes in the kkaon
max ptcl codes in the kelec
subroutine cxp2xaxsec(a, xsxp, xsxa)
subroutine ckpntotxs(p, 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
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
*Zfirst p fm *Zfirst p Zfirst p code
subroutine ckmptotxs(p, xs)
max ptcl codes in the klambda
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code knuc
subroutine ckmntotxs(p, xs)
subroutine cpipptotxs(p, xs)
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p charge
max ptcl codes in the kseethru ! subcode integer antip
*Zfirst p fm *Zfirst p mass
max ptcl codes in the kpion
max ptcl codes in the ksigma
max ptcl codes in the kbomega
subroutine cpimptotxs(p, xs)
*Zfirst p fm *Zfirst p Zfirst p Zfirst p subcode
subroutine cnbarptotxs(p, xs)