PYTHIA/pythia/pyevwt.F

   1
   2 C*********************************************************************
   3
   4       SUBROUTINE PYEVWT(WTXS)
   5
   6 C...Dummy routine, which the user can replace in order to multiply the
   7 C...standard PYTHIA differential cross-section by a process- and
   8 C...kinematics-dependent factor WTXS. For MSTP(142)=1 this corresponds
   9 C...to generation of weighted events, with weight 1/WTXS, while for
  10 C...MSTP(142)=2 it corresponds to a modification of the underlying
  11 C...physics.
  12       COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
  13       SAVE /LUDAT1/
  14       COMMON/PYINT1/MINT(400),VINT(400)
  15       COMMON/PYINT2/ISET(200),KFPR(200,2),COEF(200,20),ICOL(40,4,2)
  16       SAVE /PYINT1/,/PYINT2/
  17
  18 C...Set default weight for WTXS.
  19       WTXS=1.
  20
  21 C...Read out subprocess number.
  22       ISUB=MINT(1)
  23       ISTSB=ISET(ISUB)
  24
  25 C...Read out tau, y*, cos(theta), tau' (where defined, else =0).
  26       TAU=VINT(21)
  27       YST=VINT(22)
  28       CTH=0.
  29       IF(ISTSB.EQ.2.OR.ISTSB.EQ.4.OR.ISTSB.EQ.6) CTH=VINT(23)
  30       TAUP=0.
  31       IF(ISTSB.GE.3.AND.ISTSB.LE.5) TAUP=VINT(26)
  32
  33 C...Read out x_1, x_2, x_F, shat, that, uhat, p_T^2.
  34       X1=VINT(41)
  35       X2=VINT(42)
  36       XF=X1-X2
  37       SHAT=VINT(44)
  38       THAT=VINT(45)
  39       UHAT=VINT(46)
  40       PT2=VINT(48)
  41
  42 C...Modifications by user to be put here.
  43
  44 C...Stop program if this routine is ever called.
  45 C...You should not copy these lines to your own routine.
  46       WRITE(MSTU(11),5000)
  47       IF(RLU(0).LT.10.) STOP
  48
  49 C...Format for error printout.
  50  5000 FORMAT(1X,'Error: you did not link your PYEVWT routine ',
  51      &'correctly.'/1X,'Dummy routine in PYTHIA file called instead.'/
  52      &1X,'Execution stopped!')
  53
  54       RETURN
  55       END