00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044 #ifndef G4DATASET_HH
00045 #define G4DATASET_HH 1
00046
00047 #include <CLHEP/Units/SystemOfUnits.h>
00048
00049 #include "globals.hh"
00050 #include "G4IDataSet.hh"
00051
00052 class G4IInterpolator;
00053
00054 class G4DataSet : public G4IDataSet
00055 {
00056 public:
00057 G4DataSet(G4int argZ,
00058 G4IInterpolator* algo,
00059 G4double xUnit=CLHEP::MeV,
00060 G4double yUnit=CLHEP::barn,
00061 G4bool random=false);
00062
00063 G4DataSet(G4int argZ,
00064 G4DataVector* xData,
00065 G4DataVector* data,
00066 G4IInterpolator* algo,
00067 G4double xUnit=CLHEP::MeV,
00068 G4double yUnit=CLHEP::barn,
00069 G4bool random=false);
00070
00071 virtual ~G4DataSet();
00072
00073 virtual G4double FindValue(G4double x, G4int componentId=0) const;
00074
00075 virtual void PrintData(void) const;
00076
00077 virtual const G4IDataSet* GetComponent(G4int ) const { return 0; }
00078
00079 virtual void AddComponent(G4IDataSet* ) {}
00080
00081 virtual size_t NumberOfComponents(void) const { return 0; }
00082
00083 virtual const G4DataVector& GetEnergies(G4int ) const { return *energies; }
00084 virtual const G4DataVector& GetData(G4int ) const { return *data; }
00085 virtual void SetEnergiesData(G4DataVector* xData, G4DataVector* data, G4int componentId);
00086
00087 virtual G4bool LoadData(const G4String& fileName);
00088 virtual G4bool SaveData(const G4String& fileName) const;
00089
00090 virtual G4double RandomSelect(G4int componentId = 0) const;
00091
00092
00093 private:
00094
00095 size_t FindLowerBound(G4double energy) const;
00096 size_t FindLowerBound(G4double x, G4DataVector* values) const;
00097
00098 G4double IntegrationFunction(G4double x);
00099
00100 virtual void BuildPdf();
00101
00102 G4String FullFileName(const G4String& fileName) const;
00103
00104
00105 G4DataSet();
00106 G4DataSet(const G4DataSet& copy);
00107 G4DataSet& operator=(const G4DataSet& right);
00108
00109 G4int z;
00110
00111 G4DataVector* energies;
00112 G4DataVector* data;
00113
00114 G4IInterpolator* algorithm;
00115
00116 G4double unitEnergies;
00117 G4double unitData;
00118
00119 G4DataVector* pdf;
00120 G4bool randomSet;
00121 };
00122 #endif