C++InterfacetoTauola
phodem.F
1 C.----------------------------------------------------------------------
2 C.
3 C. PHOTOS: PHOton radiation in decays TeST program
4 C.
5 C. Purpose: Example of application of PHOTOS.
6 C.
7 C. Input Parameters: None
8 C.
9 C. Output Parameters: None
10 C.
11 C. Author(s): B. van Eijk, E. Barberio Created at: 31/05/90
12 C. Last Update: 05/06/90
13 C.
14 C.----------------------------------------------------------------------
15  PROGRAM photst
16 C IMPLICIT NONE
17  INTEGER event,nhep0
18 #include "../../include/HEPEVT.h"
19  INTEGER phlun
20  common/pholun/phlun
21 C--
22 C-- Initialise PHOTOS
23  CALL phoini
24 C--
25 C-- Loop over JETSET event until PHOTOS has generated one or more pho-
26 C-- tons. Do this for 10 JETSET events. The event record is printed
27 C-- before and after photon emission.
28  DO 20 event=1,1
29  CALL lueevt(4,91.)
30 C--
31 C-- Conversion to /HEPEVT/ standard
32  CALL luhepc(1)
33 C--
34 C-- Write event record before emission...
35  nevhep=event
36  CALL phodmp
37  nhep0=nhep
38 C--
39 C-- Generate photon(s)... Arbitrary enforced generation.
40 C-- Normally line: IF (NHEP.EQ.NHEP0) GOTO 10 must be absent!
41  10 CALL photos(1)
42  IF (nhep.EQ.nhep0) goto 10
43 C--
44 C-- 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