GeneratorInterface/AlpgenInterface/src/alpg_debug.F - master - githubmirror/cmssw

GeneratorInterface/AlpgenInterface/src/alpg_debug.F
68 строк2 KB
Giulio Eulisse
Snapshot of CMSSW_6_2_0_pre8.
26 июн 2013, 18:12
26 июн 2013, 18:12214ab73
      SUBROUTINE DBPART
C     Print the HEPEUP.
      IMPLICIT NONE
      integer iup
      INTEGER MAXNUP
      PARAMETER (MAXNUP=500)
      INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
      DOUBLE PRECISION XWGTUP,SCALUP,AQEDUP,AQCDUP,PUP,VTIMUP,SPINUP
      COMMON/HEPEUP/NUP,IDPRUP,XWGTUP,SCALUP,AQEDUP,AQCDUP,
     &              IDUP(MAXNUP),ISTUP(MAXNUP),MOTHUP(2,MAXNUP),
     &              ICOLUP(2,MAXNUP),PUP(5,MAXNUP),VTIMUP(MAXNUP),
     &              SPINUP(MAXNUP)

      WRITE(42,*) '++++++++++++++++'
      WRITE(42,*) NUP
      DO iup=1,NUP
         WRITE(42,*) 'Particle #',iup,': IDUP =',IDUP(iup),
     &        ' ISTUP =',ISTUP(iup),
     &        ' MOTHUP(1) =',MOTHUP(1,iup),
     &        ' MOTHUP(2) =',MOTHUP(2,iup),
     &        ' ICOLUP(1) =',ICOLUP(1,iup),
     &        ' ICOLUP(2) =',ICOLUP(2,iup)
      ENDDO

      END

      SUBROUTINE DBINRD
C     Print the list of INORAD'd particles
      IMPLICIT NONE
      INTEGER NMXHEP,NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
      DOUBLE PRECISION PHEP,VHEP
      PARAMETER (NMXHEP=4000)
      COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
     &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)

      INTEGER IUP,IEND,INORAD
      INTEGER MAXNUP
      PARAMETER (MAXNUP=500)
      COMMON/SHVETO/IEND,INORAD(MAXNUP)

      WRITE(43,*) '++++++++++++++++'
      DO IUP=1,NHEP
         WRITE(43,*) 'Particle #',iup,': IDHEP =',IDHEP(iup),
     &        ' ISTHEP =',ISTHEP(iup),
     &        ' INORAD =',INORAD(iup)
      ENDDO
      END

      SUBROUTINE DBCAL
C     Print the particles as they go into the calorimeter.
      IMPLICIT NONE
      INTEGER NMXHEP,NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
      DOUBLE PRECISION PHEP,VHEP
      PARAMETER (NMXHEP=4000)
      COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
     &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)

      INTEGER IUP,IEND,INORAD
      INTEGER MAXNUP
      PARAMETER (MAXNUP=500)
      COMMON/SHVETO/IEND,INORAD(MAXNUP)

      WRITE(44,*) '++++++++++++++++'
      DO IUP=1,NHEP
         WRITE(44,*) 'Particle #',iup,': IDHEP =',IDHEP(iup),
     &        ' ISTHEP =',ISTHEP(iup)
      ENDDO
      END