34 tgt.fm.
p(4) = tgt.
mass 40 pjlab.fm.
p(3) = sqrt(pj.fm.
p(4)**2-pj.
mass**2)
41 pjlab.fm.
p(4) = pj.fm.
p(4)
48 call cgeqm(pjlab, tgt, cmsp, icon)
51 *
' cms cannot be formed in chANewLund; proj and ',
83 call cfrevent(kf,
charge, ia, iz, cmsp.
mass,
a, ntp)
97 if(
a(
i).fm.
p(4) .gt. maxe)
then subroutine cerrormsg(msg, needrtn)
subroutine cibst1(init, p1, p2, po)
subroutine cgeqm(p1, p2, q, icon)
subroutine ccos2kf(code, subcode, chg, kf)
max ptcl codes in the kseethru ! subcode integer regptcl
subroutine crot3mom(p, a, n)
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
real(4), dimension(:), allocatable, save temp
*Zfirst p fm *Zfirst p Zfirst p code
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code knuc
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p charge
subroutine cmkptc(code, subcode, charge, p)
*Zfirst p fm *Zfirst p mass
*Zfirst p fm *Zfirst p Zfirst p Zfirst p subcode