JSupport/JSummaryRouter.hh

Go to the documentation of this file.

#ifndef __JSUPPORT__JSUMMARYROUTER__
#define __JSUPPORT__JSUMMARYROUTER__
 
#include "km3net-dataformat/online/JDAQ.hh"
#include "km3net-dataformat/online/JDAQSummaryslice.hh"
#include "km3net-dataformat/online/JDAQModuleIdentifier.hh"
#include "km3net-dataformat/online/JDAQPMTIdentifier.hh"
 
#include "JTools/JRouter.hh"
#include "JLang/JException.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JSUPPORT {}
namespace JPP { using namespace JSUPPORT; }
 
namespace JSUPPORT {
 
  using KM3NETDAQ::JDAQSummaryslice;
  using KM3NETDAQ::JDAQSummaryFrame;
  using KM3NETDAQ::JDAQPMTIdentifier;
  using JLANG::JValueOutOfRange;
 
 
  /**
   * Router for fast addressing of summary data in KM3NETDAQ::JDAQSummaryslice data structure
   * as a function of the optical module identifier and PMT address.
   *
   * The member method update() should be used to update the internal router for a given event.
   * The member method getRate(const JDAQPMTIdentifier&) const can subsequently be used to obtain 
   * the measured singles rate for a given PMT.
   */
  class JSummaryRouter {
  public:
    /**
     * Default constructor.
     */
    JSummaryRouter() :
      router (-1),
      summary(NULL)
    {}
 
 
    /**
     * Update router.
     *
     * \param  ps              pointer to new summary slice
     */
    void update(const JDAQSummaryslice* ps)
    {
      // reset internal router
 
      if (summary != NULL) {
        for (JDAQSummaryslice::const_iterator i = summary->begin(); i != summary->end(); ++i) {
          router.put(i->getModuleID(), router.getDefaultAddress());
        }
      }
      
      summary = ps;
 
      if (summary != NULL) {
 
        // set internal router
 
        for (JDAQSummaryslice::const_iterator i = summary->begin(); i != summary->end(); ++i) {
          router.put(i->getModuleID(), distance(summary->begin(), i));
        }
      }
    }
 
 
    /**
     * Check validity of summary data.
     *
     * \return                 true if summary data available; else false
     */
    bool is_valid() const
    {
      return summary != NULL;
    }
 
 
    /**
     * Get summary slice.
     *
     * \return                 pointer to summary slice (may be NULL)
     */
    const JDAQSummaryslice* getSummaryslice() const
    {
      return summary;
    }
 
 
    /**
     * Get address of module.
     *
     * \param  module          module
     * \return                 address
     */
    const int getAddress(const JDAQModuleIdentifier& module) const
    {
      return router.get(module.getModuleID());
    }
    
 
    /**
     * Has summary frame.
     *
     * \param  module          module
     * \return                 true if module present; else false
     */
    bool hasSummaryFrame(const JDAQModuleIdentifier& module) const
    {
      return router.has(module.getModuleID());
    }
    
 
    /**
     * Get summary frame.
     *
     * \param  module          module
     * \return                 summary frame
     */
    const JDAQSummaryFrame& getSummaryFrame(const JDAQModuleIdentifier& module) const
    {
      if (router.has(module.getModuleID()))
        return (*summary)[getAddress(module)];
      else
        THROW(JValueOutOfRange, "Module identifier " << module);
    }
    
 
    /**
     * Get summary frame.
     *
     * \param  module          module
     * \param  rate_Hz         default rate [Hz]
     * \return                 summary frame
     */
    const JDAQSummaryFrame& getSummaryFrame(const JDAQModuleIdentifier& module, const double rate_Hz) const
    {
      if (router.has(module.getModuleID())) {
 
        return (*summary)[getAddress(module)];
 
      } else {
 
        frame.setModuleIdentifier(module);
        frame.setDAQFrameStatus(JDAQFrameStatus::getInstance());     // valid status
 
        for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
          frame.setRate(i, rate_Hz);
        }
 
        return frame;
      }
    }
 
 
    /**
     * Get rate.
     *
     * \param  id              PMT identifier
     * \param  rate_Hz         default rate [Hz]
     * \return                 rate [Hz]
     */
    double getRate(const JDAQPMTIdentifier& id, const double rate_Hz) const
    {
      double R = 0.0;
 
      if (this->hasSummaryFrame(id.getModuleIdentifier())) {
        R = this->getSummaryFrame(id.getModuleIdentifier()).getRate(id.getPMTAddress());
      }
 
      if (R < rate_Hz) {
        R = rate_Hz;
      }
      
      return R;
    }
 
  private:
    JTOOLS::JRouter<int>     router;
    const JDAQSummaryslice*  summary;
    mutable JDAQSummaryFrame frame;
  };
}
 
#endif