MCNPSource.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 _MCNPSource_
#define _MCNPSource_
#include <vector>
#include <string>
#include <iosfwd>
 
#include "MCSource.hxx"
#include "Spectrum.hxx"
#include "Shape.hxx"
#include "TReference.hxx"
 
class Cell;
class Spectrum;
 
using namespace std;
 
class MathTube;
class Shape;
 
typedef Reference_ptr<Shape> Shape_ptr;
 
 
class MCNPSource : public MCSource
{
  public:
 
    MCNPSource(double NPS = 10000);
    MCNPSource(const MCNPSource &s);
    MCSource *Clone() override
    {
        return new MCNPSource(*this);   
    }
    ~MCNPSource() override; 
 
 
    void UsePreviousRunSource(string FirstSource = "") override;
    using MCSource::SetEnergy;
    void    SetEnergy(string E)
    {
        fEnergy = E;   
    }
 
    void    SetEXT(string EXT)
    {
        fEXT = EXT;   
    }
    void    SetEXT(double ext);
    void  SetCEL(int cel);
    void    SetVEC(double i = 0, double j = 0, double k = 0); 
    void    SetVEC(string vec); 
    void    SetAXS(double i = 0, double j = 0, double k = 0); 
    void    SetAXS(string axs);
    void    SetX(string x)
    {
        X = x;
    };
    void    SetY(string y)
    {
        Y = y;
    };
    void    SetZ(string z)
    {
        Z = z;
    };
    void    AddBias(string bias); 
 
    void    SetDirection(string dir)
    {
        fDirection = dir;   
    }
    void    SetRadius(string rad)
    {
        fRadius = rad;   
    }
    void    SetParticleType(int particle)
    {
        fParticleType = particle;   
    }
    int     GetParticleType()
    {
        return fParticleType;   
    }
    void    SetParticleDistribution(string ParticleDistribution)
    {
        fParticleDistribution = ParticleDistribution;   
    }
    void    SetDistribution(string D)
    {
        fDistribution = D;
    }
 
    void    SetKSRC()
    {
        fKSRC = true;
    }
    void    SetSSR()
    {
        fSSR = true;   
    }
    bool  GetSSR()
    {
        return fSSR;   
    }
    void    AddPosition(double *pos); 
    void    AddPosition(double x, double y, double z); 
    void    SetKSRCPosition(int i, double *pos)
    {
        fKSRCposition[i] = pos;   
    }
    string  Print() override; 
    void    AddTubeSource(MathTube *tube, string TubeParticleEnergy = "1e6");
 
    void    AddTubeSource(MathTube *tube, Spectrum *spectrum) ;
 
    void UseThisEnergyDistribution(Spectrum *spectrum, int SourceDistributionNum = 800);
 
    string  PrintTubeSource();
    void SetSurfaceNumber(int num)
    {
        fSurfaceNumber = num;   
    }
    
  private:
    int         fParticleType;  
    int         fSurfaceNumber; 
    string      fDirection; 
    string      fRadius;    
    string      fDistribution; 
    string      fEXT;       
    string      X;
    string      Y;
    string      Z;
    vector<string>  fVEC;   
    vector<string>  fAXS;   
    vector<string>  fBias;  
    int             fCEL;
    bool        fKSRC;      
    bool        fSSR;       
    string      fParticleDistribution; 
    vector<double *> fKSRCposition; 
 
    vector<MathTube *> fSourceTube;             
    vector<string> fTubeParticleEnergy;     
    vector<Spectrum *> fTubeParticleSpectrum; 
    vector<double> fSourceTubeRelIntensities; // maybe old
    vector<Cell *> fTubeCell; // maybe old
 
};
#endif