SerpentTallyData.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 SERPENTTALLYDATA_HXX
#define SERPENTTALLYDATA_HXX
#include <libValErr/ValErr.hxx>
#include <string>
#include <vector>
#include <map>
#include <utility>
#include <sstream>
#include <fstream>
#include <algorithm>
 
//#include "Mock.h"
#include "MUREGlobal.hxx"
#include "MiscFunction.hxx"
 
 
namespace Serpent
{
 
 
class EnergyGridResult
{
  public:
    const static unsigned NUMBER_OF_COLUMNS = 3; 
    EnergyGridResult(); 
    ~EnergyGridResult(); 
    void SetName(const string &name); 
    int GetTallyNumber(); 
    void SetTallyNumber(int number); 
    double GetBin(unsigned binNumber); 
    string Print(); 
    void Add(double value, unsigned columnIndex); 
 
    int GetSize(); 
 
  private:
    typedef vector < double > Column;
    const static unsigned VALUE_COLUMN_INDEX = 1; 
    int fTallyNumber; 
    string fName;
    vector < Column *> *fBins;
};
 
 
class BinResult
{
  public:
    static const unsigned NUMBER_OF_INDEXES = 10;
    BinResult(); 
    ValErr_t GetValErr(); 
    void SetValErr(double value, double error);  
    void AddIndex(unsigned index); 
    int GetCellBinIndex(); 
    int GetReactionBinIndex(); 
    int GetETBinIndex(); 
    int GetIndex(unsigned pos); 
    ~BinResult(); 
    string Print(); 
 
  private:
    static const unsigned ET_NUM_POS = 1;
    static const unsigned CELL_NUM_POS = 3;
    static const unsigned REACTION_NUM_POS = 6;
    vector < int > fBinIndices;
    ValErr_t *fData;
};
 
 
class TallyResult
{
  public:
    TallyResult(); 
    ~TallyResult();  
    ValErr_t Get(int CellNum = 0, int ReactionNum = 0, int ETNum = 0);  
    ValErr_t GetValErr(int index); 
    static TallyResult *Merge(vector < TallyResult *> *tallies); 
    int GetNumberOfBins(); 
    void SetName(const string &name); 
    string GetName()
    {
        return fName;   
    }
    void SetNumber(int number); 
    int GetNumber()
    {
        return fNumber;   
    }
    float *GetETValues() const;  
    int GetETSize() const; 
    void SetETSize(int size);  
    void AddBin(BinResult *bin); 
    void AllocateET();  
    void SetETBin(int index, double value); 
    string PrintBin(int index); 
 
  private:
    string fName; 
    int fNumber;  
    vector < BinResult *> fBins;
    float *fEV;
    int fEN;
 
};
 
typedef vector < TallyResult *> TallyDatabase;
typedef vector < EnergyGridResult *> EnergyGridTable;
 
 
class TallyData
{
  public:
    TallyData();    
    ~TallyData();   
 
    void Clear();  
    TallyDatabase *GetTallyDatabase()
    {
        return fTallyData;   
    }
    EnergyGridTable *GetEnergyGridTable()
    {
        return fEnergyTable;   
    }
    void AddTally(TallyResult *tally)
    {
        fTallyData -> push_back(tally);   
    }
    void AddET(EnergyGridResult *et)
    {
        fEnergyTable -> push_back(et);   
    }
    int GetNumberOfTallies()
    {
        return int(fTallyData -> size());   
    }
    int GetNumberOfET()
    {
        return int(fEnergyTable -> size());   
    }
    TallyResult *GetTally(int index)
    {
        return fTallyData -> at(index);   
    }
    EnergyGridResult *GetET(int index)
    {
        return fEnergyTable -> at(index);   
    }
 
  private:
    TallyDatabase *fTallyData;
    EnergyGridTable *fEnergyTable;
};
 
}
#endif // SERPENTTALLYDATA_HXX