4 *CMZ :- -27/07/99 13.33.03 by Mike Seymour
6 *-- Author : Ian Knowles, Bryan Webber & Kosuke Odagiri
8 C-----------------------------------------------------------------------
12 C-----------------------------------------------------------------------
14 C Prints out particle properies/decay tables in a number of formats:
16 C If (PRNDEF) ASCII to stout
18 C If (PRNTEX) LaTeX to the file HW_decays.tex
20 C Paper size and offsets as set in HWUEPR
22 C Uses the package longtable.sty
24 C Designed to be printed as landscape
26 C If (PRNWEB) HTML to the file HW_decays/index.html
28 C /PART0000001.html etc.
30 C-----------------------------------------------------------------------
32 INCLUDE 'HERWIG61.INC'
34 INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,IUNITT,IUNTW1,IUNTW2,I,NM,J,K,
42 CHARACTER*3 ASPIN(0:10)
44 CHARACTER*6 BGCOLS(5),TBCOLS(3)
46 CHARACTER*7 HWUNST,TMPNME
52 COMMON/PAPER/MMWIDE,MMLONG,MMHOFF,MMVOFF
56 DATA BGCOLS/'ffffff','0000aa','aa0000','00aa00','aa00ff'/
58 DATA TBCOLS/'ccccff','9966ff','ffff00'/
60 DATA ASPIN/' 0 ','1/2',' 1 ','3/2',' 2 ','5/2',' 3 ','7/2',
78 C Open and write out file header information for index file
96 OPEN(IUNITT,STATUS='UNKNOWN',FILE='HW_decays.tex')
100 WRITE(IUNITT,30) Z,Z,Z,MMWIDE,Z,MMLONG,Z,MMHOFF,Z,MMVOFF,
102 & Z,Z,Z,Z,Z,ZZ,Z,Z,Z,Z,Z,Z,NRES,ZZ,Z,Z,ZZ,Z,Z
106 WRITE(IUNITT,40) Z,Z,Z,MMLONG,Z,MMWIDE,Z,MMHOFF,Z,MMVOFF,
108 & Z,Z,Z,Z,Z,Z,ZZ,Z,Z,Z,Z,Z,Z,NRES,ZZ,Z,Z,Z,ZZ,Z,Z
116 OPEN(IUNTW1,STATUS='UNKNOWN',FILE='HW_decays/index.html')
118 WRITE(IUNTW1,50) BGCOLS,TBCOLS,NRES,((TBCOLS(I),I=2,3),J=1,7)
122 10 FORMAT(1H1//15X,'TABLE OF PROPERTIES OF',I4,' PARTICLES USED'/)
124 20 FORMAT(1H1//30X,'TABLE OF PROPERTIES OF',I4,' PARTICLES USED'//
126 & 5X,'Name IDPDG Mass Chg Spn Lifetime Modes ',
128 & ' Branching fractions ME codes and decay products')
130 30 FORMAT(A1,'documentclass{article}'/A1,'usepackage{longtable}'/
132 & A1,'textwidth ',I4,'mm ',A1,'textheight ',I4,'mm'/
134 & A1,'hoffset ',I4,'mm ',A1,'voffset ',I4,'mm'/
136 & A1,'pagestyle{empty}'/A1,'begin{document}'/
138 & A1,'begin{center}'/A1,'begin{longtable}{|r|c|r|r|r|r|r|r|}'/
140 & A1,'hline'/'Id HW & Name & Id PDG & Mass & Charge & Spin ',
142 & '& Lifetime & Modes ',A2/A1,'hline'/
144 & A1,'endhead'/A1,'hline'/A1,'endfoot'/A1,'hline'/
146 & A1,'multicolumn{8}{|c|}{HERWIG 6.0: Table of properties',
148 & ' of the ',I3,' particles used} ',A2/A1,'hline',A1,'hline'/
150 & 'Id HW & Name & Id PDG & Mass & Charge & Spin & ',
152 & 'Lifetime & Modes ',A2/A1,'hline'/A1,'endfirsthead')
154 40 FORMAT(A1,'documentclass{article}'/A1,'usepackage{longtable}'/
156 & A1,'textwidth ',I4,'mm ',A1,'textheight ',I4,'mm'/
158 & A1,'hoffset ',I4,'mm ',A1,'voffset ',I4,'mm'/
160 & A1,'pagestyle{empty}'/A1,'begin{document}'/A1,'begin{center}'/
162 & A1,'begin{longtable}{|r|c|r|r|r|r|r|r|c|r|ccccc|}'/
164 & A1,'hline'/'Id HW & Name & Id PDG & Mass & Charge & Spin ',
166 & '& Lifetime & Modes & B.R. & M.E. & ' /
168 & A1,'multicolumn{5}{|c|}{Decay Products} ',A2/A1,'hline'/
170 & A1,'endhead'/A1,'hline'/A1,'endfoot'/A1,'hline'/
172 & A1,'multicolumn{15}{|c|}{HERWIG 6.0: Table of properties',
174 & ' of the ',I3,' particles used} ',A2/A1,'hline',A1,'hline'/
176 & 'Id HW & Name & Id PDG & Mass & Charge & Spin & ',
178 & 'Lifetime & Modes & B.R. & M.E. & '/
180 & A1,'multicolumn{5}{|c|}{Decay Products} ',A2/A1,'hline'/
184 50 FORMAT('<!-- No Less productions -->'/'<HTML>'/'<HEAD>'/
186 & '<TITLE>HERWIG 6.0 Particle Properties</TITLE>'/'</HEAD>'/
188 & '<BODY BGCOLOR=#',A6,' TEXT=#',A6,' LINK=#',A6,
190 & ' ALINK=#',A6,' VLINK=#',A6,'>'/'<CENTER>'/
192 & '<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>',
194 & '<TR>'/'<TH COLSPAN=8 BGCOLOR=#',A6,' ALIGN="CENTER">',
196 & '<A HREF=="http://hepwww.rl.ac.uk/theory/seymour/herwig/">',
198 & 'HERWIG 6.0:</A><FONT COLOR=#',A6,'> Table of properties of',
200 & ' the ',I3,' particles used</FONT></TH>'/'<TR>'/'<TH></TH>'/
202 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Name</FONT></TH>'/
204 & '<TH BGCOLOR=#',A6,' ALIGN="CENTER"><FONT COLOR=#',A6,'>',
206 & 'Id PDG</FONT></TH>'/
208 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Mass</FONT></TH>'/
210 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Charge</FONT></TH>'/
212 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Spin</FONT></TH>'/
214 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Lifetime</FONT></TH>'/
216 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Modes</FONT></TH>'/
220 C Loop through resonances
224 C Skip particles that can't be produced or blank lines
226 IF ((VTOCDK(I).AND.VTORDK(I)).OR.
228 & (RNAME(I).EQ.' ')) GOTO 260
230 C Open and write out header information for particle file
236 WRITE(FNAMEP,'(A5,A7,A5)') 'PART_',TMPNME,'.html'
238 WRITE(FNAMEW,'(A,A17)') 'HW_decays/',FNAMEP
240 OPEN(IUNTW2,STATUS='UNKNOWN',FILE=FNAMEW)
242 WRITE(IUNTW2,60) RNAME(I),BGCOLS
244 WRITE(IUNTW2,70) TBCOLS,((TBCOLS(L),L=2,3),M=1,6)
248 60 FORMAT('<!-- No Less productions -->'/'<HTML>'/'<HEAD>'/
250 & '<TITLE>HERWIG 6.0: ',A8,' properties</TITLE>'/'</HEAD>'/
252 & '<BODY BGCOLOR=#',A6,' TEXT=#',A6,' LINK=#',A6,
254 & ' ALINK=#',A6,' VLINK=#',A6,'>'/'<CENTER>')
256 70 FORMAT('<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>'/
258 & '<TR>'/'<TH></TH>'/
260 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Name</FONT></TH>'/
262 & '<TH BGCOLOR=#',A6,' ALIGN="CENTER"><FONT COLOR=#',A6,
264 & '>Id PDG</FONT></TH>'/
266 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Mass</FONT></TH>'/
268 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Charge</FONT></TH>'/
270 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Spin</FONT></TH>'/
272 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Lifetime</FONT></TH>'/
274 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>Modes</FONT></TH>'/
278 C Trick to output charge in fractions for di/s - quarks
280 IF ((I.GE. 1.AND.I.LE. 12).OR.(I.GE.109.AND.I.LE.120).OR.
282 & (I.GE.209.AND.I.LE.218).OR.(I.GE.401.AND.I.LE.424)) THEN
292 C Write out special particles with no decay modes
294 IF (NMODES(I).EQ.0) THEN
298 IF (NPRFMT.LE.1) THEN
300 WRITE(6,80) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
302 & ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
306 WRITE(6,90) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
308 & ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
314 C Add particle to LaTeX file
318 IF (NPRFMT.LE.1) THEN
320 WRITE(IUNITT,100) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
322 & ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0,ZZ
326 WRITE(IUNITT,110) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
328 & ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0,Z,ZZ
336 C Add properties to Web index
338 WRITE(IUNTW1,120) TBCOLS(2),TBCOLS(3),I,FNAMEP,TXNAME(2,I),
340 & IDPDG(I),RMASS(I),ICHRG(I),ACHRG,
342 & ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
344 C Add properties to Web particle file
346 WRITE(IUNTW2,130) TBCOLS(2),TBCOLS(3),I,TXNAME(2,I),
348 & IDPDG(I),RMASS(I),ICHRG(I),ACHRG,
350 & ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),0
354 80 FORMAT(/1X,I3,1X,A8,' IDPDG=',I8,', M=',F8.3,', Q=',I2,', J=',
356 & A3,', T=',1P,E9.3,',',I3,' Modes')
358 90 FORMAT(/1X,I3,1X,A8,1X,I8,1X,F8.3,1X,I2,1X,A3,1X,1P,E9.3,1X,I3)
360 100 FORMAT(A1,'hline',I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,
362 & A2,'$ & ',A3,' & $',1P,E9.3,'$ & ',I3,' ',A2)
364 110 FORMAT(A1,'cline{1-8}'/
366 & I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,A2,'$ & ',A3,
368 & ' & $',1P,E9.3,'$ & ',I3,' & ',A1,'multicolumn{7}{|c|}{} ',A2)
372 & '<TD ALIGN="RIGHT" BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',I3,
376 & '<TD ALIGN="CENTER"><A HREF="',A17,'">',A37,'</A></TD>'/
378 & '<TD ALIGN="RIGHT">',I8,'</TD>'/
380 & '<TD ALIGN="RIGHT">',F8.3,'</TD>'/
382 & '<TD ALIGN="RIGHT">',I2,A2,'</TD>'/
384 & '<TD ALIGN="RIGHT">',A3,'</TD>'/
386 & '<TD ALIGN="RIGHT">',1P,E9.3,'</TD>'/
388 & '<TD ALIGN="RIGHT">',I3,'</TD>'/'</TR>')
392 & '<TD ALIGN="RIGHT" BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',I3,
396 & '<TD ALIGN="CENTER">',A37,'</TD>'/
398 & '<TD ALIGN="RIGHT">',I8,'</TD>'/
400 & '<TD ALIGN="RIGHT">',F8.3,'</TD>'/
402 & '<TD ALIGN="RIGHT">',I2,A2,'</TD>'/
404 & '<TD ALIGN="RIGHT">',A3,'</TD>'/
406 & '<TD ALIGN="RIGHT">',1P,E9.3,'</TD>'/
408 & '<TD ALIGN="RIGHT">',I3,'</TD>'/'</TR>'/'</TABLE>'/'<P>')
412 C Particle with decay modes
418 ELSEIF (VTOCDK(I)) THEN
430 C Write out properties and first decay mode
434 IF (NPRFMT.LE.1) THEN
436 WRITE(6, 80) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
438 & ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM
440 WRITE(6,140) (RNAME(IDKPRD(L,K)),L=1,5),BRFRAC(K),NME(K)
444 WRITE(6,150) I,RNAME(I),IDPDG(I),RMASS(I),ICHRG(I),
446 & ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM,BRFRAC(K),NME(K),
448 & (RNAME(IDKPRD(L,K)),L=1,5)
456 IF (NPRFMT.LE.1) THEN
458 WRITE(IUNITT,160) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
460 & ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM,ZZ,Z
462 WRITE(IUNITT,170) Z,Z,(TXNAME(1,IDKPRD(L,K)),L=1,5),Z,
464 & BRFRAC(K),Z,NME(K),ZZ
468 WRITE(IUNITT,180) Z,I,TXNAME(1,I),IDPDG(I),RMASS(I),
470 & ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM,
472 & BRFRAC(K),NME(K),(TXNAME(1,IDKPRD(L,K)),L=1,5),ZZ,Z
480 C Add properties to index
482 WRITE(IUNTW1,120) TBCOLS(2),TBCOLS(3),I,FNAMEP,TXNAME(2,I),
484 & IDPDG(I),RMASS(I),ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),
488 C Add properties to Web particle file
490 WRITE(IUNTW2,130) TBCOLS(2),TBCOLS(3),I,TXNAME(2,I),IDPDG(I),
492 & RMASS(I),ICHRG(I),ACHRG,ASPIN(INT(TWO*RSPIN(I))),RLTIM(I),NM
494 WRITE(IUNTW2,190) TBCOLS,TXNAME(2,I),
496 & ((TBCOLS(L),L=2,3),M=1,3)
498 WRITE(IUNTW2,200) TBCOLS(2),TBCOLS(3),1,BRFRAC(K),NME(K),
500 & (TXNAME(2,IDKPRD(L,K)),L=1,5)
504 140 FORMAT(5X,'BR[ -->',5(1X,A6),']=',F5.3,', ME code',I4)
506 150 FORMAT(/1X,I3,1X,A8,1X,I8,1X,F8.3,1X,I2,1X,A3,1X,1P,E9.3,1X,I3,
508 & 2X,F5.3,1X,I3,5(1X,A8))
510 160 FORMAT(A1,'hline',
512 & I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,A2,'$ & ',
514 & A3,' & $',1P,E9.3,'$ & ',I3,' ',A2/A1,'cline{2-8}')
516 170 FORMAT(' & & ',A1,'multicolumn{2}{l}{$',A1,'longrightarrow$'/
518 & 5(A37,' '),'}'/' & ',A1,'multicolumn{2}{l}{BR = ',F5.3,'} & ',
520 & A1,'multicolumn{2}{l|}{ME code = ',I3,'} ',A2)
522 180 FORMAT(A1,'hline'/
524 & I4,' & ',A37,' & $',I8,'$ & ',F8.3,' & $',I2,A2,'$ & ',
526 & A3,' & $',1P,E9.3,'$ & ',I3,' & ',F5.3,' & ',I3,
528 & 5(' & ',A37), ' ',A2/A1,'cline{2-8}')
530 190 FORMAT('<TABLE ALIGN="CENTER" BGCOLOR=#',A6,'>'/'<TR>'/
532 & '<TH COLSPAN=8 BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',A37,
534 & ' Decay Modes</FONT></TH>'/'</TR>'/'<TR>'/'<TH></TH>',
536 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>B.R.</FONT></TH>'/
538 & '<TH BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>M.E.</FONT></TH>'/
540 & '<TH BGCOLOR=#',A6,' ALIGN="CENTER" COLSPAN=5>',
542 & '<FONT COLOR=#',A6,'>Decay products</FONT></TH>'/'</TR>')
546 & '<TD ALIGN="RIGHT" BGCOLOR=#',A6,'><FONT COLOR=#',A6,'>',
550 & '<TD ALIGN="RIGHT">',F5.3,'</TD>'/
552 & '<TD ALIGN="RIGHT">',I3,'</TD>'/
554 & 5('<TD ALIGN="CENTER">',A37,'</TD>'/),'</TR>')
556 C Write out additional decay modes
558 IF (NMODES(I).GE.2) THEN
566 IF (NPRFMT.LE.1) THEN
568 WRITE(6,140) (RNAME(IDKPRD(L,K)),L=1,5),BRFRAC(K),NME(K)
572 WRITE(6,220) BRFRAC(K),NME(K),(RNAME(IDKPRD(L,K)),L=1,5)
580 IF (NPRFMT.LE.1) THEN
582 WRITE(IUNITT,170) Z,Z,(TXNAME(1,IDKPRD(L,K)),L=1,5),Z,
584 & BRFRAC(K),Z,NME(K),ZZ
588 WRITE(IUNITT,230) Z,BRFRAC(K),NME(K),
590 & (TXNAME(1,IDKPRD(L,K)),L=1,5),ZZ
596 IF (PRNWEB) WRITE(IUNTW2,200) TBCOLS(2),TBCOLS(3),J,
598 & BRFRAC(K),NME(K),(TXNAME(2,IDKPRD(L,K)),L=1,5)
602 IF (PRNTEX.AND.NPRFMT.EQ.2.AND.NMODES(I+1).EQ.0)
604 & WRITE(IUNITT,240) Z
606 220 FORMAT(54X,F5.3,1X,I3,5(1X,A8))
608 230 FORMAT(' & ',A1,'multicolumn{7}{|c|}{} & ',F5.3,' & ',I3,
610 & 5(' & ',A37),' ',A2)
612 240 FORMAT(A1,'hline')
618 C Close Web particle file
628 250 FORMAT('</TABLE>'/'</CENTER>'/'<P>'/
630 & 'Main particle <A HREF="index.html">index</A>'/
632 & '</BODY>'/'</HTML>')
636 C Close the LaTeX file
640 WRITE(IUNITT,270) Z,Z,Z
646 C Close the index file
656 270 FORMAT(A1,'end{longtable}'/A1,'end{center}'/A1,'end{document}')
658 280 FORMAT('</TABLE>'/'</CENTER>'/'</BODY>'/'</HTML>')