Material.hxx

Go to the documentation of this file.

#ifndef _MATERIAL_
#define _MATERIAL_ 

#include <math.h>
#include <vector>
#include "MureHeaders.hxx"
using namespace std;



enum ProportionUnitType{ kpUNDEFINED=-1, kpMOL=1, kpMASS, kpATCM3, kpATBCM};

enum DensityUnitType{ kdUNDEFINED=-1, kdGCM3=1, kdATCM3, kdATBCM};

//________________________________________________________________________
//
//      Material

//________________________________________________________________________

class Material 
{
 public:    

    Material();


    Material(int theNumber);
    Material(const Material &m);    
    virtual ~Material();            
    

    virtual Material* Copy();

    virtual Material* Clone(double Temperature=-1);


    void    AddNucleus(int Z, int A, double Proportion=1., ProportionUnitType ProportionUnit=kpMOL, string ModeratorCatgory="");

    void    AddNucleus(int Z, int A, int I,double Proportion, ProportionUnitType ProportionUnit=kpMOL, string ModeratorCatgory="");

    void    AddNucleus(ZAI *zai, double Proportion=1., ProportionUnitType ProportionUnit=kpMOL, string ModeratorCatgory="");

    void    AddNucleus(Nucleus_ptr nucleus, double Proportion=1., ProportionUnitType ProportionUnit=kpMOL);

    vector<Nucleus_ptr> &GetComposition(){return fComposition;} 
    Nucleus_ptr GetComposition(int i){return fComposition[i];}      
    void SetComposition(vector<Nucleus_ptr> f){fComposition = f;}   
    Nucleus_ptr         GetNucleus(int Z, int A, int I=0);      
    
    double GetProportion(int i,ProportionUnitType ProportionUnit);  

    void SetProportion(int i, double proportion,ProportionUnitType ProportionUnit);
    void ResetProportion();     


    Material* Mix(Material* M2, double part,ProportionUnitType ProportionUnit=kpMOL, Material* M3=0);

    virtual vector <Material*> Duplicate(int N,bool TrueClone=true);    
    void Proportion2AtomCM3();  
    ProportionUnitType GetProportionUnit(){return fProportionUnit;} 
    double SumProportion();     
    int GetZAIIndex(ZAI *zai);                          
    int GetZAIIndex(int Z,int A, int I);                
    void UpdateZAI(bool WithReactions=true,bool WithSF=true);   
    void AddToGlobalNucleiVector();                     
    void AddZAIDaughters(ZAI *zai);                     

    
    
    double  GetDensity(DensityUnitType DensityUnit);        

    void    SetDensity(double Density, DensityUnitType DensityUnit=kdGCM3);
    DensityUnitType GetDensityUnit(){return fDensityUnit;}          
    double  GetTemperature(){return fTemperature;}  
    void    SetTemperature(double T=293.6);         
    bool    HasModerator(){return (fModeratorName.size()>0);}       
    double  MeanMolarMass();                
    bool IsFuel(){return fFuel;}            
    bool IsModerator(){return fModerator;}  
    bool IsCladding(){return fCladding;}    
    bool IsCoolant() {return fCoolant;}     
    bool IsWaterGap() {return fWaterGap;}   
    bool IsModeratorBox() {return fModeratorBox;}   
    void SetFuel(bool flag = true){fFuel=flag;}     
    void SetModerator(bool flag = true){fModerator=flag;}   
    void SetCladding(bool flag = true){fCladding=flag;} 
    void SetCoolant(bool flag = true){fCoolant=flag;}   
    void SetWaterGap(bool flag = true){fWaterGap=flag;} 
    void SetModeratorBox(bool flag = true){fModeratorBox=flag;} 


    

    string  Print();
    string  PrintAll();     

    string  PrintModerator();
    string PrintCompositionInfo(); 

    void ForbidPrint(){fPrintable=false;}
    void EnablePrint() {fPrintable=true;}   
    
    void SetPseudoMaterial();   
    bool IsPseudoMaterial() {return fPseudoMaterial;}   
    
    bool IsPrintable(){return fPrintable;}  
    

    string PrintDensityUnit(DensityUnitType DensityUnit=kdUNDEFINED);

    string PrintProportionUnit(ProportionUnitType ProportionUnit=kpUNDEFINED);
    void DumpMaterial(string filename="material.dat");  
    void DumpCompositionInfo(string name="CompInfo.dat");   
    
    void DumpMaterial(string Filename,double Mass,double Volume);


    int     GetNumber(){return fNumber;}            
    void    SetNumber(int n){fNumber =n;}           


    void    SetDefaultXSExtension(string name);
    string  GetDefaultXSExtension(){return fDefaultXSExtension;}
    void    SetEvolution();                         
    bool    IsEvolving(){return fEvolution;}        
    bool    IsPerturbative(){return fPerturbative;} 
    void    SetPerturbative(bool flag=true){fPerturbative=flag;}    

    void SortCompositionByAmount();     
    void SortCompositionByMass();       
    void FindTotalSigmaPerNucleus();        
    Material *GetFatherOfPseudoClone(){return fFatherOfPseudoClone;}    
    void SetTemperatureEvolution(bool flag=true);   

    bool IsOutCore(){return fOutCore;}      
    Cell* SetOutCore(int TheNumber,double TheVolume);
    void SetMaterialName(string name){fMaterialName=name;}  
    string GetMaterialName(){return fMaterialName;}         
    void ResetTree(){delete fTree; fTree=0;} 
    

 protected:

    Material(Material *m);
    
    int         fNumber;            
    double      fDensity;           
    DensityUnitType fDensityUnit;   
    double      fTemperature;       
    bool        fEvolution;         
    bool        fPerturbative;      
    string      fModeratorName;     
    string      fMaterialName;      
    string      fDefaultXSExtension;
    bool        fUpdated;           
    vector<Nucleus_ptr> fComposition;   
    vector<double>   fProportion;       
    ProportionUnitType fProportionUnit; 
    static NucleiTree *fTree;           
    map<int ,map<int ,map<int,int> > >  fIndex; 
    Material *fFatherOfPseudoClone; 
    bool        fTemperatureEvolution; 
    bool        fFuel   ;   
    bool        fModerator; 
    bool        fCladding;  
    bool        fCoolant;   
    bool        fWaterGap;  
    bool        fModeratorBox;  
    
    bool        fPrintable;         
    bool        fOutCore;           
    bool        fPseudoMaterial;        

};
#endif