C...Very simple tabulation of total and partial cross section
C...parametrizations in Pythia, for gamma-p events.


      COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200) 
      COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
      COMMON/PYINT1/MINT(400),VINT(400)
      COMMON/PYINT7/SIGT(0:6,0:6,0:5)

C...Common parameters.
      MSTP(122)=0
      PARP(2)=2.
      CONV=1E3
      IOFF=1
      WRITE(6,*) '      ECM     TOTAL      VDM     ELASTIC  DIF(GAM)  ',
     &' DIF(P)   DOUBLEDIF  NONDIF'      

C...Find results for gamma-p.
      DO 150 IECM=1,20
      ECM=4.*1.3**IECM
      MSTP(14)=1
      CALL PYINIT('CMS','GAMMA','P',ECM)
      SIGTOT=CONV*SIGT(0,0,0)
      MSTP(14)=2
      CALL PYINIT('CMS','GAMMA','P',ECM)
      SIGVDM=CONV*SIGT(0,0,0)
      SIGEL=CONV*SIGT(0,0,1)
      SIGSDA=CONV*SIGT(0,0,2)
      SIGSDB=CONV*SIGT(0,0,3)
      SIGDD=CONV*SIGT(0,0,4)
      SIGND=CONV*SIGT(0,0,5)
      WRITE(6,'(F10.2,7F10.3)') ECM,SIGTOT,SIGVDM,SIGEL,SIGSDA,SIGSDB,
     &SIGDD,SIGND
  150 CONTINUE

      END