decayer.h

Go to the documentation of this file.

#ifndef CYCLUS_SRC_DECAYER_H_
#define CYCLUS_SRC_DECAYER_H_
 
#include <map>
#include <set>
 
#include "composition.h"
#include "error.h"
#include "pyne.h"
#include "use_matrix_lib.h"
 
// Undefines isnan from pyne
#ifdef isnan
  #undef isnan
#endif
 
namespace cyclus {
 
/// A map type to represent all of the parent nuclides tracked. The key
/// for this map type is the parent's Nuc number, and the value is a pair
/// that contains the corresponding decay matrix column and decay
/// constant associated with that parent.
typedef std::map< int, std::pair<int, double> > ParentMap;
 
/// A map type to represent all of the daughter nuclides tracked. The
/// key for this map type is the decay matrix column associated with the
/// parent, and the value is a vector of pairs of all the daughters for
/// that parent. Each of the daughters are represented by a pair that
/// contains the daughter's Nuc number and its branching ratio.
typedef std::map<int, std::vector<std::pair<int, double> > > DaughtersMap;
 
typedef std::vector<int> NucList;
 
/// Decayer is DEPRECATED.  Use pyne::decayers::decay.
class Decayer {
 public:
  Decayer(const CompMap& comp);
 
  ~Decayer();
 
  /// set the composition from a CompMap
  void GetResult(CompMap& comp);
 
  /// decay the material
  /// @param secs the number of seconds to decay
  void Decay(double secs);
 
  /// the number of tracked nuclides
  int n_tracked_nuclides() {
    return nuclides_tracked_.size();
  }
 
  /// the tracked nuclide at position i
  int TrackedNuclide(int i) {
    return nuclides_tracked_.at(i);
  }
 
 private:
  /// Builds the decay matrix needed for the decay calculations from
  /// the parent and daughters map variables. The resulting matrix is
  /// stored in the static variable decayMatrix.
  static void BuildDecayMatrix();
 
  /// The CompMap's parent
  static ParentMap parent_;
 
  /// The CompMap's daughters
  static DaughtersMap daughters_;
 
  /// The decay matrix
  static Matrix decay_matrix_;
 
  /// The atomic composition map
  Vector pre_vect_;
  Vector post_vect_;
 
  /// the list of tracked nuclides
  static NucList nuclides_tracked_;
 
  /// Add the nuclide to the parent/daughter maps IFF it is not in the tracked list.
  static void AddNucToMaps(int nuc);
 
  /// Add the nuclide to our list of tracked nuclides IFF it is not in the list.
  static void AddNucToList(int nuc);
 
  /// Checks if the nuclide is tracked
  static bool IsNucTracked(int nuc);
};
 
}  // namespace cyclus
 
#endif  // CYCLUS_SRC_DECAYER_H_