From 101ebbd3661e1e1431d38261873a4a2b7bd530ff Mon Sep 17 00:00:00 2001
From: adrien-matta <a.matta@surrey.ac.uk>
Date: Thu, 24 Apr 2014 15:07:01 +0100
Subject: [PATCH] * Cosmetic editing of the S1

---
 Inputs/DetectorConfiguration/e530.detector    |  4 +-
 NPSimulation/AnnularS1/AnnularS1.cc           | 52 ++++---------------
 NPSimulation/AnnularS1/AnnularS1.hh           |  9 ++--
 NPSimulation/AnnularS1/AnnularS1Scorers.cc    |  1 -
 NPSimulation/AnnularS1/AnnularS1Scorers.hh    |  2 +-
 .../ComptonTelescope/ComptonTelescope.hh      |  1 +
 NPSimulation/DummyDetector/DummyDetector.hh   |  1 +
 NPSimulation/Eurogam/Eurogam.hh               |  1 +
 NPSimulation/GASPARD/GaspardTracker.hh        |  2 +
 NPSimulation/GASPARD/GaspardTrackerAnnular.hh |  1 +
 .../GASPARD/GaspardTrackerDummyShape.hh       |  1 +
 NPSimulation/GASPARD/GaspardTrackerModule.hh  |  1 +
 NPSimulation/GASPARD/GaspardTrackerSquare.hh  |  1 +
 .../GASPARD/GaspardTrackerTrapezoid.hh        |  1 +
 NPSimulation/Helios/Helios.hh                 |  1 +
 NPSimulation/Helios/HeliosDetDummyShape.hh    |  1 +
 NPSimulation/Helios/HeliosModule.hh           |  1 +
 NPSimulation/Hyde2/Hyde2Tracker.hh            |  1 +
 NPSimulation/Hyde2/Hyde2TrackerModule.hh      |  1 +
 NPSimulation/Hyde2/Hyde2TrackerSquare1.hh     |  1 +
 NPSimulation/Hyde2/Hyde2TrackerTrapezoid1.hh  |  1 +
 NPSimulation/Hyde2/Hyde2TrackerTrapezoid2.hh  |  1 +
 NPSimulation/MUST2/MUST2Scorers.hh            |  2 +
 NPSimulation/Paris/Paris.hh                   |  4 +-
 NPSimulation/Paris/ParisCluster.hh            |  2 +
 NPSimulation/Paris/ParisModule.hh             |  2 +
 NPSimulation/Paris/ParisPhoswich.hh           |  1 +
 NPSimulation/Plastic/Plastic.hh               |  1 +
 NPSimulation/SSSD/ThinSi.hh                   |  1 +
 NPSimulation/SSSD/ThinSiScorers.hh            |  3 +-
 NPSimulation/Sharc/Sharc.hh                   |  1 +
 NPSimulation/Sharc/SharcScorers.hh            |  1 +
 NPSimulation/Shield/Shield.hh                 |  1 +
 NPSimulation/Shield/ShieldClParis.hh          |  1 +
 NPSimulation/Shield/ShieldModule.hh           |  1 +
 NPSimulation/Shield/ShieldPhParis.hh          |  1 +
 NPSimulation/Tiara/Tiara.hh                   |  1 +
 NPSimulation/Tigress/Tigress.hh               |  1 +
 NPSimulation/W1/W1.hh                         |  1 +
 NPSimulation/include/Chamber.hh               |  1 +
 .../include/EventGeneratorGammaDecay.hh       |  3 +-
 .../include/EventGeneratorIsotropic.hh        |  1 +
 .../include/EventGeneratorParticleDecay.hh    |  2 +-
 .../include/EventGeneratorTwoBodyReaction.hh  |  1 +
 NPSimulation/include/GeneralScorers.hh        |  2 +-
 NPSimulation/include/MyMagneticField.hh       |  2 +
 NPSimulation/include/Particle.hh              |  5 +-
 NPSimulation/include/ParticleStack.hh         |  3 +-
 NPSimulation/include/PhysicsList.hh           |  3 +-
 .../include/PrimaryGeneratorAction.hh         |  1 +
 NPSimulation/include/ResistiveStripScorers.hh |  1 +
 NPSimulation/include/Target.hh                |  1 +
 NPSimulation/include/VDetector.hh             |  1 +
 NPSimulation/include/VEventGenerator.hh       |  1 +
 NPSimulation/src/GeneralScorers.cc            |  2 +-
 55 files changed, 77 insertions(+), 64 deletions(-)

diff --git a/Inputs/DetectorConfiguration/e530.detector b/Inputs/DetectorConfiguration/e530.detector
index eee18b077..1f3220e4f 100644
--- a/Inputs/DetectorConfiguration/e530.detector
+++ b/Inputs/DetectorConfiguration/e530.detector
@@ -19,8 +19,8 @@
 GeneralTarget
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1
 Target
-	THICKNESS= 0
-	RADIUS= 0.1
+	THICKNESS= 4
+	RADIUS= 4
 	MATERIAL= CD2
 	ANGLE= 0
 	X= 0
diff --git a/NPSimulation/AnnularS1/AnnularS1.cc b/NPSimulation/AnnularS1/AnnularS1.cc
index 69c565a5b..3020106fe 100644
--- a/NPSimulation/AnnularS1/AnnularS1.cc
+++ b/NPSimulation/AnnularS1/AnnularS1.cc
@@ -63,24 +63,21 @@ using namespace S1SCORERS;
 
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-AnnularS1::AnnularS1()
-{
+AnnularS1::AnnularS1(){
    m_Event = new TS1Data();
 }
 
 
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-AnnularS1::~AnnularS1()
-{
+AnnularS1::~AnnularS1(){
    delete m_Event;
 }
 
 
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-void AnnularS1::AddModule(G4double PosZ)
-{
+void AnnularS1::AddModule(G4double PosZ){
    m_PosZ.push_back(PosZ);
 }
 
@@ -127,42 +124,15 @@ void AnnularS1::VolumeMaker(G4int             DetecNumber,
    G4Material* Silicon = new G4Material("Si", z = 14., a, density);
 
    // Al
-//   density = 2.702 * g / cm3;
-//   a = 26.98 * g / mole;
-//   G4Material* Aluminium = new G4Material("Aluminium", z = 13., a, density);
-
-   // Iron
-//   density = 7.874 * g / cm3;
-//   a = 55.847 * g / mole;
-//   G4Material* Iron = new G4Material("Iron", z = 26., a, density);
-
-   // CsI
-   density = 4.51 * g / cm3;
-   G4Material* CsI = new G4Material("CsI", density, ncomponents = 2);
-   CsI->AddElement(Cs , natoms = 1);
-   CsI->AddElement(I  , natoms = 1);
+   density = 2.702 * g / cm3;
+   a = 26.98 * g / mole;
+   G4Material* Aluminium = new G4Material("Aluminium", z = 13., a, density);
 
    //  Vacuum
    density = 0.000000001 * mg / cm3;
    G4Material* Vacuum = new G4Material("Vacuum", density, ncomponents = 2);
    Vacuum->AddElement(N, .7);
    Vacuum->AddElement(O, .3);
-
-   //  Mylar
-   density = 1.397 * g / cm3;
-   G4Material* Myl = new G4Material("Mylar", density, ncomponents = 3);
-   Myl->AddElement(C, natoms = 10);
-   Myl->AddElement(H, natoms = 8);
-   Myl->AddElement(O, natoms = 4);
-
-   // Havar
-   G4Material* Harvar = new G4Material("Havar", 8.3*g / cm3, 5);
-   Harvar->AddElement(Co , 42);
-   Harvar->AddElement(Cr , 20);
-   Harvar->AddElement(Ni , 13);
-   Harvar->AddElement(Fe , 19);
-   Harvar->AddElement(W  ,  1);
-
    ////////////////////////////////////////////////////////////////
    ////////////// Starting Volume Definition //////////////////////
    ////////////////////////////////////////////////////////////////
@@ -177,7 +147,6 @@ void AnnularS1::VolumeMaker(G4int             DetecNumber,
                                 0*deg, 
                                 360*deg);
 
-//   G4LogicalVolume* logicAnnularS1 = new G4LogicalVolume(solidAnnularS1, Iron, Name, 0, 0, 0);
    G4LogicalVolume* logicAnnularS1 = new G4LogicalVolume(solidAnnularS1, Vacuum, Name, 0, 0, 0);
 
    new G4PVPlacement(G4Transform3D(*rotation, position),
@@ -188,7 +157,6 @@ void AnnularS1::VolumeMaker(G4int             DetecNumber,
                                      0);
 
    logicAnnularS1->SetVisAttributes(G4VisAttributes::Invisible);
-   if (m_non_sensitive_part_visiualisation) logicAnnularS1->SetVisAttributes(G4VisAttributes(G4Colour(0.90, 0.90, 0.90)));
 
 
    // Aluminium dead layers
@@ -202,8 +170,7 @@ void AnnularS1::VolumeMaker(G4int             DetecNumber,
                                       0*deg, 
                                       360*deg); 
 
-//   G4LogicalVolume* logicAluStrip = new G4LogicalVolume(solidAluStrip, Aluminium, "logicAluStrip", 0, 0, 0);
-   G4LogicalVolume* logicAluStrip = new G4LogicalVolume(solidAluStrip, Vacuum, "logicAluStrip", 0, 0, 0);
+   G4LogicalVolume* logicAluStrip = new G4LogicalVolume(solidAluStrip, Aluminium, "logicAluStrip", 0, 0, 0);
 
    new G4PVPlacement(0, positionAluStripFront, logicAluStrip, Name + "_AluStripFront", logicAnnularS1, false, 0);
    new G4PVPlacement(0, positionAluStripBack,  logicAluStrip, Name + "_AluStripBack",  logicAnnularS1, false, 0);
@@ -229,7 +196,7 @@ void AnnularS1::VolumeMaker(G4int             DetecNumber,
 
    ///Visualisation of Silicon Strip
 //   G4VisAttributes* SiliconVisAtt = new G4VisAttributes(G4Colour(0.5, 0.5, 0.5));
-   G4VisAttributes* SiliconVisAtt = new G4VisAttributes(G4Colour(0.0, 0.0, 0.9));
+   G4VisAttributes* SiliconVisAtt = new G4VisAttributes(G4Colour(0.3, 0.3, 0.3));
    logicSilicon->SetVisAttributes(SiliconVisAtt);
 }
 
@@ -242,8 +209,7 @@ void AnnularS1::VolumeMaker(G4int             DetecNumber,
 
 // Read stream at Configfile to pick-up parameters of detector (Position,...)
 // Called in DetecorConstruction::ReadDetextorConfiguration Method
-void AnnularS1::ReadConfiguration(string Path)
-{
+void AnnularS1::ReadConfiguration(string Path){
    ifstream ConfigFile;
    ConfigFile.open(Path.c_str());
    string LineBuffer, DataBuffer;
diff --git a/NPSimulation/AnnularS1/AnnularS1.hh b/NPSimulation/AnnularS1/AnnularS1.hh
index d21ca62a6..61aa963b3 100644
--- a/NPSimulation/AnnularS1/AnnularS1.hh
+++ b/NPSimulation/AnnularS1/AnnularS1.hh
@@ -25,6 +25,7 @@
 // C++ headers
 #include <vector>
 using namespace std;
+
using namespace CLHEP;
 
 // NPTool header
 #include "VDetector.hh"
@@ -114,20 +115,16 @@ namespace ANNULARS1
 {
    // Energy/Time resolutions for the different layers
    const G4double ResoFirstStage  = 0              ;// = 52keV of Resolution   //   Unit is MeV/2.35
-//   const G4double ResoFirstStage  = 0.022          ;// = 52keV of Resolution   //   Unit is MeV/2.35
    const G4double ResoSecondStage = 0.055          ;// = 130 keV of resolution //   Unit is MeV/2.35
    const G4double ResoThirdStage  = 0              ;// = 100 keV of resolution //   Unit is MeV/2.35
-//   const G4double ResoThirdStage  = 0.043          ;// = 100 kev of resolution //   Unit is MeV/2.35
    const G4double ResoTimeGpd     = 0.212765957    ;// = 500ps                 //   Unit is  ns/2.35
    const G4double ResoTimePPAC    = 0.106382979    ;// = 250ps                 //   Unit is  ns/2.35
 
    // Geometry
    const G4double Rmin   =  46/2*mm;
    const G4double Rmax   = 100/2*mm;
-   const G4double Length = 10*mm;
-
+   
    // First stage
-//   const G4double AluStripThickness = 0.00000001*micrometer;
    const G4double AluStripThickness   = 0.4*micrometer;
    const G4double FirstStageThickness = 500*micrometer;
    const G4double FirstStageRmin      = 48/2*mm;
@@ -138,6 +135,8 @@ namespace ANNULARS1
    const G4int NbThetaStrips   = 16;
    const G4int NbThetaQuadrant = 4;
 
+   const G4double Length = AluStripThickness*2+FirstStageThickness;
+
    // Starting at the front and going in direction of third stage
    const G4double AluStripFront_PosZ = Length* -0.5       + 0.5*AluStripThickness;
    const G4double Silicon_PosZ       = AluStripFront_PosZ + 0.5*AluStripThickness   + 0.5*FirstStageThickness;
diff --git a/NPSimulation/AnnularS1/AnnularS1Scorers.cc b/NPSimulation/AnnularS1/AnnularS1Scorers.cc
index 04ef13933..1ca069598 100644
--- a/NPSimulation/AnnularS1/AnnularS1Scorers.cc
+++ b/NPSimulation/AnnularS1/AnnularS1Scorers.cc
@@ -23,7 +23,6 @@
 #include "AnnularS1Scorers.hh"
 #include "G4UnitsTable.hh"
 #include "AnnularS1.hh"
-
 using namespace ANNULARS1;
 using namespace S1SCORERS;
 
diff --git a/NPSimulation/AnnularS1/AnnularS1Scorers.hh b/NPSimulation/AnnularS1/AnnularS1Scorers.hh
index 8703693a7..98f5fe161 100644
--- a/NPSimulation/AnnularS1/AnnularS1Scorers.hh
+++ b/NPSimulation/AnnularS1/AnnularS1Scorers.hh
@@ -23,7 +23,7 @@
 
 #include "G4VPrimitiveScorer.hh"
 #include "G4THitsMap.hh"
-
+using namespace CLHEP;
 namespace S1SCORERS
 {
    // This Threshold is used in all scorers
diff --git a/NPSimulation/ComptonTelescope/ComptonTelescope.hh b/NPSimulation/ComptonTelescope/ComptonTelescope.hh
index f0a931c4b..f01b5bee4 100644
--- a/NPSimulation/ComptonTelescope/ComptonTelescope.hh
+++ b/NPSimulation/ComptonTelescope/ComptonTelescope.hh
@@ -36,6 +36,7 @@
 // C++ headers
 #include <vector>
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/DummyDetector/DummyDetector.hh b/NPSimulation/DummyDetector/DummyDetector.hh
index 74b66c39e..d99d1f501 100644
--- a/NPSimulation/DummyDetector/DummyDetector.hh
+++ b/NPSimulation/DummyDetector/DummyDetector.hh
@@ -43,6 +43,7 @@
 #include "TDUMMYDetectorData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 class DUMMYDetector : public VDetector
 {
diff --git a/NPSimulation/Eurogam/Eurogam.hh b/NPSimulation/Eurogam/Eurogam.hh
index 714c663d5..b17313567 100644
--- a/NPSimulation/Eurogam/Eurogam.hh
+++ b/NPSimulation/Eurogam/Eurogam.hh
@@ -40,6 +40,7 @@
 #include "TEurogamData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 class Eurogam : public VDetector
 {
diff --git a/NPSimulation/GASPARD/GaspardTracker.hh b/NPSimulation/GASPARD/GaspardTracker.hh
index 69175a7d5..ff730db15 100644
--- a/NPSimulation/GASPARD/GaspardTracker.hh
+++ b/NPSimulation/GASPARD/GaspardTracker.hh
@@ -30,6 +30,8 @@
 #include "GaspardTrackerModule.hh"
 
 using namespace std;
+
using namespace CLHEP;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/GASPARD/GaspardTrackerAnnular.hh b/NPSimulation/GASPARD/GaspardTrackerAnnular.hh
index e5c2882ca..bffb43bde 100644
--- a/NPSimulation/GASPARD/GaspardTrackerAnnular.hh
+++ b/NPSimulation/GASPARD/GaspardTrackerAnnular.hh
@@ -30,6 +30,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/GASPARD/GaspardTrackerDummyShape.hh b/NPSimulation/GASPARD/GaspardTrackerDummyShape.hh
index 2295cba70..65b76dc28 100644
--- a/NPSimulation/GASPARD/GaspardTrackerDummyShape.hh
+++ b/NPSimulation/GASPARD/GaspardTrackerDummyShape.hh
@@ -33,6 +33,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/GASPARD/GaspardTrackerModule.hh b/NPSimulation/GASPARD/GaspardTrackerModule.hh
index efd99b36b..72ee7ddd9 100644
--- a/NPSimulation/GASPARD/GaspardTrackerModule.hh
+++ b/NPSimulation/GASPARD/GaspardTrackerModule.hh
@@ -37,6 +37,7 @@
 #include "TGaspardTrackerData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/GASPARD/GaspardTrackerSquare.hh b/NPSimulation/GASPARD/GaspardTrackerSquare.hh
index f10b9cb9d..03651afcf 100644
--- a/NPSimulation/GASPARD/GaspardTrackerSquare.hh
+++ b/NPSimulation/GASPARD/GaspardTrackerSquare.hh
@@ -29,6 +29,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/GASPARD/GaspardTrackerTrapezoid.hh b/NPSimulation/GASPARD/GaspardTrackerTrapezoid.hh
index 9834babc3..01756747e 100644
--- a/NPSimulation/GASPARD/GaspardTrackerTrapezoid.hh
+++ b/NPSimulation/GASPARD/GaspardTrackerTrapezoid.hh
@@ -31,6 +31,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Helios/Helios.hh b/NPSimulation/Helios/Helios.hh
index 5272bbb87..3ade31aac 100644
--- a/NPSimulation/Helios/Helios.hh
+++ b/NPSimulation/Helios/Helios.hh
@@ -31,6 +31,7 @@
 #include "HeliosModule.hh"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Helios/HeliosDetDummyShape.hh b/NPSimulation/Helios/HeliosDetDummyShape.hh
index 08564b61d..7afa65fec 100644
--- a/NPSimulation/Helios/HeliosDetDummyShape.hh
+++ b/NPSimulation/Helios/HeliosDetDummyShape.hh
@@ -32,6 +32,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Helios/HeliosModule.hh b/NPSimulation/Helios/HeliosModule.hh
index b3a0b89f5..1717514d5 100644
--- a/NPSimulation/Helios/HeliosModule.hh
+++ b/NPSimulation/Helios/HeliosModule.hh
@@ -36,6 +36,7 @@
 #include "THeliosData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Hyde2/Hyde2Tracker.hh b/NPSimulation/Hyde2/Hyde2Tracker.hh
index 00c26ff5a..98d9573c0 100644
--- a/NPSimulation/Hyde2/Hyde2Tracker.hh
+++ b/NPSimulation/Hyde2/Hyde2Tracker.hh
@@ -31,6 +31,7 @@
 #include "Hyde2TrackerModule.hh"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Hyde2/Hyde2TrackerModule.hh b/NPSimulation/Hyde2/Hyde2TrackerModule.hh
index 4ea863639..0d6966c04 100644
--- a/NPSimulation/Hyde2/Hyde2TrackerModule.hh
+++ b/NPSimulation/Hyde2/Hyde2TrackerModule.hh
@@ -36,6 +36,7 @@
 #include "THyde2TrackerData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Hyde2/Hyde2TrackerSquare1.hh b/NPSimulation/Hyde2/Hyde2TrackerSquare1.hh
index ab6f53da6..8fbb4f9ee 100644
--- a/NPSimulation/Hyde2/Hyde2TrackerSquare1.hh
+++ b/NPSimulation/Hyde2/Hyde2TrackerSquare1.hh
@@ -30,6 +30,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Hyde2/Hyde2TrackerTrapezoid1.hh b/NPSimulation/Hyde2/Hyde2TrackerTrapezoid1.hh
index 595b9b115..0a2fe0cab 100644
--- a/NPSimulation/Hyde2/Hyde2TrackerTrapezoid1.hh
+++ b/NPSimulation/Hyde2/Hyde2TrackerTrapezoid1.hh
@@ -30,6 +30,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Hyde2/Hyde2TrackerTrapezoid2.hh b/NPSimulation/Hyde2/Hyde2TrackerTrapezoid2.hh
index 3f9559e2e..bb8ce7ca6 100644
--- a/NPSimulation/Hyde2/Hyde2TrackerTrapezoid2.hh
+++ b/NPSimulation/Hyde2/Hyde2TrackerTrapezoid2.hh
@@ -30,6 +30,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/MUST2/MUST2Scorers.hh b/NPSimulation/MUST2/MUST2Scorers.hh
index 46fb11198..58ba7d834 100644
--- a/NPSimulation/MUST2/MUST2Scorers.hh
+++ b/NPSimulation/MUST2/MUST2Scorers.hh
@@ -26,6 +26,8 @@
  *****************************************************************************/
 #include "G4VPrimitiveScorer.hh"
 #include "G4THitsMap.hh"
+using namespace CLHEP;
+
 namespace MUST2 {
    //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    //   This Threshold is used in the above scorer. Any energy deposit under this threshold will not create an entry.
diff --git a/NPSimulation/Paris/Paris.hh b/NPSimulation/Paris/Paris.hh
index 956f47cca..3ee1aa4f5 100644
--- a/NPSimulation/Paris/Paris.hh
+++ b/NPSimulation/Paris/Paris.hh
@@ -30,9 +30,9 @@
 // NPTool header
 #include "VDetector.hh"
 #include "ParisModule.hh"
-
 using namespace std;
-
+#include "CLHEP/Units/SystemOfUnits.h"
+using namespace CLHEP;
 
 
 class Paris : public VDetector
diff --git a/NPSimulation/Paris/ParisCluster.hh b/NPSimulation/Paris/ParisCluster.hh
index e5c9f8797..6dcb2bddf 100644
--- a/NPSimulation/Paris/ParisCluster.hh
+++ b/NPSimulation/Paris/ParisCluster.hh
@@ -31,6 +31,8 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+#include "CLHEP/Units/SystemOfUnits.h"
+using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Paris/ParisModule.hh b/NPSimulation/Paris/ParisModule.hh
index 3cf779b4c..89524b999 100644
--- a/NPSimulation/Paris/ParisModule.hh
+++ b/NPSimulation/Paris/ParisModule.hh
@@ -37,6 +37,8 @@
 #include "TParisData.h"
 
 using namespace std;
+#include"CLHEP/Units/SystemOfUnits.h"
+using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Paris/ParisPhoswich.hh b/NPSimulation/Paris/ParisPhoswich.hh
index 36da6f1d3..1d6008141 100644
--- a/NPSimulation/Paris/ParisPhoswich.hh
+++ b/NPSimulation/Paris/ParisPhoswich.hh
@@ -31,6 +31,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Plastic/Plastic.hh b/NPSimulation/Plastic/Plastic.hh
index c2c96a648..eb1231cea 100644
--- a/NPSimulation/Plastic/Plastic.hh
+++ b/NPSimulation/Plastic/Plastic.hh
@@ -41,6 +41,7 @@
 #include "TPlasticData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 class Plastic : public VDetector
 {
diff --git a/NPSimulation/SSSD/ThinSi.hh b/NPSimulation/SSSD/ThinSi.hh
index df28fda41..df8bcc5ef 100644
--- a/NPSimulation/SSSD/ThinSi.hh
+++ b/NPSimulation/SSSD/ThinSi.hh
@@ -38,6 +38,7 @@
 // NPLib 
 #include "TSSSDData.h"
 using namespace std;
+
using namespace CLHEP;
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 namespace THINSI
diff --git a/NPSimulation/SSSD/ThinSiScorers.hh b/NPSimulation/SSSD/ThinSiScorers.hh
index 61e348a19..afb115bfd 100644
--- a/NPSimulation/SSSD/ThinSiScorers.hh
+++ b/NPSimulation/SSSD/ThinSiScorers.hh
@@ -26,7 +26,8 @@
  *****************************************************************************/
 #include "G4VPrimitiveScorer.hh"
 #include "G4THitsMap.hh"
-
+#include "CLHEP/Units/SystemOfUnits.h"
+using namespace CLHEP;
 namespace THINSI {
    //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    //   This Threshold is used in the above scorer. Any energy deposit under this threshold will not create an entry.
diff --git a/NPSimulation/Sharc/Sharc.hh b/NPSimulation/Sharc/Sharc.hh
index 5249029ea..6f6ad269d 100644
--- a/NPSimulation/Sharc/Sharc.hh
+++ b/NPSimulation/Sharc/Sharc.hh
@@ -40,6 +40,7 @@
 // NPLib
 #include "TSharcData.h"
 using namespace std;
+
using namespace CLHEP;
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 namespace SHARC
diff --git a/NPSimulation/Sharc/SharcScorers.hh b/NPSimulation/Sharc/SharcScorers.hh
index 08c59c9e9..54f0fee13 100644
--- a/NPSimulation/Sharc/SharcScorers.hh
+++ b/NPSimulation/Sharc/SharcScorers.hh
@@ -29,6 +29,7 @@
 
 #include <map>
 using namespace std;
+
using namespace CLHEP;
 
 namespace SHARC {
   //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
diff --git a/NPSimulation/Shield/Shield.hh b/NPSimulation/Shield/Shield.hh
index b60a27578..5eac901e6 100644
--- a/NPSimulation/Shield/Shield.hh
+++ b/NPSimulation/Shield/Shield.hh
@@ -32,6 +32,7 @@
 #include "ShieldModule.hh"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Shield/ShieldClParis.hh b/NPSimulation/Shield/ShieldClParis.hh
index c0260d423..fe3ab4cf0 100644
--- a/NPSimulation/Shield/ShieldClParis.hh
+++ b/NPSimulation/Shield/ShieldClParis.hh
@@ -33,6 +33,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Shield/ShieldModule.hh b/NPSimulation/Shield/ShieldModule.hh
index 1499880cd..d739c61c0 100644
--- a/NPSimulation/Shield/ShieldModule.hh
+++ b/NPSimulation/Shield/ShieldModule.hh
@@ -37,6 +37,7 @@
 #include "TShieldData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Shield/ShieldPhParis.hh b/NPSimulation/Shield/ShieldPhParis.hh
index b791dce89..9a45e0f22 100644
--- a/NPSimulation/Shield/ShieldPhParis.hh
+++ b/NPSimulation/Shield/ShieldPhParis.hh
@@ -32,6 +32,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/Tiara/Tiara.hh b/NPSimulation/Tiara/Tiara.hh
index 55ebdad91..6afb431cc 100644
--- a/NPSimulation/Tiara/Tiara.hh
+++ b/NPSimulation/Tiara/Tiara.hh
@@ -42,6 +42,7 @@
 #include "TTiaraHyballData.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 namespace TIARA
diff --git a/NPSimulation/Tigress/Tigress.hh b/NPSimulation/Tigress/Tigress.hh
index 7a36d2019..1baf166b9 100644
--- a/NPSimulation/Tigress/Tigress.hh
+++ b/NPSimulation/Tigress/Tigress.hh
@@ -40,6 +40,7 @@
 // NPLib
 #include "TTigressData.h"
 using namespace std;
+
using namespace CLHEP;
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 namespace TIGRESS
diff --git a/NPSimulation/W1/W1.hh b/NPSimulation/W1/W1.hh
index 483f3cbcd..eb6c20057 100644
--- a/NPSimulation/W1/W1.hh
+++ b/NPSimulation/W1/W1.hh
@@ -36,6 +36,7 @@
 #include "TW1Data.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/include/Chamber.hh b/NPSimulation/include/Chamber.hh
index e12941cf1..4d61de321 100644
--- a/NPSimulation/include/Chamber.hh
+++ b/NPSimulation/include/Chamber.hh
@@ -37,6 +37,7 @@
 #include "VDetector.hh"
 
 using namespace std;
+
using namespace CLHEP;
 
 class Chamber : public VDetector
 {
diff --git a/NPSimulation/include/EventGeneratorGammaDecay.hh b/NPSimulation/include/EventGeneratorGammaDecay.hh
index 3cfd1d2e0..42b50178e 100644
--- a/NPSimulation/include/EventGeneratorGammaDecay.hh
+++ b/NPSimulation/include/EventGeneratorGammaDecay.hh
@@ -33,7 +33,6 @@
 #include <string>
 #include <iostream>
 using namespace std;
-
 // NPSimulation
 #include "VEventGenerator.hh"
 #include "Target.hh"
@@ -76,4 +75,4 @@ public: // Managing the different cascade
   // Read all the added cscade en instentiate every thing that is needed
   void PrepareCascade();
 };
-#endif
\ No newline at end of file
+#endif
diff --git a/NPSimulation/include/EventGeneratorIsotropic.hh b/NPSimulation/include/EventGeneratorIsotropic.hh
index c2d7c8dd2..6a7b2ab02 100644
--- a/NPSimulation/include/EventGeneratorIsotropic.hh
+++ b/NPSimulation/include/EventGeneratorIsotropic.hh
@@ -24,6 +24,7 @@
 // C++ header
 #include <string>
 using namespace std;
+
using namespace CLHEP;
 
 // G4 headers
 #include "G4Event.hh"
diff --git a/NPSimulation/include/EventGeneratorParticleDecay.hh b/NPSimulation/include/EventGeneratorParticleDecay.hh
index c38bf69b4..22ec3b032 100644
--- a/NPSimulation/include/EventGeneratorParticleDecay.hh
+++ b/NPSimulation/include/EventGeneratorParticleDecay.hh
@@ -84,4 +84,4 @@ public: // Managing the decay
         // Set everything for the decay
   void SetDecay(vector<string> DaughterName, vector<bool> shoot, vector<double> ExcitationEnergy, string CSPath , string CSName);
 };
-#endif
\ No newline at end of file
+#endif
diff --git a/NPSimulation/include/EventGeneratorTwoBodyReaction.hh b/NPSimulation/include/EventGeneratorTwoBodyReaction.hh
index 467bd7116..a43006666 100644
--- a/NPSimulation/include/EventGeneratorTwoBodyReaction.hh
+++ b/NPSimulation/include/EventGeneratorTwoBodyReaction.hh
@@ -40,6 +40,7 @@
 #include "NPReaction.h"
 
 using namespace std;
+
using namespace CLHEP;
 using namespace NPL ;
 
 
diff --git a/NPSimulation/include/GeneralScorers.hh b/NPSimulation/include/GeneralScorers.hh
index 737158743..437f5f43a 100644
--- a/NPSimulation/include/GeneralScorers.hh
+++ b/NPSimulation/include/GeneralScorers.hh
@@ -26,7 +26,7 @@
  *****************************************************************************/
 #include "G4VPrimitiveScorer.hh"
 #include "G4THitsMap.hh"
-
+using namespace CLHEP;
 namespace GENERALSCORERS
    {
          //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
diff --git a/NPSimulation/include/MyMagneticField.hh b/NPSimulation/include/MyMagneticField.hh
index 252d61951..58f0cf886 100644
--- a/NPSimulation/include/MyMagneticField.hh
+++ b/NPSimulation/include/MyMagneticField.hh
@@ -26,6 +26,8 @@
 #include "globals.hh"
 #include "G4MagneticField.hh"
 #include "G4ThreeVector.hh"
+#include "CLHEP/Units/SystemOfUnits.h"
+using namespace CLHEP;
 
 class MyMagneticField : public G4MagneticField
 {
diff --git a/NPSimulation/include/Particle.hh b/NPSimulation/include/Particle.hh
index 1fdb3c45f..ece5946ed 100644
--- a/NPSimulation/include/Particle.hh
+++ b/NPSimulation/include/Particle.hh
@@ -26,7 +26,8 @@
 // G4
 #include"G4ParticleDefinition.hh"
 #include"G4ThreeVector.hh"
-
+#include"CLHEP/Units/SystemOfUnits.h"
+using namespace CLHEP;
 class Particle{
   
 public: // Constructor and Destructor
@@ -60,4 +61,4 @@ public: // Setter and Getter
   void SetParticleMomentumDirection(G4ThreeVector);
   void SetShootStatus(bool);
 };
-#endif
\ No newline at end of file
+#endif
diff --git a/NPSimulation/include/ParticleStack.hh b/NPSimulation/include/ParticleStack.hh
index 92abb3d86..d9375bd57 100644
--- a/NPSimulation/include/ParticleStack.hh
+++ b/NPSimulation/include/ParticleStack.hh
@@ -32,6 +32,7 @@
 #include <vector>
 #include <string>
 using namespace std;
+
using namespace CLHEP;
 
 // G4
 #include "G4ParticleGun.hh"
@@ -97,4 +98,4 @@ public: // Particle management and shooting method
     // Shoot everything in the stack at the end of the event
     void ShootAllParticle(G4Event* anEvent);
 };
-#endif
\ No newline at end of file
+#endif
diff --git a/NPSimulation/include/PhysicsList.hh b/NPSimulation/include/PhysicsList.hh
index 47f37522f..23c9cb275 100644
--- a/NPSimulation/include/PhysicsList.hh
+++ b/NPSimulation/include/PhysicsList.hh
@@ -24,7 +24,8 @@
  *****************************************************************************/
 #include "G4VUserPhysicsList.hh"
 #include "globals.hh"
-
+#include "CLHEP/Units/SystemOfUnits.h"
+using namespace CLHEP;
 class PhysicsList: public G4VUserPhysicsList
 {
 public:
diff --git a/NPSimulation/include/PrimaryGeneratorAction.hh b/NPSimulation/include/PrimaryGeneratorAction.hh
index 478d60ef9..fe37ad5a6 100644
--- a/NPSimulation/include/PrimaryGeneratorAction.hh
+++ b/NPSimulation/include/PrimaryGeneratorAction.hh
@@ -34,6 +34,7 @@
 #include "DetectorConstruction.hh"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/include/ResistiveStripScorers.hh b/NPSimulation/include/ResistiveStripScorers.hh
index 55f75cba7..dcb2ee59b 100644
--- a/NPSimulation/include/ResistiveStripScorers.hh
+++ b/NPSimulation/include/ResistiveStripScorers.hh
@@ -29,6 +29,7 @@
 
 #include <map>
 using namespace std;
+
using namespace CLHEP;
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
diff --git a/NPSimulation/include/Target.hh b/NPSimulation/include/Target.hh
index 0348a77ec..c2fdf2405 100644
--- a/NPSimulation/include/Target.hh
+++ b/NPSimulation/include/Target.hh
@@ -42,6 +42,7 @@
 #include "VDetector.hh"
 
 using namespace std;
+
using namespace CLHEP;
 
 class Target : public VDetector{
 public:
diff --git a/NPSimulation/include/VDetector.hh b/NPSimulation/include/VDetector.hh
index a5ff1c7e1..befe1d462 100644
--- a/NPSimulation/include/VDetector.hh
+++ b/NPSimulation/include/VDetector.hh
@@ -50,6 +50,7 @@
 #include "TInteractionCoordinates.h"
 
 using namespace std;
+
using namespace CLHEP;
 
 
 class VDetector
diff --git a/NPSimulation/include/VEventGenerator.hh b/NPSimulation/include/VEventGenerator.hh
index 8a46fab7e..b2a9dd718 100644
--- a/NPSimulation/include/VEventGenerator.hh
+++ b/NPSimulation/include/VEventGenerator.hh
@@ -43,6 +43,7 @@
 
 using namespace CLHEP;
 using namespace std;
+
using namespace CLHEP;
 
 
 
diff --git a/NPSimulation/src/GeneralScorers.cc b/NPSimulation/src/GeneralScorers.cc
index 2c8b40421..796883071 100644
--- a/NPSimulation/src/GeneralScorers.cc
+++ b/NPSimulation/src/GeneralScorers.cc
@@ -25,7 +25,7 @@
 #include "GeneralScorers.hh"
 #include "G4UnitsTable.hh"
 using namespace GENERALSCORERS ;
-
+using namespace CLHEP;
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 //   The following function is used in many scorer. following the Detector Volume Nomenclature
-- 
GitLab