*********************************************************************** subroutine geaevopen * * open event files for gukine * * last modified: 23/7/2004 by R. Holzmann *********************************************************************** implicit none #include "geant321/gcbank.inc" #include "geant321/gclist.inc" #include "geant321/gckine.inc" #include "geaevent.inc" integer ifile, i1, ntrack, itype1, itype2, id, iw real pom1(4), pom2(4), v1(3), v2(3), e, b character*200 event_file, event_file_name do ifile=1,nevfiles event_file_name = evfilename(ifile) i1 = index(event_file_name,'/') if (i1.gt.0) then event_file = event_file_name else event_file=' ' endif open(unit=23+ifile,file=event_file,status='OLD') ipair(ifile) = 0 * ------------------------------------------------------- * check for di-electrons, in order to set flag ipair * used in gukine *-------------------------------------------------------- read(23+ifile,*) id,ntrack,e,b,iw ! read number of particles in event c if (ntrack.eq.2) then ! if 2, then check for e+e- c if (iw.le.2) then c read(23+ifile,*) pom1,itype1 c read(23+ifile,*) pom2,itype2 c else c read(23+ifile,*) pom1,v1,itype1 c read(23+ifile,*) pom2,v2,itype2 c endif c if ((itype1.eq.2 .or. itype1.eq.3.) .and. c & (itype2.eq.2 .or. itype2.eq.3.)) ipair(ifile) = 1 c endif rewind(unit=23+ifile) enddo call skipevents() ! skip events until IDEVT from /GCFLAG/ return end