ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4DNAMolecularMaterial.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4DNAMolecularMaterial.hh
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
28 //
29 // We would be very happy hearing from you, send us your feedback! :)
30 //
31 // In order for Geant4-DNA to be maintained and still open-source,
32 // article citations are crucial.
33 // If you use Geant4-DNA chemistry and you publish papers about your software,
34 // in addition to the general paper on Geant4-DNA:
35 //
36 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
37 //
38 // we would be very happy if you could please also cite the following
39 // reference papers on chemistry:
40 //
41 // J. Comput. Phys. 274 (2014) 841-882
42 // Prog. Nucl. Sci. Tec. 2 (2011) 503-508
43 
44 #ifndef G4DNAMolecularMaterial_HH
45 #define G4DNAMolecularMaterial_HH
46 
47 #include "globals.hh"
48 #include "G4ios.hh"
49 #include <map>
50 #include <vector>
51 #include "G4VStateDependent.hh"
52 
53 class G4Material;
55 
64 {
65  bool operator()(const G4Material* mat1, const G4Material* mat2) const;
66 };
67 
68 typedef std::map<const G4Material*, double, CompareMaterial> ComponentMap;
69 
96 {
97 public:
99  static void DeleteInstance();
100  void Initialize();
101  void Clear();
102 
103  virtual G4bool Notify(G4ApplicationState requestedState);
104 
105  //----------------------------------------------------------------------------
106 
123  const std::vector<double>* GetDensityTableFor(const G4Material*) const;
124 
140  const std::vector<double>* GetNumMolPerVolTableFor(const G4Material*) const;
141 
142  inline const std::vector<ComponentMap>* GetMassFractionTable() const{
143  return fpCompFractionTable;
144  }
145  inline const std::vector<ComponentMap>* GetDensityTable() const{
146  return fpCompDensityTable;
147  }
148 
149  //----------------------------------------------------------------------------
150 
152 
167 
181  const G4String&);
182 
195  void SetMolecularConfiguration(const G4String& materialName,
196  const G4String& molUserIF);
197 
198  //----------------------------------------------------------------------------
199 
209 
219  const G4Material *component,
220  G4double massFraction);
221 
222 protected:
227  virtual ~G4DNAMolecularMaterial();
228  void Create();
229  void InitializeNumMolPerVol();
230  void InitializeDensity();
231  void RecordMolecularMaterial(G4Material* parentMaterial,
232  G4Material* molecularMaterial,
233  G4double fraction);
234  void SearchMolecularMaterial(G4Material* parentMaterial,
236  double currentFraction);
237 
238  void AddMaterial(const G4Material*, double fraction);
239 
240  void PrintNotAMolecularMaterial(const char* methodName,
241  const G4Material* lookForMaterial) const;
242 
243  // Tables built for all molecular materials at initialization
244  std::vector<ComponentMap>* fpCompFractionTable;
245  std::vector<ComponentMap>* fpCompDensityTable;
246  std::vector<ComponentMap>* fpCompNumMolPerVolTable;
247 
248  mutable std::map<const G4Material*, std::vector<double>*, CompareMaterial>
250  mutable std::map<const G4Material*, std::vector<double>*, CompareMaterial>
252  mutable std::map<const G4Material*, bool, CompareMaterial> fWarningPrinted;
253 
254  std::map<int /*Material ID*/,
256 
258  size_t fNMaterials;
259 };
260 
261 #endif // G4DNAMolecularMaterial_HH