Changeset 1340 for trunk/source/processes/hadronic/cross_sections/src/G4UInelasticCrossSection.cc

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/hadronic/cross_sections/src/G4UInelasticCrossSection.cc (modified) (9 diffs)

trunk/source/processes/hadronic/cross_sections/src/G4UInelasticCrossSection.cc

@@ -78 +78 @@
 }
 
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-G4bool G4UInelasticCrossSection::IsZAApplicable(const G4DynamicParticle* dp, 
-                                                G4double Z, G4double A)
+
+G4bool
+G4UInelasticCrossSection::IsIsoApplicable(const G4DynamicParticle* dp, 
+                                          G4int Z, G4int A)
 {
   G4bool res = false;
@@ -94 +94 @@
 }
 
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-G4double G4UInelasticCrossSection::GetCrossSection(const G4DynamicParticle* dp, 
-                                                   const G4Element* elm, 
-                                                   G4double temp)
-{
-  return GetIsoZACrossSection(dp, elm->GetZ(), elm->GetN(), temp);
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-G4double G4UInelasticCrossSection::GetIsoZACrossSection(const G4DynamicParticle* dp, 
-                                                        G4double Z,
-                                                        G4double A, 
-                                                        G4double)
+
+G4double
+G4UInelasticCrossSection::GetCrossSection(const G4DynamicParticle* dp, 
+                                          const G4Element* elm, 
+                                          G4double temp)
+{
+  G4int Z = G4lrint(elm->GetZ());
+  G4int N = G4lrint(elm->GetN());
+  return GetIsoZACrossSection(dp, Z, N, temp);
+}
+
+
+G4double
+G4UInelasticCrossSection::GetZandACrossSection(const G4DynamicParticle* dp, 
+                                               G4int Z, G4int A, G4double)
 {
   G4double cross = 0.0;
   G4double ekin = dp->GetKineticEnergy();
-  G4int iz = G4int(Z);
-  if(iz > 92) iz = 92;
+  //  G4int iz = G4int(Z);
+  if(Z > 92) Z = 92;
 
     // proton and neutron
   if(fNucleon) { 
-    if(iz == 1) cross = fGheisha->GetInelasticCrossSection(dp, Z, A);
+    if(Z == 1) cross = fGheisha->GetInelasticCrossSection(dp, Z, A);
     else if(ekin > thEnergy) {
-      cross = theFac[iz]*fGlauber->GetInelasticGlauberGribov(dp, Z, A);
+      cross = theFac[Z]*fGlauber->GetInelasticGlauberGribov(dp, Z, A);
     } else {
       cross = fNucleon->GetIsoZACrossSection(dp, Z, A);
@@ -126 +126 @@
     // pions
   } else if(fUPi) {
-    if(iz == 1) cross = fGheisha->GetInelasticCrossSection(dp, Z, A);
+    if(Z == 1) cross = fGheisha->GetInelasticCrossSection(dp, Z, A);
     else if(ekin > thEnergy) {
-      cross = theFac[iz]*fGlauber->GetInelasticGlauberGribov(dp, Z, A);
+      cross = theFac[Z]*fGlauber->GetInelasticGlauberGribov(dp, Z, A);
     } else {
       cross = fUPi->GetInelasticCrossSection(dp, Z, A);
@@ -136 +136 @@
   } else {
     if(hasGlauber && ekin > thEnergy) {
-      cross = theFac[iz]*fGlauber->GetInelasticGlauberGribov(dp, Z, A);
+      cross = theFac[Z]*fGlauber->GetInelasticGlauberGribov(dp, Z, A);
     } else if(fGheisha->IsApplicable(dp, Z, A)){
       cross = fGheisha->GetInelasticCrossSection(dp, Z, A);
@@ -153 +153 @@
 }
 
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 void G4UInelasticCrossSection::BuildPhysicsTable(const G4ParticleDefinition& p)
@@ -173 +172 @@
 
 
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
 void G4UInelasticCrossSection::Initialise(const G4ParticleDefinition* p) 
 {
@@ -182 +179 @@
 
   G4NistManager* nist = G4NistManager::Instance();
-  G4double A = nist->GetAtomicMassAmu(2);
-
-  if(fGlauber->IsZAApplicable(&dp, 2.0, A)) {
+  G4int A = G4lrint(nist->GetAtomicMassAmu(2));
+
+  if(fGlauber->IsIsoApplicable(&dp, 2, A)) {
     hasGlauber = true;
 
@@ -193 +190 @@
 
       G4double Z = G4double(iz);
-      A = nist->GetAtomicMassAmu(iz);
-      csup = fGlauber->GetInelasticGlauberGribov(&dp, Z, A);
+      A = G4lrint(nist->GetAtomicMassAmu(iz));
+      csup = fGlauber->GetInelasticGlauberGribov(&dp, iz, A);
 
       // proton and neutron
       if(fNucleon) { 
-        csdn = fNucleon->GetIsoZACrossSection(&dp, Z, A);
+        csdn = fNucleon->GetIsoZACrossSection(&dp, iz, A);
 
         // pions
       } else if(fUPi) {
-        csdn = fUPi->GetInelasticCrossSection(&dp, Z, A);
+        csdn = fUPi->GetInelasticCrossSection(&dp, iz, A);
 
         // other
-      } else if(fGheisha->IsApplicable(&dp, Z, A)) {
-        csdn = fGheisha->GetInelasticCrossSection(&dp, Z, A);
+      } else if(fGheisha->IsApplicable(&dp, iz, A)) {
+        csdn = fGheisha->GetInelasticCrossSection(&dp, iz, A);
 
       } else {
@@ -218 +215 @@
 }
 
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-