PHOS/shaker/shinit.f

   1 *CMZ :          17/07/98  15.49.05  by  Federico Carminati
   2 *-- Author :
   3       SUBROUTINE SHINIT
   4 c       =================
   5
   6 c       Program Initialization
   7 c       eta/pi0, rho/pi0, omega/pi0, phi/pi0 from asymptotic value
   8 c       & mT-scaling
   9 c       [E.g.: V.Hedberg, LUNFD6/(NFFL-7073)/1987;
  10 c       M.Bourquin and J.M.Gaillard, Nucl. Phys. B114 (1976),334]
  11 c       p/pi and K/pi from Tevatron
  12 c       [T.Alexopoulos et al.: Phys. Rev. Lett. 64 (1990), 991]
  13 c       psi and Drell-Yan are rates, independent of dN/dy, for one unit
  14 c       of rapidity, for central events.
  15 c       Cross section data from J. Schukraft, priv.comm.
  16 c       Drell-Yan rate is for M > 1 GeV
  17 c       psi and Drell-Yan rates include BR --> e+ e-
  18
  19 *KEEP,SHRUNP.
  20       COMMON /SHRUNP/ VMAJ,IMIN,NRUN,NEVTOT
  21 *KEEP,SHPHYP.
  22       COMMON /SHPHYP/ JWEI,NDNDY,YLIM,PTLIM,JWEAK,JPI0,JETA,JPIC,JPRO,
  23      +                  JKAC,JKA0,JRHO,JOME,JPHI,JPSI,JDRY
  24 *KEEP,SHRATS.
  25       COMMON /SHRATS/ RETAPI,RPROPI,RKACPI,RRHOPI,ROMEPI,RPHIPI
  26 *KEEP,SHGENE.
  27       COMMON /SHGENE/ IEVT,NPI0,NETA,NPIC,NPRO,NKAC,NKA0,NRHO,NOME,
  28      +                  NPHI,NPSI,NDRY
  29 *KEEP,SHPSDY.
  30       COMMON /SHPSDY/ RNPSI,RNDRY
  31 *KEEP,LUDAT1.
  32       COMMON /LUDAT1/ MSTU(200),PARU(200),MSTJ(200),PARJ(200)
  33       SAVE /LUDAT1/
  34 *KEND.
  35
  36 c       LUND initialization
  37
  38       MSTU(11) = 3                      ! output file
  39       MSTU(4)  = 200000         ! /LUJETS/ size
  40
  41 c       SHAKER initialization
  42
  43       VMAJ = 0.0        ! Major version number
  44       IMIN = 5  ! Minor version number
  45
  46 c       particle ratios initialization
  47
  48       RETAPI = 0.17     ! eta over pi0 ratio
  49       RPROPI = 0.074    ! p over pi ratio
  50       RKACPI = 0.112    ! K over pi ratio
  51       RRHOPI = 0.15     ! rho over pi0 ratio
  52       ROMEPI = 0.14     ! omega over pi0 ratio
  53       RPHIPI = 0.016    ! phi over pi0 ratio
  54
  55       RNPSI  = 0.0018 ! J/psi rate * BR
  56       RNDRY  = 0.0005 ! Drell-Yan rate * BR for M > 1 GeV
  57
  58
  59
  60       CALL SHRATE
  61       IF (JWEAK.EQ.0) CALL SHWDIS
  62       CALL SHRNDV
  63       CALL SHSDEC(1)    ! Select decays?
  64
  65       RETURN
  66       END