Used correct x for arbuzov.f3LL

michel-decay/f-and-g/f-and-g.py

@@ -400,10 +400,11 @@ def collintegrand(Ee, o):
 
     xe = 2*Ee/Mmu
     logx = np.log(1. - xe + (Mel/Mmu)**2)
+    xa = 2*Ee*Mmu / (Mel**2 + Mmu**2)
     if o == 0:
-        return (2./Mmu) * pref3loop * (arbuzov.f3LL(xe) - 8 * logx**3)
+        return (2./Mmu) * pref3loop * (arbuzov.f3LL(xa) - 8 * logx**3)
     else:
-        return (2./Mmu) * pref3loop * (arbuzov.g3LL(xe) + 8 * logx**3)
+        return (2./Mmu) * pref3loop * (arbuzov.g3LL(xa) + 8 * logx**3)
 
 
 f3LLsub = dointegral(lambda e: collintegrand(e, 0))

[ulrich_y]

This now agrees with @gurgone_a results michel_nnloll.txt

andreaLO = np.loadtxt('muonLO.txt')
andreaLO[20:3707,1]/flo[:,1]
andreaLO[20:3707,2]/glo[:,1]

andreaNLO = np.loadtxt('muonNLO.txt')
andreaNLO[20:3707,1]/fnlo[:,1]
andreaNLO[20:3707,2]/gnlo[:,1]

andreaNNLO = np.loadtxt('muonNNLO.txt')
andreaNNLO[20:3707,1]/fnnlo[:,1]
andreaNNLO[20:3707,2]/gnnlo[:,1]

andreaYFS = np.loadtxt('muonYFS.txt')
andreaYFS[20:3707,1]/fLLYFS[:,1]
andreaYFS[20:3707,2]/gLLYFS[:,1]

andreacoll = np.loadtxt('muonLL3.txt')
andreacoll[20:3707,1]/f3LLsub[:,1]
andreacoll[20:3707,2]/g3LLsub[:,1]

andreafull = np.loadtxt('muonTOT.txt')
andreafull[20:3707,1]/FFULL[:,1]
andreafull[20:3707,2]/GFULL[:,1]