RodEvolutionControl.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 _RODEVOLUTIONCONTROL_HXX_
#define _RODEVOLUTIONCONTROL_HXX_
 
#include <fstream>
#include <iostream>
#include <sstream>
#include <vector>
#include <string>
 
#include "EvolutionControl.hxx"
 
using namespace std;
//________________________________________________________________________
//
//  RodEvolutionControl
//  @author PTO
//  @version 1.0
//________________________________________________________________________
 
class RodEvolutionControl : public EvolutionControl
{
  public:
 
    RodEvolutionControl();
    RodEvolutionControl(const RodEvolutionControl &ev); 
 
    RodEvolutionControl *Clone() override
    {
        return new RodEvolutionControl(*this);   
    }
 
    ~RodEvolutionControl() override = default;  
 
    void ControlAfterEachMCRun() override;  
    void ControlBeforeEachSubStep() override;   
    void ControlAfterEachEvolutionStep() override;  
 
    void FindNextControlRodPosition();      
    void FindNextCoefs(double value, vector < double > &V, double &Slope, double &Intersept, int &LastMCRunNumber);
 
    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() override;
 
    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