NucleiTree.hxx

Go to the documentation of this file.

/*
  This file is part of MURE,
  Copyright (C) 2007-2021 MURE developers.
 
  MURE is free software: you can redistribute it and/or modify
  it under the terms of the GNU Lesser General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  MURE 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 Lesser General Public License for more details.
 
  You should have received a copy of the GNU Lesser General Public License
  along with MURE.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef _NucleiTree_
#define _NucleiTree_
 
#include <iostream>
#include <fstream>
#include <vector>
#include <map>
#include <string>
 
#include "StringLine.hxx"
#include "FPDistribution.hxx"
 
class FPDistribution;
class ZAI;
 
using namespace std;
 
 
class NucleiTree
{
  public:
 
    NucleiTree(bool WithReactions = true, bool WithSpontaneousFis = true);
    ~NucleiTree(); 
 
 
    void ExtractZAI(ZAI *&zai, double MinHalfLife = 0, int MaxRecursionDepth = 10000, bool FromDecay = true);
    void ExtractWantedZAI(ZAI *zai,  bool FromDecay = true);
    bool IsZAIInTree(ZAI *zai);
    bool IsZAIInTree(int Z, int A, int I);
    ZAI *GetZAIInTree(ZAI *zai);
    ZAI *GetZAIInTree(int Z, int A, int I);
    ZAI *ProvideZAI(int Z, int A, int I); 
    void ClearUnUsed(); 
    string LongestTree(); 
    void DumpTree(); 
    void WriteAsciiTree(string FileName = "asciitree.dat"); 
    string PrintTreeStats(); 
    int TreeLength(ZAI *zai, int &nb); 
    void CheckConsistency(ZAI *&zai); 
    void BuildTrees(); 
    void SetZAIMass(ZAI *zai); 
    bool IsTreesBuilt()
    {
        return fTreesBuilt;   
    }
    void FissionProductSelection(ZAI *zai, double yield = 0); 
    bool IsWantedNuclei()
    {
        return fIsWantedNuclei;   
    }
 
  private:
    void FindReactions(); 
    void FindDecay(); 
    void SetZAIFissionEnergy(ZAI *zai); 
    void FillZAI(string line, ZAI *zai); 
    void AddFissionProductsToZAI(FPDistribution *FP, ZAI *zai, double E, bool replace = false);
 
 
    void RenormalizeBRS(ZAI *zai, double ought_to_be);
 
    string GetDecay(string DecayModes, double &BR, int &Iso, int &StartPos);
 
    string                  fDecayDataFile; 
    vector<ZAI *>            fZaiVector;       
    bool                    fDelete;        
    map<int, map<int, map<int, ZAI *> > > fChart;   
    int                     fDecayRecursionDepth;   
    int                     fReactionRecursionDepth;
    bool                    fTreesBuilt;        
    bool                    fWithReactions;     
    bool                    fWithSF;            
 
    map<int, map<int, bool> >   fWantedNuclei;  
    bool                    fIsWantedNuclei;    
 
};
 
#endif