Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/analysis/hdf5/src/G4Hdf5AnalysisManager.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /analysis/hdf5/src/G4Hdf5AnalysisManager.cc (Version 11.3.0) and /analysis/hdf5/src/G4Hdf5AnalysisManager.cc (Version 11.0.p2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26                                                    26 
 27 // Author: Ivana Hrivnacova, 20/07/2017 (ivana     27 // Author: Ivana Hrivnacova, 20/07/2017 (ivana@ipno.in2p3.fr)
 28                                                    28 
 29 #include "G4Hdf5AnalysisManager.hh"                29 #include "G4Hdf5AnalysisManager.hh"
 30 #include "G4Hdf5FileManager.hh"                    30 #include "G4Hdf5FileManager.hh"
 31 #include "G4Hdf5NtupleFileManager.hh"              31 #include "G4Hdf5NtupleFileManager.hh"
 32 #include "G4AnalysisManagerState.hh"               32 #include "G4AnalysisManagerState.hh"
 33 #include "G4AnalysisUtilities.hh"                  33 #include "G4AnalysisUtilities.hh"
 34 #include "G4ThreadLocalSingleton.hh"               34 #include "G4ThreadLocalSingleton.hh"
 35 #include "G4Threading.hh"                          35 #include "G4Threading.hh"
 36 #include "G4AutoLock.hh"                           36 #include "G4AutoLock.hh"
 37                                                    37 
 38 using namespace G4Analysis;                        38 using namespace G4Analysis;
 39                                                    39 
 40 // mutex in a file scope                           40 // mutex in a file scope
 41                                                    41 
 42 namespace {                                        42 namespace {
 43   //Mutex to lock master manager when opening      43   //Mutex to lock master manager when opening a file
 44   G4Mutex openFileMutex = G4MUTEX_INITIALIZER;     44   G4Mutex openFileMutex = G4MUTEX_INITIALIZER;
 45   //Mutex to lock master manager when closing      45   //Mutex to lock master manager when closing a file
 46   G4Mutex closeFileMutex = G4MUTEX_INITIALIZER     46   G4Mutex closeFileMutex = G4MUTEX_INITIALIZER;
 47 }                                                  47 }
 48                                                    48 
 49 // G4Hdf5AnalysisManager* G4Hdf5AnalysisManage     49 // G4Hdf5AnalysisManager* G4Hdf5AnalysisManager::fgMasterInstance = nullptr;
 50 // G4ThreadLocal G4bool G4Hdf5AnalysisManager:     50 // G4ThreadLocal G4bool G4Hdf5AnalysisManager::fgIsInstance = false;
 51                                                    51 
 52 //____________________________________________     52 //_____________________________________________________________________________
 53 G4Hdf5AnalysisManager* G4Hdf5AnalysisManager::     53 G4Hdf5AnalysisManager* G4Hdf5AnalysisManager::Instance()
 54 {                                                  54 {
 55   static G4ThreadLocalSingleton<G4Hdf5Analysis     55   static G4ThreadLocalSingleton<G4Hdf5AnalysisManager> instance;
 56   fgIsInstance = true;                             56   fgIsInstance = true;
 57   return instance.Instance();                      57   return instance.Instance();
 58 }                                                  58 }
 59                                                    59 
 60 //____________________________________________     60 //_____________________________________________________________________________
 61 G4bool G4Hdf5AnalysisManager::IsInstance()         61 G4bool G4Hdf5AnalysisManager::IsInstance()
 62 {                                                  62 {
 63   return fgIsInstance;                             63   return fgIsInstance;
 64 }                                                  64 }
 65                                                    65 
 66 //____________________________________________     66 //_____________________________________________________________________________
 67 G4Hdf5AnalysisManager::G4Hdf5AnalysisManager()     67 G4Hdf5AnalysisManager::G4Hdf5AnalysisManager()
 68  : G4ToolsAnalysisManager("Hdf5")                  68  : G4ToolsAnalysisManager("Hdf5")
 69 {                                                  69 {
 70 #ifdef G4MULTITHREADED                             70 #ifdef G4MULTITHREADED
 71 #ifndef H5_HAVE_THREADSAFE                         71 #ifndef H5_HAVE_THREADSAFE
 72     G4Exception("G4Hdf5AnalysisManager::G4Hdf5     72     G4Exception("G4Hdf5AnalysisManager::G4Hdf5AnalysisManager",
 73                 "Analysis_F001", FatalExceptio     73                 "Analysis_F001", FatalException,
 74                 "Your HDF5 lib is not built wi     74                 "Your HDF5 lib is not built with H5_HAVE_THREADSAFE.");
 75 #endif                                             75 #endif
 76 #endif                                             76 #endif
 77                                                    77 
                                                   >>  78   if ( ! G4Threading::IsWorkerThread() ) fgMasterInstance = this;
                                                   >>  79 
 78   // File manager                                  80   // File manager
 79   auto fileManager = std::make_shared<G4Hdf5Fi <<  81   fFileManager = std::make_shared<G4Hdf5FileManager>(fState);
 80   SetFileManager(fileManager);                 <<  82   SetFileManager(fFileManager);
 81                                                    83 
 82   // Ntuple file manager                           84   // Ntuple file manager
 83   fNtupleFileManager = std::make_shared<G4Hdf5     85   fNtupleFileManager = std::make_shared<G4Hdf5NtupleFileManager>(fState);
 84   SetNtupleFileManager(fNtupleFileManager);    <<  86   fNtupleFileManager->SetFileManager(fFileManager);
 85   fNtupleFileManager->SetFileManager(fileManag << 
 86   fNtupleFileManager->SetBookingManager(fNtupl     87   fNtupleFileManager->SetBookingManager(fNtupleBookingManager);
 87 }                                                  88 }
 88                                                    89 
 89 //____________________________________________     90 //_____________________________________________________________________________
 90 G4Hdf5AnalysisManager::~G4Hdf5AnalysisManager(     91 G4Hdf5AnalysisManager::~G4Hdf5AnalysisManager()
 91 {                                                  92 {
                                                   >>  93   if ( fState.GetIsMaster() ) fgMasterInstance = nullptr;
 92   fgIsInstance = false;                            94   fgIsInstance = false;
 93 }                                                  95 }
 94                                                    96 
 95 //                                                 97 //
 96 // protected methods                               98 // protected methods
 97 //                                                 99 //
 98                                                   100 
 99 //____________________________________________    101 //_____________________________________________________________________________
100 G4bool G4Hdf5AnalysisManager::OpenFileImpl(con    102 G4bool G4Hdf5AnalysisManager::OpenFileImpl(const G4String& fileName)
101 {                                                 103 {
                                                   >> 104   // Create ntuple manager(s)
                                                   >> 105   // and set it to base class which takes then their ownership
                                                   >> 106   SetNtupleManager(fNtupleFileManager->CreateNtupleManager());
                                                   >> 107 
                                                   >> 108   auto result = true;
                                                   >> 109 
102   G4AutoLock lock(&openFileMutex);                110   G4AutoLock lock(&openFileMutex);
103   auto result = G4ToolsAnalysisManager::OpenFi << 111   result &= fFileManager->OpenFile(fileName);
                                                   >> 112   result &= fNtupleFileManager->ActionAtOpenFile(fFileManager->GetFullFileName());
104   lock.unlock();                                  113   lock.unlock();
105                                                   114 
106   return result;                                  115   return result;
107 }                                                 116 }
108                                                   117 
109 //____________________________________________    118 //_____________________________________________________________________________
                                                   >> 119 G4bool G4Hdf5AnalysisManager::WriteImpl()
                                                   >> 120 {
                                                   >> 121   auto result = true;
                                                   >> 122 
                                                   >> 123   Message(kVL4, "write", "files");
                                                   >> 124 
                                                   >> 125   if ( G4Threading::IsWorkerThread() )  {
                                                   >> 126     result &= G4ToolsAnalysisManager::Merge();
                                                   >> 127   }
                                                   >> 128   else {
                                                   >> 129     // Open all files registered with objects
                                                   >> 130     fFileManager->OpenFiles();
                                                   >> 131 
                                                   >> 132     // Write all histograms/profile on master
                                                   >> 133     result &= G4ToolsAnalysisManager::WriteImpl();
                                                   >> 134   }
                                                   >> 135 
                                                   >> 136   // Write ASCII if activated
                                                   >> 137   if ( IsAscii() ) {
                                                   >> 138     result &= WriteAscii(fFileManager->GetFileName());
                                                   >> 139   }
                                                   >> 140 
                                                   >> 141   Message(kVL3, "write", "files", "", result);
                                                   >> 142 
                                                   >> 143   return result;
                                                   >> 144 }
                                                   >> 145 
                                                   >> 146 //_____________________________________________________________________________
110 G4bool G4Hdf5AnalysisManager::CloseFileImpl(G4    147 G4bool G4Hdf5AnalysisManager::CloseFileImpl(G4bool reset)
111 {                                                 148 {
                                                   >> 149   auto result = true;
                                                   >> 150 
112   G4AutoLock lock(&closeFileMutex);               151   G4AutoLock lock(&closeFileMutex);
113   auto result = G4ToolsAnalysisManager::CloseF << 152   result &= fFileManager->CloseFiles();
                                                   >> 153 
                                                   >> 154   if ( reset ) {
                                                   >> 155     // reset data
                                                   >> 156     result = ResetImpl();
                                                   >> 157     if ( ! result ) {
                                                   >> 158       Warn("Resetting data failed", fkClass, "OpenDirectory");
                                                   >> 159     }
                                                   >> 160   }
                                                   >> 161   result &= fNtupleFileManager->ActionAtCloseFile(reset);
                                                   >> 162 
114   lock.unlock();                                  163   lock.unlock();
                                                   >> 164 
                                                   >> 165   return result;
                                                   >> 166 }
                                                   >> 167 
                                                   >> 168 //_____________________________________________________________________________
                                                   >> 169 G4bool G4Hdf5AnalysisManager::ResetImpl()
                                                   >> 170 {
                                                   >> 171   auto result = true;
                                                   >> 172 
                                                   >> 173   result &= G4ToolsAnalysisManager::ResetImpl();
                                                   >> 174   if ( fNtupleFileManager != nullptr ) {
                                                   >> 175     result &= fNtupleFileManager->Reset();
                                                   >> 176   }
115                                                   177 
116   return result;                                  178   return result;
117 }                                                 179 }
118                                                   180