Move FairModule::fAllSensitiveVolumes to FairRunSim

Instead of using a thread local, let's use something that
is tied to the "simulaion session" (as noted in some places).

FairRunSim has all the Modules anyway, so let it also have
the list of sensitive Volumes.

Author: ChristianTackeGSI

fairroot/base/sim/FairDetector.cxx

@@ -14,6 +14,7 @@
 
 #include "FairGeoNode.h"   // for FairGeoNode
 #include "FairRootManager.h"
+#include "FairRunSim.h"
 #include "FairVolume.h"   // for FairVolume
 
 #include <TFolder.h>       // for TFolder
@@ -97,7 +98,7 @@ void FairDetector::Initialize()
     FairGeoNode* fN;
     TString cutName;
     TString copysign = "#";
-    for (auto aVol : fAllSensitiveVolumes) {
+    for (auto aVol : FairRunSim::Instance()->fAllSensitiveVolumes) {
         cutName = aVol->GetName();
         Ssiz_t pos = cutName.Index(copysign, 1);
         if (pos > 1) {

fairroot/base/sim/FairMCApplication.cxx

@@ -886,7 +886,7 @@ void FairMCApplication::InitGeometry()
     fMCEventHeader->SetRunID(runId);
 
     // Fill sensitive volumes in fVolMap
-    for (auto fv : FairModule::fAllSensitiveVolumes) {
+    for (auto fv : fRun->fAllSensitiveVolumes) {
         if (!fv) {
             continue;
         }

fairroot/base/sim/FairModule.cxx

@@ -20,7 +20,7 @@
 #include "FairGeoNode.h"        // for FairGeoNode
 #include "FairGeoParSet.h"      // for FairBaseParSet
 #include "FairMCApplication.h"
-#include "FairRun.h"          // for FairRun
+#include "FairRunSim.h"       // for FairRun, FairRunSim
 #include "FairRuntimeDb.h"    // for FairRuntimeDb
 #include "FairVolume.h"       // for FairVolume
 #include "FairVolumeList.h"   // for FairVolumeList
@@ -47,8 +47,6 @@
 #include <map>
 #include <memory>
 
-thread_local std::vector<FairVolume*> FairModule::fAllSensitiveVolumes;
-
 void FairModule::ConstructGeometry()
 {
     LOG(warn)
@@ -214,11 +212,12 @@ void FairModule::ProcessNodes(TList* aList)
         vList = new FairVolumeList();
     }
 
+    FairRunSim& run = *(FairRunSim::Instance());
     TListIter iter(aList);
     FairGeoNode* node = nullptr;
     FairGeoNode* MotherNode = nullptr;
     FairVolume* volume = nullptr;
-    FairRuntimeDb* rtdb = FairRun::Instance()->GetRuntimeDb();
+    FairRuntimeDb* rtdb = run.GetRuntimeDb();
     FairGeoParSet* par = static_cast<FairGeoParSet*>(rtdb->getContainer("FairGeoParSet"));
     TObjArray* fNodes = par->GetGeoNodes();
     while ((node = static_cast<FairGeoNode*>(iter.Next()))) {
@@ -240,7 +239,7 @@ void FairModule::ProcessNodes(TList* aList)
         if (node->isSensitive() && fActive) {
             volume->setModId(fModId);
             volume->SetModule(this);
-            fAllSensitiveVolumes.push_back(volume);
+            run.fAllSensitiveVolumes.push_back(volume);
             aVol = dynamic_cast<FairGeoVolume*>(node);
             fNodes->AddLast(aVol);
             fNbOfSensitiveVol++;
@@ -260,7 +259,7 @@ void FairModule::AddSensitiveVolume(TGeoVolume* v)
         vList->addVolume(volume);
         volume->setModId(fModId);
         volume->SetModule(this);
-        fAllSensitiveVolumes.push_back(volume);
+        FairRunSim::Instance()->fAllSensitiveVolumes.push_back(volume);
         fNbOfSensitiveVol++;
     }
 }

fairroot/base/sim/FairModule.h

@@ -22,11 +22,9 @@
 #include <TList.h>       // for TList (ptr only), TListIter
 #include <TNamed.h>      // for TNamed
 #include <TObjArray.h>   // for TObjArray
-#include <TRefArray.h>   // for TRefArray
 #include <TString.h>     // for TString, operator!=
 #include <TVirtualMC.h>
 #include <string>
-#include <vector>
 
 class FairVolumeList;
 class FairVolume;
@@ -139,8 +137,6 @@ class FairModule : public TNamed
     static thread_local inline FairVolumeList* vList{nullptr};   //!
     /**total number of volumes in a simulaion session*/
     static thread_local inline Int_t fNbOfVolumes{0};   //!
-    /**list of all sensitive volumes in  a simulaion session*/
-    static thread_local std::vector<FairVolume*> fAllSensitiveVolumes;   //!
 
     TString fMotherVolumeName{""};   //!
     FairVolume* getFairVolume(FairGeoNode* fNode);

fairroot/base/steer/FairRunSim.h

@@ -1,5 +1,5 @@
 /********************************************************************************
- * Copyright (C) 2014-2023 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH  *
+ * Copyright (C) 2014-2024 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH  *
  *                                                                              *
  *              This software is distributed under the terms of the             *
  *              GNU Lesser General Public Licence (LGPL) version 3,             *
@@ -21,12 +21,14 @@
 #include <functional>
 #include <memory>
 #include <utility>
+#include <vector>
 
 class FairField;
 class FairMCEventHeader;
 class FairMesh;
 class FairModule;
 class FairPrimaryGenerator;
+class FairVolume;
 
 /**
  * \brief Configure the Simulation session
@@ -37,6 +39,10 @@ class FairPrimaryGenerator;
  */
 class FairRunSim : public FairRun
 {
+    friend class FairModule;
+    friend class FairDetector;
+    friend class FairMCApplication;
+
   public:
     /** default ctor*/
     FairRunSim(Bool_t isMaster = kTRUE);
@@ -223,6 +229,9 @@ class FairRunSim : public FairRun
 
     bool fWasMT{false};   //!                              /** Actual MT mode used */
 
+    /**list of all sensitive volumes in  a simulaion session*/
+    std::vector<FairVolume*> fAllSensitiveVolumes;   //!
+
   protected:
     Int_t count{0};                     //!< Internal counter
     FairMCApplication* fApp{nullptr};   //!< Main VMC application