[807] | 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 | |
---|