Add an option to activate Geant4 command-line scoring (#79)

dded an option to activate Geant4 command-line scoring

- The activation should be done via calling
   'TG4RunConfiguration::SetUseOfG4Scoring()' in user g4Config.C
- The scoring (mesh, scorers etc.) can be then defined via
   Geant UI commands defined in '/score' directory in g4config.in
- The scoring data are then automatically saved for each mesh
   in a 'mesh_name.txt' file at the end of run

Author: ihrivnac

source/run/include/TG4RunConfiguration.h

@@ -102,6 +102,7 @@ class TG4RunConfiguration
   void SetMTApplication(Bool_t mtApplication);
   void SetParameter(const TString& name, Double_t value);
   void SetSpecialCutsOld();
+  void SetUseOfG4Scoring();
 
   // get methods
   TString GetUserGeometry() const;
@@ -110,6 +111,7 @@ class TG4RunConfiguration
   Bool_t IsSpecialControls() const;
   Bool_t IsSpecialCuts() const;
   Bool_t IsSpecialCutsOld() const;
+  Bool_t IsUseOfG4Scoring() const;
   Bool_t IsMTApplication() const;
 
  protected:
@@ -122,6 +124,7 @@ class TG4RunConfiguration
   Bool_t fSpecialControls;          ///< option for special controls
   Bool_t fSpecialCuts;              ///< option for special cuts
   Bool_t fSpecialCutsOld;           ///< option for special cuts old
+  Bool_t fUseOfG4Scoring;           ///< option to activate G4 commmand-line scoring
   G4UImessenger* fAGDDMessenger;    //!< XML messenger
   G4UImessenger* fGDMLMessenger;    //!< XML messenger
 
@@ -141,10 +144,21 @@ class TG4RunConfiguration
 
 // inline functions
 
+/// Activate G4 commmand-line scoring
+inline void TG4RunConfiguration::SetUseOfG4Scoring()
+{
+  fUseOfG4Scoring = true;
+}
+
 /// Return physics list selection
 inline TString TG4RunConfiguration::GetPhysicsListSelection() const
 { 
   return fPhysicsListSelection;
 }
 
+inline Bool_t TG4RunConfiguration::IsUseOfG4Scoring() const
+{
+  return fUseOfG4Scoring;
+}
+
 #endif // TG4V_RUN_CONFIGURATION_H

source/run/src/TG4RunConfiguration.cxx

@@ -49,6 +49,7 @@ TG4RunConfiguration::TG4RunConfiguration(const TString& userGeometry,
     fSpecialControls(false),
     fSpecialCuts(false),
     fSpecialCutsOld(false),
+    fUseOfG4Scoring(false),
     fAGDDMessenger(0),
     fGDMLMessenger(0),
     fParameters()

source/run/src/TG4RunManager.cxx

@@ -45,12 +45,16 @@
 #include <G4RunManager.hh>
 #endif
 
+#include <G4ScoringManager.hh>
+#include <G4VScoringMesh.hh>
+
 #include <G4UIExecutive.hh>
 #include <G4UImanager.hh>
 #include <G4UIsession.hh>
 #include <G4Version.hh>
 #include <Randomize.hh>
 
+
 #ifdef USE_G4ROOT
 #include <TG4RootNavMgr.h>
 #endif
@@ -131,6 +135,11 @@ TG4RunManager::TG4RunManager(
 
     // create and configure G4 run manager
     ConfigureRunManager();
+
+    if (runConfiguration->IsUseOfG4Scoring()) {
+      // activate G4 command-line scoring
+      G4ScoringManager::GetScoringManager();
+    }
   }
   else {
     // Get G4 worker run manager
@@ -601,6 +610,16 @@ Bool_t TG4RunManager::FinishRun()
   G4bool result = !TG4SDServices::Instance()->GetIsStopRun();
   TG4SDServices::Instance()->SetIsStopRun(false);
 
+  if (fRunConfiguration->IsUseOfG4Scoring()) {
+    // Dump all scoring meshes in file
+    auto g4ScoringManager = G4ScoringManager::GetScoringManager();
+    for (std::size_t i = 0; i < g4ScoringManager->GetNumberOfMesh(); ++i) {
+      auto meshName = g4ScoringManager->GetMesh(i)->GetWorldName();
+      auto fileName = meshName + ".txt";
+      g4ScoringManager->DumpAllQuantitiesToFile(meshName, fileName);
+    }
+  }
+
   return result;
 }