An RGE running algorithm. More...

#include <Runner.h>

Inheritance diagram for Runner:

Detailed Description

An RGE running algorithm.

Author: HEPfit Collaboration

Copyright: GNU General Public License

Renormalization group evolution of the relevant SM and THDM parameters

Definition at line 30 of file Runner.h.

Public Member Functions
double	computeThValue ()
	Empty function. More...

virtual double	RGERunner (double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order, double Rpeps, double tNLOuni)

	Runner (const StandardModel &SM_i)
	Runner constructor. More...

	~Runner ()
	Runner destructor. More...

Public Member Functions inherited from ThObservable
double	getBinMax ()
	A get method to get the maximum value of the bin. More...

double	getBinMin ()
	A get method to get the minimum value of the bin. More...

const StandardModel &	getModel ()
	A get method to get the model. More...

const std::vector< std::string >	getParametersForObservable ()
	A get method to get the parameters for the specific observable. More...

void	setBinMax (double max)
	A set method to set the maximum value of the bin. More...

void	setBinMin (double min)
	A set method to set the minimum value of the bin. More...

void	setParametersForObservable (std::vector< std::string > parametersForObservable_i)
	A set method to get the parameters for the specific observable. More...

	ThObservable (const StandardModel &SM_i)
	Constructor. More...

	ThObservable (const ThObservable &orig)
	The copy constructor. More...

virtual	~ThObservable ()
	The default destructor. More...

Public Attributes
lambda1 *	mylambda1

lambda2 *	mylambda2

lambda3 *	mylambda3

lambda4 *	mylambda4

lambda5 *	mylambda5

m11_2 *	mym11_2

m22_2 *	mym22_2

const THDM *	myTHDM

Additional Inherited Members
Protected Attributes inherited from ThObservable
double	max
	the bin maximum. More...

double	min
	The bin minimum. More...

std::vector< std::string >	parametersForObservable
	a vector of parameter namesfor the specific observable More...

const StandardModel &	SM
	A reference to an object of StandardMode class. More...

Constructor & Destructor Documentation

◆ Runner()

Runner::Runner ( const StandardModel & SM_i )

Runner constructor.

Definition at line 15 of file Runner.cpp.

                                        : ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
{
    mym11_2=new m11_2(SM_i);
    mym22_2=new m22_2(SM_i);
    mylambda1=new lambda1(SM_i);
    mylambda2=new lambda2(SM_i);
    mylambda3=new lambda3(SM_i);
    mylambda4=new lambda4(SM_i);
    mylambda5=new lambda5(SM_i);
}

◆ ~Runner()

Runner::~Runner ( )

Runner destructor.

Definition at line 26 of file Runner.cpp.

{
  delete mym11_2;
  delete mym22_2;
  delete mylambda1;
  delete mylambda2;
  delete mylambda3;
  delete mylambda4;
  delete mylambda5;
};

Member Function Documentation

◆ computeThValue()

double Runner::computeThValue ( )

virtual

Empty function.

Implements ThObservable.

Definition at line 37 of file Runner.cpp.

{
    return 0.;
}

◆ RGERunner()

double Runner::RGERunner	(	double	InitialValues[],
		unsigned long int	NumberOfRGEs,
		double	Q1,
		double	Q2,
		int	order,
		double	Rpeps,
		double	tNLOuni
	)

virtual

Definition at line 545 of file Runner.cpp.

{
    //Define which stepping function should be used
    const gsl_odeiv2_step_type * T = gsl_odeiv2_step_rk4;
 
    //Allocate space for the stepping function
    gsl_odeiv2_step * s = gsl_odeiv2_step_alloc(T, NumberOfRGEs);
 
    //Define the absolute (A) and relative (R) error on y at each step.
    //The real error will be compared to the following error estimate:
    //  A + R * |y_i|
    gsl_odeiv2_control * c = gsl_odeiv2_control_y_new(1e-6, 0.0);
 
    //Allocate space for the evolutor
    gsl_odeiv2_evolve * e = gsl_odeiv2_evolve_alloc(NumberOfRGEs);
 
    //Definition of the RGE system (the Jacobian is not necessary for the RK4 method; it's an empty function here)
    gsl_odeiv2_system RGEsystem = {RGEs, Jacobian, NumberOfRGEs, &order};
 
    //Set starting and end point as natural logarithmic scales (conversion from decadic log scale)
    double t1 = Q1*log(10.0);
    double t2 = Q2*log(10.0);
    double tNLOuni = NLOuniscale*log(10.0);
 
    //Set initial step size
    double InitialStepSize = 1e-6;
 
    //Run!
    while (t1 < t2)
    {
        int status = gsl_odeiv2_evolve_apply (e, c, s, &RGEsystem, &t1, t2, &InitialStepSize, InitialValues);
        if(status != GSL_SUCCESS) break;
 
        //intermediate checks if appropriate
        if(RGEcheck(InitialValues,t1,Rpeps,tNLOuni) != 0) break;
    }
 
    gsl_odeiv2_evolve_free (e);
    gsl_odeiv2_control_free (c);
    gsl_odeiv2_step_free (s);
 
    //Return the decadic log scale at which the evolution stopped
    return t1/log(10.0);
}

Member Data Documentation

◆ mylambda1

lambda1* Runner::mylambda1

Definition at line 53 of file Runner.h.

◆ mylambda2

lambda2* Runner::mylambda2

Definition at line 54 of file Runner.h.

◆ mylambda3

lambda3* Runner::mylambda3

Definition at line 55 of file Runner.h.

◆ mylambda4

lambda4* Runner::mylambda4

Definition at line 56 of file Runner.h.

◆ mylambda5

lambda5* Runner::mylambda5

Definition at line 57 of file Runner.h.

◆ mym11_2

m11_2* Runner::mym11_2

Definition at line 51 of file Runner.h.

◆ mym22_2

m22_2* Runner::mym22_2

Definition at line 52 of file Runner.h.

◆ myTHDM

const THDM* Runner::myTHDM

Definition at line 50 of file Runner.h.

The documentation for this class was generated from the following files: