      REAL FUNCTION broken1()
*********************************************************
*                                                       *
* This file was generated by HUWFUN.                    *
*                                                       *
*********************************************************
*
*     Ntuple Id:      51
*     Ntuple Title:   Atlfast51
*     Creation:       03/07/2002 17.59.26
*
*********************************************************
*
      LOGICAL         CHAIN
      CHARACTER*128   CFILE
      INTEGER         IDNEVT,NCHEVT,ICHEVT
      REAL            OBS(13)
*
      COMMON /PAWIDN/ IDNEVT,OBS
      COMMON /PAWCHN/ CHAIN, NCHEVT, ICHEVT
      COMMON /PAWCHC/ CFILE
*
*--   Ntuple Variable Declarations
*
      REAL PXELE(12),PYELE(12),PZELE(12),EEELE(12),PXMUO(12),PYMUO(12)
     + ,PZMUO(12),EEMUO(12),PXPHO(12),PYPHO(12),PZPHO(12),EEPHO(12)
     + ,PXMUX(12),PYMUX(12),PZMUX(12),EEMUX(12),PXJET(20),PYJET(20)
     + ,PZJET(20),EEJET(20),PTcalo(20),PTbjet(20),PTujet(20),PXJETB(20)
     + ,PYJETB(20),PZJETB(20),EEJETB(20),PXPART(40),PYPART(40)
     + ,PZPART(40),EEPART(40),PXMISS,PYMISS,PXNUE,PYNUE
      INTEGER NELE,KFELE(12),NMUO,KFMUO(12),NPHO,KFPHO(12),NMUX
     + ,KFMUX(12),NJET,KFJET(20),NJETB,KFJETB(20),NPART,KFPART(40),ISUB
     + ,JETB,JETC,JETTAU
*
      COMMON /PAWCR4/ NELE,KFELE,PXELE,PYELE,PZELE,EEELE,NMUO,KFMUO
     + ,PXMUO,PYMUO,PZMUO,EEMUO,NPHO,KFPHO,PXPHO,PYPHO,PZPHO,EEPHO,NMUX
     + ,KFMUX,PXMUX,PYMUX,PZMUX,EEMUX,NJET,KFJET,PXJET,PYJET,PZJET,EEJET
     + ,PTcalo,PTbjet,PTujet,NJETB,KFJETB,PXJETB,PYJETB,PZJETB,EEJETB
     + ,NPART,KFPART,PXPART,PYPART,PZPART,EEPART,ISUB,JETB,JETC,JETTAU
     + ,PXMISS,PYMISS,PXNUE,PYNUE
*

*
*--   Enter user code here
*
      dyevt = 1.
*
      call hfill(1110, float(nmuo), 0., 1.)
      call hfill(1101, float(njet), 0., 1.)
      if (nmuo.eq.2)
      	    temu = eemuo(1) + eemuo(2)
            tpxmu = pxmuo(1) + pxmuo(2)
            tpymu = pymuo(1) + pymuo(2)
            tpzmu = pzmuo(1) + pzmuo(2)
            tpmu = sqrt(tpxmu**2 + tpymu**2 + tpzmu**2)
            ommuos = sqrt(temu**2 - tpmu**2)
            call hfill(1102, ommuos, 0., 1.)      
      endif
      if (njet.eq.2) then
      	    tejet = ejet(1) + ejet(2)
            tpxjet = pxjet(1) + pxjet(2)
            tpyjet = pyjet(1) + pyjet(2)
            tpzjet = pzjet(1) + pzjet(2)
            tpjet = sqrt(tpxjet**2 + tpyjet**2 + tpzjet**2)
            omjets = sqrt(tejet**2 - tpjet**2)
            call hfill(1103, omjets, 0., 1.)      
      endif
      END
