(* ::Package:: *)

(************************************************************************)
(* This file was generated automatically by the Mathematica front end.  *)
(* It contains Initialization cells from a Notebook file, which         *)
(* typically will have the same name as this file except ending in      *)
(* ".nb" instead of ".m".                                               *)
(*                                                                      *)
(* This file is intended to be loaded into the Mathematica kernel using *)
(* the package loading commands Get or Needs.  Doing so is equivalent   *)
(* to using the Evaluate Initialization Cells menu command in the front *)
(* end.                                                                 *)
(*                                                                      *)
(* DO NOT EDIT THIS FILE.  This entire file is regenerated              *)
(* automatically each time the parent Notebook file is saved in the     *)
(* Mathematica front end.  Any changes you make to this file will be    *)
(* overwritten.                                                         *)
(************************************************************************)



(* ::Input:: *)
With[{dir=If[DirectoryName[$InputFileName]!="",DirectoryName[$InputFileName],
NotebookDirectory[]]},
If[FreeQ[$Path,dir],PrependTo[$Path,dir]]]


(* ::Input:: *)
BeginPackage["SMrunning`",{"SMPoleMatching`","DifferentialEquations`InterpolatingFunctionAnatomy`"}];


Print["SMrunning, v0.92, Thu 24 Oct 2013 17:44:45 F.Bezrukov <fedor@ms2.inr.ac.ru>"];


(* ::Input:: *)
Off[General::spell1];


(* ::Input:: *)
beta1loop::usage="Beta functions for the SM in one loop, including xi";


(* ::Input:: *)
beta2loop::usage="Beta functions for the SM in two loops.  xi beta function is wrong (but unused)!";


(* ::Input:: *)
betaSteinhauser::usage="2 loop everything + 3 loop gauge beta functions from 1201.5868";


beta3higgs::usage="";
beta3higgsprime::usage="";
beta3loop::usage="3-loop beta functions from Paper by Chetyrkin&Zoller arXiv:1205.xxxx";


(* ::Input:: *)
mP::usage="Numerical value of the (reduced) Planck mass";


icondsFromSolution::usage="";


(* ::Input:: *)
solveSMRG::usage="solveSMRG[mH,bfuncs,iconds] generates the solution of the RG equations for the given mH, beta functions. iconds is the name of the function giving initial conditions, eg. icondsPole1";


(* ::Input:: *)
solveSMRGdown::usage="solveSMRG[mH,bfuncs,iconds] generates the solution of the RG equations for the given mH, beta functions. iconds is the name of the function giving initial conditions, eg. icondsPole1";


(* ::Input:: *)
solvexiRG::usage="";


(* ::Input:: *)
lambdaMin::usage="lambdaMin[bfuncs_,iconds_,mH_?NumericQ] finds the minimal value of lambda for the RG evolution";


(* ::Input:: *)
muOfLambdaMin::usage="muOfLambdaMin[bfuncs_,iconds_,mH_?NumericQ] gives the value of log[mu] where the minimal value of lambda in RG evolution is obtained";


(* ::Input:: *)
referencemH::usage="referencemH[bfuncs_,iconds_] searches for the 'reference point' higgs mass, where \[Lambda][t]==\!\(\*SubscriptBox[\(\[Beta]\), \(\[Lambda]\)]\)[t]==0";


(* ::Input:: *)
g1::usage="gprime constant used in substitution rules";
g2::usage="g constant used in substitution rules";
g3::usage="gs constant used in substitution rules";
yt::usage="yt constant used in substitution rules";
\[Lambda]::usage="lambda constant used in substitution rules";
v2::usage="";
\[Alpha]0::usage="";
\[Alpha]1::usage="";
(* These moved to SMPoleMatching;
g10::usage="gprime constant used in substitution rules";
g20::usage="g constant used in substitution rules";
g30::usage="gs constant used in substitution rules";
yt0::usage="yt constant used in substitution rules";
lambda0::usage="lambda constant used in substitution rules";
mu0::usage="";
*)
\[Xi]::usage="xi constant used in substitution rules";
\[Xi]0::usage="";
\[Alpha]00::usage="";
\[Alpha]10::usage="";
v20::usage="";
\[Gamma]::usage="";
\[Beta]g1::usage="gprime constant used in substitution rules";
\[Beta]g2::usage="g constant used in substitution rules";
\[Beta]g3::usage="gs constant used in substitution rules";
\[Beta]yt::usage="yt constant used in substitution rules";
\[Beta]\[Lambda]::usage="lambda constant used in substitution rules";
\[Beta]\[Xi]::usage="";
\[Beta]v2::usage="";
\[Beta]\[Alpha]0::usage="";
\[Beta]\[Alpha]1::usage="";


(* ::Input:: *)
tmax::usage="";
tmaxxi::usage="";


vPHItree1loopSMnoScalarRG::usage="vPHItree1loopSMnoScalarRG[\[Phi]_,\[Mu]_,sol_]";


oneLoopPotentialMin::usage="oneLoopPotentialMin[beta_,icondsmH_,mH_?NumericQ]";


oneLoopPotentialMinMu::usage="oneLoopPotentialMinMu[beta_,icondsmH_,mH_?NumericQ]";


potentialminmH::usage="potentialminmH[bfuncs_,icondsmh_]";


(* ::Input:: *)
Begin["`Private`"];


(* ::Input:: *)
Off[General::spell1];


(* ::Input:: *)
beta1loop=With[
{b1=41/6,
b2=-19/6,
b3=-7},
{(* b\[Xi][yt_,g1_,g2_,g3_,\[Lambda]_,\[Xi]_,t_]->1/(4\[Pi])^2 (\[Xi]+1/6)(12*\[Lambda]+6yt^2-9/2*g2^2-9/6*g1^2), *)
\[Gamma]->-(1/(4\[Pi])^2)(9/4 g2^2+3/4 g1^2-3yt^2),
\[Beta]\[Lambda]->1/(4\[Pi])^2 (24\[Lambda]^2-6yt^4+3/8 (2g2^4+(g2^2+g1^2)^2)+(-9 g2^2-3g1^2+12yt^2)\[Lambda]),
\[Beta]yt->yt/(4\[Pi])^2 (-(9/4) g2^2-17/12 g1^2-8g3^2+9/2 yt^2),
\[Beta]g1->1/(4\[Pi])^2 g1^3 *b1,
\[Beta]g2->1/(4\[Pi])^2 g2^3 *b2,
\[Beta]g3->1/(4\[Pi])^2 g3^3 *b3}
];


(* ::Input:: *)
beta2loop=With[
{b1=41/6,
b2=-19/6,
b3=-7,
d1=17/6,
d2=3/2,
d3=2,
B11=199/18,
B12=9/2,
B13=44/3,
B21=3/2,
B22=35/6,
B23=12,
B31=11/6,
B32=9/2,
B33=-26},
{
\[Gamma]->-(1/(4\[Pi])^2)(9/4 g2^2+3/4 g1^2-3yt^2)-1/(4\[Pi])^4 (271/32 g2^4-9/16 g1^2 g2^2-431/96 g1^4-5/2 (9/4 g2^2+17/12 g1^2+8g3^2)yt^2+27/4 yt^4-6\[Lambda]^2),
\[Beta]\[Lambda]->1/(4\[Pi])^2 (24\[Lambda]^2-6*yt^4+3/8*(2g2^4+(g2^2+g1^2)^2)+(-9g2^2-3g1^2+12yt^2)\[Lambda])+
1/(4\[Pi])^4 (1/48*(915g2^6-289g2^4g1^2-559g2^2g1^4-379g1^6)+30yt^6-yt^4(8/3g1^2+32g3^2+3\[Lambda])+
\[Lambda](-73/8g2^4+39/4g2^2g1^2+629/24g1^4+108g2^2\[Lambda]+36g1^2\[Lambda]-312\[Lambda]^2)+yt^2(-9/4g2^4+21/2g2^2g1^2-19/4g1^4+\[Lambda](45/2g2^2+85/6g1^2+80g3^2-144\[Lambda]))),
\[Beta]yt->yt/(4\[Pi])^2*(-(9/4) g2^2-17/12 g1^2-8g3^2+9/2 yt^2)+
yt/(4\[Pi])^4*(-23/4g2^4-3/4g2^2 g1^2 +1187/216 g1^4+9g2^2g3^2+
19/9g1^2g3^2-108g3^4+(225/16g2^2+131/16g1^2+36g3^2)yt^2+6(-2yt^4-2yt^2 \[Lambda] +\[Lambda]^2)),
\[Beta]g1->1/(4\[Pi])^2 g1^3 *b1+1/(4\[Pi])^4 g1^3 (B11*g1^2+B12*g2^2+B13*g3^2 -d1*yt^2),
\[Beta]g2->1/(4\[Pi])^2 g2^3 *b2+1/(4\[Pi])^4 g2^3 (B21*g1^2+B22*g2^2+B23*g3^2 -d2*yt^2),
\[Beta]g3->1/(4\[Pi])^2 g3^3 *b3+1/(4\[Pi])^4 g3^3 (B31*g1^2+B32*g2^2+B33*g3^2 -d3*yt^2)
}
];


stbeta1=ToExpression["\\left(\\frac{\\alpha_1}{\\pi}\\right)^2(
    \\frac{1}{40} + \\frac{n_G}{3} 
    + \\frac{\\alpha_1}{\\pi} \\left(\\frac{9}{800} + \\frac{19n_G}{240}\\right)
    + \\frac{\\alpha_2}{\\pi} \\left(\\frac{9}{160} + \\frac{3n_G}{80}\\right)
    + \\frac{\\alpha_3}{\\pi} \\frac{11n_G}{60}
    + \\left(\\frac{\\alpha_1}{\\pi}\\right)^2 
    \\left(  \\frac{489}{512000} -\\frac{29 n_G}{2400} - \\frac{209 n_G^2}{8640} \\right)
    + \\frac{\\alpha_1}{\\pi}\\frac{\\alpha_2}{\\pi} 
    \\left(  \\frac{783}{51200} - \\frac{7 n_G}{6400} \\right)
    - \\frac{\\alpha_1}{\\pi}\\frac{\\alpha_3}{\\pi} \\frac{137n_G}{14400}
    + \\left(\\frac{\\alpha_2}{\\pi}\\right)^2 
    \\left( \\frac{3401}{20480} + \\frac{83n_G}{1920} - \\frac{11n_G^2}{960}  \\right) 
    - \\frac{\\alpha_2}{\\pi}\\frac{\\alpha_3}{\\pi} \\frac{n_G}{320}
    + \\left(\\frac{\\alpha_3}{\\pi}\\right)^2
    \\left(\\frac{275n_G}{576} - \\frac{121n_G^2}{2160}\\right)
    + n_t \\frac{\\alpha_t}{\\pi} (
      - \\frac{17}{160}
      - \\frac{\\alpha_1}{\\pi} \\frac{2827}{51200}
      - \\frac{\\alpha_2}{\\pi} \\frac{471}{2048}
      - \\frac{\\alpha_3}{\\pi} \\frac{29}{320}
      + \\frac{\\alpha_t}{\\pi}
      \\left( \\frac{339}{5120} + \\frac{303 n_t}{2560}\\right)
    )
    + \\frac{\\lambda}{4\\pi^2} \\left(
      \\frac{\\alpha_1}{\\pi}   \\frac{27}{3200} 
      + \\frac{\\alpha_2}{\\pi} \\frac{9}{640} 
      - \\frac{\\lambda}{4\\pi^2} \\frac{9}{320}
      \\right)
    )",TeXForm]


stbeta2=ToExpression["\\left(\\frac{\\alpha_2}{\\pi}\\right)^2 (
    -\\frac{43}{24} + \\frac{n_G}{3} 
    + \\frac{\\alpha_1}{\\pi} \\left(\\frac{3}{160} + \\frac{n_G}{80}\\right)
    + \\frac{\\alpha_2}{\\pi} \\left(-\\frac{259}{96} + \\frac{49n_G}{48}\\right)
    + \\frac{\\alpha_3}{\\pi} \\frac{n_G}{4}
    + \\left(\\frac{\\alpha_1}{\\pi}\\right)^2 
    \\left(  \\frac{163}{102400} -\\frac{7 n_G}{960} - \\frac{11 n_G^2}{2880} \\right)
    + \\frac{\\alpha_1}{\\pi}\\frac{\\alpha_2}{\\pi} 
    \\left(  \\frac{561}{10240} + \\frac{13 n_G}{1280} \\right)
    - \\frac{\\alpha_1}{\\pi}\\frac{\\alpha_3}{\\pi} \\frac{n_G}{960}
    + \\left(\\frac{\\alpha_2}{\\pi}\\right)^2 
    \\left(  -\\frac{667111}{110592} + \\frac{1603 n_G}{432} - \\frac{415
      n_G^2}{1728} \\right) 
    + \\frac{\\alpha_2}{\\pi}\\frac{\\alpha_3}{\\pi} \\frac{13 n_G}{64}
    + \\left(\\frac{\\alpha_3}{\\pi}\\right)^2
    \\left(\\frac{125n_G}{192} - \\frac{11n_G^2}{144}\\right)
    + n_t \\frac{\\alpha_t}{\\pi} (
      - \\frac{3}{32}
      - \\frac{\\alpha_1}{\\pi} \\frac{593}{10240}
      - \\frac{\\alpha_2}{\\pi} \\frac{729}{2048}
      - \\frac{\\alpha_3}{\\pi} \\frac{7}{64}
      + \\frac{\\alpha_t}{\\pi}
      \\left( \\frac{57}{1024} + \\frac{45 n_t}{512}\\right)
    )
    + \\frac{\\lambda}{4\\pi^2} \\left(
      \\frac{\\alpha_1}{\\pi}   \\frac{3}{640} 
      + \\frac{\\alpha_2}{\\pi} \\frac{3}{128} 
      - \\frac{\\lambda}{4\\pi^2} \\frac{3}{64}
      \\right)
    )",TeXForm]


stbeta3=ToExpression["  \\left(\\frac{\\alpha_3}{\\pi}\\right)^2 (
    -\\frac{11}{4} + \\frac{n_G}{3} 
    + \\frac{\\alpha_1}{\\pi} \\frac{11n_G}{480}
    + \\frac{\\alpha_2}{\\pi} \\frac{3n_G}{32}
    + \\frac{\\alpha_3}{\\pi} \\left(-\\frac{51}{8} + \\frac{19n_G}{12}\\right)
    + \\left(\\frac{\\alpha_1}{\\pi}\\right)^2 
    \\left(  -\\frac{13 n_G}{7680} - \\frac{121 n_G^2}{17280}  \\right)
    - \\frac{\\alpha_1}{\\pi}\\frac{\\alpha_2}{\\pi} \\frac{n_G}{2560}
    + \\frac{\\alpha_1}{\\pi}\\frac{\\alpha_3}{\\pi} \\frac{77 n_G}{2880}
    + \\left(\\frac{\\alpha_2}{\\pi}\\right)^2 
    \\left(  \\frac{241 n_G}{1536} - \\frac{11 n_G^2}{384}  \\right) 
    + \\frac{\\alpha_2}{\\pi}\\frac{\\alpha_3}{\\pi} \\frac{7 n_G}{64}
    + \\left(\\frac{\\alpha_3}{\\pi}\\right)^2
    \\left(-\\frac{2857}{128} + \\frac{5033n_G}{576} - \\frac{325n_G^2}{864}\\right)
    + n_t \\frac{\\alpha_t}{\\pi} (
      -  \\frac{1}{8}
      - \\frac{\\alpha_1}{\\pi} \\frac{101}{2560}
      - \\frac{\\alpha_2}{\\pi} \\frac{93}{512}
      - \\frac{\\alpha_3}{\\pi} \\frac{5}{8}
      + \\frac{\\alpha_t}{\\pi}
      \\left( \\frac{9}{128} + \\frac{21 n_t}{128}\\right)
    )
    )",TeXForm]


betaSteinhauser=Join[{\[Beta]g1->3/5 (4\[Pi]^2)/g1 stbeta1,\[Beta]g2->(4\[Pi]^2)/g2 stbeta2,\[Beta]g3->(4\[Pi]^2)/g3 stbeta3}/.{Subscript[\[Alpha], 1]->5/3 g1^2/(4\[Pi]),Subscript[\[Alpha], 2]->g2^2/(4\[Pi]),Subscript[\[Alpha], 3]->g3^2/(4\[Pi]),Subscript[\[Alpha], t]->yt^2/(4\[Pi]),Subscript[n, G]->3,Subscript[n, t]->1},{\[Beta]\[Lambda]->(\[Beta]\[Lambda]/.beta2loop),\[Beta]yt->(\[Beta]yt/.beta2loop)}]


belal3[yt_,g3_,la_]:=h*(12*la^2+6*la*yt^2-3*yt^4)+h^2*(-156*la^3+40*g3^2*la*yt^2-72*la^2*yt^2-16*g3^2*yt^4-(3*la*yt^4)/2+15*yt^6)+h^3*(3588*la^4+(1820*g3^4*la*yt^2)/3-1224*g3^2*la^2*yt^2+873*la^3*yt^2-32*g3^4*la*nf*yt^2-(626*g3^4*yt^4)/3+895*g3^2*la*yt^4+(1719*la^2*yt^4)/2+20*g3^4*nf*yt^4-38*g3^2*yt^6+(117*la*yt^6)/8-(1599*yt^8)/8+2016*la^4*z3-48*g3^4*la*yt^2*z3+1152*g3^2*la^2*yt^2*z3+32*g3^4*yt^4*z3-1296*g3^2*la*yt^4*z3+756*la^2*yt^4*z3+240*g3^2*yt^6*z3-198*la*yt^6*z3-36*yt^8*z3)


betaytl1 = (+yt^2*(+3/4 + 1/2*nc) + gs^2*(-3*cf))*h*yt; 
betaytl2 = (+lambda^2*(+3) + yt^2*lambda*(-6) + yt^4*(+3/4 - 9/4*nc) + 
     gs^2*yt^2*(+6*cf - 5/2*tr + 5/4*nc*ca) + 
     gs^4*(-3/2*cf^2 - 97/6*ca*cf + 10/3*nf*tr*cf))*h^2*yt;
betaytl3 = (+lambda^3*(-18) + yt^2*lambda^2*(+285/8 - 45/4*nc)
     + yt^4*lambda*(+63/2 + 45/2*nc) + 
     yt^6*(-345/32 + 9/4*z3 + 107/32*nc + 3/2*nc*z3 + 39/16*nc^2) + 
     gs^2*yt^2*lambda*(+6*cf) + 
     gs^2*yt^4*(-57/2*cf + 81/8*tr - 81/16*nc*ca)
     + gs^4*
      yt^2*(+471/16*cf^2 + 319/8*tr*cf + 717/16*ca*cf - 
        189/16*ca*tr - 33/4*nf*tr*cf - 2*nf*ca*cf - 27*z3*cf^2 - 
        18*z3*tr*cf - 27/2*z3*ca*cf + 189/8*nc*tr^2 + 
        189/16*nc*ca*cf) + 
     gs^6*(-129/2*cf^3 + 129*tr^2*cf + 387/4*ca*cf^2 - 
        3724/27*ca^2*cf + 46*nf*tr*cf^2 + 556/27*nf*ca*tr*cf + 
        140/27*nf^2*tr^2*cf - 48*z3*nf*tr*cf^2 + 48*z3*nf*ca*tr*cf))*h^3*yt;


beta3higgs=Append[DeleteCases[ betaSteinhauser,\[Beta]\[Lambda]->_],
\[Beta]\[Lambda]->(\[Beta]\[Lambda]/.beta2loop)+(2h^3 Coefficient[belal3[yt,g3,\[Lambda]],h^3]/.{h->1/(16\[Pi]^2),nf->6,z3->Zeta[3]})]


beta3higgsprime=Append[DeleteCases[ betaSteinhauser,\[Beta]\[Lambda]->_],
\[Beta]\[Lambda]->(\[Beta]\[Lambda]/.beta2loop)+(2h^3 Coefficient[belal3[yt,g3,\[Lambda]],h^3]/.{h->1/(16\[Pi]^2),nf->6,z3->Zeta[3]})+2 (-3)/(16\[Pi]^2) (-1384.7 (g3^2/(16\[Pi]^2))^3)yt^4 ]


beta3loop=Join[DeleteCases[DeleteCases[ betaSteinhauser,\[Beta]\[Lambda]->_],\[Beta]yt->_],
{\[Beta]\[Lambda]->(\[Beta]\[Lambda]/.beta2loop)+(2h^3 Coefficient[belal3[yt,g3,\[Lambda]],h^3]/.{h->1/(16\[Pi]^2),nf->6,z3->Zeta[3]})+2 (-3)/(16\[Pi]^2) (-1384.7 (g3^2/(16\[Pi]^2))^3)yt^4,
\[Beta]yt->(\[Beta]yt/.beta2loop)+2(betaytl3)/.{h->1/(16\[Pi]^2),nc->3,nf->6,z3->Zeta[3],ca -> 3, cf -> 4/3, tr -> 1/2, dR -> 3,gs->g3,lambda->\[Lambda]}} ]


(* ::Input:: *)
mP=2.4 10^18;
Log[mP/171.2]
tmax=50.
tmaxxi=12.


icondsFromSolution[sol_,muswitch_?NumericQ]:=With[{tswitch=Log[muswitch/mu0]/.sol},
Append[{g10->g1[tswitch],
g20->g2[tswitch],
g30->g3[tswitch],
yt0->yt[tswitch],
lambda0->\[Lambda][tswitch]}/.sol,
mu0->muswitch]]


(* ::Input:: *)
solveSMRG[bfuncs_,iconds_]:=Append[First[
NDSolve[{
g1'[t]==\[Beta]g1,
g2'[t]==\[Beta]g2,
g3'[t]==\[Beta]g3,
yt'[t]== \[Beta]yt,
\[Lambda]'[t]== \[Beta]\[Lambda],
g1[0]==g10,
g2[0]==g20,
g3[0]==g30,
yt[0]==yt0,
\[Lambda][0]==lambda0}/.Join[bfuncs/.{g1->g1[t],g2->g2[t],g3->g3[t],yt->yt[t],\[Lambda]->\[Lambda][t]},iconds],
{g1,g2,g3,yt,\[Lambda]},{t,0,tmax}]],
mu0->(mu0/.iconds)]


(* ::Input:: *)
solveSMRGdown[bfuncs_,iconds_]:=Append[First[
NDSolve[{
g1'[t]==\[Beta]g1,
g2'[t]==\[Beta]g2,
g3'[t]==\[Beta]g3,
yt'[t]== \[Beta]yt,
\[Lambda]'[t]== \[Beta]\[Lambda],
g1[0]==g10,
g2[0]==g20,
g3[0]==g30,
yt[0]==yt0,
\[Lambda][0]==lambda0}/.Join[bfuncs/.{g1->g1[t],g2->g2[t],g3->g3[t],yt->yt[t],\[Lambda]->\[Lambda][t]},iconds],
{g1,g2,g3,yt,\[Lambda]},{t,-1,+1}]],
mu0->(mu0/.iconds)]


(* ::Input:: *)
solveSMRGuptomin[bfuncs_,iconds_]:=Append[First[
NDSolve[{
g1'[t]==\[Beta]g1,
g2'[t]==\[Beta]g2,
g3'[t]==\[Beta]g3,
yt'[t]== \[Beta]yt,
\[Lambda]'[t]== \[Beta]\[Lambda],
g1[0]==g10,
g2[0]==g20,
g3[0]==g30,
yt[0]==yt0,
\[Lambda][0]==lambda0}/.Join[bfuncs/.{g1->g1[t],g2->g2[t],g3->g3[t],yt->yt[t],\[Lambda]->\[Lambda][t]},iconds],
{g1,g2,g3,yt,\[Lambda]},{t,0,tmax},
Method->{"EventLocator","Event"->\[Lambda]'[t]}]],
mu0->(mu0/.iconds)]


(* ::Input:: *)
solvexiRG[bfuncs_,iconds_]:=Append[First[
NDSolve[{
g1'[t]==\[Beta]g1,
g2'[t]==\[Beta]g2,
g3'[t]==\[Beta]g3,
yt'[t]== \[Beta]yt,
\[Lambda]'[t]== \[Beta]\[Lambda],
\[Xi]'[t]== \[Beta]\[Xi],
(* v2'[t]== \[Beta]v2, *)
\[Alpha]0'[t]== \[Beta]\[Alpha]0,
\[Alpha]1'[t]== \[Beta]\[Alpha]1,
g1[0]==g10,
g2[0]==g20,
g3[0]==g30,
yt[0]==yt0,
\[Lambda][0]==lambda0,
\[Xi][0]==\[Xi]0,
(* v2[0]==v20, *)
\[Alpha]0[0]==0,
\[Alpha]1[0]==0}/.Join[bfuncs/.{g1->g1[t],g2->g2[t],g3->g3[t],yt->yt[t],\[Lambda]->\[Lambda][t],\[Xi]->\[Xi][t],\[Alpha]0->\[Alpha]0[t],\[Alpha]1->\[Alpha]1[t]},iconds],
{g1,g2,g3,yt,\[Lambda],\[Xi],\[Alpha]0,\[Alpha]1(*,v2*)},{t,-tmaxxi,tmaxxi}]],
mu0->(mu0/.iconds)]


v1loopSMnoScalar[H_,G_,T_,W_,Z_,mu_]:=With[{\[Kappa]=16\[Pi]^2},
(1/\[Kappa] (-3 T^2 (mulog[T]-3/2)+(3W^2)/2 (mulog[W]-5/6)+(3Z^2)/4 (mulog[Z]-5/6)))]/.{mulog[m_]->Log[m/mu^2]}


v1loopSMnoScalar[\[Phi]_,\[Mu]_]:=With[{H=mH2tree[\[Phi]]+mH21loop[\[Phi],\[Mu]],G=mG2tree[\[Phi]]+mH21loop[\[Phi],\[Mu]]/3,T=1/2 yt^2 \[Phi]^2,W=1/4 g2^2 \[Phi]^2,Z=1/4 (g1^2+g2^2)\[Phi]^2},v1loopSMnoScalar[H,G,T,W,Z,\[Mu]]]


v1loopSMnoScalarRG[\[Phi]_,\[Mu]_,sol_]:=v1loopSMnoScalar[\[Phi],\[Mu]]/.{\[Lambda]->\[Lambda][tmu],yt->yt[tmu],g1->g1[tmu],g2->g2[tmu]}/.{tmu->Log[\[Mu]/mu0]}/.sol


vPHItree1loopSMnoScalarRG[\[Phi]_,\[Mu]_,sol_]:=\[Lambda]/4+v1loopSMnoScalar[\[Phi],\[Mu]]/\[Phi]^4/.{\[Lambda]->\[Lambda][tmu],yt->yt[tmu],g1->g1[tmu],g2->g2[tmu]}/.{tmu->Log[\[Mu]/mu0]}/.sol


oneLoopPotentialMin[beta_,icondsmH_,mH_?NumericQ]:=FindMinValue[vPHItree1loopSMnoScalarRG[Exp[t],Exp[t],solveSMRG[beta2loop,icondsmH[mH]]],{t,Log[mP],Log[0.1mP],Log[100mP]}]


oneLoopPotentialMinMu[beta_,icondsmH_,mH_?NumericQ]:=Exp[First[FindArgMin[vPHItree1loopSMnoScalarRG[Exp[t],Exp[t],solveSMRG[beta2loop,icondsmH[mH]]],{t,Log[mP],Log[0.1mP],Log[100mP]}]]]


potentialminmH[bfuncs_,icondsmh_]:=
mH/.FindRoot[oneLoopPotentialMin[bfuncs,icondsmh,mH],{mH,120,130}]


(* ::Input:: *)
lambdaMin[bfuncs_,iconds_List]:=With[{\[Lambda]sol=\[Lambda]/.solveSMRGuptomin[bfuncs,iconds]},
\[Lambda]sol[InterpolatingFunctionDomain[\[Lambda]sol][[1,2]]]]


(* ::Input:: *)
muOfLambdaMin[bfuncs_,icondsmh_,mH_?NumericQ]:=Module[{ic=icondsmh[mH],\[Lambda]sol},
\[Lambda]sol=\[Lambda]/.solveSMRGuptomin[bfuncs,ic];
mu0 Exp[InterpolatingFunctionDomain[\[Lambda]sol][[1,2]]]/.ic]


(* ::Input:: *)
referencemH[bfuncs_,icondsmh_]:=Module[{mh},mh/.FindRoot[lambdaMin[bfuncs,icondsmh[mh]],{mh,120,140}]]


(* ::Input:: *)
End[];


(* ::Input:: *)
On[General::spell1];


(* ::Input:: *)
EndPackage[]