| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 // Author: Luciano Pandola
28 //
29 // History:
30 // -----------
31 // 18 Mar 2010 L. Pandola 1st implementati
32 // 09 Mar 2012 L. Pandola Add public metho
33 // the absolute and
34 // sections indepen
35 //
36 // -------------------------------------------
37 //
38 // Class description:
39 // This class is a container for cross section
40 // calculated by Penelope models (ionisation,
41 // PhysicsTables/PhysicsVectors of
42 // a) the "hard quantities" (above the thresho
43 // 1-st order (= stopping XS), 2-nd order
44 // b) the "soft quantities" (below threshold),
45 // 1-st order (= stopping XS), 2-nd order
46 // c) total hard cross sections for individual
47 // vs. energy. Two versions are available, one
48 // (good for sampling) and one with absolute v
49 //
50 // The interface *always* uses energy and cros
51 // log(energy) and log(XS) are used.
52 //
53 // One instance per each cut-material couple s
54 // calling class.
55 //
56 // Public method to retrieve hard cross sectio
57 // total cross section and hard shell cross se
58 //
59 // Notice: all quantities stored here are *per
60 //
61 // -------------------------------------------
62
63 #ifndef G4PENELOPECROSSSECTION_HH
64 #define G4PENELOPECROSSSECTION_HH 1
65
66 #include "globals.hh"
67
68 class G4PhysicsTable;
69 class G4DataVector;
70
71 class G4PenelopeCrossSection
72 {
73
74 public:
75 //constructor: one has to give the number of
76 //(= dimension of the energy grid) and the n
77 //default).
78 explicit G4PenelopeCrossSection(size_t nOfEn
79 //
80 ~G4PenelopeCrossSection();
81
82 //! Returns total cross section at the given
83 G4double GetTotalCrossSection(G4double energ
84 //! Returns hard cross section at the given
85 G4double GetHardCrossSection(G4double energy
86 //! Returns the total stopping power due to
87 G4double GetSoftStoppingPower(G4double energ
88 //! Returns the hard cross section for the g
89 G4double GetShellCrossSection(size_t shellID
90 //! Returns the hard cross section for the g
91 G4double GetNormalizedShellCrossSection(size
92
93 size_t GetNumberOfShells() const {return fNu
94
95 //!
96 //! Public interface for the master thread
97 //!
98 void AddCrossSectionPoint(size_t binNumber,
99 G4double energy,G4
100 G4double XH2,
101 G4double XS0, G4double XS1, G4double
102 void AddShellCrossSectionPoint(size_t binNum
103 size_t shellID,G4double energy,
104 void NormalizeShellCrossSections();
105
106 G4PenelopeCrossSection & operator=(const G4P
107 G4PenelopeCrossSection(const G4PenelopeCross
108
109 private:
110 //all tables are log. XS vs. log E
111
112 //XS0, XS1, XS2 in Penelope nomenclature
113 G4PhysicsTable* fSoftCrossSections;
114
115 //XH0, XH1, XH2 in Penelope nomenclature
116 G4PhysicsTable* fHardCrossSections;
117
118 //XS for individual shells
119 G4PhysicsTable* fShellCrossSections;
120 G4PhysicsTable* fShellNormalizedCrossSection
121
122 size_t fNumberOfEnergyPoints;
123 size_t fNumberOfShells;
124 G4bool fIsNormalized;
125 };
126
127 #endif
128
129