added DeltaChord flag

FullSimLight/fullSimLight.cc

@@ -35,7 +35,9 @@
 
 #include <nlohmann/json.hpp>
 #include <fstream>
-
+#include "HistoMessenger.hh"
+#include "HistoManager.hh"
+#include "G4StepLimiterPhysics.hh"
 
 static const std::string fullSimLightShareDir=FULLSIMLIGHTSHAREDIR;
 static std::string  parMacroFileName   = fullSimLightShareDir+"/macro.g4";
@@ -167,6 +169,11 @@ int main(int argc, char** argv) {
     const FSLPhysListFactory *phyListFactory = FSLPhysListFactory::GetInstance();
     G4VModularPhysicsList *physList = phyListFactory->GetPhysList(parPhysListName);
     G4bool activateRegions = phyListFactory->GetActivateRegionsFlag();
+
+
+  physList->RegisterPhysics(new G4StepLimiterPhysics());	///added
+
+
     
     // register the final version of the physics list in the run manager
     runManager->SetUserInitialization(physList);
@@ -175,6 +182,8 @@ int main(int argc, char** argv) {
     FSLDetectorConstruction* detector = new FSLDetectorConstruction;
     ///define here deltachord///////
     detector->SetDeltaChord(parDeltaChord);
+    G4cout << "Main cc Delta Chord Length: " << parDeltaChord << G4endl;
+
     //////////////////////////////////////////////////////////////////
     if (parRunOverlapCheck) detector->SetRunOverlapCheck(true);
     if (activateRegions)    detector->SetAddRegions(true);
@@ -354,7 +363,7 @@ void GetInputArguments(int argc, char** argv) {
     }
     while (true) {
         int c, optidx = 0;
-        c = getopt_long(argc, argv, "P:pm:f:g:c:oh", options, &optidx);
+        c = getopt_long(argc, argv, "P:pm:f:g:c:ohd:", options, &optidx);
         if (c == -1)
             break;
         //