demo2.F

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file demo2.F

C***********************************************************************
C...
      PROGRAM demo

C...Demonstration job for LEPTO 6.5
      COMMON /leptou/ cut(14),lst(40),parl(30),x,y,w2,q2,u
      common/lujets/n,k(4000,5),p(4000,5),v(4000,5)
      DATA nprint/10/,nevent/100/

C...Settings
      cut(5)=10.
      cut(7)=25.
      
C...Init
      CALL timex(t1)
      CALL linit(0,-11,-27.5,820.0,4)
      CALL timex(t2)

      DO 500 ne=1,nevent
  100 CALL lepto
      IF(lst(21).NE.0) THEN
        WRITE(6,*) ' Warning: LST(21)=',lst(21),' event skipped'
        goto 100
      ENDIF
      IF(ne.LE.nprint) THEN
        CALL lulist(1)
        WRITE(6,*) ' x,y,Q2,W= ',x,y,q2,sqrt(w2)
      ENDIF
  500 CONTINUE
  600 CONTINUE
      CALL timex(t3)
      WRITE(6,1000) parl(23),parl(24),nevent,t2-t1,t3-t2,(t3-t2)/nevent
 1000 FORMAT(/,' Cross section from numerical integration:',g12.3,' pb',
     &       /,' Cross section from MonteCarlo simulation:',g12.3,' pb',
     &       /,' Number of events generated:              ',i12,
     &       /,' Time for initialization:         ',f12.3,' s',
     &       /,' Time for event generation,     total:    ',f12.3,' s',
     &       /,'                            per event:    ',f12.3,' s')
      END