chAcol.f

Go to the documentation of this file.

!         hadoron-A collision.  generic.
      subroutine chacol(pj, ia, iz, a, ntp)
      implicit none
#include  "Zcondc.h"
#include  "Zcode.h"
#include  "Zptcl.h"
#include  "Zevhnp.h"
#include  "Zevhnv.h"
#include  "Zspecial.h"
!
      type(ptcl):: pj   ! input  projectile ptcl
      integer ia    ! input. mass no. of target
      integer iz    ! input. charge no. of target
      type(ptcl):: a(*)  !  output. produced ptcls
      integer ntp   ! number of produced ptcls

      integer nc
      integer retry
      data nc/0/
      save nc
      retry = 1
      do while( retry .ne. 0) 
         call chacolx(pj, ia, iz, a, ntp)
         if(modifyx .ne. 0) then
            call ccheckx(pj, ia, iz, a, ntp, retry)
            nc = nc + 1
            if(nc .lt. 20) then
               write(0,*)
     *         '************* X modification is requested !!!!!!!!!'
            endif
         else
            retry = 0
         endif
      enddo
      end
      subroutine ccheckx(pj, ia, iz, a, ntp, retry)
      implicit none
#include  "Zcondc.h"
#include  "Zcode.h"
#include  "Zptcl.h"
#include  "Zevhnp.h"
#include  "Zevhnv.h"
#include  "Zspecial.h"
!
      type(ptcl):: pj   ! input  projectile ptcl
      integer ia    ! input. mass no. of target
      integer iz    ! input. charge no. of target
      type(ptcl):: a(*)  !  output. produced ptcls
      integer ntp   ! number of produced ptcls
      integer retry  !  output. 0 if event is acceptable

      integer i
      real*8 u, ePimax, xPimax, eNmax, xNmax, xxx

      if(pj%code .eq. knuc) then
!         find max nucleon
         epimax = 0.
         enmax = 0.
         do i = 1, ntp
            if( a(i)%code .eq. knuc  )then
               enmax = max(enmax, a(i)%fm%p(4) )
            elseif( a(i)%code .eq. kpion  )then
               epimax = max(epimax, a(i)%fm%p(4) )
            endif
         enddo
         xpimax = epimax/pj%fm%p(4)
         xnmax = enmax/pj%fm%p(4)
         call rndc(u)
         if(xnmax .gt. 0.97) then
            retry = 0
         elseif(modifyx .eq. 1) then
!            drop  high Pi x; pw=3
            if(u .lt. (1.0 - xpimax)**modifyxpw1) then
               retry = 0
            else
               retry = 1
            endif
         elseif(modifyx .eq. 2) then
!              drop high Nx case; pw=3
            if( u .lt. ((2.-xnmax)/2.)**modifyxpw1 ) then
               retry = 0
            else
               retry = 1
            endif
         elseif(modifyx .eq. 3) then
!               drop large xPimax*xNmax ; pw =3
            xxx = xpimax* xnmax
            if(xxx .lt. 0.005) then
               retry = 0
            elseif( xxx .gt. 0.08 ) then
               retry = 1
            elseif(u .lt. ((0.08-xxx)/(0.08-0.005))**modifyxpw1 ) then
               retry = 0
            else
               retry = 1
            endif
         elseif(modifyx .eq. 4) then
            xxx = xpimax/xnmax
            if(xxx .gt. 2.5) then
               retry = 1
            elseif(xxx .lt. 0.2) then
               retry = 0
            elseif(u .lt. ((2.5-xxx)/(2.5-0.2))**modifyxpw1) then
               retry = 0
            else
               retry = 1
            endif
         elseif(modifyx .eq. 5) then
            xxx=xpimax/xnmax
            if(xxx .gt. 2.)  then
               retry = 0
            elseif( u .gt.  (1.-xxx/2.)**modifyxpw1) then
               retry = 0
            else
               retry = 1
            endif
         elseif(modifyx .eq. 6 ) then
!             drop small xxx
            xxx = xpimax* xnmax
            if(xxx .gt. 0.05) then
               retry = 0
            elseif( u .lt. (1.-xxx/0.05)**modifyxpw1 ) then
               retry = 0
            else
               retry = 1
            endif
         else
            retry = 0
         endif
      else
         retry = 0
      endif
      end


      subroutine chacolx(pj, ia, iz, a, ntp)
!      if only.. is omitted "both" in here and cheavyInt.f
!      at link time of large program, a message something like
!      TargetXs for x86-64 is missing
!      comes out and stops. 
#if !defined (MacIFC)
      use modxsecmedia
#endif
      implicit none
#include  "Zcondc.h"
#include  "Zcode.h"
#include  "Zptcl.h"
#include  "Zevhnp.h"
#include  "Zevhnv.h"
#include  "Zair.h"

!
      type(ptcl):: pj   ! input  projectile ptcl
      integer ia    ! input. mass no. of target
      integer iz    ! input. charge no. of target
      type(ptcl):: a(*)  !  output. produced ptcls
      integer ntp   ! number of produced ptcls

#if defined (MacIFC)
#include "Zworkaround.h"
#endif


!
!c      if(IntModel .eq. 'int1') then
!c        if(pj.fm.p(4)  .lt. Elund) then
!            h-A collision  Fermi momentum is considered inside
      if( activemdl .eq. 'phits') then
         call cphits(pj, ia, iz, targetxs, a, ntp)
      elseif( activemdl .eq. 'dpmjet3') then
         call cdpmjet(pj, ia, iz, a, ntp)
      elseif( activemdl .eq. 'jam') then
!                  TargetXs is the xs on 'ia'. This is used 
!            in jam to get bmax.  better to include it in 
!            chAcol argment in future ?
         call cjamevent(pj, ia, iz, targetxs, a, ntp)
      elseif( activemdl .eq. 'qgsjet2') then
         call cqgsjet(pj, ia, iz, a, ntp)
      elseif( activemdl == 'epos') then
         call ceposgenoneevent(pj, ia, iz, a, ntp)
      elseif( activemdl == 'sibyll') then
         call csibyllevent(pj, ia, iz, a, ntp)
      elseif( activemdl .eq. 'fritiof1.6' .or. 
     *        activemdl .eq. 'nucrin') then
            call chalund(pj, ia, iz, a, ntp)
      elseif( activemdl .eq. 'ad-hoc') then
         call chacoladhoc(pj, ia, iz, a, ntp)
!      elseif( ActiveMdl .eq. 'fritiof7.02') then
!         if(pj.code .eq. kkaon .and. LundPara(5) .ne. 0) then
!            call chAcolAdhoc(pj, ia, iz, a, ntp)
!         elseif(pj.code .eq. keta) then
!            call chAcolAdhoc(pj, ia, iz, a, ntp)
!         else
!            call chANewLund(pj, ia, iz, a, ntp)
!         endif
      elseif(activemdl .eq. 'gheisha') then
!         if(pj.fm.p(4) .lt. Elund) then
         call chagheisha(pj, ia, iz, a, ntp)
!         elseif(pj.fm.p(4)  .lt. Elund2) then
!            call chAcolAdhoc(pj, ia, iz, a, ntp)
!         elseif(pj.fm.p(4)  .lt. Elund3) then
!            if(pj.code .eq. kkaon .and. LundPara(5) .ne. 0) then
!               call chAcolAdhoc(pj, ia, iz, a, ntp)
!            elseif(pj.code .eq. keta) then
!               call chAcolAdhoc(pj, ia, iz, a, ntp)
!            else
!               call chANewLund(pj, ia, iz, a, ntp)
!            endif
      elseif(activemdl .eq. 'incdpm3') then
!             dpmjet inclusive treatment
         call cincdpm3(pj, ia, iz, a, ntp)
      elseif(activemdl .eq. 'byenergy') then
         if( pj%fm%p(4)-pj%mass .gt. 4.0) then
            call chacoladhoc(pj, ia, iz, a, ntp)
         else
            call chalund(pj, ia, iz, a, ntp)
         endif
      else
         call cerrormsg(activemdl, 1)
         call cerrormsg('above model not supported',0)
      endif

      end