diff --git a/framework/include/vectorpostprocessors/LineMaterialSamplerBase.h b/framework/include/vectorpostprocessors/LineMaterialSamplerBase.h
index a19583259121..82688dfd7b3c 100644
--- a/framework/include/vectorpostprocessors/LineMaterialSamplerBase.h
+++ b/framework/include/vectorpostprocessors/LineMaterialSamplerBase.h
@@ -172,6 +172,7 @@ LineMaterialSamplerBase<T>::execute()
     if (!hasBlocks(elem->subdomain_id()))
       continue;
 
+    _subproblem.setCurrentSubdomainID(elem, _tid);
     _subproblem.prepare(elem, _tid);
     _subproblem.reinitElem(elem, _tid);
 
@@ -205,4 +206,3 @@ LineMaterialSamplerBase<T>::finalize()
 {
   SamplerBase::finalize();
 }
-
diff --git a/test/tests/vectorpostprocessors/line_material_sampler/line_material_real_sampler.i b/test/tests/vectorpostprocessors/line_material_sampler/line_material_real_sampler.i
index 0d59ddc3394a..106a30bfcb37 100644
--- a/test/tests/vectorpostprocessors/line_material_sampler/line_material_real_sampler.i
+++ b/test/tests/vectorpostprocessors/line_material_sampler/line_material_real_sampler.i
@@ -1,13 +1,30 @@
 [Mesh]
-  type = GeneratedMesh
-  parallel_type = replicated # Until RayTracing.C is fixed
-  dim = 2
-  xmin = 0
-  xmax = 1
-  ymin = 0
-  ymax = 1
-  nx = 4
-  ny = 4
+  [genmesh]
+    type = GeneratedMeshGenerator
+    dim = 2
+    xmin = 0
+    xmax = 1
+    ymin = 0
+    ymax = 1
+    nx = 4
+    ny = 4
+  []
+  [mesh0]
+    type = SubdomainBoundingBoxGenerator
+    input = genmesh
+    block_id = 0
+    location = INSIDE
+    bottom_left = '0 0 0'
+    top_right = '1 0.5 0'
+  []
+  [mesh01]
+    type = SubdomainBoundingBoxGenerator
+    input = mesh0
+    block_id = 1
+    location = INSIDE
+    bottom_left = '0 0.5 0'
+    top_right = '1 1 0'
+  []
 []
 
 [Variables]
@@ -71,7 +88,7 @@
 [Materials]
   [./mat]
     type = MTMaterial
-    block = 0
+    block = '0 1'
   [../]
 []
 
diff --git a/test/tests/vectorpostprocessors/line_material_sampler/tests b/test/tests/vectorpostprocessors/line_material_sampler/tests
index f2635ea431f1..dc61ee73975e 100644
--- a/test/tests/vectorpostprocessors/line_material_sampler/tests
+++ b/test/tests/vectorpostprocessors/line_material_sampler/tests
@@ -4,9 +4,10 @@
     input = 'line_material_real_sampler.i'
     csvdiff = 'out_mat_0001.csv'
     compiler = '!INTEL'
+    mesh_mode = 'REPLICATED'
 
     requirement = 'The system shall support the ability to sample a scalar material along an arbitrary line through the mesh domain.'
     design = 'LineMaterialRealSampler.md'
-    issues = '#4462'
+    issues = '#4462 #16959'
   [../]
 []