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-patch-01 $ 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