Geant4  10.02.p01
G4DataSet.hh
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 // $Id: G4DataSet.hh 70904 2013-06-07 10:34:25Z gcosmo $
28 //
29 // Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch)
30 //
31 // History:
32 // -----------
33 // 31 Jul 2001 MGP Created as G4EMDataSet
34 // 31 Jul 2008 MGP Revised and renamed to G4DataSet
35 //
36 // -------------------------------------------------------------------
37 
38 // Class description:
39 // Leaf data set
40 // Applies a Composite design pattern for data library management
41 
42 // -------------------------------------------------------------------
43 
44 #ifndef G4DATASET_HH
45 #define G4DATASET_HH 1
46 
47 #include <CLHEP/Units/SystemOfUnits.h>
48 
49 #include "globals.hh"
50 #include "G4IDataSet.hh"
51 
52 class G4IInterpolator;
53 
54 class G4DataSet : public G4IDataSet
55 {
56 public:
57  G4DataSet(G4int argZ,
58  G4IInterpolator* algo,
59  G4double xUnit=CLHEP::MeV,
60  G4double yUnit=CLHEP::barn,
61  G4bool random=false);
62 
63  G4DataSet(G4int argZ,
64  G4DataVector* xData,
66  G4IInterpolator* algo,
67  G4double xUnit=CLHEP::MeV,
68  G4double yUnit=CLHEP::barn,
69  G4bool random=false);
70 
71  virtual ~G4DataSet();
72 
73  virtual G4double FindValue(G4double x, G4int componentId=0) const;
74 
75  virtual void PrintData(void) const;
76 
77  virtual const G4IDataSet* GetComponent(G4int /* componentId */) const { return 0; }
78 
79  virtual void AddComponent(G4IDataSet* /* dataSet */) {}
80 
81  virtual size_t NumberOfComponents(void) const { return 0; }
82 
83  virtual const G4DataVector& GetEnergies(G4int /* componentId */) const { return *energies; }
84  virtual const G4DataVector& GetData(G4int /* componentId */) const { return *data; }
85  virtual void SetEnergiesData(G4DataVector* xData, G4DataVector* data, G4int componentId);
86 
87  virtual G4bool LoadData(const G4String& fileName);
88  virtual G4bool SaveData(const G4String& fileName) const;
89 
90  virtual G4double RandomSelect(G4int componentId = 0) const;
91 
92 
93 private:
94 
95  size_t FindLowerBound(G4double energy) const;
96  size_t FindLowerBound(G4double x, G4DataVector* values) const;
97 
99 
100  virtual void BuildPdf();
101 
102  G4String FullFileName(const G4String& fileName) const;
103 
104  // Hide copy constructor and assignment operator
105  G4DataSet();
106  G4DataSet(const G4DataSet& copy);
108 
110 
111  G4DataVector* energies; // Owned pointer
112  G4DataVector* data; // Owned pointer
113 
114  G4IInterpolator* algorithm; // Owned pointer
115 
118 
121 };
122 #endif /* G4DATASET_HH */
G4DataVector * energies
Definition: G4DataSet.hh:111
G4double IntegrationFunction(G4double x)
Definition: G4DataSet.cc:411
static const double MeV
Definition: G4SIunits.hh:211
size_t FindLowerBound(G4double energy) const
Definition: G4DataSet.cc:303
G4DataVector * data
Definition: G4DataSet.hh:112
G4bool randomSet
Definition: G4DataSet.hh:120
virtual void AddComponent(G4IDataSet *)
Definition: G4DataSet.hh:79
virtual size_t NumberOfComponents(void) const
Definition: G4DataSet.hh:81
G4double unitData
Definition: G4DataSet.hh:117
G4IInterpolator * algorithm
Definition: G4DataSet.hh:114
int G4int
Definition: G4Types.hh:78
virtual void PrintData(void) const
Definition: G4DataSet.cc:131
virtual G4bool LoadData(const G4String &fileName)
Definition: G4DataSet.cc:176
virtual G4bool SaveData(const G4String &fileName) const
Definition: G4DataSet.cc:229
virtual const G4DataVector & GetData(G4int) const
Definition: G4DataSet.hh:84
G4DataSet & operator=(const G4DataSet &right)
virtual void BuildPdf()
Definition: G4DataSet.cc:353
bool G4bool
Definition: G4Types.hh:79
virtual G4double RandomSelect(G4int componentId=0) const
Definition: G4DataSet.cc:383
G4double unitEnergies
Definition: G4DataSet.hh:116
virtual ~G4DataSet()
Definition: G4DataSet.cc:106
virtual const G4IDataSet * GetComponent(G4int) const
Definition: G4DataSet.hh:77
virtual G4double FindValue(G4double x, G4int componentId=0) const
Definition: G4DataSet.cc:114
G4int z
Definition: G4DataSet.hh:109
G4String FullFileName(const G4String &fileName) const
Definition: G4DataSet.cc:337
G4double energy(const ThreeVector &p, const G4double m)
const G4double x[NPOINTSGL]
G4DataVector * pdf
Definition: G4DataSet.hh:119
static const double barn
Definition: G4SIunits.hh:104
virtual const G4DataVector & GetEnergies(G4int) const
Definition: G4DataSet.hh:83
virtual void SetEnergiesData(G4DataVector *xData, G4DataVector *data, G4int componentId)
Definition: G4DataSet.cc:151
double G4double
Definition: G4Types.hh:76