GeneralTHDMRunner Class Reference

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

#include <GeneralTHDMRunner.h>

Detailed Description

An RGE running algorithm for the GeneralTHDM parameters.

Author

HEPfit Collaboration

Copyright

GNU General Public License

Renormalization group evolution of the relevant SM and GeneralTHDM parameters

Definition at line 21 of file GeneralTHDMRunner.h.

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

Public Attributes
const GeneralTHDM *	myGTHDM

Constructor & Destructor Documentation

GeneralTHDMRunner()

GeneralTHDMRunner::GeneralTHDMRunner

(

const StandardModel &

SM_i

)

GeneralTHDMRunner constructor.

Definition at line 13 of file GeneralTHDMRunner.cpp.

                                                              : myGTHDM(static_cast<const GeneralTHDM*> (&SM_i))
{}

~GeneralTHDMRunner()

GeneralTHDMRunner::~GeneralTHDMRunner ( )

virtual

Runner destructor.

Definition at line 16 of file GeneralTHDMRunner.cpp.

{};

Member Function Documentation

RGEGeneralTHDMRunner()

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

virtual

Definition at line 260 of file GeneralTHDMRunner.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 = {RGEsGTHDM, JacobianGTHDM, 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(RGEcheckGTHDM(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

myGTHDM

const GeneralTHDM* GeneralTHDMRunner::myGTHDM

Definition at line 36 of file GeneralTHDMRunner.h.

---

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

• GeneralTHDMRunner.h
• GeneralTHDMRunner.cpp