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 | |
---|