Merge branch 'master' of https://git.cern.ch/reps/AliRoot

[u/mrichter/AliRoot.git] / TEvtGen / EvtGenModels / begevtgenstorex.F

C--------------------------------------------------------------------------
C     
C     Environment:
C     This software is part of the EvtGen package developed jointly
C     for the BaBar and CLEO collaborations.  If you use all or part
C     of it, please give an appropriate acknowledgement.
C     
C     Copyright Information: See EvtGen/COPYRIGHT
C     Copyright (C) 1998      Caltech, UCSB
C     
C     Module: begevtgenstorex.F
C     
C     Description:
C     
C     Modification history:
C     
C     DJL/RYD     August 11, 1998         Module created
C     
C------------------------------------------------------------------------
      subroutine begevtgenstorex(entry,eventnum,numparticle,istat,
     +     partnum,mother,daugfirst,dauglast,
     +     px,py,pz,e,m,x,y,z,t)
      implicit none
*     
*     routine to fill the stdhep common blocks from
*     evtgen (C++). This routine allows the C++ program not to
*     have to mess with common blocks.
*     
*     Anders Ryd,  Dec 96   Created.
*     
*     
#include "EvtGenModels/stdhep.inc"
      
      logical qedrad
      integer ph_nmxhep ! this is parameter nmxhep in photos/photos_make
*                       ! Renamed here to avoid name conflict in stdhep.inc
      parameter (ph_nmxhep=10000)
      common / phoqed / qedrad(ph_nmxhep)

      integer entry,eventnum,numparticle,istat,partnum
      integer mother,daugfirst,dauglast
      integer stdhepnum, i
      double precision px,py,pz,e,m,x,y,z,t

      stdhepnum=partnum

      nevhep=eventnum
      nhep=numparticle
      isthep(entry)=istat
      idhep(entry)=stdhepnum
      jmohep(1,entry)=mother
      jmohep(2,entry)=0
      jdahep(1,entry)=daugfirst
      jdahep(2,entry)=dauglast
      phep(1,entry)=px
      phep(2,entry)=py
      phep(3,entry)=pz
      phep(4,entry)=e
      phep(5,entry)=m
      vhep(1,entry)=x
      vhep(2,entry)=y
      vhep(3,entry)=z
      vhep(4,entry)=t

      qedrad(entry)=.true.
      if (daugfirst.gt.0 .and. dauglast.gt.0) THEN
         do i=daugfirst, dauglast
            qedrad(i) = .true.
         end do
      end if

      return

      end