qgs01init.f

Go to the documentation of this file.

      subroutine qgs01init
      implicit none
#include "Zcode.h"
#include "Zptcl.h"
#include "Zmass.h"
!
!
      type(ptcl)::pj
      integer  ntp
      integer  ia, iz
      type(ptcl)::a(*)
      type(ptcl)::b(5)
      integer k, i
      integer MONIOU
      COMMON /area43/ moniou
            
      INTEGER          NSP
      COMMON /area12/  nsp
      INTEGER          ICH(95000)
      DOUBLE PRECISION ESP(4,95000)
      COMMON /area14/  esp,ich
      integer  icp, iap, iat, kicp
      integer  j
      real*8  E0
      moniou = 0

      call psasetc 
!      call PSASET    !  withDiffcode version  
      call xxaset
      call psaini 
!      call QGSPSAINI !  withDiffcode version
      return

!     ****************
      entry qgs01event(pj, ia, iz, a, ntp)
!     ****************

      if( pj%code .ne. kgnuc ) then
         iap = 1               !set proj. mass number (1-for hadron)
      else
         iap = pj%subcode
      endif
      call ccoscode2qgs( pj, kicp )
      e0 = pj%fm%p(4)/iap    !to get energy per nucleon
      iat = ia
      icp=iabs(kicp)/2+1   !primary particle class
                            ! (1- pion, 2 - nucleon, 3 - kaon)

!        icp:  projectile class
!        iap:  projectile mass number
!        iat:  target mass number

      call xxaini(e0, icp, iap, iat)
      call psconf

      ntp = 0
      DO  j = 1, nsp
!
!   0 pi0    1 pi+   -1 pi-
!   2 p  -2 pbar   3 n    -3 nbar  
!   4 K+   -4 K-  5  K0s   -5 K0L
!   6 eta    10 Lambda  -10 LambdaBar
!
         ntp = ntp + 1
         a(ntp)%fm%p(1)=esp(4, j)
         a(ntp)%fm%p(2)=esp(3, j)
         a(ntp)%fm%p(3)=esp(2, j)
         a(ntp)%fm%p(4)=esp(1, j)
         call cqgscode2cos(ich(j), a(ntp))
      enddo
      call crot3mom(pj, a, ntp)
      end
      real*8 function  qsran(X)
      real*8  X                 !  not used
      real*8 u
      call rndc(u)
      qsran = u
      end