getImpactParam.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	getimpactparam (bout)
subroutine	geteposb (bout)
subroutine	getsibyllb (bout)

Function/Subroutine Documentation

geteposb()

subroutine geteposb

(

real(8), intent(out)

bout

)

Definition at line 41 of file getImpactParam.f.

Referenced by getimpactparam().

      implicit none
      real(8),intent(out):: bout
!---------------------------------------------------------------------------
!                   epos event common block
!---------------------------------------------------------------------------

      real        phievt,bimevt,pmxevt,egyevt
     *,xbjevt,qsqevt,zppevt,zptevt
      integer     nevt,kolevt,koievt,npjevt
     *,ntgevt,npnevt,nppevt,ntnevt,ntpevt,jpnevt,jppevt,jtnevt,jtpevt
     *,nglevt,minfra,maxfra
      common/cevt/phievt,nevt,bimevt,kolevt,koievt,pmxevt,egyevt,npjevt
     *,ntgevt,npnevt,nppevt,ntnevt,ntpevt,jpnevt,jppevt,jtnevt,jtpevt
     *,xbjevt,qsqevt,nglevt,zppevt,zptevt,minfra,maxfra

      bout = bimevt

Here is the caller graph for this function:

getimpactparam()

subroutine getimpactparam

(

real(8), intent(out)

bout

)

Definition at line 3 of file getImpactParam.f.

References geteposb(), and getsibyllb().

Referenced by gencol().

      implicit none
#include  "Zptcl.h"
#include  "Zcode.h"
#include  "Zevhnp.h"
#include  "Zevhnv.h"

      real(8),intent(out):: bout  

      real(8):: rproj, rtarg,bimpac
      integer::       nwtsam,nwasam,nwbsam,nwtacc,nwaacc,nwbacc
      COMMON /dtglcp/ rproj,rtarg,bimpac,
     &                nwtsam,nwasam,nwbsam,nwtacc,nwaacc,nwbacc
      integer,save::count=0

!       for qgsjet-II/I
      integer,parameter:: iapmax=209
      real(8):: xa, xb, b
      common /qgarr7/  xa(iapmax,3),xb(iapmax,3),b

      if( activemdl == "dpmjet3" ) then
         bout = bimpac
      elseif( activemdl(1:6) == "qgsjet" ) then
         bout = b
      elseif( activemdl == "epos" ) then
         call geteposb(bout) 
      elseif (activemdl == "sibyll" ) then
         call getsibyllb(bout)
      else
         if( count < 10 ) then
            count =count + 1
            write(0,*)
     *      ' IntModel =', activemdl, ' not ready for getImpactParam' 
            write(0,*) 'so  -1.0 is returned'
         endif
         bout = -1.0
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

getsibyllb()

subroutine getsibyllb

(

real(8), intent(out)

bout

)

Definition at line 59 of file getImpactParam.f.

References parameter().

Referenced by getimpactparam().

      implicit none
      real(8),intent(out):: bout

      integer:: iamax, ntry, na
      real:: b, bmax
      parameter(iamax=56)
      COMMON /s_cncm0/ b, bmax, ntry, na
      bout = b

Here is the call graph for this function:

Here is the caller graph for this function: