Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4ITModelManager.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27// Author: Mathieu Karamitros (kara (AT) cenbg . in2p3 . fr)
28//
29// History:
30// -----------
31// 10 Oct 2011 M.Karamitros created
32//
33// -------------------------------------------------------------------
34
35#include "G4ITModelManager.hh"
36#include "G4VITStepModel.hh"
37#include "G4ITType.hh"
38#include "G4UnitsTable.hh"
39
40
42
43= default;
44
46
48{
49 std::sort(fModelInfoList.begin(), fModelInfoList.end(),
50 [](const ModelInfo& lhs, const ModelInfo& rhs) {
51 return lhs.fStartingTime < rhs.fStartingTime;
52 });
53
54 for (const auto& modelInfo : fModelInfoList)
55 {
56 modelInfo.fpModel->Initialize();
57 }
58
59 fIsInitialized = true;
60}
61
63 const G4double startingTime,
64 const G4double endTime)
65{
66 assert(pModel != nullptr);
68 {
69 G4ExceptionDescription exceptionDescription;
70 exceptionDescription
71 << "You are trying to insert a new model after initializing the model manager.";
72 G4Exception("G4ITModelManager::SetModel", "ITModelManager001",
73 FatalErrorInArgument, exceptionDescription);
74 }
75
76 fModelInfoList.emplace_back(ModelInfo({ startingTime, endTime, std::unique_ptr<G4VITStepModel>(pModel) }));
77}
78
79std::vector<G4VITStepModel*> G4ITModelManager::GetActiveModels(G4double globalTime) const
80{
81 std::vector<G4VITStepModel*> activeModels;
82
83 for (const auto& modelInfo : fModelInfoList)
84 {
85 if (modelInfo.fStartingTime < globalTime && modelInfo.fEndTime > globalTime)
86 {
87 activeModels.push_back(modelInfo.fpModel.get());
88 }
89 }
90
91 return activeModels;
92}
93
@ FatalErrorInArgument
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
std::ostringstream G4ExceptionDescription
double G4double
Definition G4Types.hh:83
std::vector< G4VITStepModel * > GetActiveModels(G4double globalTime) const
std::vector< ModelInfo > fModelInfoList
void SetModel(G4VITStepModel *aModel, G4double startingTime, G4double endTime=DBL_MAX)