ReadXSFile.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 _ReadXSFile_
#define _ReadXSFile_
 
#include <fstream>               // for ifstream, std
#include <string>                // for string
#include <vector>                // for vector
#include <libValErr/ValErr.hxx>  // for ValErr_t
 
using namespace std;
 
#ifdef MCNP5
typedef double Real; 
#else
typedef float Real; 
#endif
 
enum ReactionId
{
    kTotal = 1,
    k_Elastic = 2,
    kTotalInelastic = 4,
    kN_2Nd = 11,
    kN_2N = 16, kN_3N = 17,
    k_Fission = 18, kN_Fission = 19, kN_NFission = 20, kN_2NFission = 21,
    kN_Nalpha = 22, kN_N3alpha = 23, kN_2Nalpha = 24, kN_3Nalpha = 25,
    kN_Np = 28,
    kN_N2alpha = 29, kN_2N2alpha = 30,
    kN_Nd = 32,
    kN_Nt = 33, kN_NHe3, kN_Nd2alpha, kN_Nt2alpha,
    kN_4N = 37,
    k_3NFission = 38,
    kN_2Np = 41, kN_3Np = 42,
    kN_N2p = 44, kN_Npalpha = 45,
    kN_Np2first = 51, kN_Np2second = 52,
    kN_Np40 = 90, kN_Np2continuum = 91,
    k_Absorption = 101,
    kN_gamma = 102,
    kN_p = 103, kN_d = 104, kN_t = 105, kN_He3 = 106,
    kN_alpha = 107, kN_2alpha = 108, kN_3alpha = 109,
    kN_2p = 111, kN_palpha = 112, kN_t2alpha = 113,
    kN_d2alpha = 114, kN_pd = 115, kN_pt = 116, kN_dalpha = 117
};
 
 
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, double emin = 0, double emax = 1e40); 
    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); 
    int   GetStartEndEnergyBin(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