4 real(8):: eelec, prom, eprot, egam, sigmax, fac
9 INTEGER NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
10 DOUBLE PRECISION PHEP,VHEP
11 COMMON /poevt1/ nevhep,nhep,isthep(nmxhep),idhep(nmxhep),
12 & jmohep(2,nmxhep),jdahep(2,nmxhep),phep(5,nmxhep),
15 INTEGER IMPART,IPHIST,ICOLOR
16 COMMON /poevt2/ impart(nmxhep),iphist(2,nmxhep),icolor(2,nmxhep)
21 call pho_init(-1, irej)
28 p1(3) = sqrt(eprot**2 - prom**2)
39 call pho_setpar(1, 2212, 0, 0.
d0)
40 call pho_setpar(2, 22, 0, 0.
d0)
42 call pho_event(-1, p1, p2,sigmax, irej)
44 write(0,*)
' sigmax=',sigmax,
' ireg=',irej
49 call pho_event(3, p1, p2, fac, irej)
50 write(0,*)
' fac=',fac,
' ireg=',irej
52 write(0,*)
' nhep=',nhep
54 if( isthep(
j) == 1 )
then 55 write(0,
'(2i5, 1p, 5g13.3)')
56 * isthep(
j), idhep(
j), phep(:,
j)
61 call pho_event(-2, p1, p2, sigmax,irej)
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0