HERWIG/src/hwueec.f

   1
   2 CDECK  ID>, HWUEEC.
   3
   4 *CMZ :-        -26/04/91  14.22.30  by  Federico Carminati
   5
   6 *-- Author :    Bryan Webber and Ian Knowles
   7
   8 C-----------------------------------------------------------------------
   9
  10       SUBROUTINE HWUEEC(IL)
  11
  12 C-----------------------------------------------------------------------
  13
  14 C     Loads cross-section coefficients, for kinematically open channels,
  15
  16 C     in llbar-->qqbar; lepton label IL=1-6: e,nu_e,mu,nu_mu,tau,nu_tau.
  17
  18 C-----------------------------------------------------------------------
  19
  20       INCLUDE 'HERWIG61.INC'
  21
  22       DOUBLE PRECISION Q2
  23
  24       INTEGER IL,JL,IQ
  25
  26       Q2=EMSCA**2
  27
  28       JL=IL+10
  29
  30       MAXFL=0
  31
  32       TQWT=0.
  33
  34       DO 10 IQ=1,NFLAV
  35
  36       IF (EMSCA.GT.2.*RMASS(IQ)) THEN
  37
  38          MAXFL=MAXFL+1
  39
  40          MAPQ(MAXFL)=IQ
  41
  42          CALL HWUCFF(JL,IQ,Q2,CLQ(1,MAXFL))
  43
  44          TQWT=TQWT+CLQ(1,MAXFL)
  45
  46       ENDIF
  47
  48   10  CONTINUE
  49
  50       IF (MAXFL.EQ.0) CALL HWWARN('HWUEEC',100,*999)
  51
  52   999 END
  53
  54 CDECK  ID>, HWUEMV.
  55
  56 *CMZ :-        -30/06/94  19.31.08  by  Mike Seymour
  57
  58 *-- Author :    Mike Seymour
  59
  60 C-----------------------------------------------------------------------
  61
  62       SUBROUTINE HWUEMV(N,IFROM,ITO)
  63
  64 C-----------------------------------------------------------------------
  65
  66 C     MOVE A BLOCK OF ENTRIES IN THE EVENT RECORD
  67
  68 C     N ENTRIES IN HEPEVT STARTING AT IFROM ARE MOVED TO AFTER ITO
  69
  70 C-----------------------------------------------------------------------
  71
  72       INCLUDE 'HERWIG61.INC'
  73
  74       INTEGER N,IFROM,ITO,IMAP(0:NMXHEP),LFROM,LTO,I,IEDT(NMXHEP),IHEP,
  75
  76      $ JHEP,KHEP
  77
  78       COMMON /HWUMAP/IMAP
  79
  80       LFROM=IFROM
  81
  82       LTO=ITO
  83
  84       DO 100 I=1,N
  85
  86  100  IEDT(I)=LTO
  87
  88       CALL HWUEDT(-N,IEDT)
  89
  90       DO 300 I=1,N
  91
  92         IHEP=LTO+I
  93
  94         JHEP=IMAP(LFROM+I-1)
  95
  96         ISTHEP(IHEP)=ISTHEP(JHEP)
  97
  98         IDHW(IHEP)=IDHW(JHEP)
  99
 100         IDHEP(IHEP)=IDHEP(JHEP)
 101
 102         JMOHEP(1,IHEP)=JMOHEP(1,JHEP)
 103
 104         JMOHEP(2,IHEP)=JMOHEP(2,JHEP)
 105
 106         JDAHEP(1,IHEP)=JDAHEP(1,JHEP)
 107
 108         JDAHEP(2,IHEP)=JDAHEP(2,JHEP)
 109
 110         CALL HWVEQU(5,PHEP(1,JHEP),PHEP(1,IHEP))
 111
 112         CALL HWVEQU(4,VHEP(1,JHEP),VHEP(1,IHEP))
 113
 114         DO 200 KHEP=1,NHEP
 115
 116           IF (JMOHEP(1,KHEP).EQ.JHEP) JMOHEP(1,KHEP)=IHEP
 117
 118           IF (JMOHEP(2,KHEP).EQ.JHEP) JMOHEP(2,KHEP)=IHEP
 119
 120           IF (JDAHEP(1,KHEP).EQ.JHEP) JDAHEP(1,KHEP)=IHEP
 121
 122           IF (JDAHEP(2,KHEP).EQ.JHEP) JDAHEP(2,KHEP)=IHEP
 123
 124  200    CONTINUE
 125
 126         IEDT(I)=JHEP
 127
 128  300  CONTINUE
 129
 130       CALL HWUEDT(N,IEDT)
 131
 132  999  END
 133
 134 CDECK  ID>, HWUEPR.
 135
 136 *CMZ :-        -27/07/99  13.33.03  by  Mike Seymour
 137
 138 *-- Author :    Ian Knowles, Bryan Webber & Kosuke Odagiri
 139
 140 C-----------------------------------------------------------------------
 141
 142       SUBROUTINE HWUEPR
 143
 144 C-----------------------------------------------------------------------
 145
 146 C     Prints out event data in a number of possible formats:
 147
 148 C     If (PRNDEF) ASCII to stout
 149
 150 C     If (PRNTEX) LaTeX to the file HWEV_*******.tex
 151
 152 C                 Please check paper size and offsets given in mm
 153
 154 C                 Uses the package longtable.sty
 155
 156 C                 If (PRVTX>OR.NPRFMT.EQ.2) designed to be printed
 157
 158 C                 as landscape
 159
 160 C     If (PRNWEB) HTML  to the file HWEV_*******.html
 161
 162 C                 Call HWUDPR to create particle property files in
 163
 164 C                 the subdirectory HW_decays/
 165
 166 C     ******* gives the event number 0000001 etc.
 167
 168 C-----------------------------------------------------------------------
 169
 170       INCLUDE 'HERWIG61.INC'
 171
 172       INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,I,IST,IS,ID,MS,J,K,IUNITW,
 173
 174      & IUNITT
 175
 176       CHARACTER*1 Z
 177
 178       CHARACTER*2 ZZ
 179
 180       CHARACTER*6 BGCOLS(5),TBCOLS(3),THEAD(17,3)
 181
 182       CHARACTER*7 HWUNST,TMPNME
 183
 184       CHARACTER*16 FNAMET
 185
 186       CHARACTER*17 FNAMEW
 187
 188       CHARACTER*27 FNAMEP
 189
 190       CHARACTER*28 TITLE(11),SECTXT
 191
 192       LOGICAL FIRST(11),NEWSEC
 193
 194       COMMON/PAPER/MMWIDE,MMLONG,MMHOFF,MMVOFF
 195
 196       EXTERNAL HWUNST
 197
 198 C
 199
 200       DATA BGCOLS/'ffffff','0000aa','aa0000','00aa00','aa00ff'/
 201
 202       DATA TBCOLS/'ccccff','9966ff','ffff00'/
 203
 204       DATA THEAD/ 17*'9966ff',17*'ffff00',
 205
 206      &            'IHEP  ','  ID  ',' IDPDG',' IST  ',' MO1  ',' MO2  ',
 207
 208      &            ' DA1  ',' DA2  ',' P-X  ',' P-Y  ',' P-Z  ','ENERGY',
 209
 210      &            ' MASS ',' V-X  ',' V-Y  ',' V-Z  ',' V-C*T'/
 211
 212       DATA TITLE/'     ---INITIAL STATE---    ',
 213
 214      &           '    ---HARD SUBPROCESS---   ',
 215
 216      &           '    ---PARTON SHOWERS---    ',
 217
 218      &           '    ---GLUON SPLITTING---   ',
 219
 220      &           '   ---CLUSTER FORMATION---  ',
 221
 222      &           '    ---CLUSTER DECAYS---    ',
 223
 224      &           ' ---STRONG HADRON DECAYS--- ',
 225
 226      &           ' ---HEAVY PARTICLE DECAYS---',
 227
 228      &           '  ---H/W/Z BOSON DECAYS---  ',
 229
 230      &           ' ---SOFT UNDERLYING EVENT---',
 231
 232      &           '  ---MULTIPLE SCATTERING--- '/
 233
 234       Z=CHAR(92)
 235
 236       ZZ=Z//Z
 237
 238 C
 239
 240       IUNITT=50
 241
 242       IUNITW=51
 243
 244 C Write out any required file header information
 245
 246       TMPNME=HWUNST(NEVHEP)
 247
 248       IF (PRNTEX) THEN
 249
 250         WRITE(FNAMET,'(A5,A7,A4)') 'HWEV_',TMPNME,'.tex'
 251
 252         OPEN(IUNITT,STATUS='UNKNOWN',FILE=FNAMET)
 253
 254         IF (PRVTX.OR.NPRFMT.EQ.2) THEN
 255
 256           WRITE(IUNITT,10) Z,Z,Z,MMLONG,Z,MMWIDE,Z,MMVOFF,Z,MMHOFF,Z,Z,Z
 257
 258         ELSE
 259
 260           WRITE(IUNITT,10) Z,Z,Z,MMWIDE,Z,MMLONG,Z,MMHOFF,Z,MMVOFF,Z,Z,Z
 261
 262         ENDIF
 263
 264       ENDIF
 265
 266       IF (PRNWEB) THEN
 267
 268         WRITE(FNAMEW,'(A5,A7,A5)') 'HWEV_',TMPNME,'.html'
 269
 270         OPEN(IUNITW,STATUS='UNKNOWN',FILE=FNAMEW)
 271
 272         WRITE(IUNITW,20) BGCOLS
 273
 274       ENDIF
 275
 276    10 FORMAT(A1,'documentclass{article}'/A1,'usepackage{longtable}'/
 277
 278      & A1,'textwidth ',I4,'mm ',A1,'textheight ',I4,'mm'/
 279
 280      & A1,'hoffset   ',I4,'mm ',A1,'voffset    ',I4,'mm'/
 281
 282      & A1,'pagestyle{empty}'/A1,'begin{document}'/A1,'begin{center}')
 283
 284    20 FORMAT('<!-- No Less productions -->'/'<HTML>'/'<HEAD>'/
 285
 286      & '<TITLE>HERWIG Event Record</TITLE>'/'</HEAD>'/
 287
 288      & '<BODY BGCOLOR=#',A6,' TEXT=#',A6,' LINK=#',A6,
 289
 290      & ' ALINK=#',A6,' VLINK=#',A6,'>')
 291
 292 C Write out event header details and set up tables
 293
 294       IF (PRNDEF) THEN
 295
 296         WRITE(6,30) NEVHEP,PBEAM1,PART1,PBEAM2,PART2,IPROC,NRN,
 297
 298      &   ISTAT,IERROR,EVWGT
 299
 300       ENDIF
 301
 302       IF (PRNTEX) THEN
 303
 304         WRITE(IUNITT,40) Z,Z,Z,ISTAT,ZZ,Z,IPROC,PBEAM1,PBEAM2,NRN(1),
 305
 306      &   IERROR,ZZ,Z,Z,NEVHEP,TXNAME(1,IDHW(1)),TXNAME(1,IDHW(2)),
 307
 308      &   NRN(2),EVWGT,ZZ,Z,Z,Z
 309
 310         IF (PRVTX) THEN
 311
 312           WRITE(IUNITT,50) Z,Z,Z,Z,Z
 313
 314         ELSE
 315
 316           WRITE(IUNITT,60) Z,Z,Z,Z,Z
 317
 318         ENDIF
 319
 320       ENDIF
 321
 322       IF (PRNWEB) THEN
 323
 324         WRITE(IUNITW,70) TBCOLS(1),TBCOLS(2),(TBCOLS(2),TBCOLS(3),
 325
 326      &   I=1,4),ISTAT,TBCOLS(2),TBCOLS(3),IPROC,PBEAM1,PBEAM2,NRN(1),
 327
 328      &   TBCOLS(2),TBCOLS(3),IERROR
 329
 330         WRITE(IUNITW,71) TBCOLS(2),TBCOLS(3),NEVHEP,TXNAME(2,IDHW(1)),
 331
 332      &   TXNAME(2,IDHW(2)),NRN(2),TBCOLS(2),TBCOLS(3),EVWGT,TBCOLS(1)
 333
 334       ENDIF
 335
 336    30 FORMAT(///1X,'EVENT ',I7,':',F8.2,' GEV/C ',A8,' ON ',F8.2,
 337
 338      & ' GEV/C ',A8,' PROCESS:',I6/1X,'SEEDS: ',I11,' & ',I11,
 339
 340      & '   STATUS: ',I4,' ERROR:',I4,'  WEIGHT: ',1P,E11.4/)
 341
 342    40 FORMAT(A1,'begin{tabular}{|l|r|c|c|r|l|c|}'/A1,'hline'/
 343
 344      & A1,'multicolumn{2}{|c|}{HERWIG 6.0} & Beam 1: & Beam 2: & ',
 345
 346      & 'Seeds: & Status: & ',I4, ' ',A2/A1,'hline'/'Process: & ',I6,
 347
 348      & ' & ',F8.2,'~GeV/c & ',F8.2,'~GeV/c',' & ',I11,' & Error: & ',
 349
 350      & I4,' ',A2/A1,'cline{1-2} ',A1,'cline{6-7}'/'Event: & ',I7,' & ',
 351
 352      & A37,' & ',A37,' & ',I11,' & Weight: & ',1P,E11.4,' ',A2/A1,
 353
 354      & 'hline'/A1,'end{tabular}'/A1,'vskip 5mm')
 355
 356    50 FORMAT(A1,'begin{longtable}{|r|c|r|r|r|r|r|r|r|r|r|r|r|r|r|r|r|}'/
 357
 358      & A1,'hline'/A1,'endhead'/A1,'hline'/A1,'endfoot')
 359
 360    60 FORMAT(A1,'begin{longtable}{|r|c|r|r|r|r|r|r|r|r|r|r|r|}'/
 361
 362      & A1,'hline'/A1,'endhead'/A1,'hline'/A1,'endfoot')
 363
 364    70 FORMAT(/'<CENTER>'/'<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>'/
 365
 366      & '<TR>'/'<TH BGCOLOR=#',A6,' COLSPAN=2>',
 367
 368      & '<A HREF="http://hepwww.rl.ac.uk/theory/seymour/herwig/">',
 369
 370      & 'HERWIG 6.0</A></TH>'/
 371
 372      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Beam 1:</FONT></TH>'/
 373
 374      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Beam 2:</FONT></TH>'/
 375
 376      & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Seeds:</FONT></TH>'/
 377
 378      & '<TH BGCOLOR=#',A6,' ALIGN="LEFTT"><FONT COLOR=#',A6,
 379
 380      & '>Status:</FONT></TH>'/'<TD ALIGN="RIGHT">',I4,'</TD>'/'</TR>'/
 381
 382      & '<TR>'/
 383
 384      & '<TH BGCOLOR=#',A6,' ALIGN="LEFTT"><FONT COLOR=#',A6,
 385
 386      & '>Process:</Th>'/'<TD>',I6,'</TD>'/
 387
 388      & '<TD>',F8.2,' GeV/c</TD>'/'<TD>',F8.2,' GeV/c</TD>'/
 389
 390      & '<TD ALIGN="RIGHT">',I11,'</TD>'/
 391
 392      & '<TH BGCOLOR=#',A6,' ALIGN="LEFT"><FONT COLOR=#',A6,
 393
 394      & '>Error:</FONT></TH>'/'<TD ALIGN="RIGHT">',I4,'</TD>'/'</TR>')
 395
 396    71 FORMAT('<TR>'/
 397
 398      & '<TH BGCOLOR=#',A6,' ALIGN="LEFT"><FONT COLOR=#',A6,
 399
 400      & '>Event:</Th>'/'<TD ALIGN="RIGHT">',I7,'</TD>'/
 401
 402      & '<TD ALIGN="CENTER">',A37,'</TD>'/
 403
 404      & '<TD ALIGN="CENTER">',A37,'</TD>'/
 405
 406      & '<TD ALIGN="RIGHT">',I11,'</TD>'/
 407
 408      & '<TH BGCOLOR=#',A6,' ALIGN="LEFT"><FONT COLOR=#',A6,
 409
 410      & '>Weight:</FONT></TH>'/'<TD>',1P,E11.4,'</TD>'/'</TR>'/
 411
 412      & '</TABLE>'//'<P>'/
 413
 414      & '<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>')
 415
 416 C Initialize control flags
 417
 418       DO 80 I=1,11
 419
 420    80 FIRST(I)=.TRUE.
 421
 422 C Loop through event record
 423
 424       DO 410 I=1,NHEP
 425
 426       NEWSEC=.FALSE.
 427
 428 C First find start of new sections
 429
 430       IST=ISTHEP(I)
 431
 432       IS=IST/10
 433
 434       ID=IDHW(I)
 435
 436       IF (IST.EQ.101) THEN
 437
 438         NEWSEC=.TRUE.
 439
 440         SECTXT=TITLE(1)
 441
 442       ELSEIF (FIRST(2).AND.IS.EQ.12) THEN
 443
 444         NEWSEC=.TRUE.
 445
 446         SECTXT=TITLE(2)
 447
 448         FIRST(2)=.FALSE.
 449
 450       ELSEIF (FIRST(3).AND.IS.EQ.14) THEN
 451
 452         NEWSEC=.TRUE.
 453
 454         SECTXT=TITLE(3)
 455
 456         FIRST(3)=.FALSE.
 457
 458         FIRST(8)=.TRUE.
 459
 460         FIRST(9)=.TRUE.
 461
 462         FIRST(11)=.TRUE.
 463
 464       ELSEIF (FIRST(4).AND.IST.GE.158.AND.IST.NE.160
 465
 466      &                .AND.IST.LE.162) THEN
 467
 468         NEWSEC=.TRUE.
 469
 470         SECTXT=TITLE(4)
 471
 472         FIRST(4)=.FALSE.
 473
 474       ELSEIF (FIRST(5).AND.(IS.EQ.16.OR.IS.EQ.18)
 475
 476      &                .AND.IST.GT.162) THEN
 477
 478         NEWSEC=.TRUE.
 479
 480         SECTXT=TITLE(5)
 481
 482         FIRST(5)=.FALSE.
 483
 484       ELSEIF (IS.EQ.19.OR.IST.EQ.1.OR.IST.EQ.200) THEN
 485
 486         MS=ISTHEP(JMOHEP(1,I))/10
 487
 488         IF (MS.EQ.15.OR.MS.EQ.16.OR.MS.EQ.18) THEN
 489
 490           IF (FIRST(6)) THEN
 491
 492             NEWSEC=.TRUE.
 493
 494             SECTXT=TITLE(6)
 495
 496             FIRST(6)=.FALSE.
 497
 498           ENDIF
 499
 500         ELSEIF (FIRST(7).AND.(.NOT.FIRST(6))) THEN
 501
 502           NEWSEC=.TRUE.
 503
 504           SECTXT=TITLE(7)
 505
 506           FIRST(7)=.FALSE.
 507
 508         ENDIF
 509
 510       ELSEIF (FIRST(8).AND.(IST.EQ.125.OR.IST.EQ.155.OR.
 511
 512      &        (IST.EQ.123.AND.ISTHEP(JMOHEP(1,I)).EQ.199))) THEN
 513
 514         NEWSEC=.TRUE.
 515
 516         SECTXT=TITLE(8)
 517
 518         FIRST(3)=.TRUE.
 519
 520         FIRST(4)=.TRUE.
 521
 522         FIRST(5)=.TRUE.
 523
 524         FIRST(6)=.TRUE.
 525
 526         FIRST(7)=.TRUE.
 527
 528         FIRST(8)=.FALSE.
 529
 530       ELSEIF (FIRST(9).AND.(IST.EQ.123.OR.IST.EQ.124)) THEN
 531
 532         MS=ABS(IDHEP(JMOHEP(1,I)))
 533
 534         IF (MS.EQ.23.OR.MS.EQ.24.OR.MS.EQ.25) THEN
 535
 536           NEWSEC=.TRUE.
 537
 538           SECTXT=TITLE(9)
 539
 540           FIRST(3)=.TRUE.
 541
 542           FIRST(4)=.TRUE.
 543
 544           FIRST(5)=.TRUE.
 545
 546           FIRST(6)=.TRUE.
 547
 548           FIRST(7)=.TRUE.
 549
 550           FIRST(8)=.TRUE.
 551
 552           FIRST(9)=.FALSE.
 553
 554         ENDIF
 555
 556       ELSEIF (IST.EQ.170) THEN
 557
 558         NEWSEC=.TRUE.
 559
 560         SECTXT=TITLE(10)
 561
 562       ELSEIF (FIRST(11).AND.(ID.EQ.71.OR.ID.EQ.72)) THEN
 563
 564         NEWSEC=.TRUE.
 565
 566         SECTXT=TITLE(11)
 567
 568         FIRST(3)=.TRUE.
 569
 570         FIRST(11)=.FALSE.
 571
 572       ENDIF
 573
 574 C Prit out section heading
 575
 576       IF (NEWSEC) THEN
 577
 578         IF (PRVTX) THEN
 579
 580           IF (PRNDEF) THEN
 581
 582             IF (NPRFMT.EQ.1) THEN
 583
 584               WRITE(6, 90) SECTXT,(THEAD(J,3),J=1,17)
 585
 586             ELSE
 587
 588               WRITE(6,100) SECTXT,(THEAD(J,3),J=1,17)
 589
 590             ENDIF
 591
 592           ENDIF
 593
 594           IF (PRNTEX) WRITE(IUNITT,110) Z,Z,SECTXT,ZZ,Z,
 595
 596      &     (Z,THEAD(J,3),J=1,17),ZZ,Z
 597
 598           IF (PRNWEB) WRITE(IUNITW,120) TBCOLS(2),TBCOLS(3),
 599
 600      &     SECTXT,((THEAD(K,J),J=1,3),K=1,17)
 601
 602    90     FORMAT(/46X,A28//1X,A4,2X,A6,3X,A6,5A4,3(2X,A6),A6,3X,A5,
 603
 604      &     4(4X,A6))
 605
 606   100     FORMAT(/58X,A28//1X,A4,2X,A6,3X,A6,5A4,3(6X,A6),5X,A6,8X,A5,
 607
 608      &     4X,A6,2(5X,A6),6X,A6)
 609
 610   110     FORMAT(A1,'hline'/A1,'multicolumn{17}{|c|}{',A28,'} ',A2/A1,
 611
 612      &     'hline'/16(A1,'multicolumn{1}{|c|}{',A6,'} & '),
 613
 614      &     A1,'multicolumn{1}{|c|}{',A6,'} ',A2/A1,'hline')
 615
 616   120     FORMAT('<TR><TH COLSPAN=17 BGCOLOR=#',A6,'>',
 617
 618      &     '<FONT COLOR=#',A6,'>',A28,'</FONT></TH></TR>'/
 619
 620      &     '<TR>',17(/,1X,'<TH BGCOLOR=#',A6,'>
 621
 622      &     <FONT COLOR=',A6,'>',A6,'</FONT></TH>'),'</TR>')
 623
 624         ELSE
 625
 626           IF (PRNDEF) THEN
 627
 628             IF (NPRFMT.EQ.1) THEN
 629
 630               WRITE(6,130) SECTXT,(THEAD(J,3),J=1,13)
 631
 632             ELSE
 633
 634               WRITE(6,140) SECTXT,(THEAD(J,3),J=1,13)
 635
 636             ENDIF
 637
 638           END IF
 639
 640           IF (PRNTEX) WRITE(IUNITT,150) Z,Z,SECTXT,ZZ,Z,
 641
 642      &     (Z,THEAD(J,3),J=1,13),ZZ,Z
 643
 644           IF (PRNWEB) WRITE(IUNITW,160) TBCOLS(2),TBCOLS(3),
 645
 646      &     SECTXT,((THEAD(K,J),J=1,3),K=1,13)
 647
 648   130     FORMAT(/26X,A28//1X,A4,2X,A6,3X,A6,5A4,3(2X,A6),A6,3X,A5)
 649
 650   140     FORMAT(/36X,A28//1X,A4,2X,A6,3X,A6,5A4,3(6X,A6),5X,A6,8X,A5)
 651
 652   150     FORMAT(A1,'hline'/A1,'multicolumn{13}{|c|}{',A28,'} ',A2/A1,
 653
 654      &     'hline'/12(A1,'multicolumn{1}{|c|}{',A6,'} & '),
 655
 656      &     A1,'multicolumn{1}{|c|}{',A6,'} ',A2/A1,'hline')
 657
 658   160     FORMAT('<TR><TH COLSPAN=13 BGCOLOR=#',A6,'>',
 659
 660      &     '<FONT COLOR=#',A6,'>',A28,'</FONT></TH></TR>'/
 661
 662      &     '<TR>',13(/'<TH BGCOLOR=#',A6,'>',
 663
 664      &     '<FONT COLOR=#',A6,'>',A6,'</FONT></TH>'),'</TR>')
 665
 666         ENDIF
 667
 668       ENDIF
 669
 670 C Now print out the data line
 671
 672       IF (PRVTX) THEN
 673
 674 C     Include vertex information
 675
 676         IF (PRNDEF) THEN
 677
 678           IF (PRNDEC) THEN
 679
 680             IF (NPRFMT.EQ.1) THEN
 681
 682               WRITE(6,190) I,RNAME(IDHW(I)),IDHEP(I),IST,
 683
 684      &         JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 685
 686      &         (PHEP(J,I),J=1,5),(VHEP(J,I),J=1,4)
 687
 688             ELSE
 689
 690               WRITE(6,200) I,RNAME(IDHW(I)),IDHEP(I),IST,
 691
 692      &         JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 693
 694      &         (PHEP(J,I),J=1,5),(VHEP(J,I),J=1,4)
 695
 696             ENDIF
 697
 698           ELSE
 699
 700             IF (NPRFMT.EQ.1) THEN
 701
 702               WRITE(6,210) I,RNAME(IDHW(I)),IDHEP(I),IST,
 703
 704      &         JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 705
 706      &         (PHEP(J,I),J=1,5),(VHEP(J,I),J=1,4)
 707
 708             ELSE
 709
 710               WRITE(6,220) I,RNAME(IDHW(I)),IDHEP(I),IST,
 711
 712      &         JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 713
 714      &         (PHEP(J,I),J=1,5),(VHEP(J,I),J=1,4)
 715
 716             ENDIF
 717
 718           ENDIF
 719
 720         ENDIF
 721
 722         IF (PRNTEX) WRITE(IUNITT,230) I,TXNAME(1,IDHW(I)),IDHEP(I),
 723
 724      &   IST,JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 725
 726      &   (PHEP(J,I),J=1,5),(VHEP(J,I),J=1,4),ZZ
 727
 728         IF (PRNWEB) THEN
 729
 730           WRITE(IUNITW,240) TBCOLS(2),TBCOLS(3),I,I
 731
 732           IF (IDHEP(I).EQ.0.OR.IDHEP(I).EQ.91) THEN
 733
 734             WRITE(IUNITW,250) TXNAME(2,IDHW(I)),IDHEP(I),IST
 735
 736           ELSE
 737
 738             TMPNME=HWUNST(IDHW(I))
 739
 740             WRITE(FNAMEP,'(A15,A7,A5)')
 741
 742      &       'HW_decays/PART_',TMPNME,'.html'
 743
 744             WRITE(IUNITW,260) FNAMEP,TXNAME(2,IDHW(I)),IDHEP(I),IST
 745
 746           ENDIF
 747
 748           DO 170 J=1,2
 749
 750             IF (JMOHEP(J,I).NE.0) THEN
 751
 752               WRITE(IUNITW,270) JMOHEP(J,I),JMOHEP(J,I)
 753
 754             ELSE
 755
 756               WRITE(IUNITW,280) JMOHEP(J,I)
 757
 758             ENDIF
 759
 760   170     CONTINUE
 761
 762           DO 180 J=1,2
 763
 764             IF (JDAHEP(J,I).NE.0) THEN
 765
 766               WRITE(IUNITW,270) JDAHEP(J,I),JDAHEP(J,I)
 767
 768             ELSE
 769
 770               WRITE(IUNITW,280) JDAHEP(J,I)
 771
 772             ENDIF
 773
 774   180     CONTINUE
 775
 776           IF (NPRFMT.EQ.1) THEN
 777
 778             WRITE(IUNITW,290) (PHEP(J,I),J=1,5),(VHEP(J,I),J=1,4)
 779
 780           ELSE
 781
 782             WRITE(IUNITW,300) (PHEP(J,I),J=1,5),(VHEP(J,I),J=1,4)
 783
 784           ENDIF
 785
 786         ENDIF
 787
 788   190   FORMAT(1X,I4,1X,A8,I8,5I4,   2F8.2,2F7.1,F8.2,1P,4E10.3)
 789
 790   200   FORMAT(1X,I4,1X,A8,I8,5I4,   5F12.5,1P,4E11.4)
 791
 792   210   FORMAT(1X,Z4,1X,A8,I8,I4,4Z4,2F8.2,2F7.1,F8.2,1P,4E10.3)
 793
 794   220   FORMAT(1X,Z4,1X,A8,I8,I4,4Z4,5F12.5,1P,4E11.4)
 795
 796   230   FORMAT(I4,' & ',A37,' & $',I8,'$',5(' & ',I4),
 797
 798      &   5(' & $',F8.2,'$'),4(' & $',1P,E11.3,'$'),' ',A2)
 799
 800   240   FORMAT('<TR>'/'<TD BGCOLOR=#',A6,' ALIGN="RIGHT">',
 801
 802      &   '<FONT COLOR=#',A6,'><A NAME="',I4,'">',I4,'</A></FONT></TD>'/)
 803
 804   250   FORMAT('<TD ALIGN="CENTER">',A37,'</TD>'/'<TD ALIGN="RIGHT">',
 805
 806      &   I8,'</TD>'/'<TD ALIGN="RIGHT">',I4,'</TD>')
 807
 808   260   FORMAT('<TD ALIGN="CENTER"><A HREF="',A27,'">',A37,'</A></TD>'/
 809
 810      &   '<TD ALIGN="RIGHT">',I8,'</TD>'/
 811
 812      &   '<TD ALIGN="RIGHT">',I4,'</TD>')
 813
 814   270   FORMAT(/'<TD ALIGN="RIGHT"><A HREF="#',I4,'">',I4,'</A></TD>')
 815
 816   280   FORMAT(/'<TD ALIGN="RIGHT">',I4,'</TD>')
 817
 818   290   FORMAT(5(/'<TD ALIGN="RIGHT">',F8.2,'</TD>'),1P,
 819
 820      &   4(/'<TD ALIGN="RIGHT">',E10.3,'</TD>')/'</TR>')
 821
 822   300   FORMAT(5(/'<TD ALIGN="RIGHT">',F12.5,'</TD>'),1P,
 823
 824      &   4(/'<TD ALIGN="RIGHT">',E11.4,'</TD>')/'</TR>')
 825
 826       ELSE
 827
 828 C     Do not include vertex information
 829
 830         IF (PRNDEF) THEN
 831
 832           IF (PRNDEC) THEN
 833
 834             IF (NPRFMT.EQ.1) THEN
 835
 836               WRITE(6,330) I,RNAME(IDHW(I)),IDHEP(I),IST,
 837
 838      &         JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 839
 840      &         (PHEP(J,I),J=1,5)
 841
 842             ELSE
 843
 844               WRITE(6,340) I,RNAME(IDHW(I)),IDHEP(I),IST,
 845
 846      &         JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 847
 848      &         (PHEP(J,I),J=1,5)
 849
 850             ENDIF
 851
 852           ELSE
 853
 854             IF (NPRFMT.EQ.1) THEN
 855
 856               WRITE(6,350) I,RNAME(IDHW(I)),IDHEP(I),IST,
 857
 858      &          JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 859
 860      &          (PHEP(J,I),J=1,5)
 861
 862             ELSE
 863
 864               WRITE(6,360) I,RNAME(IDHW(I)),IDHEP(I),IST,
 865
 866      &         JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 867
 868      &         (PHEP(J,I),J=1,5)
 869
 870             ENDIF
 871
 872           ENDIF
 873
 874         ENDIF
 875
 876         IF (PRNTEX) THEN
 877
 878           IF (NPRFMT.EQ.1) THEN
 879
 880             WRITE(IUNITT,370) I,TXNAME(1,IDHW(I)),IDHEP(I),
 881
 882      &       IST,JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 883
 884      &       (PHEP(J,I),J=1,5),ZZ
 885
 886           ELSE
 887
 888             WRITE(IUNITT,380) I,TXNAME(1,IDHW(I)),IDHEP(I),
 889
 890      &       IST,JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I),
 891
 892      &       (PHEP(J,I),J=1,5),ZZ
 893
 894           ENDIF
 895
 896         ENDIF
 897
 898         IF (PRNWEB) THEN
 899
 900           WRITE(IUNITW,240) TBCOLS(2),TBCOLS(3),I,I
 901
 902           IF (IDHEP(I).EQ.0.OR.IDHEP(I).EQ.91) THEN
 903
 904             WRITE(IUNITW,250) TXNAME(2,IDHW(I)),IDHEP(I),IST
 905
 906           ELSE
 907
 908             TMPNME = HWUNST(IDHW(I))
 909
 910             WRITE(FNAMEP,'(A15,A7,A5)')
 911
 912      &       'HW_decays/PART_',TMPNME,'.html'
 913
 914             WRITE(IUNITW,260) FNAMEP,TXNAME(2,IDHW(I)),IDHEP(I),IST
 915
 916           ENDIF
 917
 918           DO 310 J=1,2
 919
 920             IF (JMOHEP(J,I).NE.0) THEN
 921
 922               WRITE(IUNITW,270) JMOHEP(J,I),JMOHEP(J,I)
 923
 924             ELSE
 925
 926               WRITE(IUNITW,280) JMOHEP(J,I)
 927
 928             ENDIF
 929
 930   310     CONTINUE
 931
 932           DO 320 J=1,2
 933
 934             IF (JDAHEP(J,I).NE.0) THEN
 935
 936               WRITE(IUNITW,270) JDAHEP(J,I),JDAHEP(J,I)
 937
 938             ELSE
 939
 940               WRITE(IUNITW,280) JDAHEP(J,I)
 941
 942             ENDIF
 943
 944   320     CONTINUE
 945
 946           IF (NPRFMT.EQ.1) THEN
 947
 948             WRITE(IUNITW,390) (PHEP(J,I),J=1,5)
 949
 950           ELSE
 951
 952             WRITE(IUNITW,400) (PHEP(J,I),J=1,5)
 953
 954           ENDIF
 955
 956         ENDIF
 957
 958   330   FORMAT(1X,I4,1X,A8,I8,5I4   ,2F8.2,2F7.1,F8.2)
 959
 960   340   FORMAT(1X,I4,1X,A8,I8,5I4   ,5F12.5)
 961
 962   350   FORMAT(1X,Z4,1X,A8,I8,I4,4Z4,2F8.2,2F7.1,F8.2)
 963
 964   360   FORMAT(1X,Z4,1X,A8,I8,I4,4Z4,5F12.5)
 965
 966   370   FORMAT(I4,' & ',A37,' & $',I8,'$',5(' & ',I4),
 967
 968      &   5(' & $',F8.2,'$'),' ',A2)
 969
 970   380   FORMAT(I4,' & ',A37,' & $',I8,'$',5(' & ',I4),
 971
 972      &   5(' & $',F12.5,'$'),' ',A2)
 973
 974   390   FORMAT(5(/'<TD ALIGN="RIGHT">',F8.2,'</TD>')/'</TR>')
 975
 976   400   FORMAT(5(/'<TD ALIGN="RIGHT">',F12.5,'</TD>')/'</TR>')
 977
 978       ENDIF
 979
 980   410 CONTINUE
 981
 982 C Close the files
 983
 984       IF (PRNTEX) THEN
 985
 986         WRITE(IUNITT,420) Z,Z,Z
 987
 988   420   FORMAT(A1,'end{longtable}'/A1,'end{center}'/A1,'end{document}')
 989
 990         CLOSE(IUNITT)
 991
 992       ENDIF
 993
 994       IF (PRNWEB) THEN
 995
 996         WRITE(IUNITW,430)
 997
 998   430   FORMAT('</TABLE>'/'</CENTER>'/'</BODY>'/'</HTML>')
 999
1000         CLOSE(IUNITW)
1001
1002       ENDIF
1003
1004       RETURN
1005
1006       END