TEvtGen/PHOTOS/photos_set.F

   1
   2
   3       SUBROUTINE PHOTOS_SET
   4 C.----------------------------------------------------------------------
   5 C.
   6 C.    Setter for PHOTOS:
   7 C.
   8 C.    Purpose:  Copies /PH_HEPEVT/ into /HEPEVT/
   9 C.
  10 C.
  11 C.    Input Parameters:   None
  12 C.
  13 C.    Output Parameters:  None
  14 C.
  15 C.    Author(s):  Z. Was                          Created at:  21/07/98
  16 C.                                                Last Update: 21/07/98
  17 C.
  18 C.----------------------------------------------------------------------
  19       IMPLICIT NONE
  20       INTEGER  d_h_nmxhep         ! maximum number of particles
  21       PARAMETER ( d_h_nmxhep=4000)
  22       REAL*8  d_h_phep,  d_h_vhep ! to be real*4 or *8  depending on host
  23       INTEGER d_h_nevhep,d_h_nhep,d_h_isthep,d_h_idhep,d_h_jmohep,
  24      $        d_h_jdahep
  25       COMMON /hepevt/
  26      $      d_h_nevhep,               ! serial number
  27      $      d_h_nhep,                 ! number of particles
  28      $      d_h_isthep(d_h_nmxhep),   ! status code
  29      $      d_h_idhep(d_h_nmxhep),    ! particle ident KF
  30      $      d_h_jmohep(2,d_h_nmxhep), ! parent particles
  31      $      d_h_jdahep(2,d_h_nmxhep), ! childreen particles
  32      $      d_h_phep(5,d_h_nmxhep),   ! four-momentum, mass [GeV]
  33      $      d_h_vhep(4,d_h_nmxhep)    ! vertex [mm]
  34 * ----------------------------------------------------------------------
  35 !      LOGICAL d_h_qedrad
  36 !      COMMON /phoqed/
  37 !     $     d_h_qedrad(d_h_nmxhep)    ! Photos flag
  38       INTEGER NMXHEP
  39       PARAMETER (NMXHEP=10000)
  40       INTEGER IDHEP,ISTHEP,JDAHEP,JMOHEP,NEVHEP,NHEP
  41       REAL*8 PHEP,VHEP
  42       COMMON/PH_HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
  43      &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
  44       LOGICAL QEDRAD
  45       COMMON/PHOQED/QEDRAD(NMXHEP)
  46       INTEGER K,L
  47
  48          d_h_nevhep=  nevhep             ! serial number
  49          d_h_nhep  =  nhep               ! number of particles
  50          DO K=1,nhep
  51            d_h_isthep(k)    =isthep(k)   ! status code
  52            d_h_idhep(k)     =idhep(k)    ! particle ident KF
  53            d_h_jmohep(1,k)  =jmohep(1,k) ! parent particles
  54            d_h_jdahep(1,k)  =jdahep(1,k) ! childreen particles
  55            d_h_jmohep(2,k)  =jmohep(2,k) ! parent particles
  56            d_h_jdahep(2,k)  =jdahep(2,k) ! childreen particles
  57            DO l=1,4
  58            d_h_phep(l,k)    =phep(l,k)   ! four-momentum, mass [GeV]
  59            d_h_vhep(l,k)    =vhep(l,k)   ! vertex [mm]
  60            ENDDO
  61            d_h_phep(5,k)    =phep(5,k)   ! four-momentum, mass [GeV]
  62          ENDDO
  63       END