EvolutionControl.hxx

Go to the documentation of this file.

#ifndef _EVOLUTIONCONTROL_HXX_
#define _EVOLUTIONCONTROL_HXX_

#include <fstream>
#include <iostream>
#include <sstream>
#include <vector>
#include <string>
#include "MureHeaders.hxx"
using namespace std;
//________________________________________________________________________
//
//                      EvolutionControl

//________________________________________________________________________

class EvolutionControl
{
 public:

    EvolutionControl();
    EvolutionControl(const EvolutionControl &ev); 
    
    virtual EvolutionControl* Clone(){return new EvolutionControl(*this);} 

    virtual ~EvolutionControl(){}           
    
    virtual void ControlAtEachMCNPStep();           
    virtual void ControlAtEachRKStep();             

// ADDED FOR ADS SIMULATION - BY NT
    virtual void ControlBeforeMCNPStep();           
// ADDED FOR ADS SIMULATION - BY NT

    virtual void ControlAfterEndOfRKIntegration();  

    void SetTimeCutting(vector<double>T);       
    vector<double> & GetMCNPRunTimes(){return fMCNPRunTimes;}   

    
    void SetWantedKeff(double keff=1.){fWantedKeff=keff;}                   
    void SetWantedHNProportion(double HNproportion=0.01){fWantedHNProportion=HNproportion;} 
    double GetWantedHNProportion() {return fWantedHNProportion;}            
    void SetLockToInitialKeff(bool flag=true){fLockToInitialKeff = flag;}   
    void SetConstantPower(bool flag=true){fConstantPower=flag;}             
    bool IsConstantPower(){return fConstantPower;}                          
    void SetFuelReprocessing(FuelReprocessing* processing){fFuelReprocessing=processing;}
    FuelReprocessing* GetFuelReprocessing(){return fFuelReprocessing;}


    
    virtual void CalculateCoeffsAtStep(int s){}                     
    virtual void SetRunFuelTemperatureCoeff(bool flag=true){}       
    virtual void SetRunVoidCoeff(bool flag=true) {}                 
    virtual void SetRunModeratorTemperatureCoeff(bool flag=true){}  
    virtual void SetRunPoisonReactivityWorth(bool flag=true){}      
    virtual void SetCriticalitySource(MCNPSource* s){}              
    virtual void LaunchReactivityCoeffsCalculations(){}             
    virtual void ControlKeff(double Time){}                         
    void RunMCNPCriticality(string FileName);   

    
    virtual double EstimatedKeff();                 
    virtual void FindKeffSlope();                   
    virtual void FindAbsorptionSlope();             
        
    virtual void FindFissionSlope();                
    virtual void FindNuFissionSlope();              
    virtual void ExtrapolateRates();                


    
    virtual void WriteControlParameters(){} 
    virtual void WriteCurrentRates();               


    void PrintKeffWarning();


    void FindNextCoefs(double value, vector<double> &V, double &Slope,double &Intersept,int &LastMCNPRunNumber);

    
 protected:
    
    virtual void InitVector();      
    vector<double>fMCNPRunTimes;    
    
    bool    fConstantPower;             
    double  fWantedKeff;                
    double fWantedHNProportion;         
    vector<double>  fMCNPTime;          
    
    vector<double>  fMCNPKeff;          
    vector<double>  fMCNPDeltaKeff;     

    vector<double>  fKeffWarningTime;   
    vector<double>  fKeffWarningGap;    
    vector<double>  fKeffWarning3sigma; 
    vector<int> fKeffWarningMCNPNum;    
    
    double  fFitKeffSlope;                  
    double  fFitKeffIntersept;              
    
    vector<double>  fMCNPAbs;               
    vector<double>  fMCNPDeltaAbs;          
    double  fFitAbsSlope;                   
    double  fFitAbsIntercept;               
    
    vector<double>  fMCNPFiss;              
    vector<double>  fMCNPDeltaFiss;         
    double  fFitFissSlope;                  
    double  fFitFissIntercept;              
    
    vector<double>  fMCNPNuFiss;            
    vector<double>  fMCNPDeltaNuFiss;       
    double  fFitNuFissSlope;                
    double  fFitNuFissIntercept;            
    
    double fOldKeff;                

    bool fLockToInitialKeff ;       
    double fSumOfFission ;      
    double fSumOfCapture ;      
    double fSumOfN2N ;          
    double fTNF;                
        
    FuelReprocessing    *fFuelReprocessing; 
 
};

#endif