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