runAndEvent/RE05/pythia_main.f

0001 *
0002       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0003       IMPLICIT INTEGER(I-N)
0004       INTEGER PYK,PYCHGE,PYCOMP
0005 *
0006       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0007       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
0008       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
0009       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0010 *
0011       EXTERNAL PYDATA
0012       CHARACTER FRAME*4,BEAM*10,TARGET*10
0013       REAL WIN
0014       REAL PDT(500,5)
0015 *
0016       NEVNT  = 100
0017 *
0018       MSEL = 0
0019       MSUB(102) = 1
0020       MSUB(123) = 1
0021       MSUB(124) = 1
0022 *
0023       PMAS(6,1) = 172.9
0024       PMAS(25,1) = 125.
0025 *
0026 *     turn off all Higgs decays except the ZZ
0027 *
0028       IH = PYCOMP(25)
0029       DO IDC=MDCY(IH,2),MDCY(IH,2)+MDCY(IH,3)-1
0030         IF (KFDP(IDC,1).NE.23.AND.MDME(IDC,1).EQ.1) MDME(IDC,1)=0
0031       ENDDO
0032 *
0033 *     turn off all Z decays except mumu or ee
0034 *
0035       IZ = PYCOMP(23)
0036       DO IDC=MDCY(IZ,2),MDCY(IZ,2)+MDCY(IZ,3)-1
0037         IF (MDME(IDC,1).EQ.1) THEN
0038           IF ((IABS(KFDP(IDC,1)).NE.13)
0039      >   .AND.(IABS(KFDP(IDC,1)).NE.11)) THEN
0040             MDME(IDC,1)=0
0041           ENDIF
0042         ENDIF
0043       ENDDO
0044 ***********************************************************
0045       CALL PYINIT('CMS','p','p',14000D0)
0046 *
0047       DO IEVT = 1, NEVNT
0048 *
0049        CALL PYEVNT
0050        CALL PYEDIT(11)
0051        CALL PYEDIT(12)
0052        CALL PYEDIT(15)
0053 *
0054        CALL PYHEPC(1)
0055 *
0056        CALL HEP2G4
0057 *
0058       ENDDO
0059 ***********************************************************
0060       STOP
0061       END
0062 ***********************************************************
0063       SUBROUTINE HEP2G4
0064 *
0065 * Output /HEPEVT/ event structure to G4HEPEvtInterface
0066 *
0067 * M.Asai (asai@kekvax.kek.jp)  --  24/09/96
0068 *
0069 ***********************************************************
0070       PARAMETER (NMXHEP=4000)
0071       COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
0072      >JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
0073       DOUBLE PRECISION PHEP,VHEP
0074 *
0075       WRITE(6,*) NHEP
0076       DO IHEP=1,NHEP
0077        WRITE(6,10)
0078      >  ISTHEP(IHEP),IDHEP(IHEP),JDAHEP(1,IHEP),JDAHEP(2,IHEP),
0079      >  PHEP(1,IHEP),PHEP(2,IHEP),PHEP(3,IHEP),PHEP(5,IHEP)
0080 10    FORMAT(I4,I10,I5,I5,4(1X,D15.8))
0081       ENDDO
0082 *
0083       RETURN
0084       END