Fix the energy scaling of the generated photons

ECP-WarpX · Jan 8, 2025 · c50b855 · c50b855
1 parent 2dcd95a
commit c50b855
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/PhotonCreationFunc.H b/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/PhotonCreationFunc.H
@@ -168,10 +168,11 @@ public:
                 auto& uz1 = soa_products_data[0].m_rdata[PIdx::uz][product_index];
 
                 // Normalize out the electron velocity and multiply by the photon momentum, E/c
+                // Also, the photon momentum is normalized by m_e
                 auto u1 = std::sqrt(ux1*ux1 + uy1*uy1 + uz1*uz1);
-                ux1 *= (p_product_data[i]/PhysConst::c)/u1;
-                uy1 *= (p_product_data[i]/PhysConst::c)/u1;
-                uz1 *= (p_product_data[i]/PhysConst::c)/u1;
+                ux1 *= (p_product_data[i]/PhysConst::c)/u1/PhysConst::m_e;
+                uy1 *= (p_product_data[i]/PhysConst::c)/u1/PhysConst::m_e;
+                uz1 *= (p_product_data[i]/PhysConst::c)/u1/PhysConst::m_e;
             }
         });