| 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$
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