| 1 | *
|
|---|
| 2 | COMMON /LUDAT2/ KCHG(500,3),PMAS(500,4),PARF(2000),VCKM(4,4)
|
|---|
| 3 | COMMON /LUDAT3/ MDCY(500,3),MDME(2000,2),BRAT(2000),KFDP(2000,5)
|
|---|
| 4 | COMMON /PYSUBS/ MSEL,MSUB(200),KFIN(2,-40:40),CKIN(200)
|
|---|
| 5 | COMMON /PYPARS/ MSTP(200),PARP(200),MSTI(200),PARI(200)
|
|---|
| 6 | *********** COMMON /LUJETS/ N,K(4000,5),P(4000,5),V(4000,5)
|
|---|
| 7 | EXTERNAL LUDATA,PYDATA
|
|---|
| 8 | CHARACTER FRAME*4,BEAM*10,TARGET*10
|
|---|
| 9 | REAL WIN
|
|---|
| 10 | INTEGER LUCOMP
|
|---|
| 11 | REAL PDT(500,5)
|
|---|
| 12 | *
|
|---|
| 13 | FRAME = 'CMS'
|
|---|
| 14 | BEAM = 'P'
|
|---|
| 15 | TARGET = 'P'
|
|---|
| 16 | WIN = 14000.
|
|---|
| 17 | NEVNT = 3
|
|---|
| 18 | *
|
|---|
| 19 | MSEL = 0
|
|---|
| 20 | MSUB(102) = 1
|
|---|
| 21 | MSUB(123) = 1
|
|---|
| 22 | MSUB(124) = 1
|
|---|
| 23 | *
|
|---|
| 24 | PMAS(6,1) = 176.
|
|---|
| 25 | PMAS(25,1) = 500.
|
|---|
| 26 | CKIN(1) = 470.
|
|---|
| 27 | CKIN(2) = 530.
|
|---|
| 28 | *
|
|---|
| 29 | * turn off all Higgs decays except the ZZ
|
|---|
| 30 | *
|
|---|
| 31 | IH = LUCOMP(25)
|
|---|
| 32 | DO IDC=MDCY(IH,2),MDCY(IH,2)+MDCY(IH,3)-1
|
|---|
| 33 | IF (KFDP(IDC,1).NE.23.AND.MDME(IDC,1).EQ.1) MDME(IDC,1)=0
|
|---|
| 34 | ENDDO
|
|---|
| 35 | *
|
|---|
| 36 | * turn off all Z decays except mumu or ee
|
|---|
| 37 | *
|
|---|
| 38 | IZ = LUCOMP(23)
|
|---|
| 39 | DO IDC=MDCY(IZ,2),MDCY(IZ,2)+MDCY(IZ,3)-1
|
|---|
| 40 | IF (MDME(IDC,1).EQ.1) THEN
|
|---|
| 41 | IF ((IABS(KFDP(IDC,1)).NE.13)
|
|---|
| 42 | > .AND.(IABS(KFDP(IDC,1)).NE.11)) THEN
|
|---|
| 43 | MDME(IDC,1)=0
|
|---|
| 44 | ENDIF
|
|---|
| 45 | ENDIF
|
|---|
| 46 | ENDDO
|
|---|
| 47 | ***********************************************************
|
|---|
| 48 | CALL PYINIT(FRAME,BEAM,TARGET,WIN)
|
|---|
| 49 | *
|
|---|
| 50 | DO IEVT = 1, NEVNT
|
|---|
| 51 | *
|
|---|
| 52 | CALL PYEVNT
|
|---|
| 53 | CALL LUEDIT(11)
|
|---|
| 54 | CALL LUEDIT(12)
|
|---|
| 55 | CALL LUEDIT(15)
|
|---|
| 56 | *
|
|---|
| 57 | CALL LUHEPC(1)
|
|---|
| 58 | *
|
|---|
| 59 | CALL HEP2G4
|
|---|
| 60 | *
|
|---|
| 61 | ENDDO
|
|---|
| 62 | ***********************************************************
|
|---|
| 63 | STOP
|
|---|
| 64 | END
|
|---|
| 65 | ***********************************************************
|
|---|
| 66 | SUBROUTINE HEP2G4
|
|---|
| 67 | *
|
|---|
| 68 | * Output /HEPEVT/ event structure to G4HEPEvtInterface
|
|---|
| 69 | *
|
|---|
| 70 | * M.Asai (asai@kekvax.kek.jp) -- 24/09/96
|
|---|
| 71 | *
|
|---|
| 72 | ***********************************************************
|
|---|
| 73 | PARAMETER (NMXHEP=2000)
|
|---|
| 74 | COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
|
|---|
| 75 | >JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
|
|---|
| 76 | DOUBLE PRECISION PHEP,VHEP
|
|---|
| 77 | *
|
|---|
| 78 | WRITE(6,*) NHEP
|
|---|
| 79 | DO IHEP=1,NHEP
|
|---|
| 80 | WRITE(6,10)
|
|---|
| 81 | > ISTHEP(IHEP),IDHEP(IHEP),JDAHEP(1,IHEP),JDAHEP(2,IHEP),
|
|---|
| 82 | > PHEP(1,IHEP),PHEP(2,IHEP),PHEP(3,IHEP),PHEP(5,IHEP)
|
|---|
| 83 | 10 FORMAT(4I5,4(1X,D15.8))
|
|---|
| 84 | ENDDO
|
|---|
| 85 | *
|
|---|
| 86 | RETURN
|
|---|
| 87 | END
|
|---|
| 88 |
|
|---|
