NucleiTree.hxx

Go to the documentation of this file.

#ifndef _NucleiTree_
#define _NucleiTree_ 

#include <math.h>
#include <iostream>
#include <fstream>
#include <vector>
#include <map>
#include <string>
//#include "MureHeaders.hxx"
#include "StringLine.hxx"
#include "FPDistribution.hxx"
using namespace std;


class NucleiTree
{
 public:

    NucleiTree(string DecayDataFile="../data/chart.JEF3T", string MassDataFile="../data/Mass.dat", 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);
    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 RenormalizeBRS(ZAI* zai, double ought_to_be);

    string GetDecay(string DecayModes, double &BR,int &Iso, int &StartPos);
    
    string                  fDecayDataFile; 
    vector<ZAI*>            fNuclei;        
    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