TEvtGen/PHOTOS/photos_get.F

   1
   2       SUBROUTINE  PHOTOS_GET
   3 C.----------------------------------------------------------------------
   4 C.
   5 C.    Getter for PHOTOS:
   6 C.
   7 C.    Purpose:  Copies /HEPEVT/ into /PH_HEPEVT/
   8 C.
   9 C.
  10 C.    Input Parameters:   None
  11 C.
  12 C.    Output Parameters:  None
  13 C.
  14 C.    Author(s):  Z. Was                          Created at:  21/07/98
  15 C.                                                Last Update: 21/07/98
  16 C.
  17 C.----------------------------------------------------------------------
  18
  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          nevhep=  d_h_nevhep             ! serial number
  48          nhep  =  d_h_nhep               ! number of particles
  49          DO K=1,nhep
  50            isthep(k)    =d_h_isthep(k)   ! status code
  51            idhep(k)     =d_h_idhep(k)    ! particle ident KF
  52            jmohep(1,k)  =d_h_jmohep(1,k) ! parent particles
  53            jdahep(1,k)  =d_h_jdahep(1,k) ! childreen particles
  54            jmohep(2,k)  =d_h_jmohep(2,k) ! parent particles
  55            jdahep(2,k)  =d_h_jdahep(2,k) ! childreen particles
  56            DO l=1,4
  57            phep(l,k)    =d_h_phep(l,k)   ! four-momentum, mass [GeV]
  58            vhep(l,k)    =d_h_vhep(l,k)   ! vertex [mm]
  59            ENDDO
  60            phep(5,k)    =d_h_phep(5,k)   ! four-momentum, mass [GeV]
  61          ENDDO
  62       END