SUBROUTINE HERA C...Initialize parameters in Ariadne and Jetset CALL ARTUNE('4.10') C...Call a user supplied routine setting C...the parameters and switches in LEPTO CALL LEPSET C...Call a user supplied routine setting C...the parameters and switches in Ariadne CALL ARISET C...Initialize Ariadne to run with LEPTO CALL ARINIT('LEPTO') C...Initialize LEPTO for HERA CALL LINIT(0,11,-26.5,820.0,4) C...Loop over a number of events DO 100 IEVE=1,10 C...Generate an event with LEPTO CALL LEPTO C...Apply the Dipole Cascade CALL AREXEC C...Call a user supplied analysis routine CALL HERANA 100 CONTINUE END SUBROUTINE LEPSET COMMON /LEPTOU/ CUT(14),LST(40),PARL(30),X,Y,W2,XQ2,U SAVE /LEPTOU/ COMMON /PYPARS/ MSTP(200),PARP(200),MSTI(200),PARI(200) SAVE /PYPARS/ C...Use structure funtions from Pythia LST(15)=0 C...Or rather use the LDC version of pystfu giving the structure C...functions evolved with LDC (Make a dummy call to pystfu to ensure C...correct version of pystfu is linked.) IF (MSTP(51).EQ.940801) CALL PYSTFU MSTP(51)=0 C...Switch off hadronization LST(7)=0 C...Suppress printouts from LEPTO LST(3)=0 C...Set x, Q2 and W2 ranges in Lepto CUT(1) = 0.00008 CUT(2) = 0.5 CUT(5) = 5.0 CUT(6) = 1280.0 CUT(7) = 1000.0 RETURN END SUBROUTINE ARISET COMMON /ARDAT1/ PARA(40),MSTA(40) SAVE /ARDAT1/ COMMON /ARHIDE/ PHAR(400),MHAR(400) SAVE /ARHIDE/ COMMON /LDCDAT/ MSTLDC(100),PARLDC(100) C...Master switch for using LDC in Ariadne MSTA(32)=-32 C...Some parameters in LDCMC PARLDC(4)=1.0 PARLDC(6)=0.2 RETURN END SUBROUTINE HERANA C...Trivial analysis - print the event CALL LULIST(1) RETURN END SUBROUTINE DUMMY9 C...This is a dummy routine to ensure the relevant block data routines C...are linked from the JETSET/PYTHIA and LEPTO libraries. This routine C...should never be called. CALL LUDATA CALL PYDATA CALL LEPTOD RETURN END