CDECK ID>, HWUEEC. *CMZ :- -26/04/91 14.22.30 by Federico Carminati *-- Author : Bryan Webber and Ian Knowles C----------------------------------------------------------------------- SUBROUTINE HWUEEC(IL) C----------------------------------------------------------------------- C Loads cross-section coefficients, for kinematically open channels, C in llbar-->qqbar; lepton label IL=1-6: e,nu_e,mu,nu_mu,tau,nu_tau. C----------------------------------------------------------------------- INCLUDE 'HERWIG61.INC' DOUBLE PRECISION Q2 INTEGER IL,JL,IQ Q2=EMSCA**2 JL=IL+10 MAXFL=0 TQWT=0. DO 10 IQ=1,NFLAV IF (EMSCA.GT.2.*RMASS(IQ)) THEN MAXFL=MAXFL+1 MAPQ(MAXFL)=IQ CALL HWUCFF(JL,IQ,Q2,CLQ(1,MAXFL)) TQWT=TQWT+CLQ(1,MAXFL) ENDIF 10 CONTINUE IF (MAXFL.EQ.0) CALL HWWARN('HWUEEC',100,*999) 999 END CDECK ID>, HWUEMV. *CMZ :- -30/06/94 19.31.08 by Mike Seymour *-- Author : Mike Seymour C----------------------------------------------------------------------- SUBROUTINE HWUEMV(N,IFROM,ITO) C----------------------------------------------------------------------- C MOVE A BLOCK OF ENTRIES IN THE EVENT RECORD C N ENTRIES IN HEPEVT STARTING AT IFROM ARE MOVED TO AFTER ITO C----------------------------------------------------------------------- INCLUDE 'HERWIG61.INC' INTEGER N,IFROM,ITO,IMAP(0:NMXHEP),LFROM,LTO,I,IEDT(NMXHEP),IHEP, $ JHEP,KHEP COMMON /HWUMAP/IMAP LFROM=IFROM LTO=ITO DO 100 I=1,N 100 IEDT(I)=LTO CALL HWUEDT(-N,IEDT) DO 300 I=1,N IHEP=LTO+I JHEP=IMAP(LFROM+I-1) ISTHEP(IHEP)=ISTHEP(JHEP) IDHW(IHEP)=IDHW(JHEP) IDHEP(IHEP)=IDHEP(JHEP) JMOHEP(1,IHEP)=JMOHEP(1,JHEP) JMOHEP(2,IHEP)=JMOHEP(2,JHEP) JDAHEP(1,IHEP)=JDAHEP(1,JHEP) JDAHEP(2,IHEP)=JDAHEP(2,JHEP) CALL HWVEQU(5,PHEP(1,JHEP),PHEP(1,IHEP)) CALL HWVEQU(4,VHEP(1,JHEP),VHEP(1,IHEP)) DO 200 KHEP=1,NHEP IF (JMOHEP(1,KHEP).EQ.JHEP) JMOHEP(1,KHEP)=IHEP IF (JMOHEP(2,KHEP).EQ.JHEP) JMOHEP(2,KHEP)=IHEP IF (JDAHEP(1,KHEP).EQ.JHEP) JDAHEP(1,KHEP)=IHEP IF (JDAHEP(2,KHEP).EQ.JHEP) JDAHEP(2,KHEP)=IHEP 200 CONTINUE IEDT(I)=JHEP 300 CONTINUE CALL HWUEDT(N,IEDT) 999 END CDECK ID>, HWUEPR. *CMZ :- -27/07/99 13.33.03 by Mike Seymour *-- Author : Ian Knowles, Bryan Webber & Kosuke Odagiri C----------------------------------------------------------------------- SUBROUTINE HWUEPR C----------------------------------------------------------------------- C Prints out event data in a number of possible formats: C If (PRNDEF) ASCII to stout C If (PRNTEX) LaTeX to the file HWEV_*******.tex C Please check paper size and offsets given in mm C Uses the package longtable.sty C If (PRVTX>OR.NPRFMT.EQ.2) designed to be printed C as landscape C If (PRNWEB) HTML to the file HWEV_*******.html C Call HWUDPR to create particle property files in C the subdirectory HW_decays/ C ******* gives the event number 0000001 etc. C----------------------------------------------------------------------- INCLUDE 'HERWIG61.INC' INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,I,IST,IS,ID,MS,J,K,IUNITW, & IUNITT CHARACTER*1 Z CHARACTER*2 ZZ CHARACTER*6 BGCOLS(5),TBCOLS(3),THEAD(17,3) CHARACTER*7 HWUNST,TMPNME CHARACTER*16 FNAMET CHARACTER*17 FNAMEW CHARACTER*27 FNAMEP CHARACTER*28 TITLE(11),SECTXT LOGICAL FIRST(11),NEWSEC COMMON/PAPER/MMWIDE,MMLONG,MMHOFF,MMVOFF EXTERNAL HWUNST C DATA BGCOLS/'ffffff','0000aa','aa0000','00aa00','aa00ff'/ DATA TBCOLS/'ccccff','9966ff','ffff00'/ DATA THEAD/ 17*'9966ff',17*'ffff00', & 'IHEP ',' ID ',' IDPDG',' IST ',' MO1 ',' MO2 ', & ' DA1 ',' DA2 ',' P-X ',' P-Y ',' P-Z ','ENERGY', & ' MASS ',' V-X ',' V-Y ',' V-Z ',' V-C*T'/ DATA TITLE/' ---INITIAL STATE--- ', & ' ---HARD SUBPROCESS--- ', & ' ---PARTON SHOWERS--- ', & ' ---GLUON SPLITTING--- ', & ' ---CLUSTER FORMATION--- ', & ' ---CLUSTER DECAYS--- ', & ' ---STRONG HADRON DECAYS--- ', & ' ---HEAVY PARTICLE DECAYS---', & ' ---H/W/Z BOSON DECAYS--- ', & ' ---SOFT UNDERLYING EVENT---', & ' ---MULTIPLE SCATTERING--- '/ Z=CHAR(92) ZZ=Z//Z C IUNITT=50 IUNITW=51 C Write out any required file header information TMPNME=HWUNST(NEVHEP) IF (PRNTEX) THEN WRITE(FNAMET,'(A5,A7,A4)') 'HWEV_',TMPNME,'.tex' OPEN(IUNITT,STATUS='UNKNOWN',FILE=FNAMET) IF (PRVTX.OR.NPRFMT.EQ.2) THEN WRITE(IUNITT,10) Z,Z,Z,MMLONG,Z,MMWIDE,Z,MMVOFF,Z,MMHOFF,Z,Z,Z ELSE WRITE(IUNITT,10) Z,Z,Z,MMWIDE,Z,MMLONG,Z,MMHOFF,Z,MMVOFF,Z,Z,Z ENDIF ENDIF IF (PRNWEB) THEN WRITE(FNAMEW,'(A5,A7,A5)') 'HWEV_',TMPNME,'.html' OPEN(IUNITW,STATUS='UNKNOWN',FILE=FNAMEW) WRITE(IUNITW,20) BGCOLS ENDIF 10 FORMAT(A1,'documentclass{article}'/A1,'usepackage{longtable}'/ & A1,'textwidth ',I4,'mm ',A1,'textheight ',I4,'mm'/ & A1,'hoffset ',I4,'mm ',A1,'voffset ',I4,'mm'/ & A1,'pagestyle{empty}'/A1,'begin{document}'/A1,'begin{center}') 20 FORMAT(''/''/'
'/ & '', & '', & 'HERWIG 6.0 | '/ & 'Beam 1: | '/ & 'Beam 2: | '/ & 'Seeds: | '/ & 'Status: | '/'',I4,' | '/'|
---|---|---|---|---|---|---|
Process: | '/'',I6,' | '/ & '',F8.2,' GeV/c | '/'',F8.2,' GeV/c | '/ & '',I11,' | '/ & 'Error: | '/'',I4,' | '/'
Event: | '/'',I7,' | '/ & '',A37,' | '/ & '',A37,' | '/ & '',I11,' | '/ & 'Weight: | '/'',1P,E11.4,' | '/'
'/ & '
', & '',A28,' | ||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
& ',A6,' | '),'||||||||||||||||
', & '',A28,' | ||||||||||||||||
', & '',A6,' | '),'||||||||||||||||
', & '',I4,' | '/) 250 FORMAT('',A37,' | '/'', & I8,' | '/'',I4,' | ') 260 FORMAT('',A37,' | '/ & '',I8,' | '/ & '',I4,' | ') 270 FORMAT(/'',I4,' | ') 280 FORMAT(/'',I4,' | ') 290 FORMAT(5(/'',F8.2,' | '),1P, & 4(/'',E10.3,' | ')/'',F12.5,' | '),1P, & 4(/'',E11.4,' | ')/'') ELSE C Do not include vertex information IF (PRNDEF) THEN IF (PRNDEC) THEN IF (NPRFMT.EQ.1) THEN WRITE(6,330) I,RNAME(IDHW(I)),IDHEP(I),IST, & JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I), & (PHEP(J,I),J=1,5) ELSE WRITE(6,340) I,RNAME(IDHW(I)),IDHEP(I),IST, & JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I), & (PHEP(J,I),J=1,5) ENDIF ELSE IF (NPRFMT.EQ.1) THEN WRITE(6,350) I,RNAME(IDHW(I)),IDHEP(I),IST, & JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I), & (PHEP(J,I),J=1,5) ELSE WRITE(6,360) I,RNAME(IDHW(I)),IDHEP(I),IST, & JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I), & (PHEP(J,I),J=1,5) ENDIF ENDIF ENDIF IF (PRNTEX) THEN IF (NPRFMT.EQ.1) THEN WRITE(IUNITT,370) I,TXNAME(1,IDHW(I)),IDHEP(I), & IST,JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I), & (PHEP(J,I),J=1,5),ZZ ELSE WRITE(IUNITT,380) I,TXNAME(1,IDHW(I)),IDHEP(I), & IST,JMOHEP(1,I),JMOHEP(2,I),JDAHEP(1,I),JDAHEP(2,I), & (PHEP(J,I),J=1,5),ZZ ENDIF ENDIF IF (PRNWEB) THEN WRITE(IUNITW,240) TBCOLS(2),TBCOLS(3),I,I IF (IDHEP(I).EQ.0.OR.IDHEP(I).EQ.91) THEN WRITE(IUNITW,250) TXNAME(2,IDHW(I)),IDHEP(I),IST ELSE TMPNME = HWUNST(IDHW(I)) WRITE(FNAMEP,'(A15,A7,A5)') & 'HW_decays/PART_',TMPNME,'.html' WRITE(IUNITW,260) FNAMEP,TXNAME(2,IDHW(I)),IDHEP(I),IST ENDIF DO 310 J=1,2 IF (JMOHEP(J,I).NE.0) THEN WRITE(IUNITW,270) JMOHEP(J,I),JMOHEP(J,I) ELSE WRITE(IUNITW,280) JMOHEP(J,I) ENDIF 310 CONTINUE DO 320 J=1,2 IF (JDAHEP(J,I).NE.0) THEN WRITE(IUNITW,270) JDAHEP(J,I),JDAHEP(J,I) ELSE WRITE(IUNITW,280) JDAHEP(J,I) ENDIF 320 CONTINUE IF (NPRFMT.EQ.1) THEN WRITE(IUNITW,390) (PHEP(J,I),J=1,5) ELSE WRITE(IUNITW,400) (PHEP(J,I),J=1,5) ENDIF ENDIF 330 FORMAT(1X,I4,1X,A8,I8,5I4 ,2F8.2,2F7.1,F8.2) 340 FORMAT(1X,I4,1X,A8,I8,5I4 ,5F12.5) 350 FORMAT(1X,Z4,1X,A8,I8,I4,4Z4,2F8.2,2F7.1,F8.2) 360 FORMAT(1X,Z4,1X,A8,I8,I4,4Z4,5F12.5) 370 FORMAT(I4,' & ',A37,' & $',I8,'$',5(' & ',I4), & 5(' & $',F8.2,'$'),' ',A2) 380 FORMAT(I4,' & ',A37,' & $',I8,'$',5(' & ',I4), & 5(' & $',F12.5,'$'),' ',A2) 390 FORMAT(5(/'',F8.2,' | ')/'') 400 FORMAT(5(/'',F12.5,' | ')/'') ENDIF 410 CONTINUE C Close the files IF (PRNTEX) THEN WRITE(IUNITT,420) Z,Z,Z 420 FORMAT(A1,'end{longtable}'/A1,'end{center}'/A1,'end{document}') CLOSE(IUNITT) ENDIF IF (PRNWEB) THEN WRITE(IUNITW,430) 430 FORMAT('