PHOS/shaker/shrate.f

   1 *CMZ :          17/07/98  15.49.05  by  Federico Carminati
   2 *-- Author :
   3       SUBROUTINE SHRATE
   4 c       =================
   5
   6 c       Compute rates for various particles.
   7 c       !!! eta --> gamma gamma only !!!
   8
   9 *KEEP,SHPHYP.
  10       COMMON /SHPHYP/ JWEI,NDNDY,YLIM,PTLIM,JWEAK,JPI0,JETA,JPIC,JPRO,
  11      +                  JKAC,JKA0,JRHO,JOME,JPHI,JPSI,JDRY
  12 *KEEP,SHRATS.
  13       COMMON /SHRATS/ RETAPI,RPROPI,RKACPI,RRHOPI,ROMEPI,RPHIPI
  14 *KEEP,SHPRAT.
  15       COMMON /SHPRAT/ PI0R,ETAR,RHOR,OMER,PHIR,PSIR,DRYR
  16 *KEEP,SHGENE.
  17       COMMON /SHGENE/ IEVT,NPI0,NETA,NPIC,NPRO,NKAC,NKA0,NRHO,NOME,
  18      +                  NPHI,NPSI,NDRY
  19 *KEEP,SHPSDY.
  20       COMMON /SHPSDY/ RNPSI,RNDRY
  21 *KEND.
  22
  23       CHA = NDNDY*2*YLIM
  24       NCHA = CHA
  25       NPIC = JPIC*NCHA/(1+JPRO*RPROPI+JKAC*RKACPI+JKA0*RKACPI)
  26       NPRO = JPRO*RPROPI*NPIC
  27       NKAC = JKAC*RKACPI*NPIC
  28       NKA0 = JKA0*RKACPI*NPIC
  29
  30       IF (JPIC.EQ.0) NPI0=JPI0*NCHA/2
  31       IF (JPIC.EQ.1) NPI0=JPI0*NPIC/2
  32       NETA=JETA*RETAPI*0.389*NPI0  !eta/pi0 * npi0 * BR[eta --> gamma gamma]
  33
  34 c       Initialize /SHPRAT/ rates...
  35
  36 c       ... to half of charged multiplicity for pi0
  37
  38       PI0R = FLOAT(NCHA)/2.
  39
  40 c       ... to pi0 * eta/pi0 * BR --> gamma gamma for eta
  41
  42       ETAR = PI0R*RETAPI*.389
  43
  44 c       ... to pi0 * vect/pi0 * BR --> e+e- for vector mesons
  45
  46       RHOR = PI0R*RRHOPI*4.44E-5
  47       OMER = PI0R*ROMEPI*7.07E-5
  48       PHIR = PI0R*RPHIPI*3.11E-4
  49
  50 c       ... to calculated rate for J/psi and Drell-Yan
  51
  52       PSIR = RNPSI
  53       DRYR = RNDRY
  54
  55
  56       RETURN
  57       END