pythiaeRHIC/pythia/pymemx.f

0001
0002 C*********************************************************************
0003
0004 C...PYMEMX
0005 C...Generates maximum ME weight in some initial-state showers.
0006 C...Inparameter MECOR: kind of hard scattering process
0007 C...Outparameter WTFF: maximum weight for fermion -> fermion
0008 C...             WTGF: maximum weight for gluon/photon -> fermion
0009 C...             WTFG: maximum weight for fermion -> gluon/photon
0010 C...             WTGG: maximum weight for gluon -> gluon
0011
0012       SUBROUTINE PYMEMX(MECOR,WTFF,WTGF,WTFG,WTGG)
0013
0014 C...Double precision and integer declarations.
0015       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0016       IMPLICIT INTEGER(I-N)
0017       INTEGER PYK,PYCHGE,PYCOMP
0018 C...Commonblocks.
0019       COMMON/PYJETS/N,NPAD,K(4000,5),P(4000,5),V(4000,5)
0020       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0021       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0022       COMMON/PYINT1/MINT(400),VINT(400)
0023       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
0024       SAVE /PYJETS/,/PYDAT1/,/PYPARS/,/PYINT1/,/PYINT2/
0025
0026 C...Default maximum weight.
0027       WTFF=1D0
0028       WTGF=1D0
0029       WTFG=1D0
0030       WTGG=1D0
0031
0032 C...Select maximum weight by process.
0033       IF(MECOR.EQ.1) THEN
0034         WTFF=1D0
0035         WTGF=3D0
0036       ELSEIF(MECOR.EQ.2) THEN
0037         WTFG=1D0
0038         WTGG=1D0
0039       ENDIF
0040
0041       RETURN
0042       END