AlpgenInterface/src/alpg_debug.F

0001       SUBROUTINE DBPART
0002 C     Print the HEPEUP.
0003       IMPLICIT NONE
0004       integer iup
0005       INTEGER MAXNUP
0006       PARAMETER (MAXNUP=500)
0007       INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
0008       DOUBLE PRECISION XWGTUP,SCALUP,AQEDUP,AQCDUP,PUP,VTIMUP,SPINUP
0009       COMMON/HEPEUP/NUP,IDPRUP,XWGTUP,SCALUP,AQEDUP,AQCDUP,
0010      &              IDUP(MAXNUP),ISTUP(MAXNUP),MOTHUP(2,MAXNUP),
0011      &              ICOLUP(2,MAXNUP),PUP(5,MAXNUP),VTIMUP(MAXNUP),
0012      &              SPINUP(MAXNUP)
0013
0014       WRITE(42,*) '++++++++++++++++'
0015       WRITE(42,*) NUP
0016       DO iup=1,NUP
0017          WRITE(42,*) 'Particle #',iup,': IDUP =',IDUP(iup),
0018      &        ' ISTUP =',ISTUP(iup),
0019      &        ' MOTHUP(1) =',MOTHUP(1,iup),
0020      &        ' MOTHUP(2) =',MOTHUP(2,iup),
0021      &        ' ICOLUP(1) =',ICOLUP(1,iup),
0022      &        ' ICOLUP(2) =',ICOLUP(2,iup)
0023       ENDDO
0024
0025       END
0026
0027       SUBROUTINE DBINRD
0028 C     Print the list of INORAD'd particles
0029       IMPLICIT NONE
0030       INTEGER NMXHEP,NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
0031       DOUBLE PRECISION PHEP,VHEP
0032       PARAMETER (NMXHEP=4000)
0033       COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
0034      &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
0035
0036       INTEGER IUP,IEND,INORAD
0037       INTEGER MAXNUP
0038       PARAMETER (MAXNUP=500)
0039       COMMON/SHVETO/IEND,INORAD(MAXNUP)
0040
0041       WRITE(43,*) '++++++++++++++++'
0042       DO IUP=1,NHEP
0043          WRITE(43,*) 'Particle #',iup,': IDHEP =',IDHEP(iup),
0044      &        ' ISTHEP =',ISTHEP(iup),
0045      &        ' INORAD =',INORAD(iup)
0046       ENDDO
0047       END
0048
0049       SUBROUTINE DBCAL
0050 C     Print the particles as they go into the calorimeter.
0051       IMPLICIT NONE
0052       INTEGER NMXHEP,NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
0053       DOUBLE PRECISION PHEP,VHEP
0054       PARAMETER (NMXHEP=4000)
0055       COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
0056      &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
0057
0058       INTEGER IUP,IEND,INORAD
0059       INTEGER MAXNUP
0060       PARAMETER (MAXNUP=500)
0061       COMMON/SHVETO/IEND,INORAD(MAXNUP)
0062
0063       WRITE(44,*) '++++++++++++++++'
0064       DO IUP=1,NHEP
0065          WRITE(44,*) 'Particle #',iup,': IDHEP =',IDHEP(iup),
0066      &        ' ISTHEP =',ISTHEP(iup)
0067       ENDDO
0068       END