16 CALL simple_ee_to_z_to_tautau_event
19 WRITE(*,*)
"##############" 20 WRITE(*,*)
"PHODMP: BEFORE" 21 WRITE(*,*)
"##############" 28 WRITE(*,*)
"##############" 29 WRITE(*,*)
"PHODMP: AFTER" 30 WRITE(*,*)
"##############" 34 SUBROUTINE simple_ee_to_z_to_tautau_event
42 CALL add_particle( 11, 6, 1.7763568394002505d-15, -3.55658944257
44 CALL add_particle( -11, 6, -1.7763568394002505d-15, 3.54883522047
46 CALL add_particle( 23, 5, 0d0, 0d0,
48 CALL add_particle( 15, 2, -2.3191595992562256d+01, -2.63105009206
50 CALL add_particle( -15, 2, 2.3191595992562256d+01, 2.63105009206
52 CALL add_particle( 16, 1, -1.2566536214715378d+00, -1.79702511383
54 CALL add_particle(-211, 1, -2.1935073012334062d+01, -2.45136240172
56 CALL add_particle( -16, 1, 8.4364531743909055d+00, 8.32028308316
58 CALL add_particle( 211, 1, 1.4755273459419701d+01, 1.79903660479
62 SUBROUTINE even_simpler_z_to_ee_event
70 CALL add_particle( 23, 2, 0.000000000000000d+00, 0.0000000000000
72 CALL add_particle( 11, 1, -3.018114335125028d+01, -1.2040719054541
74 CALL add_particle(-11, 1, 3.018114335125028d+01, 1.2040719054541
78 SUBROUTINE add_particle(ID,STATUS,PX,PY,PZ,E,M,MOTHER1,MOTHER2,DAUGHTER1,DAUGHTER2)
79 INTEGER ID,STATUS,MOTHER1,MOTHER2,DAUGHTER1,DAUGHTER2
91 parameter(nmxhep=10000)
93 INTEGER nevhep,nhep,isthep,idhep,jmohep,
122 jmohep(1,nhep)=mother1
123 jmohep(2,nhep)=mother2
124 jdahep(1,nhep)=daughter1
125 jdahep(2,nhep)=daughter2
144 DOUBLE PRECISION SUMVEC(5)
148 parameter(nmxhep=10000)
150 INTEGER nevhep,nhep,isthep,idhep,jmohep,
169 WRITE(phlun,9010) nevhep
175 IF (jdahep(1,i).EQ.0)
THEN 177 20 sumvec(j)=sumvec(j)+phep(j,i)
178 IF (jmohep(2,i).EQ.0)
THEN 179 WRITE(phlun,9030) i,idhep(i),jmohep(1,i),(phep(j,i),j=1,5)
181 WRITE(phlun,9040) i,idhep(i),jmohep(1,i),jmohep(2,i),(phep
185 IF (jmohep(2,i).EQ.0)
THEN 186 WRITE(phlun,9050) i,idhep(i),jmohep(1,i),jdahep(1,i),
187 & jdahep(2,i),(phep(j,i),j=1,5)
189 WRITE(phlun,9060) i,idhep(i),jmohep(1,i),jmohep(2,i),
190 & jdahep(1,i),jdahep(2,i),(phep(j,i),j=1,5)
194 sumvec(5)=sqrt(sumvec(4)**2-sumvec(1)**2-sumvec(2)**2-
196 WRITE(phlun,9070) (sumvec(j),j=1,5)
198 9000
FORMAT(1h0,80(
'='))
199 9010
FORMAT(1h ,29x,
'Event No.:',i10)
200 9020
FORMAT(1h0,1x,
'Nr',3x,
'Type',3x,
'Parent(s)',2x,
'Daughter(s)',6x,
201 &
'Px',7x,
'Py',7x,
'Pz',7x,
'E',4x,
'Inv. M.')
202 9030
FORMAT(1h ,i4,i7,3x,i4,9x,
'Stable',2x,5f9.2)
203 9040
FORMAT(1h ,i4,i7,i4,
' - ',i4,5x,
'Stable',2x,5f9.2)
204 9050
FORMAT(1h ,i4,i7,3x,i4,6x,i4,
' - ',i4,5f9.2)
205 9060
FORMAT(1h ,i4,i7,i4,
' - ',i4,2x,i4,
' - ',i4,5f9.2)
206 9070
FORMAT(1h0,23x,
'Vector Sum: ', 5f9.2)
207 9080
FORMAT(1h0,6x,
'Particle Parameters')