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

Timestamp:

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

Author:

garnier

Message:

update processes

File:

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

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

@@ -50 +50 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
-G4UInelasticCrossSection::G4UInelasticCrossSection(const G4ParticleDefinition* p) 
-{
+G4UInelasticCrossSection::G4UInelasticCrossSection(const G4ParticleDefinition*) 
+{
+  verboseLevel = 0;
   hasGlauber = false;
-  thEnergy   = 100.*GeV;
+  thEnergy   = 90.*GeV;
   fGlauber   = new G4GlauberGribovCrossSection();
   fGheisha   = G4HadronCrossSections::Instance();
   fNucleon   = 0;
   fUPi       = 0;
-  if(p == G4Proton::Proton() || p == G4Neutron::Neutron())
-    fNucleon = new G4NucleonNuclearCrossSection();
-  else if(p == G4PionPlus::PionPlus() || p == G4PionMinus::PionMinus())
-    fUPi = new G4UPiNuclearCrossSection();
-  verboseLevel = 0;
-  Initialise(p);
 }
 
@@ -117 +112 @@
   G4double cross = 0.0;
   G4double ekin = dp->GetKineticEnergy();
-  G4int iz = G4int(Z + 0.5);
+  G4int iz = G4int(Z);
   if(iz > 92) iz = 92;
 
@@ -160 +155 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
-void G4UInelasticCrossSection::BuildPhysicsTable(const G4ParticleDefinition&)
-{
+void G4UInelasticCrossSection::BuildPhysicsTable(const G4ParticleDefinition& p)
+{
+  if(&p == G4Proton::Proton() || &p == G4Neutron::Neutron()) {
+    fNucleon = new G4NucleonNuclearCrossSection();
+  } else if(&p == G4PionPlus::PionPlus() || &p == G4PionMinus::PionMinus()) {
+    fUPi = new G4UPiNuclearCrossSection();
+  }
+  Initialise(&p);
 }