pymemx.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pymemx.f

 
C*********************************************************************
 
C...PYMEMX
C...Generates maximum ME weight in some initial-state showers.
C...Inparameter MECOR: kind of hard scattering process
C...Outparameter WTFF: maximum weight for fermion -> fermion
C...             WTGF: maximum weight for gluon/photon -> fermion
C...             WTFG: maximum weight for fermion -> gluon/photon
C...             WTGG: maximum weight for gluon -> gluon
 
      SUBROUTINE pymemx(MECOR,WTFF,WTGF,WTFG,WTGG)
 
C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
C...Commonblocks.
      common/pyjets/n,npad,k(4000,5),p(4000,5),v(4000,5)
      common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
      common/pypars/mstp(200),parp(200),msti(200),pari(200)
      common/pyint1/mint(400),vint(400)
      common/pyint2/iset(500),kfpr(500,2),coef(500,20),icol(40,4,2)
      SAVE /pyjets/,/pydat1/,/pypars/,/pyint1/,/pyint2/
 
C...Default maximum weight.
      wtff=1d0
      wtgf=1d0
      wtfg=1d0
      wtgg=1d0
 
C...Select maximum weight by process.
      IF(mecor.EQ.1) THEN
        wtff=1d0
        wtgf=3d0
      ELSEIF(mecor.EQ.2) THEN
        wtfg=1d0
        wtgg=1d0
      ENDIF
 
      RETURN
      END