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 // G4ProductionCuts << 27 // 26 // 28 // Class description: << 27 // $Id: G4ProductionCuts.hh 70369 2013-05-29 14:59:24Z gcosmo $ >> 28 // GEANT4 tag $Name: geant4-09-04-ref-00 $ 29 // 29 // 30 // A G4ProductionCuts object must be created a << 30 // 31 // value desired for a given geometrical regio << 31 // ------------------------------------------------------------ 32 // Production cyts are applicable for gamma, p << 32 // GEANT 4 class header file 33 << 33 // 34 // Author: H.Kurashige, 17 September 2002 - Fi << 34 // 35 // ------------------------------------------- << 35 // Class Description 36 #ifndef G4ProductionCuts_hh << 36 // This class is 37 #define G4ProductionCuts_hh 1 << 37 // >> 38 // ------------------------------------------------------------ >> 39 // First Implementation 17 Sep. 2002 H.Kurahige >> 40 // Add cuts for proton 28 Jul. 2009 H.Kurashige >> 41 // ------------------------------------------------------------ 38 42 39 #include <vector> << 43 #ifndef G4ProductionCuts_h >> 44 #define G4ProductionCuts_h 1 40 45 41 #include "globals.hh" 46 #include "globals.hh" 42 #include "G4ios.hh" 47 #include "G4ios.hh" >> 48 #include <vector> 43 #include "G4ParticleDefinition.hh" 49 #include "G4ParticleDefinition.hh" 44 50 45 enum G4ProductionCutsIndex 51 enum G4ProductionCutsIndex 46 { 52 { 47 idxG4GammaCut = 0, << 53 idxG4GammaCut =0, 48 idxG4ElectronCut, 54 idxG4ElectronCut, 49 idxG4PositronCut, 55 idxG4PositronCut, >> 56 50 idxG4ProtonCut, // for proton 57 idxG4ProtonCut, // for proton 51 58 52 NumberOfG4CutIndex 59 NumberOfG4CutIndex 53 }; 60 }; 54 61 55 class G4ProductionCuts 62 class G4ProductionCuts 56 { 63 { >> 64 public: // with description >> 65 // constructor >> 66 G4ProductionCuts(); >> 67 >> 68 // copy constructor >> 69 G4ProductionCuts(const G4ProductionCuts &right); >> 70 >> 71 G4ProductionCuts & operator=(const G4ProductionCuts &right); >> 72 >> 73 public: >> 74 // destructor >> 75 virtual ~G4ProductionCuts(); >> 76 >> 77 // equal opperators >> 78 G4int operator==(const G4ProductionCuts &right) const; >> 79 G4int operator!=(const G4ProductionCuts &right) const; >> 80 >> 81 public: // with description >> 82 // Set Cuts methods >> 83 void SetProductionCut(G4double cut, G4int index = -1); >> 84 void SetProductionCut(G4double cut, G4ParticleDefinition* ptcl); >> 85 void SetProductionCut(G4double cut, const G4String& pName); >> 86 // Set the productionCut in range with an index to particle type >> 87 // if index is omitted, the value is applied to all particles >> 88 >> 89 G4double GetProductionCut(G4int index) const; >> 90 // Get the productionCut in range with an index to particle type >> 91 >> 92 G4double GetProductionCut(const G4String& name) const; >> 93 // Get the productionCut in range with a name of particle type >> 94 >> 95 void SetProductionCuts(std::vector<G4double>&); >> 96 // Set the vector of production cuts in range for all particles >> 97 >> 98 const std::vector<G4double>& GetProductionCuts() const; >> 99 // Get the vector of production cuts in range for all particles >> 100 >> 101 G4bool IsModified() const; >> 102 // return true if any cut value has been modified >> 103 // after last calculation of PhysicsTable >> 104 >> 105 void PhysicsTableUpdated(); >> 106 // inform end of calculation of PhysicsTable to ProductionCut >> 107 57 public: 108 public: >> 109 static G4int GetIndex(const G4String& name); >> 110 static G4int GetIndex(const G4ParticleDefinition* ptcl); 58 111 59 G4ProductionCuts(); << 112 protected: 60 // Constructor << 113 std::vector<G4double> fRangeCuts; >> 114 G4bool isModified; 61 115 62 G4ProductionCuts(const G4ProductionCuts& r << 116 private: 63 // Copy constructor << 117 static G4ThreadLocal G4ParticleDefinition* gammaDef; >> 118 static G4ThreadLocal G4ParticleDefinition* electDef; >> 119 static G4ThreadLocal G4ParticleDefinition* positDef; 64 120 65 virtual ~G4ProductionCuts(); << 121 static G4ThreadLocal G4ParticleDefinition* protonDef; // for proton 66 // Destructor << 67 122 68 G4ProductionCuts& operator=(const G4Produc << 123 }; 69 // Assignment operator << 70 << 71 G4bool operator==(const G4ProductionCuts& << 72 G4bool operator!=(const G4ProductionCuts& << 73 // Equality operators << 74 << 75 void SetProductionCut(G4double cut, G4int << 76 void SetProductionCut(G4double cut, G4Part << 77 void SetProductionCut(G4double cut, const << 78 // Set the production cut in range for a << 79 // and index, G4ParticleDefinition* or p << 80 // particle is not a photon, e-, e+ or p << 81 // no effect. << 82 void SetProductionCut(G4double cut); << 83 // Set the production cut in range for p << 84 124 85 G4double GetProductionCut(G4int index) con << 86 // Get the production cut in range with << 87 125 88 G4double GetProductionCut(const G4String& << 126 inline 89 // Get the production cut in range with << 127 void G4ProductionCuts::SetProductionCut(G4double cut, G4int index) 90 << 128 { 91 void SetProductionCuts(std::vector<G4doubl << 129 if (index<0) { 92 // Set the vector of production cuts in << 130 for(G4int i = 0; i < NumberOfG4CutIndex; i++) { >> 131 fRangeCuts[i] = cut; >> 132 } >> 133 isModified = true; >> 134 >> 135 } else if (index < NumberOfG4CutIndex) { >> 136 fRangeCuts[index] = cut; >> 137 isModified = true; >> 138 } >> 139 } 93 140 94 const std::vector<G4double>& GetProduction << 141 inline 95 // Get the vector of production cuts in << 142 void G4ProductionCuts::SetProductionCut(G4double cut, G4ParticleDefinition* ptcl) >> 143 { >> 144 G4int idx = -1; >> 145 if(ptcl) idx = GetIndex(ptcl); >> 146 if(idx>=0) SetProductionCut(cut,idx); >> 147 } 96 148 97 G4bool IsModified() const; << 149 inline 98 // Return true if any cut value has been << 150 void G4ProductionCuts::SetProductionCut(G4double cut, const G4String& pName) 99 // after last calculation of Physics Tab << 151 { >> 152 G4int idx = GetIndex(pName); >> 153 if(idx>=0) SetProductionCut(cut,idx); >> 154 } 100 155 101 void PhysicsTableUpdated(); << 156 inline 102 // Inform end of calculation of Physics << 157 G4double G4ProductionCuts::GetProductionCut(G4int index) const 103 << 158 { 104 static G4int GetIndex(const G4String& name << 159 G4double cut=-1.0; 105 static G4int GetIndex(const G4ParticleDefi << 160 if ( (index>=0) && (index<NumberOfG4CutIndex) ) { >> 161 cut = fRangeCuts[index]; >> 162 } >> 163 return cut; >> 164 } 106 165 107 protected: << 166 inline >> 167 G4double G4ProductionCuts::GetProductionCut(const G4String& name) const >> 168 { >> 169 return GetProductionCut(GetIndex(name)); >> 170 } 108 171 109 std::vector<G4double> fRangeCuts; << 172 110 G4bool isModified = true; << 173 inline 111 }; << 174 const std::vector<G4double>& G4ProductionCuts::GetProductionCuts() const >> 175 { >> 176 return fRangeCuts; >> 177 } >> 178 >> 179 inline >> 180 G4bool G4ProductionCuts::IsModified() const >> 181 { >> 182 return isModified; >> 183 } >> 184 >> 185 inline >> 186 void G4ProductionCuts::PhysicsTableUpdated() >> 187 { >> 188 isModified = false; >> 189 } 112 190 113 #endif 191 #endif >> 192 >> 193 >> 194 >> 195 >> 196 >> 197 >> 198 >> 199 >> 200 >> 201 >> 202 114 203