TEvtGen/PHOTOS/phocin.F

   1       SUBROUTINE PHOCIN(ISEEDF,ISEEDS)
   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 ISEEDF, ISEEDS
  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
  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
  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