MureTally.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 _MURETALLY_HXX_
#define _MURETALLY_HXX_
 
#include <vector>
#include <algorithm>
#include <iosfwd>
#include <string>
 
using namespace std;
#include "Material.hxx"
#include "Reaction.hxx"
#include "MureTimeEnergyGrid.hxx"
#include "MCSource.hxx"
 
class Cell;
class Material;
class MureTallyBin;
class PinCell;
class Reaction;
 
 
class TallyMultiplicator
{
  public:
    TallyMultiplicator(Material *M, int reaction, double constant = 1);
    TallyMultiplicator(Material *M, Reaction *R, double constant = 1);
    TallyMultiplicator(int MatNumber, int reaction, double constant = 1);
 
    TallyMultiplicator(int MatNumber, Reaction *R, double constant = 1);
 
    TallyMultiplicator(const TallyMultiplicator &fm); 
    ~TallyMultiplicator(); 
 
    int GetMaterialNumber()
    {
        return fMatNumber;   
    }
    double GetConstant()
    {
        return fConstant;   
    }
    Reaction *GetReaction()
    {
        return fReaction;   
    }
 
  private:
    void Init(int MatNumber, Reaction *R, double constant = 1); 
    double              fConstant; 
    int                 fMatNumber; 
    Reaction            *fReaction;  
};
 
 
 
 
 
class MureTally
{
  public:
 
    typedef MureTimeEnergyGrid::BinRange ETBinRange;
 
    MureTally();
    MureTally(const MureTally &t); 
    virtual MureTally *Clone() = 0;
    virtual ~MureTally(); 
 
    int     GetNumber()
    {
        return fNumber;   
    }
    void    SetNumber(int num)
    {
        fNumber = num;   
    }
 
    bool    NeedVolume(); 
    virtual vector < double >   GetBinVolume(int i) = 0; 
    virtual void    SetBinVolume(int i, double V, int k = 0) = 0; 
    virtual string Print() = 0;
    virtual string PrintMultiplicator() = 0;
    void    SetComment(string str)
    {
        fComment = str;   
    }
    void    AddEnergy(int NE, double *E);
    void    AddEnergy(vector < double > E);
    void    AddEnergy(int NE, double Emin, double Emax, bool Linear = false);
 
    void    AddEnergy(double Emin, double Emax, int N_by_Decade = 10);
 
    string PrintDebugET();
 
    void    AddMultiGroupEnergy(double Emin, double Emax, double decad_mult = 1);
 
    void    AddTime(int NT, double *T);
 
    void    AddTime(int NT, double Tmin, double Tmax, bool Linear = false);
 
    void    AddTime(double Tmin, double Tmax, int N_by_Decade = 10);
 
    void Add(Cell *C, bool SetCellTally = true); 
    void Add(PinCell *C, int i, bool SetCellTally = true); 
    void Add(Shape_ptr S); 
    void Add(int universe);
    void Add(MureTallyBin *theTallyBin);
 
    void Add(TallyMultiplicator *Multiplicator)
    {
        fMultiplicator.push_back(Multiplicator);   
    }
    void AddMultiplicator(Material *M, int reaction, double constant = 1); 
    void AddMultiplicator(Material *M, Reaction *R, double constant = 1); 
    void AddMultiplicator(int M, int R, double constant); 
    vector < TallyMultiplicator *> &GetMultiplicator()
    {
        return fMultiplicator;   
    }
 
    MureTallyBin    *GetBin(int i)
    {
        return fBin[i];   
    }
    vector < MureTallyBin *>      &GetBin()
    {
        return fBin;   
    }
 
    int IsInTallyBin(Cell *C); 
    int IsInTallyBin(int num); 
    int IsInTallyBin(Shape_ptr S); 
    void SetPrintable(bool flag = true)
    {
        fIsPrintable = flag;   
    }
    bool IsPrintable()
    {
        return fIsPrintable;   
    }
    int GetType()
    {
        return fType;   
    }
    void SetType(int type)
    {
        fType = type;   
    }
    void SkipVolumeCalculation()
    {
        fSkipVolumeCalc = true;   
    }
    bool IsVolumeCalculationSkept()
    {
        return fSkipVolumeCalc;   
    }
    void SetPinCellTally(bool flag)
    {
        fIsPinCellTally = flag;   
    }
    bool IsPinCellTally()
    {
        return fIsPinCellTally;   
    }
 
  protected:
      
    bool IsMaterialUsed(int num);
    virtual MureTimeEnergyGrid *GetNewTimeEnergyGrid() = 0; 
 
    int fType;      
    int fNumber;    
    vector < TallyMultiplicator *>  fMultiplicator; 
    MureTimeEnergyGrid *fTimeEnergyGrid;        
    vector < MureTallyBin *> fBin;  
    bool fIsPrintable;  
    string fComment;    
    bool fSkipVolumeCalc;
    bool fIsPinCellTally;
 
};
 
 
 
 
 
#endif