diff --git a/NPSimulation/src/GaspardTrackerDummyShape.cc b/NPSimulation/src/GaspardTrackerDummyShape.cc
index 1867e685e4793e3754e509def0597614f2514212..f8376f23f4c8dddf8649bde345015e890ab4dcc4 100644
--- a/NPSimulation/src/GaspardTrackerDummyShape.cc
+++ b/NPSimulation/src/GaspardTrackerDummyShape.cc
@@ -561,8 +561,6 @@ void GaspardTrackerDummyShape::ConstructDetector(G4LogicalVolume* world)
          MMv = m_X1_Y128[i] - m_X1_Y1[i];
          MMv = MMv.unit();
 
-         G4ThreeVector MMscal = MMu.dot(MMv);
-
          MMw = MMu.cross(MMv);
 //         if (MMw.z() > 0) MMw = MMv.cross(MMu) ;
          MMw = MMw.unit();
diff --git a/NPSimulation/src/GaspardTrackerSquare.cc b/NPSimulation/src/GaspardTrackerSquare.cc
index 71e021a8e7654bc974534bc8036ef05129d243cb..a5ac645e3fe55cae3534e7fb7bafbf05024fa283 100644
--- a/NPSimulation/src/GaspardTrackerSquare.cc
+++ b/NPSimulation/src/GaspardTrackerSquare.cc
@@ -743,7 +743,6 @@ void GaspardTrackerSquare::ConstructDetector(G4LogicalVolume* world)
          MMv = MMv.unit()                ;
          G4cout << "Norm MMv X = " << MMv(0) << " , Y = " << MMv(1) << " , Z = " << MMv(2) << G4endl;
 
-         G4ThreeVector MMscal = MMu.dot(MMv);
          G4cout << "Norm MMu.MMv X = " << MMv(0) << " , Y = " << MMv(1) << " , Z = " << MMv(2) << G4endl;
 
          MMw = MMu.cross(MMv)                  ;
diff --git a/NPSimulation/src/GaspardTrackerTrapezoid.cc b/NPSimulation/src/GaspardTrackerTrapezoid.cc
index 1a501934aee355b17b4a2ea32127491e28cef858..0d107928803d0a2e36b4d04a9f9bc0266afa4852 100644
--- a/NPSimulation/src/GaspardTrackerTrapezoid.cc
+++ b/NPSimulation/src/GaspardTrackerTrapezoid.cc
@@ -635,9 +635,6 @@ void GaspardTrackerTrapezoid::ConstructDetector(G4LogicalVolume* world)
          MMv = MMv.unit()                ;
          G4cout << "Norm MMv: " << MMv << G4endl;
 
-         G4ThreeVector MMscal = MMu.dot(MMv);
-         G4cout << "Norm MMu.MMv: " << MMscal << G4endl;
-
          MMw = MMu.cross(MMv)                  ;
 //         if (MMw.z() > 0) MMw = MMv.cross(MMu) ;
          MMw = MMw.unit()                      ;
diff --git a/NPSimulation/src/ParisCluster.cc b/NPSimulation/src/ParisCluster.cc
index 464b611fe9489e03c7f70e49c1b6ed76d9a8da9c..044741537c6fcdda059451be41848161bc3f1fe6 100644
--- a/NPSimulation/src/ParisCluster.cc
+++ b/NPSimulation/src/ParisCluster.cc
@@ -564,8 +564,6 @@ void ParisCluster::ConstructDetector(G4LogicalVolume* world)
          MMv = m_X1_Y128[i] - m_X1_Y1[i];
          MMv = MMv.unit();
 
-         G4ThreeVector MMscal = MMu.dot(MMv);
-
          MMw = MMu.cross(MMv);
 //         if (MMw.z() > 0) MMw = MMv.cross(MMu) ;
          MMw = MMw.unit();
diff --git a/NPSimulation/src/ParisPhoswich.cc b/NPSimulation/src/ParisPhoswich.cc
index cd7a6a4f4342f650dab49819f8af4c21ee3ff4ab..32fe9b6cfe1e9c748eebebbe0e96017ad5d1f969 100644
--- a/NPSimulation/src/ParisPhoswich.cc
+++ b/NPSimulation/src/ParisPhoswich.cc
@@ -475,8 +475,6 @@ void ParisPhoswich::ConstructDetector(G4LogicalVolume* world)
          MMv = m_X1_Y128[i] - m_X1_Y1[i];
          MMv = MMv.unit();
 
-         G4ThreeVector MMscal = MMu.dot(MMv);
-
          MMw = MMu.cross(MMv);
 //         if (MMw.z() > 0) MMw = MMv.cross(MMu) ;
          MMw = MMw.unit();
diff --git a/NPSimulation/src/ShieldClParis.cc b/NPSimulation/src/ShieldClParis.cc
index 2b77daa1872f600c4d441b9d0551885fcc843bc7..f3fe6a464dddf0450840d05d4fa8e398faea9632 100644
--- a/NPSimulation/src/ShieldClParis.cc
+++ b/NPSimulation/src/ShieldClParis.cc
@@ -463,8 +463,6 @@ void ShieldClParis::ConstructDetector(G4LogicalVolume* world)
          //MMv = -0.5 * (m_X1_Y1[i] + m_X128_Y128[i] - m_X1_Y128[i] - m_X128_Y1[i]);
          MMv = MMv.unit();
 
-         G4ThreeVector MMscal = MMu.dot(MMv);
-
          MMw = MMu.cross(MMv);
 //         if (MMw.z() > 0) MMw = MMv.cross(MMu) ;
          MMw = MMw.unit();
diff --git a/NPSimulation/src/ShieldPhParis.cc b/NPSimulation/src/ShieldPhParis.cc
index 0f6b140dcf54dbd26f18d4ed95396aea0ea61cc1..8d8edc55db700584e415c66ac2189649a034cdba 100644
--- a/NPSimulation/src/ShieldPhParis.cc
+++ b/NPSimulation/src/ShieldPhParis.cc
@@ -473,8 +473,6 @@ void ShieldPhParis::ConstructDetector(G4LogicalVolume* world)
          //MMv = -0.5 * (m_X1_Y1[i] + m_X128_Y128[i] - m_X1_Y128[i] - m_X128_Y1[i]);
          MMv = MMv.unit();
 
-         G4ThreeVector MMscal = MMu.dot(MMv);
-
          MMw = MMu.cross(MMv);
 //         if (MMw.z() > 0) MMw = MMv.cross(MMu) ;
          MMw = MMw.unit();