MCNPTally.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 _TALLY_HXX_
#define _TALLY_HXX_
 
#include <vector>
#include <algorithm>
#include <iosfwd>
#include <string>
 
using namespace std;
 
#include "MureTally.hxx"
#include "Shape.hxx"
#include "Reaction.hxx"
#include "Cell.hxx"
#include "MureGroupBin.hxx"
#include "MureSimpleBin.hxx"
#include "MureLatticeBin.hxx"
#include "MCNPTimeEnergyGrid.hxx"
#include "MCSource.hxx"
#include "MureTimeEnergyGrid.hxx"
 
class Cell;
class MureTallyBin;
class PinCell;
 
namespace MCNP
{
 
enum TallyType {kSurfaceCurrent = 1, kSurfaceFlux, kCellFlux = 4,
                kEnergyDeposition = 6, kFissionEnergyDeposition
               };
 
 
 
class SimpleBin : public MureSimpleBin
{
  public:
    SimpleBin(): MureSimpleBin() {} 
    SimpleBin(const SimpleBin &b) = default; 
    SimpleBin *Clone() override
    {
        return new SimpleBin(*this);   
    }
    SimpleBin(Cell *C): MureSimpleBin(C) {} 
    SimpleBin(PinCell *C, int layeridx): MureSimpleBin(C, layeridx) {} 
    SimpleBin(Shape_ptr S): MureSimpleBin(S) {} 
    SimpleBin(int universe): MureSimpleBin(universe) {} 
 
 
    string  Print(int BeforeRparenthesis = 0) override; 
};
 
 
class GroupBin : public MureGroupBin
{
  public:
    GroupBin(): MureGroupBin() {} 
    GroupBin(const GroupBin &b) = default; 
    GroupBin *Clone() override
    {
        return new GroupBin(*this);   
    }
    GroupBin(Cell *C): MureGroupBin(C) {} 
    GroupBin(PinCell *C, int layeridx): MureGroupBin(C, layeridx) {} 
    GroupBin(Shape_ptr S): MureGroupBin(S) {} 
    GroupBin(MureTallyBin *theTallyBin): MureGroupBin(theTallyBin) {} 
    GroupBin(int universe): MureGroupBin(universe) {} 
 
    string Print(int BeforeRparenthesis = 0) override;
};
 
 
class LatticeBin: public MureLatticeBin
{
  public:
    LatticeBin(const LatticeBin &b) = default; 
    LatticeBin *Clone() override
    {
        return new LatticeBin(*this);   
    }
    LatticeBin(Cell *C, string pos = ""): MureLatticeBin(C, pos) {} 
    LatticeBin(PinCell *C, int layeridx, string pos = ""): MureLatticeBin(C, layeridx, pos) {} 
    LatticeBin(MureTallyBin *theTallyBin): MureLatticeBin(theTallyBin) {}; 
    LatticeBin(int universe): MureLatticeBin(universe) {} 
 
    string Print(int BeforeRparenthesis = 0) override;
 
};
 
 
 
 
 
class Tally: public MureTally
{
  public:
    Tally(int type = kCellFlux, const char *particle = "N");
    Tally(int type, int particle);
    Tally(const Tally &t); 
    Tally *Clone() override
    {
        return new Tally(*this);
    }
    ~Tally() override = default; 
 
    void AddFluenceToDoseConversion(bool ELinear = false, bool CFLinear = false);
 
    string Print() override; 
    string PrintMultiplicator() override; 
 
    void SetTallySpecialTreatment(string Treatment)
    {
        fTallySpecialTreatment = Treatment;
    }
 
    vector < double > GetBinVolume(int i) override; 
    double GetBinSurface(int i); 
    void SetBinSurface(int i, double S); 
    void SetBinVolume(int i, double V, int k = 0) override; 
 
    string  GetParticle()
    {
        return fParticle;   
    }
    void SetParticleN()
    {
        fParticle = "N";   
    }
    void SetParticleP()
    {
        fParticle = "P";   
    }
    void SetParticleNP()
    {
        fParticle = "N, P";   
    }
    void SetParticleNE()
    {
        fParticle = "N, E";   
    }
    void SetParticleNPE()
    {
        fParticle = "N, P, E";   
    }
    void SetParentheses(bool flag = true)
    {
        fParentheses = flag;   
    }
 
    vector < double > GetEnergyBinningVector()
    {
        return fTimeEnergyGrid -> GetEnergyBins();   
    }
 
  protected:
    TimeEnergyGrid *GetNewTimeEnergyGrid() override; 
 
  private:
    string  fParticle;      
    bool    fParentheses;   
    bool    fStar;  
    vector  < double >  fEnergy;    
    vector  < double >  fConversionFactor;
    bool    fELinearFTD;    
    bool    fCFLinearFTD;   
    string  fTallySpecialTreatment;
};
 
 
}
 
#endif