Changeset 962 for trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4PreCompoundDeuteron.hh

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/include/G4PreCompoundDeuteron.hh (modified) (3 diffs)

trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4PreCompoundDeuteron.hh

@@ -24 +24 @@
 // ********************************************************************
 //
+// by V. Lara
 //
-// $Id: G4PreCompoundDeuteron.hh,v 1.7 2007/10/01 10:41:59 ahoward Exp $
-// GEANT4 tag $Name:  $
-//
-// by V. Lara
+//J. M. Quesada (July 08) 
 
 #ifndef G4PreCompoundDeuteron_h
@@ -36 +34 @@
 #include "G4ReactionProduct.hh"
 #include "G4Deuteron.hh"
-
 #include "G4DeuteronCoulombBarrier.hh"
+#include "G4PreCompoundParameters.hh"
 
 
@@ -66 +64 @@
 
 
-  G4ReactionProduct * GetReactionProduct() const
-  {
-    G4ReactionProduct * theReactionProduct =
-      new G4ReactionProduct(G4Deuteron::DeuteronDefinition());
-    theReactionProduct->SetMomentum(GetMomentum().vect());
-    theReactionProduct->SetTotalEnergy(GetMomentum().e());
-#ifdef PRECOMPOUND_TEST
-    theReactionProduct->SetCreatorModel("G4PrecompoundModel");
-#endif
-    return theReactionProduct;
-  }   
+  G4ReactionProduct * GetReactionProduct() const;
+
   
 private:
 
-// added Rj method according to literature and JMQ - formula from Jose Quesada
-  virtual G4double GetRj(const G4int NumberParticles, const G4int NumberCharged)
-  {
-    G4double rj = 1.0;
-    G4double denominator = NumberParticles*(NumberParticles-1);
-    if(denominator != 0) rj = 2.0*static_cast<G4double>(NumberCharged*(NumberParticles-NumberCharged))/static_cast<G4double>(denominator); //JMQ re-corrected 23/8/07
-    //    return static_cast<G4double>(NumberCharged*(NumberCharged-1))/static_cast<G4double>(NumberParticles*(NumberParticles-1));
+  virtual G4double GetRj(const G4int NumberParticles, const G4int NumberCharged);
 
-    return rj;
-  }
+  virtual G4double CrossSection(const  G4double K) ; 
 
+  virtual G4double FactorialFactor(const G4double N, const G4double P);
 
+  virtual G4double CoalescenceFactor(const G4double A);
 
-  virtual G4double GetAlpha()
-  {
-    G4double C = 0.0;
-    G4double aZ = GetZ() + GetRestZ();
-    if (aZ >= 70) 
-      {
-        C = 0.10;
-      } 
-    else 
-      {
-        C = ((((0.15417e-06*aZ) - 0.29875e-04)*aZ + 0.21071e-02)*aZ - 0.66612e-01)*aZ + 0.98375; 
-      }
-    return 1.0 + C/2.0;
-  }
-    
-  virtual G4double GetBeta()
-  {
-    return -GetCoulombBarrier();
-  }
+  G4double GetOpt0(const G4double K);
+  G4double GetOpt12(const G4double K);
+  G4double GetOpt34(const G4double K);
 
-  virtual G4double FactorialFactor(const G4double N, const G4double P)
-  {
-    return (N-1.0)*(N-2.0)*(P-1.0)*P/2.0;
-  }
+  G4double GetAlpha();
   
-  virtual G4double CoalescenceFactor(const G4double A)
-  {
-    return 16.0/A;
-  }    
-private:
-  
-  G4DeuteronCoulombBarrier theDeuteronCoulombBarrier;
+  G4double GetBeta();
+
+//data members
+
+      G4DeuteronCoulombBarrier theDeuteronCoulombBarrier;
+      G4double ResidualA;
+      G4double ResidualZ; 
+      G4double theA;
+      G4double theZ;
+      G4double ResidualAthrd;
+      G4double FragmentA;
+      G4double FragmentAthrd;
 
 };
+#endif
 
-#endif
- 
-