diff --git a/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/DistanceCalculatorSaggingOff.cxx b/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/DistanceCalculatorSaggingOff.cxx
index 637f0725d065f0f9d11c5201ecee7c4b3809cfe6..4c73adce874c793f41544b717e10fb768aa1b48e 100644
--- a/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/DistanceCalculatorSaggingOff.cxx
+++ b/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/DistanceCalculatorSaggingOff.cxx
@@ -49,10 +49,6 @@ namespace LArWheelCalculator_Impl
 #else // parameterized sine
     double cos_a, sin_a;
     lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
-    // std::cout << sin_a << " " << cos_a << std::endl;
-    // double cos_a_old, sin_a_old;
-    // lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
-    // std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
 #endif
     // determination of the nearest quarter-wave number
     int nqwave = (z < 0.) ? 0 : int(z / lwc()->m_QuarterWaveLength);
@@ -165,10 +161,6 @@ namespace LArWheelCalculator_Impl
 #else // parameterized sine
     double cos_a, sin_a;
     lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
-    // std::cout << sin_a << " " << cos_a << std::endl;
-    // double cos_a_old, sin_a_old;
-    // lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
-    // std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
 #endif
 
     bool sqw = false;
@@ -260,10 +252,6 @@ namespace LArWheelCalculator_Impl
 #else // parameterized sine
     double cos_a, sin_a;
     lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
-    // std::cout << sin_a << " " << cos_a << std::endl;
-    // double cos_a_old, sin_a_old;
-    // lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
-    // std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
 #endif
 
     int nqwave;
@@ -352,10 +340,6 @@ namespace LArWheelCalculator_Impl
 #else // parameterized sine
     double cos_a, sin_a;
     lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
-    // std::cout << sin_a << " " << cos_a << std::endl;
-    // double cos_a_old, sin_a_old;
-    // lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
-    // std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
 #endif
 
     bool sqw = false;
@@ -451,10 +435,6 @@ namespace LArWheelCalculator_Impl
 #else
     double cos_a, sin_a;
     lwc()->m_vsincos_par.eval(P.y(), sin_a, cos_a);
-    // std::cout << sin_a << " " << cos_a << std::endl;
-    // double cos_a_old, sin_a_old;
-    // lwc()->parameterized_sincos(P.y(), sin_a_old, cos_a_old);
-    // std::cout << "old " << sin_a_old << " " << cos_a_old << std::endl;
 #endif
 
     // determination of the nearest quarter-wave number
diff --git a/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/sincos_poly.cxx b/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/sincos_poly.cxx
index b0df77c991241602d4fc0dd5378b6f4e641688ea..86cd40c3f6009734b06e997b14bd09b3be7da383 100644
--- a/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/sincos_poly.cxx
+++ b/DetectorDescription/GeoModel/GeoSpecialShapes/src/LArWheelCalculator_Impl/sincos_poly.cxx
@@ -190,8 +190,10 @@ void LArWheelCalculator::fill_sincos_parameterization()
     double sin_a_v, cos_a_v;
     parameterized_sincos(r, sin_a, cos_a);
     m_vsincos_par.eval(r, sin_a_v, cos_a_v);
+#if DEBUGPRINT
     std::cout << "default: " << r << " " << sin_a << " " << cos_a << std::endl;
     std::cout << "vec: " << r << " " << sin_a_v << " " << cos_a_v << std::endl;
+#endif
     double ds = fabs(scalpha.sn - sin_a);
     if(ds > dsin){
       dsin = ds;
@@ -256,4 +258,4 @@ void LArWheelCalculator::fill_sincos_parameterization()
   std::cout << "Approximation is " << timeSinCos/timePoly << " faster " << std::endl;
 #endif
 
-}
\ No newline at end of file
+}