| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 <<
26 ////////////////////////////////////////////// <<
27 // 22 //
28 // class G4MaterialPropertyVector <<
29 // 23 //
30 // Class description: << 24 // $Id: G4MaterialPropertyVector.hh,v 1.9 2003/06/16 16:56:20 gunter Exp $
>> 25 // GEANT4 tag $Name: geant4-07-00-cand-01 $
>> 26 //
>> 27 //
>> 28 ////////////////////////////////////////////////////////////////////////
>> 29 // G4MaterialPropertyVector Class Definition
>> 30 ////////////////////////////////////////////////////////////////////////
31 // 31 //
32 // A one-to-one mapping from Photon Energy to <<
33 <<
34 // File: G4MaterialPropertyVector.hh 32 // File: G4MaterialPropertyVector.hh
35 // 33 //
36 // Version: 2.0 << 34 // Description: A one-to-one mapping from Photon Momentum to some
>> 35 // optical property
>> 36 // Version: 1.0
37 // Created: 1996-02-08 37 // Created: 1996-02-08
38 // Author: Juliet Armstrong 38 // Author: Juliet Armstrong
39 // Updated: 2011-10-13 by Peter Gumplinger << 39 // Updated: 1999-10-29 add method and class descriptors
40 // remove the class: simply typed << 40 // 1997-03-25 by Peter Gumplinger
>> 41 // > value.h -> templates.hh
>> 42 // mail: gum@triumf.ca
41 // 43 //
42 ////////////////////////////////////////////// 44 ////////////////////////////////////////////////////////////////////////
43 45
44 #ifndef G4MaterialPropertyVector_h << 46 #ifndef G4MaterialPropertyVector_h
45 #define G4MaterialPropertyVector_h 1 47 #define G4MaterialPropertyVector_h 1
46 48
47 #include "G4PhysicsFreeVector.hh" << 49 /////////////
48 using G4MaterialPropertyVector = G4PhysicsFree << 50 // Includes
>> 51 /////////////
>> 52
>> 53 #include "G4MPVEntry.hh"
>> 54 #include <vector>
>> 55 #include <functional>
>> 56
>> 57 // Class Description:
>> 58 // A one-to-one mapping from Photon Momentum to some optical property.
>> 59 // Class Description - End:
>> 60
>> 61 /////////////////////
>> 62 // Class Definition
>> 63 /////////////////////
>> 64
>> 65 class G4MaterialPropertyVector {
>> 66
>> 67 struct MPVEntry_less
>> 68 : public std::binary_function<G4MPVEntry*, G4MPVEntry*, G4bool>
>> 69 {
>> 70 G4bool operator()(G4MPVEntry* x, G4MPVEntry* y) { return *x < *y; }
>> 71 };
>> 72
>> 73 public: // Without description
>> 74
>> 75 //////////////
>> 76 // Operators
>> 77 //////////////
>> 78
>> 79 G4bool operator ++();
>> 80 G4MaterialPropertyVector&
>> 81 operator =(const G4MaterialPropertyVector &right);
>> 82
>> 83 /////////////////
>> 84 // Constructors
>> 85 /////////////////
>> 86
>> 87 G4MaterialPropertyVector() : MPV(0)
>> 88 {
>> 89 CurrentEntry = -1;
>> 90 NumEntries = 0;
>> 91 };
>> 92
>> 93 public: // With description
>> 94
>> 95 G4MaterialPropertyVector(G4double *PhotonMomenta,
>> 96 G4double *PropertyValues,
>> 97 G4int NumElements);
>> 98 // Constructor of G4MaterialPropertyVector object.
>> 99
>> 100 public: // Without description
>> 101
>> 102 G4MaterialPropertyVector(const G4MaterialPropertyVector &right);
>> 103
>> 104 ///////////////
>> 105 // Destructor
>> 106 ///////////////
>> 107
>> 108 ~G4MaterialPropertyVector();
>> 109
>> 110 ////////////
>> 111 // Methods
>> 112 ////////////
>> 113
>> 114 public: // With description
>> 115
>> 116 void ResetIterator();
>> 117
>> 118 void AddElement(G4double aPhotonMomentum,
>> 119 G4double aPropertyValue);
>> 120 // Add a new element (pair of numbers) to the G4MaterialPropertyVector.
>> 121 void RemoveElement(G4double aPhotonMomentum);
>> 122 // Remove the element with given x-value.
>> 123
>> 124 G4double GetProperty(G4double aPhotonMomentum) const;
>> 125 // Returns the y-value for given x-value (with interpolation).
>> 126 G4double GetPhotonMomentum(G4double aProperty) const;
>> 127 // Returns the x-value for given y-value (with interpolation).
>> 128 // NOTE: Assumes that the y-value is an increasing function of
>> 129 // the x-value. Returns the x-value corresponding to the
>> 130 // y-value passed in. If several x-values correspond to
>> 131 // the y-value passed in, the method returns the first
>> 132 // x-value in the vector that corresponds to that value.
>> 133 // For use with G4MaterialPropertyVector iterator: return
>> 134 // property (or Photon momentum) at current point of iterator.
>> 135
>> 136 G4double GetProperty() const;
>> 137 G4double GetPhotonMomentum() const;
>> 138
>> 139 G4double GetMaxProperty() const;
>> 140 G4double GetMinProperty() const;
>> 141 G4double GetMaxPhotonMomentum() const;
>> 142 G4double GetMinPhotonMomentum() const;
>> 143
>> 144 //////////
>> 145 // Tests
>> 146 //////////
>> 147
>> 148 void DumpVector();
>> 149
>> 150 private:
>> 151
>> 152 /////////////////////
>> 153 // Helper Functions
>> 154 /////////////////////
>> 155
>> 156 G4MPVEntry GetEntry(G4int i) const;
>> 157
>> 158 void GetAdjacentBins(G4double aPhotonMomentum,
>> 159 G4int *left,G4int *right) const;
>> 160
>> 161 /////////////////////////
>> 162 // Private Data Members
>> 163 /////////////////////////
>> 164
>> 165 std::vector<G4MPVEntry*> MPV;
>> 166 G4int NumEntries;
>> 167 G4int CurrentEntry;
>> 168 };
>> 169
>> 170 ///////////////////
>> 171 // Inline methods
>> 172 ///////////////////
>> 173
>> 174 inline
>> 175 G4double G4MaterialPropertyVector::GetMaxProperty() const
>> 176 {
>> 177 return MPV.back()->GetProperty();
>> 178 }
>> 179
>> 180 inline
>> 181 G4double G4MaterialPropertyVector::GetMinProperty() const
>> 182 {
>> 183 return MPV.front()->GetProperty();
>> 184 }
>> 185
>> 186 inline
>> 187 G4double G4MaterialPropertyVector::GetMaxPhotonMomentum() const
>> 188 {
>> 189 return MPV.back()->GetPhotonMomentum();
>> 190 }
>> 191
>> 192 inline
>> 193 G4double G4MaterialPropertyVector::GetMinPhotonMomentum() const
>> 194 {
>> 195 return MPV.front()->GetPhotonMomentum();
>> 196 }
49 197
50 #endif /* G4MaterialPropertyVector_h */ 198 #endif /* G4MaterialPropertyVector_h */
51 199