RodEvolutionControl.hxx

Go to the documentation of this file.

#ifndef _RODEVOLUTIONCONTROL_HXX_
#define _RODEVOLUTIONCONTROL_HXX_

#include <fstream>
#include <iostream>
#include <sstream>
#include <vector>
#include <string>
#include "MureHeaders.hxx"
using namespace std;
//________________________________________________________________________
//
//  RodEvolutionControl
//  @author PTO
//  @version 1.0
//________________________________________________________________________

class RodEvolutionControl : public EvolutionControl
{
 public:

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

    virtual ~RodEvolutionControl(){}            
    
    virtual void ControlAtEachMCNPStep();           
    virtual void ControlAtEachRKStep();             
    
    virtual void ControlAfterEndOfRKIntegration();  
    
    void FindNextControlRodPosition();      
    void FindNextCoefs(double value, vector<double> &V, double &Slope,double &Intersept,int &LastMCNPRunNumber);

    void FindControlRodWorthParameters(double H, double keff_x0,double keff_xH,double keff_xH4);

    double RodDifferentialWorth(double x);

    double RodWorth(double x);

    double IntegrateRodWorth(double a,double b);
    bool Dichotom(double DeltaRho, double &FinalPos);   
    void FindRodWorthWeigthSlope(double CurrentWorth);  
    
 protected:
    void InitVector();
    
    double  fRodH;                  
    double  fRho0;                  
    double  fRho_rod_cte;           
    double  fRho_rod_sin;           
    double  fRodLength;             
    int     fRodNormalMove;         
    int     fRodAxis;               
    double  fOldDeltaRho;           

    double fDLength;                
    vector<double> fGlobalWorth;    
    double  fGlobalWorthSlope;      
    double  fGlobalWorthIntercept;  
    double  fGlobalWorthWeight;     
};

#endif