JOscProbInterpolator8D.hh

Go to the documentation of this file.

#ifndef __JOSCPROB__JOSCPROBINTERPOLATOR8D__
#define __JOSCPROB__JOSCPROBINTERPOLATOR8D__

#include <string>

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

#include "Jeep/JMessage.hh"

#include "JLang/JException.hh"

#include "JTools/JDistance.hh"

#include "JOscProb/JOscProbTable_t.hh"
#include "JOscProb/JOscProbTable.hh"
#include "JOscProb/JOscProbToolkit.hh"
#include "JOscProb/JOscChannel.hh"
#include "JOscProb/JOscParameters.hh"


namespace JOSCPROB {};
namespace JPP { using namespace JOSCPROB; }

namespace JOSCPROB {

  using JLANG::JValueOutOfRange;
  using JLANG::JNullPointerException;

  using JOSCPROB::JOscProbTable8D;
  

  /**
   * Auxiliary class for interpolating an oscillation probability table in 8 dimensions in terms of:\n
   * 1. \f$sin^2(\theta_{12})\f$
   * 2. \f$\Delta m_{21}^2\f$
   * 3. \f$sin^2(\theta_{13})\f$
   * 4. \f$\Delta m_{31}^2\f$
   * 5. \f$sin^2(\theta_{23})\f$
   * 6. \f$\delta_{\text{CP}}\f$
   * 7. \f$log_{10}(E)\f$
   * 8. \f$cos(\theta)\f$
   */
  template<class JDistance_t = JTOOLS::JDistance<typename JOscProbTable8D_t::JFunction2D_t::argument_type> >
  struct JOscProbInterpolator8D :
    public JOscParameters
  {
    typedef JOscProbTable8D<JDistance_t>                      table_type;
    typedef typename table_type::data_type                     data_type;
    
    
    /**
     * Default constructor.
     */
    JOscProbInterpolator8D() :
      pTable(),
      JOscParameters()
    {}


    /**
     * Constructor.
     *
     * \param  pTable              pointer to 8D oscillation probability table
     * \param  parameters          oscillation parameters            
     */
    JOscProbInterpolator8D(const table_type*     pTable,
                           const JOscParameters& parameters) :
      pTable        (pTable),
      JOscParameters(parameters)
    {}


    /**
     * Get oscillation probability for given oscillation channel.
     *
     * \param  channel            oscillation channel
     * \param  log10E             logarithmic neutrino energy
     * \param  costh              cosine zenith angle
     * \return                    oscillation probability
     */
    virtual double operator()(const JOscChannel& channel,
                              const double       log10E,
                              const double       costh) const
    {
      using namespace std;
      using namespace JPP;
      
      const JOscChannel* p = find(getOscChannel, getOscChannel + NUMBER_OF_OSCCHANNELS, channel);

      if (pTable != NULL) {
        
        if (p != end(getOscChannel) && is_valid()) {

          const size_t index = distance(getOscChannel, p);
        
          const data_type& probabilities = (*pTable)(this->sinsqTh12, this->dM21sq,
                                                     this->sinsqTh13, this->dM31sq,
                                                     this->sinsqTh23, this->deltaCP,
                                                     log10E,          costh);
          return probabilities[index];
        
        } else {

          if (is_valid()) {
            THROW(JValueOutOfRange, "JOscProbInterpolator8D::operator(): Invalid oscillation channel "    << channel    << endl);
          } else {
            THROW(JValueOutOfRange, "JOscProbInterpolator8D::operator(): Invalid oscillation parameters " << static_cast<const JOscParameters&>(*this) << endl);          
          }
        }
        
      } else {

        THROW(JNullPointerException, "JOscProbInterpolator8D::operator(): Unspecified oscillation probability table." << endl);
      }
    }

    
  private:

    const table_type* pTable; //!< Pointer to 8D oscillation probability table
  };
}


#endif