qgs01init.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	qgs01init
real *8 function	qsran (X)

Function/Subroutine Documentation

qgs01init()

subroutine qgs01init

(

)

Definition at line 2 of file qgs01init.f.

References a, b, crot3mom(), kgnuc, psaini(), psasetc(), psconf(), xxaini(), and xxaset().

Referenced by init().

      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)

Here is the call graph for this function:

Here is the caller graph for this function:

qsran()

real*8 function qsran

(

real*8

X

)

Definition at line 74 of file qgs01init.f.

References rndc().

Referenced by crossc_kk(), ixxdef(), pscajet(), psconf(), pscs(), psgea(), pshot(), psrec(), psshar(), psvdef(), pszsim(), xxddfr(), xxdec2(), xxdec3(), xxdpr(), xxdtg(), xxfragm(), xxgener(), and xxstr().

      real*8  x                 !  not used
      real*8 u
      call rndc(u)
      qsran = u

Here is the call graph for this function:

Here is the caller graph for this function: