Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 /// \file radiobiology/src/Manager.cc 28 /// \brief Implementation of the RadioBio::Man 29 30 #include "Manager.hh" 31 32 #include "G4AccumulableManager.hh" 33 34 #include <map> 35 #include <mutex> 36 37 std::mutex init_mutex; 38 39 namespace RadioBio 40 { 41 42 //....oooOO0OOooo........oooOO0OOooo........oo 43 44 #define width 15L 45 46 // To create the instance the DetectorConstruc 47 // be passed (to account for correct voxelizat 48 Manager* Manager::CreateInstance() 49 { 50 if (fInstance) { 51 G4Exception("RadioBioManager::createInstan 52 "Creating another, new, instan 53 delete fInstance; 54 } 55 56 fInstance = new Manager(); 57 return fInstance; 58 } 59 60 //....oooOO0OOooo........oooOO0OOooo........oo 61 62 Manager* Manager::GetInstance() 63 { 64 return fInstance; 65 } 66 67 //....oooOO0OOooo........oooOO0OOooo........oo 68 69 // Destructor deletes all the quantities 70 Manager::~Manager() 71 { 72 for (auto q : fQuantities) 73 delete q.second; 74 } 75 76 //....oooOO0OOooo........oooOO0OOooo........oo 77 78 void Manager::InitializeAll() 79 { 80 std::lock_guard<std::mutex> lock(init_mutex) 81 for (auto const& q : fQuantities) 82 (q.second)->Initialize(); 83 } 84 85 //....oooOO0OOooo........oooOO0OOooo........oo 86 87 void Manager::ComputeAll() 88 { 89 for (auto const& q : fQuantities) 90 (q.second)->Compute(); 91 } 92 93 //....oooOO0OOooo........oooOO0OOooo........oo 94 95 void Manager::ResetAll() 96 { 97 for (auto const& q : fQuantities) 98 (q.second)->Reset(); 99 } 100 101 //....oooOO0OOooo........oooOO0OOooo........oo 102 103 void Manager::StoreAll() 104 { 105 for (auto const& q : fQuantities) 106 (q.second)->Store(); 107 } 108 109 //....oooOO0OOooo........oooOO0OOooo........oo 110 111 VRadiobiologicalQuantity* Manager::GetQuantity 112 { 113 return fQuantities.find(str)->second; 114 } 115 116 //....oooOO0OOooo........oooOO0OOooo........oo 117 118 void Manager::PrintParameters() 119 { 120 G4cout << "********************************* 121 << "*** right now registered quantiti 122 for (auto const& q : fQuantities) 123 G4cout << "*** " << q.first; 124 G4cout << "*** but their calculation might b 125 << "*** active. Ask for parameters of 126 << "********************************* 127 } 128 129 //....oooOO0OOooo........oooOO0OOooo........oo 130 131 void Manager::DigestAccumulables() 132 { 133 for (auto q : fQuantities) { 134 // Hook in the accumulable manager the one 135 G4VAccumulable* GenAcc = G4AccumulableMana 136 137 if (!GenAcc) { 138 G4Exception("RadioBioManager::AddFromAcc 139 } 140 141 // If calculation is not set enabled, exit 142 if (!q.second->IsCalculationEnabled()) con 143 144 // Add from the accumulable. 145 q.second->AddFromAccumulable(GenAcc); 146 } 147 } 148 149 //....oooOO0OOooo........oooOO0OOooo........oo 150 151 bool Manager::Register(VRadiobiologicalQuantit 152 { 153 if (q == nullptr) { 154 G4Exception("RadioBioManager::Register", " 155 "Asking to register a quantity 156 return false; 157 } 158 159 if (fQuantities.find(name) != fQuantities.en 160 G4Exception("RadioBioManager::Register", " 161 "Registering two radiobiologic 162 return false; 163 } 164 fQuantities[name] = q; 165 return true; 166 } 167 168 //....oooOO0OOooo........oooOO0OOooo........oo 169 170 } // namespace RadioBio