| Geant4 Cross Reference |
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 /// \file scavenger/include/ParserChemReaction.hh
27 /// \brief Definition of the scavenger::ParserChemReaction class
28
29 #ifndef SCAVENGER_ParserChemReaction_h
30 #define SCAVENGER_ParserChemReaction_h 1
31
32 #include "globals.hh"
33
34 #include <array>
35 #include <vector>
36
37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
38 namespace scavenger
39 {
40
41 /// Parser to read user files defining chemical reactions and
42 /// scavengers (reaction with background)
43
44 class ParserChemReaction
45 {
46 public:
47 ParserChemReaction() = default;
48
49 ~ParserChemReaction() = default;
50
51 void ReadReactionFile(const G4String& fileName);
52
53 inline auto GetListReactant1() { return fListReactant1; }
54
55 inline auto GetListReactant2() { return fListReactant2; }
56
57 inline auto GetListProduct() { return fListProduct; }
58
59 inline auto GetListRate() { return fListRate; }
60
61 private:
62 void AddReaction(const G4String& reactionString, const G4String& type);
63
64 void ReadReservoir(const G4String& reservoirString);
65
66 G4double GetScavengerConcentration(const G4String& name);
67
68 static void ReadReaction(const G4String& reactionString, std::vector<G4String>& reactant,
69 std::vector<G4String>& product, G4double& reactionRate);
70
71 void ImplementReaction(const G4String& reactant1, const G4String& reactant2,
72 const std::vector<G4String>& product, const G4double& reactionRate,
73 const G4String& type);
74
75 static void ReplaceString(G4String& aString, const G4String& from, const G4String& to);
76
77 std::array<std::vector<G4String>, 5> fListReactant1;
78 std::array<std::vector<G4String>, 5> fListReactant2;
79 std::array<std::vector<std::vector<G4String>>, 5> fListProduct;
80 std::array<std::vector<G4double>, 5> fListRate;
81 std::map<G4String, G4double> fReservoirConcentrationMap;
82 };
83
84 } // namespace scavenger
85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
86
87 #endif
88