phodem.F

00001 C.----------------------------------------------------------------------
00002 C.
00003 C.    PHOTOS:   PHOton radiation in decays TeST program
00004 C.
00005 C.    Purpose:  Example of application of PHOTOS.
00006 C.
00007 C.    Input Parameters:   None
00008 C.
00009 C.    Output Parameters:  None
00010 C.
00011 C.    Author(s):  B. van Eijk, E. Barberio        Created at:  31/05/90
00012 C.                                                Last Update: 05/06/90
00013 C.
00014 C.----------------------------------------------------------------------
00015       PROGRAM PHOTST
00016 C      IMPLICIT NONE
00017       INTEGER EVENT,NHEP0
00018 #include "../../include/HEPEVT.h"
00019       INTEGER PHLUN
00020       COMMON/PHOLUN/PHLUN
00021 C--
00022 C--   Initialise PHOTOS
00023       CALL PHOINI
00024 C--
00025 C--   Loop over JETSET event until PHOTOS has generated one or more pho-
00026 C--   tons.  Do this  for 10 JETSET events.  The event record is printed
00027 C--   before and after photon emission.
00028       DO 20 EVENT=1,1 
00029         CALL LUEEVT(4,91.)
00030 C--
00031 C--   Conversion to /HEPEVT/ standard
00032         CALL LUHEPC(1)
00033 C--
00034 C--   Write event record before emission...
00035         NEVHEP=EVENT
00036         CALL PHODMP
00037         NHEP0=NHEP
00038 C--
00039 C--   Generate photon(s)... Arbitrary enforced generation.
00040 C--   Normally line: IF (NHEP.EQ.NHEP0) GOTO 10 must be absent!
00041    10   CALL PHOTOS(1)
00042         IF (NHEP.EQ.NHEP0) GOTO 10
00043 C--
00044 C--   Write event record...
00045       WRITE(PHLUN,9050)
00046       WRITE(PHLUN,9040)
00047         CALL PHODMP
00048    20 CONTINUE
00049       WRITE(PHLUN,9000)
00050       WRITE(PHLUN,9010)
00051       WRITE(PHLUN,9020)
00052       WRITE(PHLUN,9030)
00053       WRITE(PHLUN,9020)
00054       WRITE(PHLUN,9010)
00055       STOP
00056  9000 FORMAT(1H1)
00057  9010 FORMAT(1H ,80('*'))
00058  9020 FORMAT(1H ,'*',78X,'*')
00059  9030 FORMAT(1H ,'**** PHOTOS Test Run has successfully ended',32X,
00060      &' ****')
00061  9040 FORMAT(1H ,26X,'=== after PHOTOS: ===')
00062  9050 FORMAT(1H0,80('='))
00063       END