

(* Numerics with c1susy3l.m *)

(* This file should be in the directory of c1susy3l.m
 To run this script use the command: 
 math < c1susy3l_numerics.m

*)


(* Input Values taken from ModMhMax-Scenario Run File QCC1h2mepsModMhMax3R.m
 
 Scenario input values at 1 TeV: 
 mtildet = 370.0 GeV
 At= 1500 GeV
 MA = 1000 GeV
 tanbeta = 20 

*)

Get[Directory[]<>"/c1susy3l.m"];

NumberRule = {
   muR -> 173.3(*= MtOS*),
   mt -> 158.8239618252126775`18.20091602519064(* TSIL value*),
   Mz -> 91.19999999999999,
   Mw -> 80.399,
   Cw -> Mw/Mz,
   Sw -> Sqrt[1 - (Mw/Mz)^2],
   alpha -> -0.0526868498,
   tanb -> 19.8085859,
   muSUSY -> 197.733698,
   mgl -> 860.33885,
   mst1 -> 340.0029964389096,
   mst2 -> 1044.670924691974,
   thetat -> -1.3673755590330332,
   msq -> 1041.5719769827467,
   ca -> 3,
   cf -> 4/3,
   na -> 8,
   d33 -> 5/6,
   nq -> 5,
   Cthetat -> Cos[thetat],
   Sthetat -> Sin[thetat],
   Calpha -> Cos[alpha],
   Salpha -> Sin[alpha],
   Sbeta -> Sin[ArcTan[tanb]],
   Cbeta -> Cos[ArcTan[tanb]],
   Dm1glst1 -> mgl - mst1,
   Dm2st1st2 -> mst1^2 - mst2^2,
   Dm2sqst1 -> msq^2 - mst1^2,
   Dm2glst2 -> (mgl^2 - mst2^2),
   Dm2sqgl -> (msq^2 - mgl^2),
   Dm2st1t -> mst1^2 - mt^2,
   lmmgl -> Log[muR^2/mgl^2],
   lmmst1 -> Log[muR^2/mst1^2],
   lmmst2 -> Log[muR^2/mst2^2],
   lmmt -> Log[muR^2/mt^2],
   lmmsq -> Log[muR^2/msq^2]};
NumberRule = Sort[(#[[1]] -> (#[[2]] /. NumberRule) & /@ NumberRule)];

Print["Input values:"];
Print[NumberRule];


Print["1-Loop: ",{"C1exact","C1h1","C1h2","C1h3"}];
Print[(Coefficient[#, api5] & /@ {C1exact, C1h1, C1h2, C1h3}) /. NumberRule];


Print["2-Loop: ",{"C1exact","C1h1","C1h2","C1h3"}];
Print[(Coefficient[#, api5^2] & /@ {C1exact, C1h1, C1h2, C1h3}) /. NumberRule];

