Merge branch 'fix_particle_sorting' into nodal_aux_3

1 files changed, 19 insertions, 1 deletions

diff --git a/Source/Particles/PhotonParticleContainer.cpp b/Source/Particles/PhotonParticleContainer.cpp
index 4a75ec9f3..3c70a957f 100644
--- a/Source/Particles/PhotonParticleContainer.cpp
+++ b/Source/Particles/PhotonParticleContainer.cpp
@@ -21,7 +21,25 @@ using namespace amrex;
 PhotonParticleContainer::PhotonParticleContainer (AmrCore* amr_core, int ispecies,
                                                   const std::string& name)
     : PhysicalParticleContainer(amr_core, ispecies, name)
-{}
+{
+
+    ParmParse pp(species_name);
+
+#ifdef WARPX_QED
+        //IF do_qed is enabled, find out if Breit Wheeler process is enabled
+        if(do_qed)
+            pp.query("do_qed_breit_wheeler", do_qed_breit_wheeler);
+
+        //Check for processes which do not make sense for photons
+        bool test_quantum_sync;
+        pp.query("do_qed_quantum_sync", test_quantum_sync);
+        AMREX_ALWAYS_ASSERT_WITH_MESSAGE(
+        test_quantum_sync == 0,
+        "ERROR: do_qed_quantum_sync can be 1 for species NOT listed in particles.photon_species only!");
+        //_________________________________________________________
+#endif
+
+}
 
 void PhotonParticleContainer::InitData()
 {