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 // $Id: G4PenelopeCrossSection.hh 74626 2013-10-17 07:00:59Z gcosmo $ 26 // 27 // 27 // Author: Luciano Pandola 28 // Author: Luciano Pandola 28 // 29 // 29 // History: 30 // History: 30 // ----------- 31 // ----------- 31 // 18 Mar 2010 L. Pandola 1st implementati 32 // 18 Mar 2010 L. Pandola 1st implementation. 32 // 09 Mar 2012 L. Pandola Add public metho 33 // 09 Mar 2012 L. Pandola Add public method (and machinery) to return 33 // the absolute and 34 // the absolute and the normalized shell cross 34 // sections indepen 35 // sections independently. 35 // 36 // 36 // ------------------------------------------- 37 // ------------------------------------------------------------------- 37 // 38 // 38 // Class description: 39 // Class description: 39 // This class is a container for cross section 40 // This class is a container for cross sections and transport momenta 40 // calculated by Penelope models (ionisation, 41 // calculated by Penelope models (ionisation, bremsstrahlung). It stores 41 // PhysicsTables/PhysicsVectors of 42 // PhysicsTables/PhysicsVectors of 42 // a) the "hard quantities" (above the thresho 43 // a) the "hard quantities" (above the threshold), 0-th order (cross section) 43 // 1-st order (= stopping XS), 2-nd order 44 // 1-st order (= stopping XS), 2-nd order (= straggling XS) 44 // b) the "soft quantities" (below threshold), 45 // b) the "soft quantities" (below threshold), 0-th order (cross section) 45 // 1-st order (= stopping XS), 2-nd order 46 // 1-st order (= stopping XS), 2-nd order (= straggling XS) 46 // c) total hard cross sections for individual 47 // c) total hard cross sections for individual oscillators 47 // vs. energy. Two versions are available, one 48 // vs. energy. Two versions are available, one with normalized values 48 // (good for sampling) and one with absolute v 49 // (good for sampling) and one with absolute values. 49 // 50 // 50 // The interface *always* uses energy and cros 51 // The interface *always* uses energy and cross sections, while internally 51 // log(energy) and log(XS) are used. 52 // log(energy) and log(XS) are used. 52 // 53 // 53 // One instance per each cut-material couple s 54 // One instance per each cut-material couple should be created by the 54 // calling class. 55 // calling class. 55 // 56 // 56 // Public method to retrieve hard cross sectio 57 // Public method to retrieve hard cross section, soft stopping power, 57 // total cross section and hard shell cross se 58 // total cross section and hard shell cross sections. 58 // 59 // 59 // Notice: all quantities stored here are *per 60 // Notice: all quantities stored here are *per molecule* 60 // 61 // 61 // ------------------------------------------- 62 // ------------------------------------------------------------------- 62 63 63 #ifndef G4PENELOPECROSSSECTION_HH 64 #ifndef G4PENELOPECROSSSECTION_HH 64 #define G4PENELOPECROSSSECTION_HH 1 65 #define G4PENELOPECROSSSECTION_HH 1 65 66 66 #include "globals.hh" 67 #include "globals.hh" 67 68 68 class G4PhysicsTable; 69 class G4PhysicsTable; 69 class G4DataVector; 70 class G4DataVector; 70 71 71 class G4PenelopeCrossSection 72 class G4PenelopeCrossSection 72 { 73 { 73 74 74 public: 75 public: 75 //constructor: one has to give the number of 76 //constructor: one has to give the number of points in each PhysicsVector 76 //(= dimension of the energy grid) and the n 77 //(= dimension of the energy grid) and the number of shells (0 is the 77 //default). 78 //default). 78 explicit G4PenelopeCrossSection(size_t nOfEn << 79 G4PenelopeCrossSection(size_t nOfEnergyPoints,size_t nOfShells=0); 79 // 80 // 80 ~G4PenelopeCrossSection(); 81 ~G4PenelopeCrossSection(); 81 82 82 //! Returns total cross section at the given 83 //! Returns total cross section at the given energy 83 G4double GetTotalCrossSection(G4double energ 84 G4double GetTotalCrossSection(G4double energy) const; 84 //! Returns hard cross section at the given 85 //! Returns hard cross section at the given energy 85 G4double GetHardCrossSection(G4double energy 86 G4double GetHardCrossSection(G4double energy) const; 86 //! Returns the total stopping power due to 87 //! Returns the total stopping power due to soft collisions 87 G4double GetSoftStoppingPower(G4double energ 88 G4double GetSoftStoppingPower(G4double energy) const; 88 //! Returns the hard cross section for the g 89 //! Returns the hard cross section for the given shell (per molecule) 89 G4double GetShellCrossSection(size_t shellID 90 G4double GetShellCrossSection(size_t shellID,G4double energy) const; 90 //! Returns the hard cross section for the g 91 //! Returns the hard cross section for the given shell (normalized to 1) 91 G4double GetNormalizedShellCrossSection(size 92 G4double GetNormalizedShellCrossSection(size_t shellID,G4double energy) const; 92 93 93 size_t GetNumberOfShells() const {return fNu << 94 size_t GetNumberOfShells() const {return numberOfShells;}; 94 95 95 //! 96 //! 96 //! Public interface for the master thread 97 //! Public interface for the master thread 97 //! 98 //! 98 void AddCrossSectionPoint(size_t binNumber, 99 void AddCrossSectionPoint(size_t binNumber, 99 G4double energy,G4 100 G4double energy,G4double XH0, G4double XH1, 100 G4double XH2, 101 G4double XH2, 101 G4double XS0, G4double XS1, G4double 102 G4double XS0, G4double XS1, G4double XS2); 102 void AddShellCrossSectionPoint(size_t binNum 103 void AddShellCrossSectionPoint(size_t binNumber, 103 size_t shellID,G4double energy, 104 size_t shellID,G4double energy,G4double xs); 104 void NormalizeShellCrossSections(); 105 void NormalizeShellCrossSections(); 105 106 106 G4PenelopeCrossSection & operator=(const G4P << 107 private: 107 G4PenelopeCrossSection(const G4PenelopeCross << 108 G4PenelopeCrossSection & operator=(const G4PenelopeCrossSection &right); >> 109 G4PenelopeCrossSection(const G4PenelopeCrossSection&); >> 110 >> 111 G4bool isNormalized; >> 112 >> 113 size_t numberOfEnergyPoints; >> 114 size_t numberOfShells; 108 115 109 private: << 110 //all tables are log. XS vs. log E 116 //all tables are log. XS vs. log E 111 117 112 //XS0, XS1, XS2 in Penelope nomenclature 118 //XS0, XS1, XS2 in Penelope nomenclature 113 G4PhysicsTable* fSoftCrossSections; << 119 G4PhysicsTable* softCrossSections; 114 120 115 //XH0, XH1, XH2 in Penelope nomenclature 121 //XH0, XH1, XH2 in Penelope nomenclature 116 G4PhysicsTable* fHardCrossSections; << 122 G4PhysicsTable* hardCrossSections; 117 123 118 //XS for individual shells 124 //XS for individual shells 119 G4PhysicsTable* fShellCrossSections; << 125 G4PhysicsTable* shellCrossSections; 120 G4PhysicsTable* fShellNormalizedCrossSection << 126 G4PhysicsTable* shellNormalizedCrossSections; 121 127 122 size_t fNumberOfEnergyPoints; << 123 size_t fNumberOfShells; << 124 G4bool fIsNormalized; << 125 }; 128 }; 126 129 127 #endif 130 #endif 128 131 129 132