cfixModel.f

Go to the documentation of this file.

      subroutine cfixmodel( aPtcl )
      implicit none
!
!        fix hadronic interaction model (including 
!        photoproduction of hadrons, and possibly
!        muon nucear interaction).
!        and set ActiveMdl
!
#include "Zptcl.h"
#include "Zevhnp.h"
#include "Zevhnv.h"
#include "Zcode.h"

      type(ptcl)::aPtcl ! input. particle for hadronic interaction

      real*8 kepn 
      integer i
      integer :: icon

      kepn= aptcl%fm%p(4)- aptcl%mass
      if( aptcl%code .eq. kgnuc ) then
         kepn =  kepn / aptcl%subcode
      endif

      do i = 1, noofmdls
         if( kepn  .lt. inteerg(i)) then
            activemdl = modellist(i)
            goto 10
         endif
      enddo

      call cerrormsg('IntModel list shown below is strange', 1)
      call cerrormsg(intmodel, 1)
      call cpdpmjetinp
      write(0,*) ' current K%E/n=',kepn,
     *  ' code=',aptcl%code, ' subcode=',
     *  aptcl%subcode, ' charge=',aptcl%charge
      write(0,*) '  NoOfMdls=', noofmdls
      do i = 1,  noofmdls
         write(0,*) ' i=', i, ' InteErg(i)=',
     *   inteerg(i)
      enddo
      stop 12345
 10   continue
      do i = 1, noofmdls2
         if( kepn  .lt. inteerg2(i)) then
            activemdl2 = modellist2(i) ! used only for mfp calc
            goto 15
         endif
      enddo

      call cerrormsg('XsecModel list shown below is strange', 1)
      call cerrormsg(xsecmodel, 1)
      call cpdpmjetinp
      write(0,*) ' current K%E/n=',kepn,
     *  ' code=',aptcl%code, ' subcode=',
     *  aptcl%subcode, ' charge=',aptcl%charge
      write(0,*) '  NoOfMdls2=', noofmdls2
      do i = 1,  noofmdls2
         write(0,*) ' i=', i, ' InteErg2(i)=',
     *   inteerg2(i)
      enddo
!!      stop 123456
      stop 23456    !! jaxa length should be <=5
 15   continue

      if( aptcl%code == knuc .and.  aptcl%subcode == antip ) then
         if(activemdl == "phits" ) then
          !  ' for p-bar or n-bar, phits cannot be used'
            call eprescue(icon)
            if( icon /= 0 ) then
               call eprescueerr('phits cannot be used for p-bar', aptcl)
               stop
            endif
         endif
      elseif(aptcl%code == kkaon ) then
         if(activemdl == "phits" ) then
            call eprescue(icon)
            if( icon /= 0 ) then
               call eprescueerr('phits cannot be used for Kaon', aptcl)
               stop
            endif
         endif
      endif
      ! If photon or muon, they will be 
         !  managed in photon interaction routine
      ! other hadronic particles,say, lambda0, will be decayed
      ! if the model cannot cope with it.  (cseeColPossible).
 20   continue      
      end
      subroutine eprescue(icon)
      implicit none
!
!        fix hadronic interaction model (including 
!        photoproduction of hadrons, and possibly
!        muon nucear interaction).
!        and set ActiveMdl
!
#include "Zptcl.h"
#include "Zevhnp.h"
#include "Zevhnv.h"
#include "Zcode.h"
      integer,intent(out):: icon  ! 0.  rescuee could be set
                           !    1.    no resucee could be fount
!             see if jam or dpmjet3 can be used, if not error
      integer::i

      icon  = 0
      do i = 1, noofmdls
         if( modellist(i) == "dpmjet3") then
            activemdl= "dpmjet3"
            goto 20
         elseif(modellist(i) == "jam") then
            activemdl= "jam"
            goto 20
         elseif( modellist(i) == "nucrin") then
            activemdl= "nucrin"
            goto 20
         endif
      enddo
      icon = 1
 20   continue
      end

      subroutine eprescueerr(msg, aPtcl)
      implicit none
#include "Zptcl.h"
      character(*),intent(in):: msg
      type(ptcl)::aPtcl !  input. problematic  ptcl

      write(0,'(a)')  trim(msg),
     * ' so that dpmjet3 or jam  (or nucrin)  must be given in ',
     * ' the IntModel list as a rescuee: ',
     * " If you don't want to use such a rescuee as a normal ",
     * ' interaction model, you may specify it as follows:',
     * " IntModel='"//'"abc" E1 "dpmjet3" E2 "xyz"',
     * ' and set E1=E2, then, model abc is used below E1=E2. ',
     * ' When energy is >= E1, model xyz is used. So dpmjet3 will',
     * ' not be used at any energy. It will be used only as ',
     * ' rescuee.',
     * ' ',
     * 'particle info: '
      write(0,'(a, g13.4)') ' ptcl code=',aptcl%code,
     *               '  charge  =',aptcl%charge,
     *               '  subcode =',aptcl%subcode,
     *              '  Ek      =', aptcl%fm%p(4)-aptcl%mass, ' GeV'
      end
!       
      subroutine cqactivemdl(model)
      implicit none
!
!        inquire the current active model
!
#include "Zptcl.h"
#include "Zevhnv.h"

      character*16 model
      model = activemdl
      end
      subroutine cqactivemdl2(model)
      implicit none
#include "Zptcl.h"
#include "Zevhnv.h"

      character*16 model
      model = activemdl2
      end