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