C     =================
      SUBROUTINE MUOFIL
C     =================

      IMPLICIT NONE

#include "partap.inc"
#include "mfrtz.inc"
#include "zescommon.inc"
#include "zestmpcommon.inc"

      INTEGER BMuon1, BMuon2, BMnrmips
      INTEGER nomu
      REAL    mupt(2),muth(2),muph(2), BMmmtrk

C -----------------------------------
C --- reconstruct BMUON variables ---
C -----------------------------------
      ZESTMP_Bmuon1  = 0
      ZESTMP_Bmuon2  = 0
      ZES_BMmmtrk = 0.0
      CALL MBZES(Bmuon1, BMuon2, BMmmtrk)
      ZESTMP_Bmuon1  = Bmuon1
      ZESTMP_Bmuon2  = Bmuon2
C 17.12.98 AFM : Bug fix. Created
C new integer variables BMuon1 and BMuon2
      ZES_Bmuon1  = Bmuon1
      ZES_Bmuon2  = Bmuon2

      ZES_BMmmtrk = BMmmtrk

C -----------------------------------
C --- reconstruct FMUON variables ---
C -----------------------------------
      ZESTMP_FMuon1  = 0
      ZES_FMNtrk  = 0
      ZES_FMNtrp  = 0
      ZES_FMRecp  = 0.0
      ZES_FMAngr  = 0
      ZES_FMMatch = 0.0
      ZESTMP_FMuon1  = MIN(COUTAB(MFRTZ),1)
      CALL MFFASE(ZES_FMNtrk,
     &            ZES_FMNtrp,
     &            ZES_FMRecp,
     &            ZES_FMAngr,
     &            ZES_FMMatch)

C -------------------------------
C --- run CAL+CTD muon finder ---
C -------------------------------
      ZES_nomu     = 0
      ZES_mupt1    = 0.0
      ZES_muth1    = 0.0
      ZES_muph1    = 0.0
      ZES_mupt2    = 0.0
      ZES_muth2    = 0.0
      ZES_muph2    = 0.0
      ZES_BMnrmips = 0

      CALL MUONFIND(nomu,mupt,muth,muph,BMnrmips)
      ZES_nomu     = nomu
      ZES_mupt1    = mupt(1)
      ZES_muth1    = muth(1)
      ZES_muph1    = muph(1)
      ZES_mupt2    = mupt(2)
      ZES_muth2    = muth(2)
      ZES_muph2    = muph(2)

      ZES_BMnrmips = BMnrmips

      RETURN
      END