fermi_function.h

Go to the documentation of this file.

/* Author(s) :    Francois Mauger <mauger@lpccaen.in2p3.fr>
 * Creation date: 2014-06-06
 * Last modified: 2014-06-06
 *
 * License:
 * Copyright 2014 F. Mauger
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or (at
 * your option) any later version.
 *
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor,
 * Boston, MA 02110-1301, USA.
 *
 * Description:
 *
 *   Fermi function.
 *
 * History:
 *
 */

#ifndef GENBB_HELP_FERMI_FUNCTION_H
#define GENBB_HELP_FERMI_FUNCTION_H 1

namespace genbb {

  class fermi_function
  {
  public:

    enum mode_type {
      MODE_NON_RELATIVISTIC = 0x1,      
      MODE_SPHERICALLY_SYMMETRIC = 0x2, 
      MODE_SHAPE_ONLY = 0x10            
    };

    enum decay_type {
      DT_INVALID    = 0,
      DT_BETA_MINUS = 1,
      DT_BETA_PLUS  = 2
    };

    bool is_valid() const;

    bool is_beta_plus() const;

    bool is_beta_minus() const;

    fermi_function();

    fermi_function(int z_,
                   int a_,
                   decay_type beta_type_,
                   unsigned int mode_,
                   double ke_cut_ = -1.0);

    void set(int z_,
             int a_,
             decay_type beta_type_,
             unsigned int mode_,
             double ke_cut_ = -1.0);

    void reset();

    double operator()(double ke_) const;

    static double eval_non_relativistic_approx(int z_daughter_,
                                               decay_type type_,
                                               double ke_cut_,
                                               bool shape_only_,
                                               double ke_);

    static double eval_spherically_symmetric_approx(int z_daughter_,
                                                    int a_daughter_,
                                                    decay_type type_,
                                                    double ke_cut_,
                                                    bool shape_only_,
                                                    double ke_);

  private:

    unsigned int       _mode_;       
    int                _Z_;          
    int                _A_;          
    decay_type         _decay_type_; 
    double             _ke_cut_;     
  };

} // end of namespace genbb

#endif // GENBB_HELP_FERMI_FUNCTION_H

// Local Variables: --
// mode: c++ --
// End: --