From 54bab6d1a5646161fcdb0c24a4029cede4df7a1a Mon Sep 17 00:00:00 2001
From: "Pyry Runko (Puck)" <pyjorunk@student.jyu.fi>
Date: Mon, 17 Feb 2025 12:46:45 +0200
Subject: [PATCH] photon waave function continued
---
src/ipglasma.cpp | 14 +++++++++++---
1 file changed, 11 insertions(+), 3 deletions(-)
diff --git a/src/ipglasma.cpp b/src/ipglasma.cpp
index 7f9a084..97df109 100644
--- a/src/ipglasma.cpp
+++ b/src/ipglasma.cpp
@@ -671,18 +671,26 @@ double IPGlasma::SumAmplitudes()
// return longitudinal photon polarization wave function
// source arxiv.org/pdf/hep-ph/0606272
-double IPGlasma::WFGammaqqbar(double r)
+// quark flavor f: 0=u, 1=d
+double IPGlasma::WFGammaqqbar(double r, int f)
{
//Assuming r is the transverse size of the dipole
//constants: e_f (quark charge), e, N_C,
//variables: r, Q, z, epsilon
//epsilon = sqrt(z(1-z)Q^2-m^2_f
+
+ // fine structure constant set to 1/137
+ double e = cmath.sqrt(4*cmath.pi*1/137)
+ double m = 0.16e-9 //GeV
+ double e_f;
+ if (f){ e_f = -1/3};
+ else { e_f = 2/3};
+
double z = 0.5 //GeV^0
double QQ = 10 //GeV^2
- double m_f = 0.16e-9 //GeV
double Delta = 0.5 //GeV
double l = a/100 // << a^-1 (a = lattice const)
- double eps = sqrt(z*(1-z)*QQ*m_f*m_f)
+ double eps = cmath.sqrt(z*(1-z)*QQ*m_f*m_f)
return e_f*e*cmath.sqrt(NC)*2*Q*z*(1-z)*std::cyl_bessel_k(0, eps*r)0.5/cmath.pi
}
--
GitLab