C++ Interface to Tauola
phodem.F
1C.----------------------------------------------------------------------
2C.
3C. PHOTOS: PHOton radiation in decays TeST program
4C.
5C. Purpose: Example of application of PHOTOS.
6C.
7C. Input Parameters: None
8C.
9C. Output Parameters: None
10C.
11C. Author(s): B. van Eijk, E. Barberio Created at: 31/05/90
12C. Last Update: 05/06/90
13C.
14C.----------------------------------------------------------------------
15 PROGRAM photst
16C IMPLICIT NONE
17 INTEGER EVENT,NHEP0
18#include "../../include/HEPEVT.h"
19 INTEGER PHLUN
20 common/pholun/phlun
21C--
22C-- Initialise PHOTOS
23 CALL phoini
24C--
25C-- Loop over JETSET event until PHOTOS has generated one or more pho-
26C-- tons. Do this for 10 JETSET events. The event record is printed
27C-- before and after photon emission.
28 DO 20 event=1,1
29 CALL lueevt(4,91.)
30C--
31C-- Conversion to /HEPEVT/ standard
32 CALL luhepc(1)
33C--
34C-- Write event record before emission...
35 nevhep=event
36 CALL phodmp
37 nhep0=nhep
38C--
39C-- Generate photon(s)... Arbitrary enforced generation.
40C-- Normally line: IF (NHEP.EQ.NHEP0) GOTO 10 must be absent!
41 10 CALL photos(1)
42 IF (nhep.EQ.nhep0) GOTO 10
43C--
44C-- Write event record...
45 WRITE(phlun,9050)
46 WRITE(phlun,9040)
47 CALL phodmp
48 20 CONTINUE
49 WRITE(phlun,9000)
50 WRITE(phlun,9010)
51 WRITE(phlun,9020)
52 WRITE(phlun,9030)
53 WRITE(phlun,9020)
54 WRITE(phlun,9010)
55 stop
56 9000 FORMAT(1h1)
57 9010 FORMAT(1h ,80('*'))
58 9020 FORMAT(1h ,'*',78x,'*')
59 9030 FORMAT(1h ,'**** PHOTOS Test Run has successfully ended',32x,
60 &' ****')
61 9040 FORMAT(1h ,26x,'=== after PHOTOS: ===')
62 9050 FORMAT(1h0,80('='))
63 END