| Geant4 Cross Reference |
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: Mathieu Karamitros (kara (AT) cenbg 27 // Author: Mathieu Karamitros (kara (AT) cenbg . in2p3 . fr)
28 // 28 //
29 // History: 29 // History:
30 // ----------- 30 // -----------
31 // 10 Oct 2011 M.Karamitros created 31 // 10 Oct 2011 M.Karamitros created
32 // 32 //
33 // ------------------------------------------- 33 // -------------------------------------------------------------------
34 34
35 #include "G4ITModelManager.hh" 35 #include "G4ITModelManager.hh"
36 #include "G4VITStepModel.hh" <<
37 #include "G4ITType.hh" 36 #include "G4ITType.hh"
38 #include "G4UnitsTable.hh" 37 #include "G4UnitsTable.hh"
39 38
>> 39 using namespace std;
40 40
41 G4ITModelManager::G4ITModelManager() << 41 G4ITModelManager::G4ITModelManager() : fIsInitialized(FALSE)
42 << 42 {;}
43 = default; <<
44 43
45 G4ITModelManager::~G4ITModelManager() = defaul << 44 G4ITModelManager::~G4ITModelManager()
46 <<
47 void G4ITModelManager::Initialize() <<
48 { 45 {
49 std::sort(fModelInfoList.begin(), fModelIn << 46 //dtor
50 [](const ModelInfo& lhs, const M << 47 mapModels::iterator it ;
51 return lhs.fStartingTime < r <<
52 }); <<
53 48
54 for (const auto& modelInfo : fModelInfoLis << 49 for(it = fModels.begin() ; it != fModels.end() ; it ++)
55 { 50 {
56 modelInfo.fpModel->Initialize(); << 51 delete it->second;
57 } 52 }
58 << 53 fModels.clear();
59 fIsInitialized = true; <<
60 } 54 }
61 55
62 void G4ITModelManager::SetModel(G4VITStepModel << 56 G4ITModelManager::G4ITModelManager(const G4ITModelManager& right)
63 const G4double <<
64 const G4double <<
65 { 57 {
66 assert(pModel != nullptr); << 58 mapModels::const_iterator it = right.fModels.begin();
67 if (fIsInitialized) << 59
>> 60 for ( ; it != right.fModels.end() ; it++)
68 { 61 {
69 G4ExceptionDescription exceptionDescri << 62 fModels[it->first] = it->second->Clone();
70 exceptionDescription <<
71 << "You are trying to insert a new <<
72 G4Exception("G4ITModelManager::SetMode <<
73 FatalErrorInArgument, exce <<
74 } 63 }
75 64
76 fModelInfoList.emplace_back(ModelInfo({ st << 65 fIsInitialized = right.fIsInitialized;
77 } 66 }
78 67
79 std::vector<G4VITStepModel*> G4ITModelManager: << 68 G4ITModelManager& G4ITModelManager::operator=(const G4ITModelManager& rhs)
80 { 69 {
81 std::vector<G4VITStepModel*> activeModels; << 70 if (this == &rhs) return *this; // handle self assignment
>> 71 //assignment operator
>> 72 return *this;
>> 73 }
82 74
83 for (const auto& modelInfo : fModelInfoLis << 75 void G4ITModelManager::Initialize()
>> 76 {
>> 77 mapModels::iterator it = fModels.begin();
>> 78
>> 79 for( ; it != fModels.end() ; it ++)
84 { 80 {
85 if (modelInfo.fStartingTime < globalTi << 81 G4VITModel* model = it ->second;
>> 82 if(model!=0)
86 { 83 {
87 activeModels.push_back(modelInfo.f << 84 model->Initialize();
88 } 85 }
89 } 86 }
>> 87 }
90 88
91 return activeModels; << 89 void G4ITModelManager::SetModel(G4VITModel* aModel, G4double startingTime)
>> 90 {
>> 91 assert(fIsInitialized == FALSE);
>> 92 if(fIsInitialized == true)
>> 93 {
>> 94 G4ExceptionDescription exceptionDescription ;
>> 95 exceptionDescription << "You are trying to insert a new model after initialization of th model manager.";
>> 96 G4Exception("G4ITModelManager::SetModel","ITModelManager001",
>> 97 FatalErrorInArgument,exceptionDescription);
>> 98 }
>> 99 fModels[startingTime] = aModel;
92 } 100 }
93 101
>> 102 G4VITModel* G4ITModelManager::GetModel(const G4double globalTime)
>> 103 {
>> 104 if(!fModels.empty())
>> 105 {
>> 106 mapModels::reverse_iterator rit = fModels.rbegin();
>> 107 if(rit != fModels.rend())
>> 108 {
>> 109 if(globalTime > rit->first)
>> 110 {
>> 111 return rit->second;
>> 112 }
>> 113 else
>> 114 {
>> 115 mapModels::iterator it = fModels.begin();
>> 116
>> 117 if(globalTime < it->first)
>> 118 {
>> 119 G4ExceptionDescription exceptionDescription ;
>> 120 exceptionDescription << "No model was found at time ";
>> 121 exceptionDescription << G4BestUnit(globalTime,"Time");
>> 122 exceptionDescription << ". The first model is registered at time : ";
>> 123 exceptionDescription << G4BestUnit(it->first,"Time") << ". ";
>> 124 G4Exception("G4ITModelManager::GetModel","ITModelManager003",
>> 125 FatalErrorInArgument,exceptionDescription);
>> 126 }
>> 127
>> 128 it = fModels.lower_bound(globalTime);
>> 129
>> 130 if(it != fModels.end()) return it->second;
>> 131 }
>> 132 }
>> 133 }
>> 134
>> 135 G4ExceptionDescription exceptionDescription ;
>> 136 exceptionDescription << "No model was found.";
>> 137 G4Exception("G4ITModelManager::GetModel","ITModelManager004",
>> 138 FatalErrorInArgument,exceptionDescription);
>> 139 return 0;
>> 140 }
94 141