RunnerTHDMW Class Reference

An RGE running algorithm for the THDMW parameters. More...

#include <RunnerTHDMW.h>

Detailed Description

An RGE running algorithm for the THDMW parameters.

Author

HEPfit Collaboration

Copyright

GNU General Public License

Renormalization group evolution of the relevant SM and THDMW parameters

Definition at line 33 of file RunnerTHDMW.h.

Public Member Functions
virtual double	RGERunnercustodialMW (double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order, double Rpeps, double tNLOuni)
virtual double	RGERunnerMW (double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order, double Rpeps, double tNLOuni)
virtual double	RGERunnerTHDMW (double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order, double Rpeps, double tNLOuni)
	RunnerTHDMW (const StandardModel &SM_i)
	RunnerTHDMW constructor. More...
virtual	~RunnerTHDMW ()
	Runner destructor. More...

Public Attributes
const THDMW *	myTHDMW

Constructor & Destructor Documentation

RunnerTHDMW()

RunnerTHDMW::RunnerTHDMW

(

const StandardModel &

SM_i

)

RunnerTHDMW constructor.

Definition at line 13 of file RunnerTHDMW.cpp.

                                                  : myTHDMW(static_cast<const THDMW*> (&SM_i))
{}

~RunnerTHDMW()

RunnerTHDMW::~RunnerTHDMW ( )

virtual

Runner destructor.

Definition at line 16 of file RunnerTHDMW.cpp.

{};

Member Function Documentation

RGERunnercustodialMW()

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

virtual

Definition at line 927 of file RunnerTHDMW.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 = {RGEscustodialMW, JacobianTHDMW, 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(RGEcheckcustodialMW(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);
}

RGERunnerMW()

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

virtual

Definition at line 882 of file RunnerTHDMW.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 = {RGEsMW, JacobianTHDMW, 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(RGEcheckMW(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);
}

RGERunnerTHDMW()

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

virtual

Definition at line 837 of file RunnerTHDMW.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 = {RGEsTHDMW, JacobianTHDMW, 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(RGEcheckTHDMW(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

myTHDMW

const THDMW* RunnerTHDMW::myTHDMW

Definition at line 50 of file RunnerTHDMW.h.

---

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

• RunnerTHDMW.h
• RunnerTHDMW.cpp