TEvtGen/PHOTOS/phlupa.F

   1       SUBROUTINE PHLUPA(IPOINT)
   2 C.----------------------------------------------------------------------
   3 C.
   4 C.    PHLUPA:   debugging tool
   5 C.
   6 C.    Purpose:  NONE, eventually may printout content of the
   7 C.              /PHOEVT/ common
   8 C.
   9 C.    Input Parameters:   Common /PHOEVT/ and /PHNUM/
  10 C.                        latter may have number of the event.
  11 C.
  12 C.    Output Parameters:  None
  13 C.
  14 C.    Author(s):  Z. Was                          Created at:  30/05/93
  15 C.                                                Last Update: 10/08/93
  16 C.
  17 C.----------------------------------------------------------------------
  18       INTEGER NMXPHO
  19       PARAMETER (NMXPHO=10000)
  20       INTEGER IDPHO,ISTPHO,JDAPHO,JMOPHO,NEVPHO,NPHO
  21       REAL*8 PPHO,VPHO
  22       COMMON/PHOEVT/NEVPHO,NPHO,ISTPHO(NMXPHO),IDPHO(NMXPHO),
  23      &JMOPHO(2,NMXPHO),JDAPHO(2,NMXPHO),PPHO(5,NMXPHO),VPHO(4,NMXPHO)
  24       COMMON /PHNUM/ IEV
  25       INTEGER PHLUN
  26       COMMON/PHOLUN/PHLUN
  27       DIMENSION SUM(5)
  28       IF (IPOINT.LT.3000) RETURN
  29       IOUT=56
  30       IF (IEV.LT.1000) THEN
  31       DO I=1,5
  32         SUM(I)=0.0D0
  33       ENDDO
  34       WRITE(PHLUN,*) 'EVENT NR=',IEV,
  35      $            'WE ARE TESTING /PHOEVT/ at IPOINT=',IPOINT
  36       WRITE(PHLUN,10)
  37       I=1
  38       WRITE(PHLUN,20) IDPHO(I),PPHO(1,I),PPHO(2,I),PPHO(3,I),
  39      $                      PPHO(4,I),PPHO(5,I),JDAPHO(1,I),JDAPHO(2,I)
  40       I=2
  41       WRITE(PHLUN,20) IDPHO(I),PPHO(1,I),PPHO(2,I),PPHO(3,I),
  42      $                      PPHO(4,I),PPHO(5,I),JDAPHO(1,I),JDAPHO(2,I)
  43       WRITE(PHLUN,*) ' '
  44       DO I=3,NPHO
  45       WRITE(PHLUN,20) IDPHO(I),PPHO(1,I),PPHO(2,I),PPHO(3,I),
  46      $                      PPHO(4,I),PPHO(5,I),JMOPHO(1,I),JMOPHO(2,I)
  47         DO J=1,4
  48           SUM(J)=SUM(J)+PPHO(J,I)
  49         ENDDO
  50       ENDDO
  51       SUM(5)=SQRT(ABS(SUM(4)**2-SUM(1)**2-SUM(2)**2-SUM(3)**2))
  52       WRITE(PHLUN,30) SUM
  53  10   FORMAT(1X,'  ID      ','p_x      ','p_y      ','p_z      ',
  54      $                   'E        ','m        ',
  55      $                   'ID-MO_DA1','ID-MO DA2' )
  56  20   FORMAT(1X,I4,5(F9.3),2I9)
  57  30   FORMAT(1X,' SUM',5(F9.3))
  58       ENDIF
  59       END