JGen2.hh

Go to the documentation of this file.

#ifndef __JASTRONOMY__JGEN2__
#define __JASTRONOMY__JGEN2__
 
#include <vector>
#include <numeric>
 
#include "JAstronomy/JPseudoExperiment.hh"
#include "JAstronomy/JAspera.hh"
 
 
/**
 * \file
 *
 * Set of pseudo experments.
 * \author mdejong
 */
namespace JASTRONOMY {}
namespace JPP { using namespace JASTRONOMY; }
 
namespace JASTRONOMY {
 
  /**
   * Auxiliary data structure to fit signal strength using likelihood ratio for multiple pseudo experiments.
   */
  struct JGen2 :
    public std::vector<JPseudoExperiment>,
    public JPseudoExperiment_t
  {
    /**
     * Set scaling factors of signal and background strengths.
     *
     * \param  fS        signal strength
     * \param  fB        background strength
     */
    void set(const double fS, const double fB)
    {
      for (auto& i : static_cast<std::vector<JPseudoExperiment>&>(*this)) {
        i.set(fS, fB);
      }
    }
 
 
    /**
     * Get total signal.
     *
     * \return           signal
     */
    double getSignal() const
    {
      return std::accumulate(this->begin(), this->end(), 0.0, [](const double y, const JPseudoExperiment& px) { return y + px.getSignal(); } );
    }
 
 
    /**
     * Get total background.
     *
     * \return           background
     */
    double getBackground() const
    {
      return std::accumulate(this->begin(), this->end(), 0.0, [](const double y, const JPseudoExperiment& px) { return y + px.getBackground(); } );
    }
 
 
    /**
     * Get fit method.
     *
     * \return            fit
     */
    virtual JAspera& getAspera() override
    {
      return this->fit;
    }
 
 
    /**
     * Generate pseudo experiment and transfer S/N values to fit method.
     *
     * \param  out        output
     * \return            result
     */
    virtual stats_type run(JAspera& out) const
    {
      stats_type px = { 0, 0 };
      
      for (const auto& i : static_cast<const std::vector<JPseudoExperiment>&>(*this)) {
        px += i.run(out);
      }
 
      return px;
    }
 
  protected:
    JAspera  fit;                                            //!< fit
  };
}
 
#endif