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.6 2001/07/17 15:54:39 verderi Exp $ >> 25 // GEANT4 tag $Name: geant4-04-00 $ >> 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 "g4std/vector" >> 55 >> 56 // Class Description: >> 57 // A one-to-one mapping from Photon Momentum to some optical property. >> 58 // Class Description - End: >> 59 >> 60 ///////////////////// >> 61 // Class Definition >> 62 ///////////////////// >> 63 >> 64 class G4MaterialPropertyVector { >> 65 >> 66 public: // Without description >> 67 >> 68 ////////////// >> 69 // Operators >> 70 ////////////// >> 71 >> 72 G4bool operator ++(); >> 73 G4MaterialPropertyVector& >> 74 operator =(const G4MaterialPropertyVector &right); >> 75 >> 76 ///////////////// >> 77 // Constructors >> 78 ///////////////// >> 79 >> 80 G4MaterialPropertyVector() : MPV(0) >> 81 { >> 82 CurrentEntry = -1; >> 83 NumEntries = 0; >> 84 }; >> 85 >> 86 public: // With description >> 87 >> 88 G4MaterialPropertyVector(G4double *PhotonMomenta, >> 89 G4double *PropertyValues, >> 90 G4int NumElements); >> 91 // Constructor of G4MaterialPropertyVector object. >> 92 >> 93 public: // Without description >> 94 >> 95 G4MaterialPropertyVector(const G4MaterialPropertyVector &right); >> 96 >> 97 /////////////// >> 98 // Destructor >> 99 /////////////// >> 100 >> 101 ~G4MaterialPropertyVector(); >> 102 >> 103 //////////// >> 104 // Methods >> 105 //////////// >> 106 >> 107 public: // With description >> 108 >> 109 void ResetIterator(); >> 110 >> 111 void AddElement(G4double aPhotonMomentum, >> 112 G4double aPropertyValue); >> 113 // Add a new element (pair of numbers) to the G4MaterialPropertyVector. >> 114 void RemoveElement(G4double aPhotonMomentum); >> 115 // Remove the element with given x-value. >> 116 >> 117 G4double GetProperty(G4double aPhotonMomentum) const; >> 118 // Returns the y-value for given x-value (with interpolation). >> 119 G4double GetPhotonMomentum(G4double aProperty) const; >> 120 // Returns the x-value for given y-value (with interpolation). >> 121 // NOTE: Assumes that the y-value is an increasing function of >> 122 // the x-value. Returns the x-value corresponding to the >> 123 // y-value passed in. If several x-values correspond to >> 124 // the y-value passed in, the method returns the first >> 125 // x-value in the vector that corresponds to that value. >> 126 // For use with G4MaterialPropertyVector iterator: return >> 127 // property (or Photon momentum) at current point of iterator. >> 128 >> 129 G4double GetProperty() const; >> 130 G4double GetPhotonMomentum() const; >> 131 >> 132 G4double GetMaxProperty() const; >> 133 G4double GetMinProperty() const; >> 134 G4double GetMaxPhotonMomentum() const; >> 135 G4double GetMinPhotonMomentum() const; >> 136 >> 137 ////////// >> 138 // Tests >> 139 ////////// >> 140 >> 141 void DumpVector(); >> 142 >> 143 private: >> 144 >> 145 ///////////////////// >> 146 // Helper Functions >> 147 ///////////////////// >> 148 >> 149 G4MPVEntry GetEntry(G4int i) const; >> 150 >> 151 void GetAdjacentBins(G4double aPhotonMomentum, >> 152 G4int *left,G4int *right) const; >> 153 >> 154 ///////////////////////// >> 155 // Private Data Members >> 156 ///////////////////////// >> 157 >> 158 //G4RWTPtrSortedVector<G4MPVEntry> MPV; >> 159 G4std::vector<G4MPVEntry*> MPV; >> 160 G4int NumEntries; >> 161 G4int CurrentEntry; >> 162 }; >> 163 >> 164 /////////////////// >> 165 // Inline methods >> 166 /////////////////// >> 167 >> 168 inline >> 169 G4double G4MaterialPropertyVector::GetMaxProperty() const >> 170 { >> 171 return MPV.back()->GetProperty(); >> 172 } >> 173 >> 174 inline >> 175 G4double G4MaterialPropertyVector::GetMinProperty() const >> 176 { >> 177 return MPV.front()->GetProperty(); >> 178 } >> 179 >> 180 inline >> 181 G4double G4MaterialPropertyVector::GetMaxPhotonMomentum() const >> 182 { >> 183 return MPV.back()->GetPhotonMomentum(); >> 184 } >> 185 >> 186 inline >> 187 G4double G4MaterialPropertyVector::GetMinPhotonMomentum() const >> 188 { >> 189 return MPV.front()->GetPhotonMomentum(); >> 190 } 49 191 50 #endif /* G4MaterialPropertyVector_h */ 192 #endif /* G4MaterialPropertyVector_h */ 51 193