PoisonEvolutionControl.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 _POISONEVOLUTIONCONTROL_HXX_
#define _POISONEVOLUTIONCONTROL_HXX_
 
#include <fstream>
#include <iostream>
#include <sstream>
#include <vector>
#include <string>
 
#include "EvolutionControl.hxx"
#include "ControlMaterial.hxx"
 
using namespace std;
//________________________________________________________________________
//
//    PoisonEvolutionControl
//________________________________________________________________________
 
class PoisonEvolutionControl : public EvolutionControl
{
  public:
 
    PoisonEvolutionControl();
    PoisonEvolutionControl(const PoisonEvolutionControl &ev); 
    PoisonEvolutionControl *Clone() override
    {
        return new PoisonEvolutionControl(*this);   
    }
 
    ~PoisonEvolutionControl() override = default;    
 
    void ControlAfterEachMCRun() override;   
    void ControlBeforeEachSubStep() override;
 
    virtual void FindNextPoisonProportion();
 
    void ExtrapolateRates() override;        
 
    void WriteControlParameters() override;  
 
    void SetConstantPoisonDecrease(bool flag = true)
    {
        fConstantPoisonDecrease = flag;   
    }
    bool GetConstantPoisonDecrease()
    {
        return fConstantPoisonDecrease;   
    }
 
    void SetPoisonDecreasePerUnitTime(double x)
    {
        fPoisonDecreasePerUnitTime = x;   
    }
    double GetPoisonDecreasePerUnitTime()
    {
        return fPoisonDecreasePerUnitTime;   
    }
 
    void AddReactions (ControlMaterial *M)
    {
        M -> AddControlReaction( - 2);   
    }
    void WriteCurrentRates() override;    
 
  protected:
 
    void    InitVector() override;  
 
    vector  < double >  fMCPoisonAbs;       
    vector  < double >  fMCDeltaPoisonAbs;  
    double  fFitPoisonAbsSlope;             
    double  fFitPoisonAbsIntercept;         
 
    vector  < double > fNextPoisonProportion;
    vector  < double > fInitialPoisonProportion;
    bool    fConstantPoisonDecrease;    
    double  fPoisonDecreasePerUnitTime; 
};
 
#endif