ReadXSFile.hxx

Go to the documentation of this file.

#ifndef _ReadXSFile_
#define _ReadXSFile_

#include<fstream> 
#include<iostream> 
#include <string>
using namespace std;

#ifdef MCNP5
    typedef double Real; 
#else
    typedef float Real; 
#endif

enum  ENDFMT
    {
        kTOTAL=1,kELASTIC,
        kN2N=16,kN3N,
        kTOTALFISS,
        kNFISS,kNNpFISS,kN2NFISS,
        kNNpALPHA,kNNpP,kNNpD,kNNpT,
        kN3NFISS=38,
        kNNp2first=51,kNNp2second,kNNp40=90,kNNp2continuum,
        kABSORPTION=101,
        kNGAMMA,
        kNP,kND,kNT,KNHe3,kNALPHA
    };  


class ReadXSFile
{
 public:
    ReadXSFile(ifstream *inputfile,int start_record,int table_length,int record_length,int nrecord_per_entry,string FN);
    ReadXSFile(ifstream *inputfile,int start_record,int table_length,string FN);
    ReadXSFile(const ReadXSFile & XSF); 
    ~ReadXSFile(); 
    
    Real    GetInterpolXS(int MT,Real E); 
    Real    GetXS(int MT,Real E); 
    Real    GetXS(int MT,int BinE); 
    double  SumXS(int MT); 
    Real*   GetEnergy(); 
    Real    GetEnergy(int i); 
    int     GetNEnergy(){return fNEnergy;} 
    int     GetXSMinEnergyBin(int MT); 
    int     GetXSMaxEnergyBin(int MT); 
    int     FindEnergyBin(double E,int N, float *Ephi );    
    int     FindEnergyBin(Real E,int N, Real *Earray);  
    int     FindEnergyBin(Real E); 
    bool    IsXSFound(){return fXSFound;} 
    Real*   GetTotalNuBar(int &NE,Real* &E,bool &Todelete); 
    Real    GetTotalInterpolateNuBar(Real E); 
    Real*   GetPromptNuBar(int &NE,Real* &E); 
    ValErr_t SigmaPhi(int N, vector<ValErr_t> &Phi,float *Ephi, int MT); 
    ValErr_t NuSigmaFisPhi(int N, vector<ValErr_t> &Phi,float *Ephi); 

 private:
    void    ReadNXS(); 
    void    ReadJXS(); 
    void    ReadAsciiNXS(); 
    void    ReadAsciiJXS(); 
    void    AffectPointer(); 
    Real*   GetXS(int MT); 
    Real    GetTotalFissionXS(Real E); 
    Real    GetFissionXS(int binE); 
    void    XSFound(){fXSFound=true;} 
    int     GetXSMinMaxEnergyBin(int MT,int min); 
    
    ifstream *OriginalFile; 
    int     fNXS[16]; 
    int     fJXS[32]; 
    Real    *fReactionType; 
    int     fNReaction; 
    Real    *fEnergy; 
    int     fNEnergy; 
    Real    *fTotalXS; 
    Real    *fAbsorptionXS; 
    Real    *fElasticXS; 
    Real    **fXS; 
    Real    *XXS; 
    int     fNXXS; 
    Real    *fLocXS; 
    bool    fXSFound; 

};

#endif