JWeightEventToolkit.hh

Go to the documentation of this file.

#ifndef __JAANET__JWEIGHTEVENTTOOLKIT__
#define __JAANET__JWEIGHTEVENTTOOLKIT__

#include <vector>

#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/Evt.hh"

#include "JLang/JSharedPointer.hh"

#include "JAAnet/JHead.hh"
#include "JAAnet/JHeadToolkit.hh"
#include "JAAnet/JWeightEvent.hh"
#include "JAAnet/JFlux.hh"
#include "JAAnet/JFluxHelper.hh"
#include "JAAnet/JWeightDAQ.hh"
#include "JAAnet/JWeightMupage.hh"
#include "JAAnet/JWeightGSeaGen.hh"
#include "JAAnet/JWeightKM3BUU.hh"

/**
 * \author bjung
 */

namespace JAANET {

  using JLANG::JSharedPointer;
  using JLANG::JValueOutOfRange;


  /**
   * Look-up table for event weighters.
   */
  struct JEventWeighter :
    public std::vector< JSharedPointer<JWeightEvent> >
  {
    /*
     * Constructor
     */
    JEventWeighter()
    {
      this->push_back(new JWeightGSeaGen());
      this->push_back(new JWeightKM3BUU());
      this->push_back(new JWeightMupage());
      this->push_back(new JWeightDAQ());
    }

      
    /**
     * Get event weighter corresponding to given header.
     *
     * \param  header          header
     * \return                 event weighter
     */
    const JWeightEvent& operator()(const JHead& header) const
    {
      for (const_iterator i = this->begin(); i != this->end(); ++i) {
        
        if ((*i)->check(header)) {
          
          return *(*i);
        }
      }
      
      THROW(JValueOutOfRange, "JEventWeighter::operator(): No event weighter found for given header.");
    }
  };


  /**
   * Check if given event weighter contains a specifiable flux function.
   *
   * \param  weighter        event weighter
   * \return                 true if given event weighter contains a specifiable flux function; else false
   */
  inline bool hasFlux(const JWeightEvent& weighter)
  {
    return dynamic_cast<const JFluxHelper*>(&weighter) != NULL;
  }
  
  
  extern JEventWeighter getEventWeighter;               //!< Function object for mapping header to event weighter.
}

#endif