chncol.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	chncol (pj, tg, a, ntp, icon)
subroutine	chncol2 (pj, tg, a, ntp, icon)
subroutine	cdcycp (a, nin, n)
subroutine	cpikcd (a, ntp)
subroutine	cnnbdc (a, b, n)
subroutine	cddbdc (a, b, n)
subroutine	cspipm (pj, x, a)
subroutine	cskchg (pj, x, a)

Function/Subroutine Documentation

cdcycp()

subroutine cdcycp	(	type(ptcl), dimension(*)	a,
		integer	nin,
		integer	n
	)

Definition at line 158 of file chncol.f.

References cddbdc(), cnnbdc(), kddb, and knnb.

Referenced by chncol2().

!      ****************************************************
       implicit none
!----       include  '../../Zptcl.h'
#include  "Zptcl.h"
!----       include '../../Zcode.h'
#include  "Zcode.h"
       integer nin, n
       type(ptcl):: a(*)
       integer i, k, nx
       type(ptcl):: p(2)
          n=nin
          do i=1, nin
             k=a(i)%code
             if(k .eq. knnb .or. k .eq. kddb) then
                 if(k .eq. knnb) then
                     call cnnbdc(a(i), p, nx)
                 elseif(k .eq. kddb) then
                     call cddbdc(a(i), p, nx)
                 endif
!                 put n or D in the i-th pos. and append c antiptcl 
!                 at the bottome and increase n
                 a(i) = p(1)
                 a(n +1 ) = p(2)
                 n=n+1
             endif    
          enddo

Here is the call graph for this function:

Here is the caller graph for this function:

cddbdc()

subroutine cddbdc	(	type(ptcl)	a,
		type(ptcl), dimension(2)	b,
		integer	n
	)

Definition at line 259 of file chncol.f.

References c2bdcy(), cmkptc(), kd0, kd0b, kdmes, and rndc().

Referenced by cdcycp(), and cinteddb().

      implicit none
!----       include '../../Zptcl.h'
#include  "Zptcl.h"
!----       include '../../Zcode.h'
#include  "Zcode.h"
       integer n
       type(ptcl):: a, b(2) 
!
       integer ic
       real*8 u 
!
         call rndc(u)
         ic=int(u*3)-1 ! charge
         if(ic .eq. 0) then
             call cmkptc(kdmes,  kd0, ic,  b(1))
             call cmkptc(kdmes,  kd0b, -ic, b(2))
         else    
             call cmkptc(kdmes,  0, ic,  b(1))
             call cmkptc(kdmes,  0, -ic,  b(2))
         endif
         call c2bdcy(a, b(1), b(2))
         n=2

Here is the call graph for this function:

Here is the caller graph for this function:

chncol()

subroutine chncol	(	type(ptcl)	pj,
		type(ptcl)	tg,
		type(ptcl), dimension(*)	a,
		integer	ntp,
		integer	icon
	)

Definition at line 11 of file chncol.f.

References chncol2(), rndc(), and true.

Referenced by chacoladhoc().

!             pj:  type ptcl. Input. to give the incident
!     projectile hadron in the Lab. system.
!             tg:  type ptcl. Input. to give the target
!                  nucleon in the Lab. system.
!              a:  array of type ptcl. Output. to get
!                  produced particles.
!            ntp:  integer. Output. to get the total number of produced
!                  particles in a.
!           icon:  Integer. Output. if 0, particle generation is ok
!                   if non 0, you may give up generation.
!
!                 Note: target recoil is put in a(ntp-1),
!                       projectile recoil is put in a(ntp)
!
       implicit none

#include  "Zptcl.h"
#include  "Zevhnv.h"

!
       type(ptcl):: pj, tg,  a(*)
       integer ntp, icon
       integer i
       real*8 u, xmax

!       real*8 BigXRejCnst/1./, BigXRejPw/2./
!         default
        real*8 bigxrejcnst/.4/, bigxrejpw/2.2/
!         enhance
!        real*8 BigXRejCnst/.8/, BigXRejPw/2.2/
!         denhance
!       real*8 BigXRejCnst/.01/, BigXRejPw/1./
!       real*8 BigXRejCnst/0.0003/, BigXRejPw/3.5/

       do while (.true.)
          call chncol2(pj, tg, a, ntp, icon)
          if(icon .ne. 0) goto 100
!             find max secondary energy.
          xmax = 0.
          do i = 1, ntp-2
             if(xmax .lt. a(i)%fm%p(4)) xmax = a(i)%fm%p(4)
          enddo
          xmax = xmax/pj%fm%p(4)
!              if very large x appears, discard it
!              with some probabilty to adjust
!              the x-distribution.
!
!c          if(xmax .gt. 0.4) then
!c             call rndc(u)
!c             if(u .lt.  (1. - (xmax-0.40)/0.6 )**2 ) then
!c                goto 100
!c             endif
!c          else
!c             goto 100
!c          endif
          call rndc(u)
          if(u  .lt. 
     *         (  bigxrejcnst/(bigxrejcnst +
     *         (xmax/(1.0-xmax))**bigxrejpw) ) )  goto 100

       enddo
 100   continue

Here is the call graph for this function:

Here is the caller graph for this function:

chncol2()

subroutine chncol2	(	type(ptcl)	pj,
		type(ptcl)	tg,
		type(ptcl), dimension(*)	a,
		integer	ntp,
		integer	icon
	)

Definition at line 84 of file chncol.f.

References cconsvchg(), cdcycp(), cgnlp(), cibst1(), cpikcd(), crot3mom(), cs2lp(), and true.

Referenced by chncol().

!             pj:  type ptcl. Input. to give the incident
!                  projectile hadron in the Lab. system.
!             tg:  type ptcl. Input. to give the target
!                  nucleon in the Lab. system.
!              a:  array of type ptcl. Output. to get
!                  produced particles.
!            ntp:  integer. Output. to get the total number of produced
!                  particles in a.
!           icon:  Integer. Output. if 0, particle generation is ok
!                   if non 0, you may give up generation.
!
!                 Note: target recoil is put in a(ntp-1),
!                       projectile recoil is put in a(ntp)
!
       implicit none

#include  "Zptcl.h"
#include  "Zevhnv.h"
!
!
       type(ptcl):: pj, tg,  a(*)
       integer ntp, icon, jcon, nfin
       integer i, outc
       type(ptcl):: pjin
!         use pt=0 incident
      pjin = pj
      pjin%fm%p(1) = 0.
      pjin%fm%p(2) = 0.
      pjin%fm%p(3) = sqrt(pjin%fm%p(4)**2 - pjin%mass**2)
!

!        *** until loop**until succeed
      do while (.true.)
!            generate 2 leading ptcls  Zevhnv become ready for use.
          call cs2lp(pjin, tg, icon)
          if(icon .ne. 0) then
!              neglect this event; because of very low energy
!              no particle can be produced
               ntp = 0
              jcon = 0
          else
!               generation of  particles other than the leadings.
              call cgnlp(a, ntp, jcon)
          endif

          if(jcon .eq. 0 ) goto 10
       enddo   
   10 continue
!           ptcls in 'a' should have cms energy, here.
!           give final ptcl charge/ subcode for 
!           pi+-,  k0,k0~,k+,k-
      call cpikcd(a, ntp)
      outc =rtglab%charge+rpjlab%charge - pj%charge- tg%charge
      call cconsvchg(outc, a, ntp, icon)
      if(icon .eq. 0) then
!              boost  to lab system
         do   i=1, ntp
            call cibst1(i, cmsp, a(i), a(i))
         enddo
!          decay of "composit ptcls" (nn~, dd~)
         call cdcycp(a, ntp, nfin)
         ntp=nfin
!            store recoils in a
         a(ntp+1) = rtglab
         a(ntp+2) = rpjlab
         ntp = ntp +2
!     -------------- rotate 
         call crot3mom(pj, a, ntp)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cnnbdc()

subroutine cnnbdc	(	type(ptcl)	a,
		type(ptcl), dimension(2)	b,
		integer	n
	)

Definition at line 231 of file chncol.f.

References c2bdcy(), cmkptc(), kneutron, kneutronb, knuc, and rndc().

Referenced by cdcycp(), and cintennb().

      implicit none
!----       include '../../Zptcl.h'
#include  "Zptcl.h"
!----       include '../../Zcode.h'
#include  "Zcode.h"
       integer n
       type(ptcl):: a, b(2) 
!
       integer ic
       real*8 u 
!
         call rndc(u)
         ic=u*2         ! charge
         if(ic .eq. 0) then
             call cmkptc(knuc,  kneutron, ic,   b(1))
             call cmkptc(knuc,  kneutronb, -ic,  b(2))
         else    
             call cmkptc(knuc, 0, ic,  b(1))
             call cmkptc(knuc, 0, -ic, b(2))
         endif    
         call c2bdcy(a, b(1), b(2))
         n=2

Here is the call graph for this function:

Here is the caller graph for this function:

cpikcd()

subroutine cpikcd	(	type(ptcl), dimension(ntp)	a,
		integer	ntp
	)

Definition at line 188 of file chncol.f.

References cskchg(), cspipm(), kkaon, and kpion.

Referenced by chncol2().

!
!*******************************************************
      implicit none
!----      include '../../Zptcl.h'
#include  "Zptcl.h"
!----      include '../../Zcode.h'
#include  "Zcode.h"
!----      include '../Zevhnv.h'
#include  "Zevhnv.h"
      integer ntp
      type(ptcl)::  a(ntp) 
!
      integer i, k
      real*8 x
        if(ntp .eq. 1) then
!              nothing to do. charge assignment should have
!              been done in c1pion
        else
             do i=1, ntp
                x = a(i)%fm%p(3)/pjcms%fm%p(3)
                k = a(i)%code
                if(k .eq. kpion .and. a(i)%charge .ne. 0) then
                     if(x .gt. 0.) then
                         call cspipm(pjcms, x, a(i))
                     else
                         call cspipm(tgcms, -x, a(i))
                     endif
                elseif(k .eq. kkaon) then 
!                            set kaon charge
                     if(x .gt. 0.) then
                        call cskchg(pjcms, x, a(i))
                     else
                        call cskchg(tgcms, -x, a(i))
                     endif
                endif
             enddo
        endif

Here is the call graph for this function:

Here is the caller graph for this function:

cskchg()

subroutine cskchg	(	type(ptcl)	pj,
		real*8	x,
		type(ptcl)	a
	)

Definition at line 333 of file chncol.f.

References cmkptc(), d0, k0l, k0s, kkaon, regptcl, and rndc().

Referenced by cpikcd().

!     ************************************************
      implicit none

#include  "Zptcl.h"
#include  "Zcode.h"

      type(ptcl):: pj, a
      real*8 x
!
      real*8 f, u, up
       if(a%charge .ne. 0) then
!             k+/k- for p incident
           if(x .lt. .3d0) then
              f=exp(4.36d0*x)
           elseif(x .lt. .6d0) then
              f=3.7d0*exp(6.5d0*(x-.3d0))
           else
              f=27.d0*exp(11.3d0*(x-.6d0))
           endif
           if(pj%charge .eq. 1) then
              up=f/(1.d0+f)
           elseif(pj%charge .eq. -1) then
              up=1./(1.d0+f)
           else
              if(pj%subcode .eq. 0) then
                 up=0.5d0
              elseif(pj%subcode .eq. regptcl) then
                 up=1.d0/(1.d0+f)
              else
                 up=f/(1.d0+f)
              endif
           endif
           call rndc(u)
           if(u .lt. up) then
               call cmkptc(kkaon, 0, 1,  a)
           else
               call cmkptc(kkaon, 0,  -1,  a)
           endif
       else
!            k0
           call rndc(u)
           if(u .lt. .50d0) then
              call rndc(u)
              if(u .lt. 0.5) then
!                 k0 short
                 call cmkptc(kkaon, k0s, 0,  a)
              else
                 call cmkptc(kkaon, -k0s, 0,  a)
              endif
           else
              call rndc(u)
              if(u .lt. 0.5) then
                 call cmkptc(kkaon, k0l, 0, a)
              else
                 call cmkptc(kkaon, -k0l, 0, a)
              endif
           endif
       endif

Here is the call graph for this function:

Here is the caller graph for this function:

cspipm()

subroutine cspipm	(	type(ptcl)	pj,
		real*8	x,
		type(ptcl)	a
	)

Definition at line 292 of file chncol.f.

References antip, d0, regptcl, and rndc().

Referenced by cpikcd().

       implicit none

#include  "Zptcl.h"
#include  "Zcode.h"

       type(ptcl):: pj, a
       real*8 x
!
       real*8 f, up, u
!
         if(x .lt. .6d0) then
             f = 1.d0+3.33d0*x
         else
             f = 3.0d0* exp( 4.4d0*(x-.6d0))
         endif
         if(pj%charge .eq. 1) then
             up=f/(1.d0+f)
         elseif(pj%charge .eq. -1) then
             up=1.d0/(1.d0+f)
         else
             if(pj%subcode .eq. 0) then
                 up=0.5
             elseif(pj%subcode .eq. regptcl) then
                 up=1.d0/(1.d0+f)
             else
                 up=f/(1.d0+f)
             endif
         endif
         call rndc(u)
         if(u .lt. up) then
              a%charge = 1
              a%subcode = regptcl
         else
              a%charge = -1
              a%subcode = antip
         endif

Here is the call graph for this function:

Here is the caller graph for this function: