[u/mrichter/AliRoot.git] / PHOS / shaker / shinit.f

*CMZ :          17/07/98  15.49.05  by  Federico Carminati
*-- Author :
      SUBROUTINE SHINIT
c	=================

c	Program Initialization
c	eta/pi0, rho/pi0, omega/pi0, phi/pi0 from asymptotic value
c       & mT-scaling
c	[E.g.: V.Hedberg, LUNFD6/(NFFL-7073)/1987;
c       M.Bourquin and J.M.Gaillard, Nucl. Phys. B114 (1976),334]
c	p/pi and K/pi from Tevatron
c	[T.Alexopoulos et al.: Phys. Rev. Lett. 64 (1990), 991]
c	psi and Drell-Yan are rates, independent of dN/dy, for one unit
c       of rapidity, for central events.
c       Cross section data from J. Schukraft, priv.comm.
c	Drell-Yan rate is for M > 1 GeV
c	psi and Drell-Yan rates include BR --> e+ e-

*KEEP,SHRUNP.
      COMMON /SHRUNP/ VMAJ,IMIN,NRUN,NEVTOT
*KEEP,SHPHYP.
      COMMON /SHPHYP/ JWEI,NDNDY,YLIM,PTLIM,JWEAK,JPI0,JETA,JPIC,JPRO,
     +                  JKAC,JKA0,JRHO,JOME,JPHI,JPSI,JDRY
*KEEP,SHRATS.
      COMMON /SHRATS/ RETAPI,RPROPI,RKACPI,RRHOPI,ROMEPI,RPHIPI
*KEEP,SHGENE.
      COMMON /SHGENE/ IEVT,NPI0,NETA,NPIC,NPRO,NKAC,NKA0,NRHO,NOME,
     +                  NPHI,NPSI,NDRY
*KEEP,SHPSDY.
      COMMON /SHPSDY/ RNPSI,RNDRY
*KEEP,LUDAT1.
      COMMON /LUDAT1/ MSTU(200),PARU(200),MSTJ(200),PARJ(200)
      SAVE /LUDAT1/
*KEND.

c	LUND initialization

      MSTU(11) = 3			! output file
      MSTU(4)  = 200000		! /LUJETS/ size

c	SHAKER initialization

      VMAJ = 0.0	! Major version number
      IMIN = 5	! Minor version number

c	particle ratios initialization

      RETAPI = 0.17	! eta over pi0 ratio
      RPROPI = 0.074	! p over pi ratio
      RKACPI = 0.112	! K over pi ratio
      RRHOPI = 0.15	! rho over pi0 ratio
      ROMEPI = 0.14	! omega over pi0 ratio
      RPHIPI = 0.016	! phi over pi0 ratio

      RNPSI  = 0.0018 ! J/psi rate * BR
      RNDRY  = 0.0005 ! Drell-Yan rate * BR for M > 1 GeV


      CALL SHRATE	
      IF (JWEAK.EQ.0) CALL SHWDIS
      CALL SHRNDV
      CALL SHSDEC(1)	! Select decays?

      RETURN
      END
Commit	Line	Data
fe4da5cc	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