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