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;
  

  /**
   * Auxiliary class for interpolating an oscillation probability table in 8 dimensions in terms of:
   *
   * 1. \f$ \sin^{2}(\theta_{12}) \f$
   * 2. \f$ \Delta m_{21}^{2}     \left[\mathrm{eV}^{2}\right] \f$
   * 3. \f$ \sin^{2}(\theta_{13}) \f$
   * 4. \f$ \Delta m_{31}^{2}     \left[\mathrm{eV}^{2}\right] \f$
   * 5. \f$ \sin^{2}(\theta_{23}) \f$
   * 6. \f$ \delta_{\mathrm{CP}}  \left[\pi\right] \f$
   * 7. \f$ E                     \left[\mathrm{GeV}\right]  \f$
   * 8. \f$ \sin(\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  E                  neutrino energy     [GeV]
     * \param  costh              cosine zenith angle
     * \return                    oscillation probability
     */
    virtual double operator()(const JOscChannel& channel,
                              const double       E,
                              const double       costh) const
    {
      using namespace std;
      using namespace JPP;

      if (pTable != NULL && is_valid()) {
      
        const JOscChannel* p = find(getOscChannel, getOscChannel + NUMBER_OF_OSCCHANNELS, channel);
        
        if (p != end(getOscChannel)) {

          const double L     = getL(costh);
          const double sinth = (fabs(costh) < 1.0 ? sqrt((1 - costh) * (1 + costh)) : 0.0);

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

          THROW(JValueOutOfRange, "JOscProbInterpolator8D::operator(): Invalid oscillation channel " << channel << endl);
        }
        
      } else {

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

    
  private:

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


#endif