diff --git a/openmc/model/surface_composite.py b/openmc/model/surface_composite.py
index 6271154fb8c..51019b0d8e4 100644
--- a/openmc/model/surface_composite.py
+++ b/openmc/model/surface_composite.py
@@ -451,10 +451,10 @@ def __init__(self, xmin, xmax, ymin, ymax, zmin, zmax, **kwargs):
         self.zmax = openmc.ZPlane(z0=zmax, **kwargs)
 
     def __neg__(self):
-        return +self.xmin & -self.xmax & +self.ymin & -self.ymax & +self.zmin & -self.zmax
+        return -self.xmax & +self.xmin & -self.ymax & +self.ymin & -self.zmax & +self.zmin
 
     def __pos__(self):
-        return -self.xmin | +self.xmax | -self.ymin | +self.ymax | -self.zmin | +self.zmax
+        return +self.xmax | -self.xmin | +self.ymax | -self.ymin | +self.zmax | -self.zmin
 
 
 class XConeOneSided(CompositeSurface):
diff --git a/tests/regression_tests/torus/inputs_true.dat b/tests/regression_tests/torus/inputs_true.dat
index 0045d7e466b..0ccbda4008e 100644
--- a/tests/regression_tests/torus/inputs_true.dat
+++ b/tests/regression_tests/torus/inputs_true.dat
@@ -3,7 +3,7 @@
   <cell id="1" material="1" region="-2" universe="1" />
   <cell id="2" material="1" region="-3" universe="1" />
   <cell id="3" material="1" region="-1" universe="1" />
-  <cell id="4" material="2" region="4 -5 6 -7 8 -9 2 3 1" universe="1" />
+  <cell id="4" material="2" region="-5 4 -7 6 -9 8 2 3 1" universe="1" />
   <surface coeffs="0.0 0.0 0.0 3 1.5 1" id="1" type="z-torus" />
   <surface coeffs="6 0.0 0.0 3 1.5 1" id="2" type="x-torus" />
   <surface coeffs="6 0.0 0.0 6 1 0.75" id="3" type="y-torus" />