C...Simple program for ppbar -> W' -> W + Z, illustrating
C...the basics of process selection.

      COMMON/LUJETS/N,K(4000,5),P(4000,5),V(4000,5)
      COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
      COMMON/LUDAT2/KCHG(500,3),PMAS(500,4),PARF(2000),VCKM(4,4)
      COMMON/LUDAT3/MDCY(500,3),MDME(2000,2),BRAT(2000),KFDP(2000,5)
      COMMON/PYSUBS/MSEL,MSUB(200),KFIN(2,-40:40),CKIN(200) 
      COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)

C...Set process type and W' mass.
      MSEL=0
      MSUB(142)=1
      PMAS(34,1)=400.

C...Switch on W' -> WZ and off all the others.
      DO 100 ID=MDCY(34,2),MDCY(34,2)+MDCY(34,3)-1
      MDME(ID,1)=MIN(0,MDME(ID,1))
      IF(KFDP(ID,1).EQ.24.AND.KFDP(ID,2).EQ.23) MDME(ID,1)=1
  100 CONTINUE

C...Initialize and list resonance widths.
      CALL PYINIT('CMS','P','PBAR',1800.)
      CALL PYSTAT(2)

C...Generate events and list first two.
      DO 110 IEV=1,10
      CALL PYEVNT
      IF(IEV.LE.2) CALL LULIST(1)
  110 CONTINUE

C...Print cross section from run.
      CALL PYSTAT(1)

      END