MeshTally.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 _MESHTALLY_HXX_
#define _MESHTALLY_HXX_
 
#include <vector>
#include <algorithm>
#include <iosfwd>
#include <string>
 
using namespace std;
 
 
class MeshTally
{
  public:
    MeshTally();
    virtual MeshTally *Clone() = 0; 
    virtual ~MeshTally() = default;
    virtual string Print()
    {
        return "";
    }
 
};
 
 
 
class TMeshTally: public MeshTally
{
  public:
    TMeshTally();
    TMeshTally(const TMeshTally  &mt); 
    MeshTally *Clone() override
    {
        return new TMeshTally(*this);   
    }
    ~TMeshTally() override = default;
    void AddMesh(string meshtype = "rectangular", int type = 1, string particle = "N", string obs = "flux");
    void SetCorA(string cora);
    void SetCorB(string corb);
    void SetCorC(string corc);
 
    void SetDoseParameter(int H = 10, int InterpolationMethod = 1, string Units = "sievert", double NormFactor = 1.0);
 
    string Print() override; 
 
  private:
    vector < int >      fMeshType;  
    vector < int >      fType;      
    vector < string >   fMeshParticle;  
    vector < string >   fMeshObs;       
 
    vector < string >   fMeshcora;  
    vector < string >   fMeshcorb;  
    vector < string >   fMeshcorc;  
 
    int fH;                     
    int fInterpolationMethod;   
    int fDoseUnits;             
    double fNormFactor;         
};
 
 
class FMeshTally: public MeshTally
{
  public:
    FMeshTally(string meshtype = "rectangular", string particle = "N");
    FMeshTally(const FMeshTally  &mt); 
    ~FMeshTally() override = default;
    MeshTally *Clone() override
    {
        return new FMeshTally(*this);   
    }
 
    void SetOrigin(vector < double > Origin)
    {
        fOrigin = Origin;
    }
    void SetCylinderAxis(double x = 0., double y = 0., double z = 1.);
    void SetCylinderVector(double x = 1., double y = 0., double z = 0.);
 
    void SetIMesh(vector < double > IMESH)
    {
        fIMESH = IMESH;
    }
 
    void SetIFineGrid(vector < int > IINTS);
 
    void SetJMesh(vector < double > JMESH)
    {
        fJMESH = JMESH;
    }
 
    void SetJFineGrid(vector < int > JINTS);
 
    void SetKMesh(vector < double > KMESH);
 
    void SetKFineGrid(vector < int > KINTS);
 
    void SetEMesh(vector < double > E);
    void SetEFineGrid(vector < int > EINTS);
 
    string Print() override; 
 
    int GetNumber()
    {
        return fNumber;
    }
 
  private:
    string fType;   
    string fParticle;   
    int fNumber; 
    vector < double > fOrigin;
 
    vector < double > fAxis;  
    vector < double > fVector;
 
    vector < double > fIMESH; 
    vector < double > fJMESH;
    vector < double > fKMESH;
    vector < double > fEMESH;
    vector < int > fIINTS;  
    vector < int > fJINTS;  
    vector < int > fKINTS;  
    vector < int > fEINTS;  
    bool    fIfineGrid; 
    bool    fJfineGrid; 
    bool    fKfineGrid; 
};
 
#endif