pyevwt.f

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

 
C*********************************************************************
 
C...PYEVWT
C...Dummy routine, which the user can replace in order to multiply the
C...standard PYTHIA differential cross-section by a process- and
C...kinematics-dependent factor WTXS. For MSTP(142)=1 this corresponds
C...to generation of weighted events, with weight 1/WTXS, while for
C...MSTP(142)=2 it corresponds to a modification of the underlying
C...physics.
 
      SUBROUTINE pyevwt(WTXS)
 
C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
C...Commonblocks.
      common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
      common/pyint1/mint(400),vint(400)
      common/pyint2/iset(500),kfpr(500,2),coef(500,20),icol(40,4,2)
      SAVE /pydat1/,/pyint1/,/pyint2/
 
C...Set default weight for WTXS.
      wtxs=1d0
 
C...Read out subprocess number.
      isub=mint(1)
      istsb=iset(isub)
 
C...Read out tau, y*, cos(theta), tau' (where defined, else =0).
      tau=vint(21)
      yst=vint(22)
      cth=0d0
      IF(istsb.EQ.2.OR.istsb.EQ.4) cth=vint(23)
      taup=0d0
      IF(istsb.GE.3.AND.istsb.LE.5) taup=vint(26)
 
C...Read out x_1, x_2, x_F, shat, that, uhat, p_T^2.
      x1=vint(41)
      x2=vint(42)
      xf=x1-x2
      shat=vint(44)
      that=vint(45)
      uhat=vint(46)
      pt2=vint(48)
 
C...Modifications by user to be put here.
 
C...Stop program if this routine is ever called.
C...You should not copy these lines to your own routine.
      WRITE(mstu(11),5000)
      CALL pystop(4)
 
C...Format for error printout.
 5000 FORMAT(1x,'Error: you did not link your PYEVWT routine ',
     &'correctly.'/1x,'Dummy routine in PYTHIA file called instead.'/
     &1x,'Execution stopped!')
 
      RETURN
      END