JOscParameters.hh

Go to the documentation of this file.

#ifndef __JOSCPROB__JOSCPARAMETERS__
#define __JOSCPROB__JOSCPARAMETERS__

#include <istream>
#include <ostream>
#include <iomanip>

#include "Jeep/JPrint.hh"
#include "Jeep/JProperties.hh"


/**
 * \author bjung, mdejong
 */

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

namespace JOSCPROB {

  using JLANG::JEquationParameters;
  

  /**
   * Data structure for oscillation parameters.
   */
  struct JOscParameters
  {
    /**
     * Default constructor.\n
     * \n
     * Values taken from the NuFIT three-flavour global analysis best fit values in:\n
     * https://arxiv.org/abs/2111.03086?context=hep-ex\n
     * for normal ordering including the Super-Kamiokande atmospheric data.
     */
    JOscParameters() :
      sinsqTh12(0.304),
      dM21sq   (7.42e-5),
      sinsqTh13(0.02246),
      dM31sq   (2.510e-3),
      sinsqTh23(0.450),
      deltaCP  (1.28)
    {}


    /**
     * Constructor.
     * 
     * \param  sinsqTh12       Squared sine of the PMNS mixing angle between the first and second neutrino mass eigenstates [-]
     * \param  dM21sq          Squared mass difference between the first and second neutrino mass eigenstates               [eV2]
     * \param  sinsqTh13       Squared sine of the PMNS mixing angle between the first and third neutrino mass eigenstates  [-]
     * \param  dM31sq          Squared mass difference between the first and third neutrino mass eigenstates                [eV2]
     * \param  sinsqTh23       Squared sine of the PMNS mixing angle between the second and third neutrino mass eigenstates [-]
     * \param  deltaCP         PMNS phase angle                                                                             [rad]
     */
    JOscParameters(const double sinsqTh12,
                   const double dM21sq,
                   const double sinsqTh13,
                   const double dM31sq,
                   const double sinsqTh23,
                   const double deltaCP) :
      sinsqTh12(sinsqTh12),
      dM21sq   (dM21sq),
      sinsqTh13(sinsqTh13),
      dM31sq   (dM31sq),
      sinsqTh23(sinsqTh23),
      deltaCP  (deltaCP)
    {}


    /**
     * Get oscillation parameters.
     *
     * \return                 oscillation parameters
     */
    const JOscParameters& getOscParameters() const
    {
      return static_cast<const JOscParameters&>(*this);
    }


    /**
     * Set oscillation parameters.
     *
     * \param  parameters      oscillation parameters
     */
    void setOscParameters(const JOscParameters& parameters)
    {
      static_cast<JOscParameters&>(*this) = parameters;
    }


    /**
     * Check validity of oscillation parameters.
     *
     * \return                 true if valid; else false
     */
    bool is_valid() const
    {
      if (this->sinsqTh12 < 0.0 ||
          this->dM21sq    < 0.0 ||
          this->sinsqTh13 < 0.0 ||
          this->dM31sq    < 0.0 ||
          this->sinsqTh23 < 0.0) {
        return false;
      }

      return true;
    }


    /**
     * Stream input of oscillation parameters.
     *
     * \param  in              input stream
     * \param  object          oscillation parameters
     * \return                 input stream
     */
    friend inline std::istream& operator>>(std::istream& in, JOscParameters& object)
    {
      JProperties properties(object.getProperties());

      in >> properties;

      object.setProperties(properties);

      return in;
    }


    /**
     * Stream output of oscillation parameters.
     *
     * \param  out             output stream
     * \param  object          oscillation parameters
     * \return                 output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const JOscParameters& object)
    {
      return out << object.getProperties();
    }


    /**
     * Get equation parameters.
     *
     * \return                 equation parameters
     */
    static inline JEquationParameters& getEquationParameters()
    {
      static JEquationParameters equation("=", "\n\r;,", "./", "#");

      return equation;
    }
      

    /**
     * Set equation parameters.
     *
     * \param  equation        equation parameters
     */
    static inline void setEquationParameters(const JEquationParameters& equation)
    {
      getEquationParameters() = equation;
    }


    /**
     * Get properties of this class.
     *
     * \param  equation        equation parameters
     */
    JProperties getProperties(const JEquationParameters& equation = JOscParameters::getEquationParameters())
    {
      return JOscParametersHelper(*this, equation);
    }


    /**
     * Get properties of this class.
     *
     * \param  equation        equation parameters
     */
    JProperties getProperties(const JEquationParameters& equation = JOscParameters::getEquationParameters()) const
    {
      return JOscParametersHelper(*this, equation);
    }


    /**
     * Set properties of this class
     *
     * \param  properties      properties
     */
    void setProperties(const JProperties& properties)
    {
      this->sinsqTh12 = properties.getValue<double>("sinsqTh12");
      this->dM21sq    = properties.getValue<double>("dM21sq");
      this->sinsqTh13 = properties.getValue<double>("sinsqTh13");
      this->dM31sq    = properties.getValue<double>("dM31sq");
      this->sinsqTh23 = properties.getValue<double>("sinsqTh23");
      this->deltaCP   = properties.getValue<double>("deltaCP");
    }


    double sinsqTh12;  //!< Squared sine of the PMNS mixing angle between the first and second neutrino mass eigenstates [-]
    double dM21sq;     //!< Squared mass difference between the first and second neutrino mass eigenstates               [eV2]
    double sinsqTh13;  //!< Squared sine of the PMNS mixing angle between the first and third neutrino mass eigenstates  [-]
    double dM31sq;     //!< Squared mass difference between the first and third neutrino mass eigenstates                [eV2]
    double sinsqTh23;  //!< Squared sine of the PMNS mixing angle between the second and third neutrino mass eigenstates [-]
    double deltaCP;    //!< PMNS phase angle                                                                             [pi * rad]

  private:
    /**
     * Auxiliary class for I/O of oscillation parameters.
     */
    struct JOscParametersHelper :
      public JProperties
    {
      /**
       * Constructor.
       *
       * \param  object        oscillation parameters
       * \param  equation      equation    parameters 
       */
      template<class JOscParameters_t>
      JOscParametersHelper(JOscParameters_t&          object,
                           const JEquationParameters& equation) :
        JProperties(equation, 1)
      {
        this->insert(gmake_property(object.sinsqTh12));
        this->insert(gmake_property(object.dM21sq));
        this->insert(gmake_property(object.sinsqTh13));
        this->insert(gmake_property(object.dM31sq));
        this->insert(gmake_property(object.sinsqTh23));
        this->insert(gmake_property(object.deltaCP));   
      }
    };
  };
}

#endif