| Geant4 Cross Reference |
>> 1 // This code implementation is the intellectual property of
>> 2 // the GEANT4 collaboration.
1 // 3 //
2 // ******************************************* << 4 // By copying, distributing or modifying the Program (or any work
3 // * License and Disclaimer << 5 // based on the Program) you indicate your acceptance of this statement,
4 // * << 6 // and all its terms.
5 // * The Geant4 software is copyright of th << 7 //
6 // * the Geant4 Collaboration. It is provided << 8 // $Id: G4Isotope.hh,v 1.1.8.1.2.1 1999/12/07 20:49:16 gunter Exp $
7 // * conditions of the Geant4 Software License << 9 // GEANT4 tag $Name: geant4-01-00 $
8 // * LICENSE and available at http://cern.ch/ << 10 //
9 // * include a list of copyright holders. << 11 //
10 // * << 12 // ------------------------------------------------------------
11 // * Neither the authors of this software syst << 13 // GEANT 4 class header file
12 // * institutes,nor the agencies providing fin << 14 //
13 // * work make any representation or warran << 15 // For information related to this code contact:
14 // * regarding this software system or assum << 16 // CERN, CN Division, ASD group
15 // * use. Please see the license in the file << 17 //
16 // * for the full disclaimer and the limitatio << 18 // ----------------- class G4Isotope ------------------
17 // * << 19 //
18 // * This code implementation is the result << 20 // Torre Wenaus, November 1995
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 // class description <<
27 // 21 //
28 // An isotope is a chemical isotope defined by 22 // An isotope is a chemical isotope defined by its name,
29 // << 23 // Z (atomic number),
30 // << 24 // N (number of nucleons),
31 // << 25 // A (mass of a mole).
32 // <<
33 // If A is not defined it is taken from Geant4 <<
34 // 26 //
35 // The class contains as a private static memb 27 // The class contains as a private static member the table of defined
36 // isotopes (an ordered vector of isotopes). 28 // isotopes (an ordered vector of isotopes).
37 // 29 //
38 // Isotopes can be assembled into elements via 30 // Isotopes can be assembled into elements via the G4Element class.
39 // << 31
40 // 20.08.11: Add flag fm for isomer level (mma << 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
41 // 15.11.05: GetIsotope(isotopeName, G4bool wa << 33
42 // 31.03.05: A becomes optional. Taken from Ni << 34 // 17-01-97, aesthetic rearrangement, M.Maire
43 // 26.02.02: fIndexInTable renewed << 35 // 04-08-98, new method GetIsotope(isotopeName), M.Maire
44 // 14.09.01: fCountUse: nb of elements which u << 36
45 // 13.09.01: stl migration. Suppression of the << 37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
46 // 30.03.01: suppression of the warning messag <<
47 // 04.08.98: new method GetIsotope(isotopeName <<
48 // 17.01.97: aesthetic rearrangement (mma) <<
49 38
50 #ifndef G4ISOTOPE_HH 39 #ifndef G4ISOTOPE_HH
51 #define G4ISOTOPE_HH 40 #define G4ISOTOPE_HH
52 41
53 #include "G4ios.hh" 42 #include "G4ios.hh"
>> 43 #include "g4rw/tpordvec.h"
54 #include "globals.hh" 44 #include "globals.hh"
55 45
56 #include <vector> <<
57 <<
58 class G4Isotope; 46 class G4Isotope;
59 using G4IsotopeTable = std::vector<G4Isotope*> << 47 typedef G4RWTPtrOrderedVector<G4Isotope> G4IsotopeTable;
60 <<
61 class G4Isotope <<
62 { <<
63 public: // with description <<
64 // Make an isotope <<
65 G4Isotope(const G4String& name, // its name <<
66 G4int z, // atomic number <<
67 G4int n, // number of nucleons <<
68 G4double a = 0., // mass of mole <<
69 G4int mlevel = 0); // isomer level <<
70 <<
71 ~G4Isotope(); <<
72 <<
73 G4Isotope(const G4Isotope&) = delete; <<
74 G4Isotope& operator=(const G4Isotope&) = del <<
75 48
76 // Retrieval methods <<
77 const G4String& GetName() const { return fNa <<
78 49
79 // Atomic number << 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
80 G4int GetZ() const { return fZ; } <<
81 51
82 // Number of nucleous << 52 class G4Isotope
83 G4int GetN() const { return fN; } << 53 {
84 << 54 public:
85 // Atomic mass of mole in Geant4 units with <<
86 G4double GetA() const { return fA; } <<
87 <<
88 // Isomer level <<
89 G4int Getm() const { return fm; } <<
90 <<
91 static G4Isotope* GetIsotope(const G4String& <<
92 <<
93 static const G4IsotopeTable* GetIsotopeTable <<
94 <<
95 static std::size_t GetNumberOfIsotopes(); <<
96 <<
97 std::size_t GetIndex() const { return fIndex <<
98 <<
99 friend std::ostream& operator<<(std::ostream <<
100 <<
101 friend std::ostream& operator<<(std::ostream <<
102 <<
103 friend std::ostream& operator<<(std::ostream <<
104 55
105 G4bool operator==(const G4Isotope&) const; << 56 // Make an isotope
106 G4bool operator!=(const G4Isotope&) const; << 57 G4Isotope(const G4String& name, //its name
>> 58 G4int z, //atomic number
>> 59 G4int n, //number of nucleons
>> 60 G4double a); //mass of mole
>> 61
>> 62 ~G4Isotope();
>> 63
>> 64 // Retrieval methods
>> 65 G4String GetName() const {return fName;};
>> 66 G4int GetZ() const {return fZ;};
>> 67 G4int GetN() const {return fN;};
>> 68 G4double GetA() const {return fA;};
>> 69 size_t GetIndex() const {return fIndexInTable;};
>> 70
>> 71 static G4Isotope* GetIsotope(G4String name);
>> 72
>> 73 static
>> 74 const G4IsotopeTable* GetIsotopeTable() {return &theIsotopeTable;};
>> 75 static size_t GetNumberOfIsotopes() {return theIsotopeTable.length();};
>> 76
>> 77 friend
>> 78 ostream& operator<<(ostream&, G4Isotope*);
>> 79
>> 80 friend
>> 81 ostream& operator<<(ostream&, G4Isotope&);
>> 82
>> 83 friend
>> 84 ostream& operator<<(ostream&, G4IsotopeTable);
>> 85
>> 86 G4int operator==(const G4Isotope &right) const;
>> 87 G4int operator!=(const G4Isotope &right) const;
>> 88
>> 89 private:
107 90
108 void SetName(const G4String& name) { fName = << 91 G4Isotope(G4Isotope &right);
>> 92 G4Isotope & operator=(const G4Isotope &right);
109 93
110 private: 94 private:
111 G4String fName; // name of the Isotope <<
112 G4int fZ; // atomic number <<
113 G4int fN; // number of nucleons <<
114 G4double fA; // atomic mass of a mole <<
115 G4int fm; // isomer level <<
116 95
117 static G4IsotopeTable& GetIsotopeTableRef(); << 96 G4String fName; // name of the Isotope
118 << 97 G4int fZ; // atomic number
119 std::size_t fIndexInTable; // index in the << 98 G4int fN; // number of nucleons
>> 99 G4double fA; // mass of a mole
>> 100
>> 101 static
>> 102 G4IsotopeTable theIsotopeTable;
>> 103 size_t fIndexInTable; // index in the Isotope Table
120 }; 104 };
>> 105
>> 106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
>> 107
>> 108 inline
>> 109 G4Isotope* G4Isotope::GetIsotope(G4String isotopeName)
>> 110 {
>> 111 // search the isotope by its name
>> 112 for (G4int J=0 ; J<theIsotopeTable.length() ; J++)
>> 113 {
>> 114 if(theIsotopeTable[J]->GetName() == isotopeName)
>> 115 return theIsotopeTable[J];
>> 116 }
>> 117
>> 118 G4cerr << " Warning from GetIsotope(name). The isotope: " << isotopeName
>> 119 << " does not exist in the IsotopeTable. Return NULL pointer \n";
>> 120 return NULL;
>> 121 }
121 122
122 #endif 123 #endif
123 124