]>
Commit | Line | Data |
---|---|---|
da0e9ce3 | 1 | SUBROUTINE PHOCIN |
2 | C.---------------------------------------------------------------------- | |
3 | C. | |
4 | C. PHOTOS: PHOton Common INitialisation | |
5 | C. | |
6 | C. Purpose: Initialisation of parameters in common blocks. | |
7 | C. | |
8 | C. Input Parameters: None | |
9 | C. | |
10 | C. Output Parameters: Commons /PHOLUN/, /PHOPHO/, /PHOCOP/, /PHPICO/ | |
11 | C. and /PHSEED/. | |
12 | C. | |
13 | C. Author(s): B. van Eijk Created at: 26/11/89 | |
14 | C. Z. Was Last Update: 10/08/93 | |
15 | C. | |
16 | C.---------------------------------------------------------------------- | |
17 | IMPLICIT NONE | |
18 | INTEGER NMXHEP | |
19 | PARAMETER (NMXHEP=10000) | |
20 | LOGICAL QEDRAD | |
21 | COMMON/PHOQED/QEDRAD(NMXHEP) | |
22 | INTEGER PHLUN | |
23 | COMMON/PHOLUN/PHLUN | |
24 | REAL*8 ALPHA,XPHCUT | |
25 | COMMON/PHOCOP/ALPHA,XPHCUT | |
26 | REAL*8 PI,TWOPI | |
27 | COMMON/PHPICO/PI,TWOPI | |
28 | INTEGER ISEED,I97,J97 | |
29 | REAL*8 URAN,CRAN,CDRAN,CMRAN | |
30 | COMMON/PHSEED/ISEED(2),I97,J97,URAN(97),CRAN,CDRAN,CMRAN | |
31 | INTEGER PHOMES | |
32 | PARAMETER (PHOMES=10) | |
33 | INTEGER STATUS | |
34 | COMMON/PHOSTA/STATUS(PHOMES) | |
35 | LOGICAL INTERF,ISEC,IFTOP | |
36 | REAL*8 FINT,FSEC | |
37 | COMMON /PHOKEY/ FSEC,FINT,INTERF,ISEC,IFTOP | |
38 | INTEGER INIT,I | |
39 | SAVE INIT | |
40 | DATA INIT/ 0/ | |
41 | C-- | |
42 | C-- Return if already initialized... | |
43 | IF (INIT.NE.0) RETURN | |
44 | INIT=1 | |
45 | C-- | |
46 | C-- Preset switch for photon emission to 'TRUE' for each particle in | |
47 | C-- /PH_HEPEVT/, this interface is needed for KORALB and KORALZ... | |
48 | DO 10 I=1,NMXHEP | |
49 | 10 QEDRAD(I)=.TRUE. | |
50 | C-- | |
51 | C-- Logical output unit for printing of PHOTOS error messages | |
52 | PHLUN=6 | |
53 | C-- | |
54 | C-- Set cut parameter for photon radiation | |
55 | XPHCUT=0.01D0 | |
56 | C-- | |
57 | C-- Define some constants | |
58 | ALPHA=0.00729735039D0 | |
59 | PI=3.14159265358979324D0 | |
60 | TWOPI=6.28318530717958648D0 | |
61 | C-- | |
62 | C-- Default seeds Marsaglia and Zaman random number generator | |
63 | ISEED(1)=1802 | |
64 | ISEED(2)=9373 | |
65 | C-- | |
66 | C-- Iitialization for extra options | |
67 | C-- (1) | |
68 | C-- Interference weight for two body symmetric channels only. | |
69 | INTERF=.TRUE. | |
70 | C-- (2) | |
71 | C-- Second order - double photon switch | |
72 | ISEC=.TRUE. | |
73 | C-- (3) | |
74 | C-- Emision in the hard process g g (q qbar) --> t tbar | |
75 | C-- t --> W b | |
76 | IFTOP=.TRUE. | |
77 | C-- | |
78 | C-- further initialization done automatically | |
79 | IF (INTERF) THEN | |
80 | C-- best choice is if FINT=2**N where N+1 is maximal number | |
81 | C-- of charged daughters | |
82 | C-- see report on overweihted events | |
83 | FINT=2.0D0 | |
84 | ELSE | |
85 | FINT=1.0D0 | |
86 | ENDIF | |
87 | C-- Initialise status counter for warning messages | |
88 | DO 20 I=1,PHOMES | |
89 | 20 STATUS(I)=0 | |
90 | RETURN | |
91 | END |