cfixMCSmodel.f

Go to the documentation of this file.

      subroutine cfixmcsmodel( aPtcl )
      use modmcscontrol
      implicit none
!
!        fix how to treat MCS  for this e-/e+
!
#include "Zptcl.h"
#include "Zcode.h"

      type(ptcl):: aPtcl ! input. e+/e-

      integer i
      integer :: icon

      if( aptcl%code /= kelec) then
         write(0,*) 'cfixMSmodel is called non e-/e+'
         write(0,*) ' the code,subcode,charge= ', 
     *    aptcl%code, aptcl%subcode, aptcl%charge
         stop 
      endif
      kegev = aptcl%fm%p(4)- aptcl%mass
      keev = kegev*1.0d9 
      if(kegev > maxergmcs) then
         activemcs = 'Mol'
      else
         do i = 1, noofhowmcs
            if( kegev <  mcserg(i)) then
               if(kegev > minergmcs) then
                  activemcs = howmcslist(i)
               else
                  activemcs = 'Mol'
               endif
!!!!!!!!!!!!!
!!               write(0,*) 'E=',KEeV,' model=',ActiveMCS
!!!!!!!!!!!
               return           !!!!!!!!!!
            endif
         enddo
         call cerrormsg('MCSModel and energy range  strange', 1)
         do i = 1,  noofhowmcs
            write(0,*) howmcslist(i), mcserg(i)
         enddo
         stop 44444
      endif
      end