EvolutionControl.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 _EVOLUTIONCONTROL_HXX_
#define _EVOLUTIONCONTROL_HXX_
 
#include <fstream>
#include <iostream>
#include <sstream>
#include <vector>
#include <string>
 
#include "FuelReprocessing.hxx"
#include "MCSource.hxx"
 
class FuelReprocessing;
class MCSource;
 
using namespace std;
//________________________________________________________________________
//
//                      EvolutionControl
 
//________________________________________________________________________
 
class EvolutionControl
{
  public:
 
    EvolutionControl();
    EvolutionControl(const EvolutionControl &ev); 
 
    virtual EvolutionControl *Clone()
    {
        return new EvolutionControl(*this);   
    }
 
    virtual ~EvolutionControl(); 
 
    virtual void ControlBeforeFirstMCRun() {};  
    virtual void ControlAfterRestart(int RestartRun) {};    
 
 
    virtual void ControlAfterEachMCRun();
 
 
    virtual void ControlBeforeEachSubStep();
 
 
    virtual void ControlAfterEachEvolutionStep();
 
    virtual void ControlAfterLastMCRun() {};    
 
    void SetMCRunTime(vector < double > T);     
    vector < double >  &GetMCRunTimes()
    {
        return fMCRunTimes;   
    }
 
    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(MCSource *s) {}               
    virtual void LaunchReactivityCoeffsCalculations() {}                
    virtual void ControlKeff(double Time) {}                            
    void RunMCCriticality(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 &LastMCRunNumber);
 
 
  protected:
 
    virtual void InitVector();  
    vector < double > fMCRunTimes;  
 
    bool    fConstantPower;     
    double  fWantedKeff;        
    double fWantedHNProportion; 
    vector < double >   fMCTime;    
 
    vector < double >   fMCKeff;        
    vector < double >   fMCDeltaKeff;   
    double fOldKeff;                
    bool fLockToInitialKeff;        
 
    vector < double >   fKeffWarningTime;   
    vector < double >   fKeffWarningGap;    
    vector < double >   fKeffWarning3sigma; 
    vector < int >  fKeffWarningMCNum;  
 
    double  fFitKeffSlope;      
    double  fFitKeffIntersept;  
 
    vector < double >   fMCAbs;         
    vector < double >   fMCDeltaAbs;    
    double  fFitAbsSlope;           
    double  fFitAbsIntercept;       
 
    vector < double >   fMCFiss;        
    vector < double >   fMCDeltaFiss;   
    double  fFitFissSlope;          
    double  fFitFissIntercept;      
 
    vector < double >   fMCNuFiss;      
    vector < double >   fMCDeltaNuFiss; 
    double  fFitNuFissSlope;        
    double  fFitNuFissIntercept;    
 
 
    double fSumOfFission ;      
    double fSumOfCapture ;      
    double fSumOfN2N ;          
    double fTNF;                
 
    FuelReprocessing  *fFuelReprocessing;   
 
};
 
#endif