Add: contribution of dFs to expansion

Expansion checks work.
lcarpino · Nov 27, 2019 · 0da1c23 · 0da1c23
1 parent 8512c1c
commit 0da1c23
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 3 deletions.
diff --git a/ARES/Driver/Expansion.wl b/ARES/Driver/Expansion.wl
@@ -58,13 +58,16 @@ BeginPackage["ARES`Driver`Expansion`"]
           muR, xmuR, logX0, logXV,
           OrderExpansion, Ltilde,
           ExpansionOpts, alphaSOpts,
-          res0, res1, res2, res
+          res0, res1, res2, res,
+          xq, xqb
         },
 
         (* set up the event *)
         event = BuildEvent[eventConfig];
         legs = event["legs"];
         dipoles = event["dipoles"];
+        xq  = eventConfig[[1]];
+        xqb = eventConfig[[2]];
 
         (* set up the observable parametrisation for this event geometry *)
         obsSC = obs["SCParametrisation"][event];

diff --git a/ARES/Expansion/AdditiveInterface.wl b/ARES/Expansion/AdditiveInterface.wl
@@ -160,12 +160,16 @@ BeginPackage["ARES`Expansion`AdditiveInterface`"]
     Fs21ab[RpNLL11_?NumericQ, dipole_?AssociationQ, obsSC_?AssociationQ] :=
       Module[
         {
-          a, res
+          a, xa, xb, I, res
         },
 
         a = obsSC["ktpow"];
+        xa = dipole["legs"][[1]]["x"];
+        xb = dipole["legs"][[2]]["x"];
 
-        res = dipole["col"] (8 Pi beta0)/a^2 I[dipole, obsSC]
+        I = - Log[xa + xb - 1] Pi^2/6 RpNLL11;
+
+        res = dipole["col"] 2/a I
       ]
 
     Fcorrel21l[RpNLL11_?NumericQ, leg_?AssociationQ, obsSC_?AssociationQ, I_] :=