| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 // 23 //
>> 24 // $Id: G4EMDataSet.hh,v 1.4 2003/05/20 20:03:34 pia Exp $
>> 25 // GEANT4 tag $Name: geant4-05-02 $
27 // 26 //
28 // Author: Maria Grazia Pia (Maria.Grazia.Pia@ 27 // Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch)
29 // 28 //
30 // History: 29 // History:
31 // ----------- 30 // -----------
32 // 31 Jul 2001 MGP Created << 31 // 31 Jul 2001 MGP Created
33 // 32 //
34 // 15 Jul 2009 N.A.Karakatsanis New metho <<
35 // to enhanc <<
36 // ------------------------------------------- 33 // -------------------------------------------------------------------
37 34
38 // Class description: 35 // Class description:
39 // Low Energy Electromagnetic Physics 36 // Low Energy Electromagnetic Physics
40 // Data set for an electromagnetic physics pro 37 // Data set for an electromagnetic physics process
41 // A strategy pattern is used to encapsulate a 38 // A strategy pattern is used to encapsulate algorithms for data interpolation
>> 39 // Further documentation available from http://www.ge.infn.it/geant4/lowE
>> 40
42 // ------------------------------------------- 41 // -------------------------------------------------------------------
43 42
44 #ifndef G4EMDATASET_HH 43 #ifndef G4EMDATASET_HH
45 #define G4EMDATASET_HH 1 44 #define G4EMDATASET_HH 1
46 45
47 #include <CLHEP/Units/SystemOfUnits.h> <<
48 #include "globals.hh" 46 #include "globals.hh"
>> 47 #include "G4DataVector.hh"
49 #include "G4VEMDataSet.hh" 48 #include "G4VEMDataSet.hh"
50 49
51 class G4VDataSetAlgorithm; 50 class G4VDataSetAlgorithm;
52 51
53 class G4EMDataSet : public G4VEMDataSet << 52 class G4EMDataSet : public G4VEMDataSet {
54 { << 53
55 public: 54 public:
56 explicit G4EMDataSet(G4int argZ, <<
57 G4VDataSetAlgorithm* algo, <<
58 G4double xUnit=CLHEP::MeV, <<
59 G4double yUnit=CLHEP::barn, <<
60 G4bool random=false); <<
61 <<
62 explicit G4EMDataSet(G4int argZ, <<
63 G4DataVector* xData, <<
64 G4DataVector* data, <<
65 G4VDataSetAlgorithm* algo, <<
66 G4double xUnit=CLHEP::MeV, <<
67 G4double yUnit=CLHEP::barn, <<
68 G4bool random=false); <<
69 <<
70 explicit G4EMDataSet(G4int argZ, <<
71 G4DataVector* xData, <<
72 G4DataVector* data, <<
73 G4DataVector* xLogData, <<
74 G4DataVector* Logdata, <<
75 G4VDataSetAlgorithm* algo, <<
76 G4double xUnit=CLHEP::MeV, <<
77 G4double yUnit=CLHEP::barn, <<
78 G4bool random=false); <<
79 55
80 virtual ~G4EMDataSet(); << 56 G4EMDataSet(G4int Z,
>> 57 G4DataVector* points,
>> 58 G4DataVector* values,
>> 59 const G4VDataSetAlgorithm* interpolation,
>> 60 G4double unitE = MeV, G4double unitData = barn);
>> 61
>> 62 G4EMDataSet(G4int Z,
>> 63 const G4String& dataFile,
>> 64 const G4VDataSetAlgorithm* interpolation,
>> 65 G4double unitE = MeV, G4double unitData = barn);
>> 66
>> 67 ~G4EMDataSet();
81 68
82 virtual G4double FindValue(G4double x, G4int << 69 G4double FindValue(G4double e, G4int id = 0) const;
83 70
84 virtual void PrintData(void) const; << 71 virtual const G4VEMDataSet* GetComponent(G4int i) const;
85 72
86 virtual const G4VEMDataSet* GetComponent(G4i << 73 virtual void AddComponent(G4VEMDataSet* dataSet);
87 74
88 virtual void AddComponent(G4VEMDataSet* /* d << 75 virtual size_t NumberOfComponents() const;
89 76
90 virtual std::size_t NumberOfComponents(void) << 77 void PrintData() const;
91 78
92 virtual const G4DataVector& GetEnergies(G4in << 79 const G4DataVector& GetEnergies(G4int) const { return *energies; }
93 virtual const G4DataVector& GetData(G4int /* << 80 const G4DataVector& GetData(G4int) const { return *data; }
94 virtual const G4DataVector& GetLogEnergies(G <<
95 virtual const G4DataVector& GetLogData(G4int <<
96 81
97 virtual void SetEnergiesData(G4DataVector* x << 82 private:
98 virtual void SetLogEnergiesData(G4DataVector << 83
99 G4DataVector << 84 // Hide copy constructor and assignment operator
100 G4DataVector << 85 G4EMDataSet& operator=(const G4EMDataSet& right);
101 G4DataVector << 86 G4EMDataSet(const G4EMDataSet&);
102 G4int compon << 87
>> 88 void LoadData(const G4String& dataFile);
>> 89
>> 90 G4int FindBinLocation(G4double energy) const;
>> 91
>> 92 G4int z;
>> 93
>> 94 G4DataVector* energies; // Owned pointer
>> 95 G4DataVector* data; // Owned pointer
>> 96
>> 97 const G4VDataSetAlgorithm* algorithm; // Owned pointer
>> 98
>> 99 G4double unit1;
>> 100 G4double unit2;
>> 101
>> 102 size_t numberOfBins;
>> 103
>> 104 };
>> 105
>> 106 #endif
>> 107
103 108
104 109
105 virtual G4bool LoadData(const G4String& file <<
106 virtual G4bool LoadNonLogData(const G4String <<
107 110
108 virtual G4bool SaveData(const G4String& file <<
109 111
110 virtual G4double RandomSelect(G4int componen <<
111 <<
112 112
113 private: <<
114 113
115 std::size_t FindLowerBound(G4double energy) <<
116 std::size_t FindLowerBound(G4double x, G4Dat <<
117 114
118 G4double IntegrationFunction(G4double x); <<
119 virtual void BuildPdf(); <<
120 G4String FullFileName(const G4String& fileNa <<
121 115
122 // Hide copy constructor and assignment oper <<
123 explicit G4EMDataSet(); <<
124 G4EMDataSet(const G4EMDataSet& copy) = delet <<
125 G4EMDataSet& operator=(const G4EMDataSet& ri <<
126 <<
127 G4DataVector* energies; // Owned <<
128 G4DataVector* data; // Owned <<
129 G4DataVector* log_energies; // Owned <<
130 G4DataVector* log_data; // Owned <<
131 G4VDataSetAlgorithm* algorithm; // Owned <<
132 G4DataVector* pdf; <<
133 116
134 G4double unitEnergies; <<
135 G4double unitData; <<
136 117
137 G4int z; <<
138 G4bool randomSet; <<
139 }; <<
140 #endif /* G4EMDATASET_HH */ <<
141 118