Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // 26 // >> 27 // $Id: G4VCrossSectionHandler.hh 66241 2012-12-13 18:34:42Z gunter $ 27 // 28 // 28 // Author: Maria Grazia Pia (Maria.Grazia.Pia@ 29 // Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch) 29 // 30 // 30 // History: 31 // History: 31 // ----------- 32 // ----------- 32 // 16 Sep 2001 MGP Created 33 // 16 Sep 2001 MGP Created 33 // 26 Sep 2001 V.Ivanchenko Hide copy 34 // 26 Sep 2001 V.Ivanchenko Hide copy constructor and assignement operator 34 // 18 Apr 2002 V.Ivanchenko Move membe 35 // 18 Apr 2002 V.Ivanchenko Move member function ValueForMaterial to public 35 // 21 Jan 2003 V.Ivanchenko Cut per re 36 // 21 Jan 2003 V.Ivanchenko Cut per region 36 // 15 Jul 2009 N.A.Karakatsanis New method 37 // 15 Jul 2009 N.A.Karakatsanis New methods added for loading logarithmic data 37 // to enhance 38 // to enhance computing performance of interpolation 38 // 39 // 39 // ------------------------------------------- 40 // ------------------------------------------------------------------- 40 41 41 // Class description: 42 // Class description: 42 // Low Energy Electromagnetic Physics 43 // Low Energy Electromagnetic Physics 43 // Base class for cross section manager for an 44 // Base class for cross section manager for an electromagnetic physics process 44 45 45 // ------------------------------------------- 46 // ------------------------------------------------------------------- 46 47 47 #ifndef G4VCROSSSECTIONHANDLER_HH 48 #ifndef G4VCROSSSECTIONHANDLER_HH 48 #define G4VCROSSSECTIONHANDLER_HH 1 49 #define G4VCROSSSECTIONHANDLER_HH 1 49 50 50 #include <map> 51 #include <map> 51 #include <vector> 52 #include <vector> 52 #include <CLHEP/Units/SystemOfUnits.h> 53 #include <CLHEP/Units/SystemOfUnits.h> 53 54 54 #include "globals.hh" 55 #include "globals.hh" 55 #include "G4DataVector.hh" 56 #include "G4DataVector.hh" 56 #include "G4MaterialCutsCouple.hh" 57 #include "G4MaterialCutsCouple.hh" 57 58 58 class G4VDataSetAlgorithm; 59 class G4VDataSetAlgorithm; 59 class G4VEMDataSet; 60 class G4VEMDataSet; 60 class G4Material; 61 class G4Material; 61 class G4Element; 62 class G4Element; 62 63 63 class G4VCrossSectionHandler { 64 class G4VCrossSectionHandler { 64 65 65 public: 66 public: 66 explicit G4VCrossSectionHandler(); << 67 67 explicit G4VCrossSectionHandler(G4VDataSetAl << 68 G4VCrossSectionHandler(); >> 69 >> 70 G4VCrossSectionHandler(G4VDataSetAlgorithm* interpolation, 68 G4double minE = 250*CLHEP::eV, 71 G4double minE = 250*CLHEP::eV, 69 G4double maxE = 100*C 72 G4double maxE = 100*CLHEP::GeV, 70 G4int nBins = 200, 73 G4int nBins = 200, 71 G4double unitE = CLHEP::MeV, 74 G4double unitE = CLHEP::MeV, 72 G4double unitData = C 75 G4double unitData = CLHEP::barn, 73 G4int minZ = 1, G4int maxZ = 99); 76 G4int minZ = 1, G4int maxZ = 99); 74 77 75 virtual ~G4VCrossSectionHandler(); 78 virtual ~G4VCrossSectionHandler(); 76 79 77 void Initialise(G4VDataSetAlgorithm* interpo << 80 void Initialise(G4VDataSetAlgorithm* interpolation = 0, 78 G4double minE = 250*CLHEP::eV, 81 G4double minE = 250*CLHEP::eV, 79 G4double maxE = 100*CLHEP::G 82 G4double maxE = 100*CLHEP::GeV, 80 G4int numberOfBins = 200, 83 G4int numberOfBins = 200, 81 G4double unitE = CLHEP::MeV, 84 G4double unitE = CLHEP::MeV, 82 G4double unitData = CLHEP::b 85 G4double unitData = CLHEP::barn, 83 G4int minZ = 1, G4int maxZ = 99); 86 G4int minZ = 1, G4int maxZ = 99); 84 87 85 G4int SelectRandomAtom(const G4MaterialCutsC 88 G4int SelectRandomAtom(const G4MaterialCutsCouple* couple, G4double e) const; >> 89 86 const G4Element* SelectRandomElement(const G 90 const G4Element* SelectRandomElement(const G4MaterialCutsCouple* material, 87 G4double e) const; 91 G4double e) const; 88 92 89 G4int SelectRandomShell(G4int Z, G4double e) 93 G4int SelectRandomShell(G4int Z, G4double e) const; 90 G4VEMDataSet* BuildMeanFreePathForMaterials( << 94 >> 95 G4VEMDataSet* BuildMeanFreePathForMaterials(const G4DataVector* energyCuts = 0); >> 96 91 G4double FindValue(G4int Z, G4double e) cons 97 G4double FindValue(G4int Z, G4double e) const; >> 98 92 G4double FindValue(G4int Z, G4double e, G4in 99 G4double FindValue(G4int Z, G4double e, G4int shellIndex) const; >> 100 93 G4double ValueForMaterial(const G4Material* 101 G4double ValueForMaterial(const G4Material* material, G4double e) const; >> 102 94 void LoadData(const G4String& dataFile); 103 void LoadData(const G4String& dataFile); >> 104 95 void LoadNonLogData(const G4String& dataFile 105 void LoadNonLogData(const G4String& dataFile); >> 106 96 void LoadShellData(const G4String& dataFile) 107 void LoadShellData(const G4String& dataFile); >> 108 97 void PrintData() const; 109 void PrintData() const; 98 void Clear(); << 99 110 100 G4VCrossSectionHandler(const G4VCrossSection << 111 void Clear(); 101 G4VCrossSectionHandler & operator=(const G4V << 102 112 103 protected: 113 protected: 104 114 105 G4int NumberOfComponents(G4int Z) const; 115 G4int NumberOfComponents(G4int Z) const; >> 116 106 void ActiveElements(); 117 void ActiveElements(); 107 118 108 // Factory method 119 // Factory method 109 virtual std::vector<G4VEMDataSet*>* BuildCro 120 virtual std::vector<G4VEMDataSet*>* BuildCrossSectionsForMaterials(const G4DataVector& energyVector, 110 const G4DataVector* ene << 121 const G4DataVector* energyCuts = 0) = 0; 111 122 112 // Factory method 123 // Factory method 113 virtual G4VDataSetAlgorithm* CreateInterpola 124 virtual G4VDataSetAlgorithm* CreateInterpolation(); >> 125 114 const G4VDataSetAlgorithm* GetInterpolation( 126 const G4VDataSetAlgorithm* GetInterpolation() const { return interpolation; } 115 127 >> 128 116 private: 129 private: 117 G4VDataSetAlgorithm* interpolation; << 118 G4DataVector activeZ; << 119 130 120 std::map<G4int,G4VEMDataSet*,std::less<G4int << 131 // Hide copy constructor and assignment operator 121 std::vector<G4VEMDataSet*>* crossSections; << 132 G4VCrossSectionHandler(const G4VCrossSectionHandler&); >> 133 G4VCrossSectionHandler & operator=(const G4VCrossSectionHandler &right); >> 134 >> 135 G4VDataSetAlgorithm* interpolation; 122 136 123 G4double eMin; 137 G4double eMin; 124 G4double eMax; 138 G4double eMax; >> 139 G4int nBins; >> 140 125 G4double unit1; 141 G4double unit1; 126 G4double unit2; 142 G4double unit2; 127 143 128 G4int zMin; 144 G4int zMin; 129 G4int zMax; 145 G4int zMax; 130 G4int nBins; << 146 >> 147 G4DataVector activeZ; >> 148 >> 149 std::map<G4int,G4VEMDataSet*,std::less<G4int> > dataMap; >> 150 >> 151 std::vector<G4VEMDataSet*>* crossSections; 131 }; 152 }; 132 153 133 #endif 154 #endif 134 155 135 156 136 157 137 158 138 159 139 160 140 161 141 162 142 163 143 164 144 165 145 166