PoisonEvolutionControl.hxx

Go to the documentation of this file.

#ifndef _POISONEVOLUTIONCONTROL_HXX_
#define _POISONEVOLUTIONCONTROL_HXX_

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

class PoisonEvolutionControl :public EvolutionControl
{
 public:

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

    virtual ~PoisonEvolutionControl(){}         
    
    virtual void ControlAtEachMCNPStep();           
    virtual void ControlAtEachRKStep();             
    
    virtual void ControlAfterEndOfRKIntegration();  
    
    virtual void FindNextPoisonProportion();        
    
    virtual void ExtrapolateRates();                
    
    virtual void WriteControlParameters();  

    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();               
    
 protected:
    
    void InitVector();                      

    vector<double>  fMCNPPoisonAbs;         
    vector<double>  fMCNPDeltaPoisonAbs;    
    double  fFitPoisonAbsSlope;             
    double  fFitPoisonAbsIntercept;         

    vector<double> fNextPoisonProportion; 
    vector<double> fInitialPoisonProportion;
    bool fConstantPoisonDecrease;           
    double fPoisonDecreasePerUnitTime;      
};

#endif