8 #ifndef WARPX_schwinger_process_wrapper_h_
9 #define WARPX_schwinger_process_wrapper_h_
15 #include <picsar_qed/physics/schwinger/schwinger_pair_engine_core.hpp>
42 const amrex::ParticleReal Ex,
const amrex::ParticleReal Ey,
const amrex::ParticleReal Ez,
43 const amrex::ParticleReal Bx,
const amrex::ParticleReal By,
const amrex::ParticleReal Bz,
44 const amrex::ParticleReal PoissonToGaussianThreshold,
47 using namespace amrex;
48 namespace pxr_p = picsar::multi_physics::phys;
49 namespace pxr_sh = picsar::multi_physics::phys::schwinger;
51 const auto expectedPairNumber =
52 pxr_sh::expected_pair_number<amrex::Real, pxr_p::unit_system::SI>(
53 Ex, Ey, Ez, Bx, By, Bz, dV,
dt);
55 if (expectedPairNumber <= PoissonToGaussianThreshold) {
60 return numpairs > 0._rt ? numpairs : 0._rt;
#define AMREX_FORCE_INLINE
#define AMREX_GPU_HOST_DEVICE
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getSchwingerProductionNumber(const amrex::Real dV, const amrex::Real dt, const amrex::ParticleReal Ex, const amrex::ParticleReal Ey, const amrex::ParticleReal Ez, const amrex::ParticleReal Bx, const amrex::ParticleReal By, const amrex::ParticleReal Bz, const amrex::ParticleReal PoissonToGaussianThreshold, amrex::RandomEngine const &engine)
Definition: SchwingerProcessWrapper.H:41
Real RandomNormal(Real mean, Real stddev)
unsigned int RandomPoisson(Real lambda)
float dt
Definition: stencil.py:442