alpha_decay.h

Go to the documentation of this file.

/* Author(s) :    Francois Mauger <mauger@lpccaen.in2p3.fr>
 * Creation date: 2014-11-08
 * Last modified: 2014-11-10
 *
 * 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:
 *
 *   Alpha decay driver.
 *
 * History:
 *
 */

#ifndef GENBB_HELP_ALPHA_DECAY_H
#define GENBB_HELP_ALPHA_DECAY_H 1

// Standard library:
#include <string>
#include <fstream>

// Third party:
// - Bayeux/datatools:
#include <datatools/properties.h>
#include <datatools/i_tree_dump.h>
// - Bayeux/mygsl:
#include <mygsl/rng.h>

// This project:
#include <genbb_help/base_decay_driver.h>
#include <genbb_help/genbb_macros.h>
#include <genbb_help/two_body_decay.h>

namespace genbb {

  // Forward declaration:
  class nuclear_level;

  class alpha_decay  : public base_decay_driver
  {

  public:

    alpha_decay();

    virtual ~alpha_decay();

    bool has_q_alpha() const;

    void set_q_alpha(double);

    double get_q_alpha() const;

    bool has_A() const;

    void set_A(int);

    int get_A() const;

    bool has_Z() const;

    void set_Z(int);

    int get_Z() const;

    int get_Z_parent() const;

    int get_Z_daughter() const;

    int get_A_parent() const;

    int get_A_daughter() const;

    void set_mass_parent(double);

    double get_mass_parent() const;

    void set_mass_daughter(double);

    double get_mass_daughter() const;

    void set_energy_parent(double);

    double get_energy_parent() const;

    void set_energy_daughter(double);

    double get_energy_daughter() const;

    bool is_alpha_generated() const;

    void set_alpha_generated(bool);

    bool is_daughter_generated() const;

    void set_daughter_generated(bool);

    void set_levels(const nuclear_level & lstart_, const nuclear_level & lstop_);

    virtual void tree_dump(std::ostream & out_ = std::clog,
                           const std::string & title_  = "",
                           const std::string & indent_ = "",
                           bool inherit_ = false) const;

    int fire_event_kinematics(mygsl::rng & prng_,
                              double & ka_,
                              double & thetaa_,
                              double & phia_,
                              double & kr_,
                              double & thetar,
                              double & phir_) const;

    virtual void initialize(const datatools::properties & config_);

    virtual void reset();

    virtual int fill(mygsl::rng & prng_, genbb::primary_event & event_);

  protected:

    void _set_default();

    void _init();

  private:

    double _q_alpha_;         
    int    _A_;               
    int    _Z_;               
    double _mass_parent_;     
    double _energy_parent_;   
    double _mass_daughter_;   
    double _energy_daughter_; 
    bool   _alpha_generated_; 
    bool   _daughter_generated_; 

    // Working data:
    double _M_;               
    double _m1_;              
    double _m2_;              
    two_body_decay _tbd_;     

    GENBB_BDD_REGISTRATION_INTERFACE(alpha_decay)

  };

} // end of namespace genbb

/***************
 * OCD support *
 ***************/
#include <datatools/ocd_macros.h>
DOCD_CLASS_DECLARATION(genbb::alpha_decay)

#endif // GENBB_HELP_ALPHA_DECAY_H

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