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 // 27 // 28 // Author: Maria Grazia Pia (Maria.Grazia.Pia@ 28 // Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch) 29 // 29 // 30 // History: 30 // History: 31 // ----------- 31 // ----------- 32 // 16 Sep 2001 MGP Created 32 // 16 Sep 2001 MGP Created 33 // 26 Sep 2001 V.Ivanchenko Hide copy 33 // 26 Sep 2001 V.Ivanchenko Hide copy constructor and assignement operator 34 // 18 Apr 2002 V.Ivanchenko Move membe 34 // 18 Apr 2002 V.Ivanchenko Move member function ValueForMaterial to public 35 // 21 Jan 2003 V.Ivanchenko Cut per re 35 // 21 Jan 2003 V.Ivanchenko Cut per region 36 // 15 Jul 2009 N.A.Karakatsanis New method 36 // 15 Jul 2009 N.A.Karakatsanis New methods added for loading logarithmic data 37 // to enhance 37 // to enhance computing performance of interpolation 38 // 38 // 39 // ------------------------------------------- 39 // ------------------------------------------------------------------- 40 40 41 // Class description: 41 // Class description: 42 // Low Energy Electromagnetic Physics 42 // Low Energy Electromagnetic Physics 43 // Base class for cross section manager for an 43 // Base class for cross section manager for an electromagnetic physics process 44 44 45 // ------------------------------------------- 45 // ------------------------------------------------------------------- 46 46 47 #ifndef G4VCROSSSECTIONHANDLER_HH 47 #ifndef G4VCROSSSECTIONHANDLER_HH 48 #define G4VCROSSSECTIONHANDLER_HH 1 48 #define G4VCROSSSECTIONHANDLER_HH 1 49 49 50 #include <map> 50 #include <map> 51 #include <vector> 51 #include <vector> 52 #include <CLHEP/Units/SystemOfUnits.h> 52 #include <CLHEP/Units/SystemOfUnits.h> 53 53 54 #include "globals.hh" 54 #include "globals.hh" 55 #include "G4DataVector.hh" 55 #include "G4DataVector.hh" 56 #include "G4MaterialCutsCouple.hh" 56 #include "G4MaterialCutsCouple.hh" 57 57 58 class G4VDataSetAlgorithm; 58 class G4VDataSetAlgorithm; 59 class G4VEMDataSet; 59 class G4VEMDataSet; 60 class G4Material; 60 class G4Material; 61 class G4Element; 61 class G4Element; 62 62 63 class G4VCrossSectionHandler { 63 class G4VCrossSectionHandler { 64 64 65 public: 65 public: 66 explicit G4VCrossSectionHandler(); << 66 67 explicit G4VCrossSectionHandler(G4VDataSetAl << 67 G4VCrossSectionHandler(); >> 68 >> 69 G4VCrossSectionHandler(G4VDataSetAlgorithm* interpolation, 68 G4double minE = 250*CLHEP::eV, 70 G4double minE = 250*CLHEP::eV, 69 G4double maxE = 100*C 71 G4double maxE = 100*CLHEP::GeV, 70 G4int nBins = 200, 72 G4int nBins = 200, 71 G4double unitE = CLHEP::MeV, 73 G4double unitE = CLHEP::MeV, 72 G4double unitData = C 74 G4double unitData = CLHEP::barn, 73 G4int minZ = 1, G4int maxZ = 99); 75 G4int minZ = 1, G4int maxZ = 99); 74 76 75 virtual ~G4VCrossSectionHandler(); 77 virtual ~G4VCrossSectionHandler(); 76 78 77 void Initialise(G4VDataSetAlgorithm* interpo << 79 void Initialise(G4VDataSetAlgorithm* interpolation = 0, 78 G4double minE = 250*CLHEP::eV, 80 G4double minE = 250*CLHEP::eV, 79 G4double maxE = 100*CLHEP::G 81 G4double maxE = 100*CLHEP::GeV, 80 G4int numberOfBins = 200, 82 G4int numberOfBins = 200, 81 G4double unitE = CLHEP::MeV, 83 G4double unitE = CLHEP::MeV, 82 G4double unitData = CLHEP::b 84 G4double unitData = CLHEP::barn, 83 G4int minZ = 1, G4int maxZ = 99); 85 G4int minZ = 1, G4int maxZ = 99); 84 86 85 G4int SelectRandomAtom(const G4MaterialCutsC 87 G4int SelectRandomAtom(const G4MaterialCutsCouple* couple, G4double e) const; >> 88 86 const G4Element* SelectRandomElement(const G 89 const G4Element* SelectRandomElement(const G4MaterialCutsCouple* material, 87 G4double e) const; 90 G4double e) const; 88 91 89 G4int SelectRandomShell(G4int Z, G4double e) 92 G4int SelectRandomShell(G4int Z, G4double e) const; 90 G4VEMDataSet* BuildMeanFreePathForMaterials( << 93 >> 94 G4VEMDataSet* BuildMeanFreePathForMaterials(const G4DataVector* energyCuts = 0); >> 95 91 G4double FindValue(G4int Z, G4double e) cons 96 G4double FindValue(G4int Z, G4double e) const; >> 97 92 G4double FindValue(G4int Z, G4double e, G4in 98 G4double FindValue(G4int Z, G4double e, G4int shellIndex) const; >> 99 93 G4double ValueForMaterial(const G4Material* 100 G4double ValueForMaterial(const G4Material* material, G4double e) const; >> 101 94 void LoadData(const G4String& dataFile); 102 void LoadData(const G4String& dataFile); >> 103 95 void LoadNonLogData(const G4String& dataFile 104 void LoadNonLogData(const G4String& dataFile); >> 105 96 void LoadShellData(const G4String& dataFile) 106 void LoadShellData(const G4String& dataFile); >> 107 97 void PrintData() const; 108 void PrintData() const; 98 void Clear(); << 99 109 100 G4VCrossSectionHandler(const G4VCrossSection << 110 void Clear(); 101 G4VCrossSectionHandler & operator=(const G4V << 102 111 103 protected: 112 protected: 104 113 105 G4int NumberOfComponents(G4int Z) const; 114 G4int NumberOfComponents(G4int Z) const; >> 115 106 void ActiveElements(); 116 void ActiveElements(); 107 117 108 // Factory method 118 // Factory method 109 virtual std::vector<G4VEMDataSet*>* BuildCro 119 virtual std::vector<G4VEMDataSet*>* BuildCrossSectionsForMaterials(const G4DataVector& energyVector, 110 const G4DataVector* ene << 120 const G4DataVector* energyCuts = 0) = 0; 111 121 112 // Factory method 122 // Factory method 113 virtual G4VDataSetAlgorithm* CreateInterpola 123 virtual G4VDataSetAlgorithm* CreateInterpolation(); >> 124 114 const G4VDataSetAlgorithm* GetInterpolation( 125 const G4VDataSetAlgorithm* GetInterpolation() const { return interpolation; } 115 126 >> 127 116 private: 128 private: 117 G4VDataSetAlgorithm* interpolation; << 118 G4DataVector activeZ; << 119 129 120 std::map<G4int,G4VEMDataSet*,std::less<G4int << 130 // Hide copy constructor and assignment operator 121 std::vector<G4VEMDataSet*>* crossSections; << 131 G4VCrossSectionHandler(const G4VCrossSectionHandler&); >> 132 G4VCrossSectionHandler & operator=(const G4VCrossSectionHandler &right); >> 133 >> 134 G4VDataSetAlgorithm* interpolation; 122 135 123 G4double eMin; 136 G4double eMin; 124 G4double eMax; 137 G4double eMax; >> 138 G4int nBins; >> 139 125 G4double unit1; 140 G4double unit1; 126 G4double unit2; 141 G4double unit2; 127 142 128 G4int zMin; 143 G4int zMin; 129 G4int zMax; 144 G4int zMax; 130 G4int nBins; << 145 >> 146 G4DataVector activeZ; >> 147 >> 148 std::map<G4int,G4VEMDataSet*,std::less<G4int> > dataMap; >> 149 >> 150 std::vector<G4VEMDataSet*>* crossSections; 131 }; 151 }; 132 152 133 #endif 153 #endif 134 154 135 155 136 156 137 157 138 158 139 159 140 160 141 161 142 162 143 163 144 164 145 165