MCNPSource.hxx

Go to the documentation of this file.

#ifndef _MCNPSource_
#define _MCNPSource_

#include <math.h>
#include <vector>
#include <string>
#include "MureHeaders.hxx"

using namespace std;


class Tube;
class Shape;
typedef Reference_ptr<Shape> Shape_ptr;


class MCNPSource
{
 public:

    MCNPSource(double NPS=10000);   
    MCNPSource(const MCNPSource &s);
    ~MCNPSource(); 
    
    void    SetNPS(double NPS=10000){fNPS=NPS;}
    double      GetNPS(){return fNPS;} 
    void    SetEnergy(double E=2.5e6); 
    void    SetEnergy(string E){fEnergy=E;} 
    void    SetPosition(double *pos=0); 
    void    SetPosition(double x=0,double y=0,double z=0); 
    double* GetPosition(){return fPosition;} 
    void    SetSurface(int num) {fSurface=abs(num);}
    void    SetSurface(Shape_ptr TheShape) {fSurfaceShape=TheShape;}
    void    SetWeight(double wgt){fWeight=wgt;} 
    void    SetNormal(double nrm){fNormal=nrm;} 
    void    SetEXT(string EXT){fEXT=EXT;} 
    void    SetEXT(double ext);
    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    AddBias(string bias); 
    
    void    SetExtern(bool state=true){fExtern=state;}  
    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    SetKcode(int ActiveCycle=100,int InactiveCycle=10,double ExpectedKeff=1.);
    bool    GetKcode(){return fKcode;}
    int     GetNActiveCycle(){return fActiveCycle;} 
    int     GetNInactiveCycle(){return fInactiveCycle;}

    void    SetKSRC(){fKSRC=true;}
    // Change by NT
    void    SetSSR(){fSSR=true;}
    // End of change by NT
    void    AddPosition(double *pos); 
    void    AddPosition(double x,double y,double z); 
    void    SetKSRCPosition(int i,double *pos){fKSRCposition[i]=pos;} 
    string  Print(); 

    void    AddTubeSource(Tube *tube, string TubeParticleEnergy="1e6");
    
    void    AddTubeSource(Tube *tube,Spectrum* spectrum) ;
    
    void UseThisEnergyDistribution(Spectrum* spectrum,int SourceDistributionNum=800);
    
    string  PrintTubeSource();
 
 private:
    double          fNPS;       
    string      fEnergy;    
    double      *fPosition; 
    int         fSurface;   
    Shape_ptr   fSurfaceShape;  
    double      fNormal;    
    double      fWeight;    
    string      fDirection; 
    string      fRadius;    
    string      fDistribution; 
    string      fEXT;       
    vector<string>  fVEC;   
    vector<string>  fAXS;   
    vector<string>  fBias;  
    bool        fExtern;    
    bool        fKcode;     
    bool        fKSRC;      
    // change by NT
    bool        fSSR;       
    // End of change by NT
    int         fActiveCycle;
    int         fInactiveCycle; 
    double      fKeff;      
    int         fParticleType;  
    string      fParticleDistribution; 
    vector<double*> fKSRCposition;
    
    vector<Tube*> fSourceTube;              
    vector<string> fTubeParticleEnergy;     
    vector<Spectrum*> fTubeParticleSpectrum;

};
#endif