Changeset 962 for trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/src/G4VPreCompoundNucleon.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/src/G4VPreCompoundNucleon.cc (modified) (1 diff)

trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/src/G4VPreCompoundNucleon.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4VPreCompoundNucleon.cc,v 1.5 2007/07/23 09:56:40 ahoward Exp $
-// GEANT4 tag $Name:  $
+
+// $Id: G4VPreCompoundNucleon.cc,v 1.9 2008/09/22 10:18:36 ahoward Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
+
 //
 // by V. Lara
-
-#include "G4VPreCompoundNucleon.hh"
-#include "G4PreCompoundParameters.hh"
-//#include "G4EvaporationLevelDensityParameter.hh"
+//
+//J.M. Quesada (Apr. 2008) DUMMY abstract base class for nucleons.
+// it is not ihherited by anything
+// Suppressed
 
 
-G4double G4VPreCompoundNucleon::
-ProbabilityDistributionFunction(const G4double eKin, 
-                                const G4Fragment& aFragment)
-{
-  if ( !IsItPossible(aFragment) ) return 0.0;
-
-  const G4double r0 = G4PreCompoundParameters::GetAddress()->Getr0();
-
-  G4double U = aFragment.GetExcitationEnergy();
-  G4double P = aFragment.GetNumberOfParticles();
-  G4double H = aFragment.GetNumberOfHoles();
-  G4double N = P + H;
-
-  // g = (6.0/pi2)*a*A
-  //  G4EvaporationLevelDensityParameter theLDP;
-  G4double g0 = (6.0/pi2)*aFragment.GetA() * 
-    G4PreCompoundParameters::GetAddress()->GetLevelDensity();
-  //    theLDP.LevelDensityParameter(G4int(aFragment.GetA()),G4int(aFragment.GetZ()),U);
-  G4double g1 = (6.0/pi2)*GetRestA() * 
-    G4PreCompoundParameters::GetAddress()->GetLevelDensity();
-    //    theLDP.LevelDensityParameter(G4int(GetRestA()),G4int(GetRestZ()),U);
-
-
-  G4double A0 = ((P*P+H*H+P-H)/4.0 - H/2.0)/g0;
-  G4double A1 = (A0*g0 - P/2.0)/g1;
-
-  G4double E0 = std::max(0.0,U - A0);
-  if (E0 == 0.0) return 0.0;
-  G4double E1 = std::max(0.0,U - eKin - GetBindingEnergy() - A1);
-
-  G4double Probability = 2.0/(hbarc*hbarc*hbarc) * GetReducedMass() * 
-      GetRj(aFragment.GetNumberOfParticles(), aFragment.GetNumberOfCharged()) * r0 * r0 * std::pow(GetRestA(),2.0/3.0) * GetAlpha() * (eKin + GetBeta()) *
-      P*(N-1.0) * std::pow(g1*E1/(g0*E0),N-2.0)/E0 *
-      g1/(g0*g0);
-
-  return Probability;
-}
-
-