ReactionList.hxx

Go to the documentation of this file.

#ifndef _ReactionList_
#define _ReactionList_ 

#include <math.h>
#include <map>
#include <string>
#include <fstream>

#include "MureHeaders.hxx"

using namespace std;


enum ReactionId {
            kTotal=1,
            kN_2Nd=11,
            kN_2N=16,kN_3N,
            k_Fission=18,kN_Fission,kN_NFission,kN_2NFission,
            kN_Nalpha=22,kN_N3alpha,kN_2Nalpha,kN_3Nalpha,
            kN_Np=28,
            kN_N2alpha=29,kN_2N2alpha,
            kN_Nd=32,
            kN_Nt=33,kN_NHe3,kN_Nd2alpha,kN_Nt2alpha,
            kN_4N=37,
            k_3NFission=38,
            kN_2Np=41,kN_3Np,
            kN_N2p=44,kN_Npalpha, 
            kN_gamma=102,
            kN_p=103,kN_d,kN_t,kN_He3,kN_alpha,kN_2alpha,kN_3alpha,
            kN_2p=111,kN_palpha,kN_t2alpha,
            kN_d2alpha,kN_pd,kN_pt,kN_dalpha
         };

enum DecayId {kGamma=0,kBeta_minus=1,kBeta_plus,kIsomericTransition,
              kAlpha,kNeutron,kSpontaneousFission,kProton,kUnknown=10};

class ZAI;


class ReactionList
{
 public:

    ReactionList(ZAI *zai);
    ~ReactionList(); 


    void FindReactions();
    void SetReactionListDir(string str="ReactionList"); 
    void SetReactionThreshold(double sigma=1e-2){fReactionThreshold=sigma;} 
    void InitAll(bool state); 
    void InitJon();         
    void InitFissionTest(); 
    void InitFrancoTest();  
    void InitU5Only();      
    void InitPTO();         


    void ReactionAutoDetection();

    void N_fission(bool state=true){fList[k_Fission]=state;} 
    void N_gamma(bool state=true){fList[kN_gamma]=state;} 
    void N_2N(bool state=true){fList[kN_2N]=state;} 
    void N_3N(bool state=true){fList[kN_3N]=state;} 
    void N_4N(bool state=true){fList[kN_4N]=state;} 
    void N_p(bool state=true){fList[kN_p]=state;} 
    void N_d(bool state=true){fList[kN_d]=state;} 
    void N_t(bool state=true){fList[kN_t]=state;} 
    void N_He3(bool state=true){fList[kN_He3]=state;} 
    void N_alpha(bool state=true){fList[kN_alpha]=state;} 
    void N_Np(bool state=true){fList[kN_Np]=state;} 
    void N_Nd(bool state=true){fList[kN_Nd]=state;} 
    void N_Nt(bool state=true){fList[kN_Nt]=state;} 
    void N_NHe3(bool state=true){fList[kN_NHe3]=state;} 
    void N_Nalpha(bool state=true){fList[kN_Nalpha]=state;} 
    void N_2Np(bool state=true){fList[kN_2Np]=state;} 
    void N_3Np(bool state=true){fList[kN_3Np]=state;} 
    void N_2Nd(bool state=true){fList[kN_2Nd]=state;} 
    void N_2Nalpha(bool state=true){fList[kN_2Nalpha]=state;} 
    void N_3Nalpha(bool state=true){fList[kN_3Nalpha]=state;} 
    void N_2p(bool state=true){fList[kN_2p]=state;} 
    void N_2alpha(bool state=true){fList[kN_2alpha]=state;} 
    void N_3alpha(bool state=true){fList[kN_3alpha]=state;} 
    void N_N2p(bool state=true){fList[kN_N2p]=state;} 
    void N_N2alpha(bool state=true){fList[kN_N2alpha]=state;} 
    void N_N3alpha(bool state=true){fList[kN_N3alpha]=state;} 
    void N_pd(bool state=true){fList[kN_pd]=state;} 
    void N_pt(bool state=true){fList[kN_pt]=state;} 
    void N_palpha(bool state=true){fList[kN_palpha]=state;} 
    void N_dalpha(bool state=true){fList[kN_dalpha]=state;} 
    void N_d2alpha(bool state=true){fList[kN_d2alpha]=state;} 
    void N_t2alpha(bool state=true){fList[kN_t2alpha]=state;} 
    void N_Npalpha(bool state=true){fList[kN_Npalpha]=state;} 
    void N_2N2alpha(bool state=true){fList[kN_2N2alpha]=state;} 
    void N_Nd2alpha(bool state=true){fList[kN_Nd2alpha]=state;} 
    void N_Nt2alpha(bool state=true){fList[kN_Nt2alpha]=state;} 
    

    bool IsReactionAllowed(int reaction_code);                  
    int  GetReactionArray(int i){return fReactionArray[i];} 

    
    void WriteReactionList(int other_isomer);   
    void InitReactionArray(); 
    string Print(); 
    string PrintReaction(int reaction_code);

 
 protected:
 
    ZAI*            fZAI; 
    map<int ,bool>  fList; 
    static int      fN; 
    static int      *fReactionArray;
    static string   *fReactionName; 
    double  fReactionThreshold; 
    ifstream    fXsdir; 
    string      fXsdirName; 
    string      fReactionListDir; 
};                  
#endif