Changeset 962 for trunk/source/processes/hadronic/cross_sections/include/G4BGGNucleonElasticXS.hh

Timestamp:

Apr 6, 2009, 12:30:29 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/cross_sections/include/G4BGGNucleonElasticXS.hh (modified) (5 diffs)

trunk/source/processes/hadronic/cross_sections/include/G4BGGNucleonElasticXS.hh

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4BGGNucleonElasticXS.hh,v 1.1 2007/03/13 15:19:30 vnivanch Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4BGGNucleonElasticXS.hh,v 1.2 2008/12/01 16:50:23 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -60 +60 @@
 class G4GlauberGribovCrossSection;
 class G4NucleonNuclearCrossSection;
+class G4HadronNucleonXsc;
 
 class G4BGGNucleonElasticXS : public G4VCrossSectionDataSet
@@ -93 +94 @@
   void Initialise();
 
-  G4double thEnergy;  // threshold of Glauber model
-  G4double theFac[93];
+  G4double CoulombFactor(G4double kinEnergy, G4double A);
+
+  G4double fGlauberEnergy;  
+  G4double fLowEnergy;  
+  G4double theGlauberFac[93];
+  G4double theCoulombFac[93];
 
   const G4ParticleDefinition*     particle;
   G4GlauberGribovCrossSection*    fGlauber;
   G4NucleonNuclearCrossSection*   fNucleon;
+  G4HadronNucleonXsc*             fHadron;
+  G4bool                          isProton;
+  G4bool                          isInitialized;
 };
 
@@ -104 +112 @@
 
 inline
-G4bool G4BGGNucleonElasticXS::IsApplicable(const G4DynamicParticle* dp, 
-                                           const G4Element*  elm)
+G4bool G4BGGNucleonElasticXS::IsApplicable(const G4DynamicParticle*, 
+                                           const G4Element*)
 {
-  return IsZAApplicable(dp, elm->GetZ(), elm->GetN());
+  return true;
+  //  return IsZAApplicable(dp, elm->GetZ(), elm->GetN());
 }
 
@@ -113 +122 @@
 
 inline
-G4bool G4BGGNucleonElasticXS::IsZAApplicable(const G4DynamicParticle* dp, 
-                                             G4double Z, G4double/* A*/)
+G4bool G4BGGNucleonElasticXS::IsZAApplicable(const G4DynamicParticle*, 
+                                             G4double /*Z*/, G4double/* A*/)
 {
-  return (dp->GetDefinition() == particle && Z > 1.5);
+  return true;
+  //  return (dp->GetDefinition() == particle);
 }