diff --git a/Examples/Tests/particle_boundary_process/analysis_absorption.py b/Examples/Tests/particle_boundary_process/analysis_absorption.py
index fdde2622684..79e8d0e4bc6 100755
--- a/Examples/Tests/particle_boundary_process/analysis_absorption.py
+++ b/Examples/Tests/particle_boundary_process/analysis_absorption.py
@@ -1,15 +1,19 @@
 #!/usr/bin/env python3
 
-import sys
-
-import yt
-
 # This test shoots a beam of electrons at cubic embedded boundary geometry
 # At time step 40, none of the particles have hit the boundary yet. At time
 # step 60, all of them should have been absorbed by the boundary. In the
 # absence of the cube, none of the particles would have had time to exit
 # the problem domain yet.
 
+import os
+import sys
+
+import yt
+
+sys.path.insert(1, "../../../../warpx/Regression/Checksum/")
+from checksumAPI import evaluate_checksum
+
 # all particles are still there
 ds40 = yt.load("diags/diag1000040")
 np40 = ds40.index.particle_headers["electrons"].num_particles
@@ -20,3 +24,9 @@
 ds60 = yt.load(filename)
 np60 = ds60.index.particle_headers["electrons"].num_particles
 assert np60 == 0
+
+# compare checksums
+evaluate_checksum(
+    test_name=os.path.split(os.getcwd())[1],
+    output_file=sys.argv[1],
+)
diff --git a/Regression/Checksum/benchmarks_json/test_3d_particle_absorption.json b/Regression/Checksum/benchmarks_json/test_3d_particle_absorption.json
new file mode 100644
index 00000000000..ce6e2fcf79b
--- /dev/null
+++ b/Regression/Checksum/benchmarks_json/test_3d_particle_absorption.json
@@ -0,0 +1,10 @@
+{
+  "lev=0": {
+    "Bx": 202106.71291347666,
+    "By": 202106.71291347663,
+    "Bz": 3371.897999274175,
+    "Ex": 38304043178806.11,
+    "Ey": 38304043178806.11,
+    "Ez": 83057027925874.84
+  }
+}