- Timestamp:
- Nov 5, 2010, 3:45:55 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/cross_sections/include/G4BGGPionElasticXS.hh
r1337 r1340 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4BGGPionElasticXS.hh,v 1. 5 2009/11/19 11:46:13 vnivanchExp $27 // GEANT4 tag $Name: geant4-09-04-beta-01$26 // $Id: G4BGGPionElasticXS.hh,v 1.7 2010/10/20 09:07:38 stesting Exp $ 27 // GEANT4 tag $Name: hadr-cross-V09-03-12 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 55 55 #include "G4ParticleDefinition.hh" 56 56 #include "G4Element.hh" 57 #include "G4HadTmpUtil.hh" 57 58 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......59 59 60 60 class G4GlauberGribovCrossSection; … … 73 73 G4bool IsApplicable(const G4DynamicParticle*, const G4Element*); 74 74 75 // virtual 76 // G4bool IsZAApplicable(const G4DynamicParticle*, G4double Z, G4double A); 77 75 78 virtual 76 G4bool Is ZAApplicable(const G4DynamicParticle*, G4double Z, G4doubleA);79 G4bool IsIsoApplicable(const G4DynamicParticle*, G4int Z, G4int A); 77 80 78 81 virtual … … 80 83 const G4Element*, G4double aTemperature = 0.); 81 84 85 // virtual 86 // G4double GetIsoZACrossSection(const G4DynamicParticle*, G4double /*Z*/, 87 // G4double /*A*/, G4double aTemperature = 0.); 88 82 89 virtual 83 G4double Get IsoZACrossSection(const G4DynamicParticle*, G4double/*Z*/,84 G4double/*A*/, G4double aTemperature = 0.);90 G4double GetZandACrossSection(const G4DynamicParticle*, G4int /*Z*/, 91 G4int /*A*/, G4double aTemperature = 0.); 85 92 86 93 virtual … … 94 101 void Initialise(); 95 102 96 G4double CoulombFactor(G4double kinEnergy, G4double A); 103 // G4double CoulombFactor(G4double kinEnergy, G4double A); 104 G4double CoulombFactor(G4double kinEnergy, G4int A); 97 105 98 106 G4double fGlauberEnergy; … … 116 124 { 117 125 return true; 118 //return IsZAApplicable(dp, elm->GetZ(), elm->GetN());119 126 } 120 127 121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 128 //inline 129 //G4bool G4BGGPionElasticXS::IsZAApplicable(const G4DynamicParticle*, 130 // G4double /*Z*/, G4double/* A*/) 122 131 123 132 inline 124 G4bool G4BGGPionElasticXS::Is ZAApplicable(const G4DynamicParticle*,125 G4double /*Z*/, G4double/* A*/)133 G4bool G4BGGPionElasticXS::IsIsoApplicable(const G4DynamicParticle*, 134 G4int /*Z*/, G4int/* A*/) 126 135 { 127 136 return false; 128 // return (dp->GetDefinition() == particle);129 137 } 130 138 131 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......132 139 133 inline 140 inline 134 141 G4double G4BGGPionElasticXS::GetCrossSection(const G4DynamicParticle* dp, 135 142 const G4Element* elm, 136 143 G4double temp) 137 144 { 138 return GetIsoZACrossSection(dp, elm->GetZ(), elm->GetN(), temp); 145 // return GetIsoZACrossSection(dp, elm->GetZ(), elm->GetN(), temp); 146 G4int Z = G4lrint(elm->GetZ()); 147 G4int N = G4lrint(elm->GetN()); 148 return GetZandACrossSection(dp, Z, N, temp); 139 149 } 140 150 141 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......142 143 151 #endif
Note: See TracChangeset
for help on using the changeset viewer.