TEvtGen/EvtGenModels/begevtgenstorex.F

   1 C--------------------------------------------------------------------------
   2 C
   3 C     Environment:
   4 C     This software is part of the EvtGen package developed jointly
   5 C     for the BaBar and CLEO collaborations.  If you use all or part
   6 C     of it, please give an appropriate acknowledgement.
   7 C
   8 C     Copyright Information: See EvtGen/COPYRIGHT
   9 C     Copyright (C) 1998      Caltech, UCSB
  10 C
  11 C     Module: begevtgenstorex.F
  12 C
  13 C     Description:
  14 C
  15 C     Modification history:
  16 C
  17 C     DJL/RYD     August 11, 1998         Module created
  18 C
  19 C------------------------------------------------------------------------
  20       subroutine begevtgenstorex(entry,eventnum,numparticle,istat,
  21      +     partnum,mother,daugfirst,dauglast,
  22      +     px,py,pz,e,m,x,y,z,t)
  23       implicit none
  24 *
  25 *     routine to fill the stdhep common blocks from
  26 *     evtgen (C++). This routine allows the C++ program not to
  27 *     have to mess with common blocks.
  28 *
  29 *     Anders Ryd,  Dec 96   Created.
  30 *
  31 *
  32 #include "EvtGenModels/stdhep.inc"
  33
  34       logical qedrad
  35       integer ph_nmxhep ! this is parameter nmxhep in photos/photos_make
  36 *                       ! Renamed here to avoid name conflict in stdhep.inc
  37       parameter (ph_nmxhep=10000)
  38       common / phoqed / qedrad(ph_nmxhep)
  39
  40       integer entry,eventnum,numparticle,istat,partnum
  41       integer mother,daugfirst,dauglast
  42       integer stdhepnum, i
  43       double precision px,py,pz,e,m,x,y,z,t
  44
  45       stdhepnum=partnum
  46
  47       nevhep=eventnum
  48       nhep=numparticle
  49       isthep(entry)=istat
  50       idhep(entry)=stdhepnum
  51       jmohep(1,entry)=mother
  52       jmohep(2,entry)=0
  53       jdahep(1,entry)=daugfirst
  54       jdahep(2,entry)=dauglast
  55       phep(1,entry)=px
  56       phep(2,entry)=py
  57       phep(3,entry)=pz
  58       phep(4,entry)=e
  59       phep(5,entry)=m
  60       vhep(1,entry)=x
  61       vhep(2,entry)=y
  62       vhep(3,entry)=z
  63       vhep(4,entry)=t
  64
  65       qedrad(entry)=.true.
  66       if (daugfirst.gt.0 .and. dauglast.gt.0) THEN
  67          do i=daugfirst, dauglast
  68             qedrad(i) = .true.
  69          end do
  70       end if
  71
  72       return
  73
  74       end
  75
  76