Changeset 962 for trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4PreCompoundHe3.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/G4PreCompoundHe3.hh (modified) (3 diffs)

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

@@ -24 +24 @@
 // ********************************************************************
 //
+// by V. Lara
 //
-// $Id: G4PreCompoundHe3.hh,v 1.6 2007/10/01 10:42:00 ahoward Exp $
-// GEANT4 tag $Name:  $
-//
-// by V. Lara
+//J. M. Quesada (July 08) 
 
 #ifndef G4PreCompoundHe3_h
@@ -36 +34 @@
 #include "G4ReactionProduct.hh"
 #include "G4He3.hh"
-
 #include "G4He3CoulombBarrier.hh"
-
+#include "G4PreCompoundParameters.hh"
 
 class G4PreCompoundHe3 : public G4PreCompoundIon
@@ -66 +63 @@
 
 
-  G4ReactionProduct * GetReactionProduct() const
-  {
-    G4ReactionProduct * theReactionProduct =
-      new G4ReactionProduct(G4He3::He3Definition());
-    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
-  virtual G4double GetRj(const G4int NumberParticles, const G4int NumberCharged)
-  {
-    G4double rj = 1.0;
-    G4double denominator = NumberParticles*(NumberParticles-1)*(NumberParticles-2);
-    if(denominator != 0) rj = 3.0*static_cast<G4double>(NumberCharged*(NumberCharged-1)*(NumberParticles-NumberCharged))/static_cast<G4double>(denominator); //JMQ 23/8/07
+  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);
+
+  G4double GetOpt0(const G4double K);
+  G4double GetOpt12(const G4double K);
+  G4double GetOpt34(const G4double K);
+
+  G4double GetAlpha();
+  
+  G4double GetBeta();
+
+//data members
+
+      G4He3CoulombBarrier theHe3CoulombBarrier;
+        G4double ResidualA;
+      G4double ResidualZ; 
+      G4double theA;
+      G4double theZ;
+      G4double ResidualAthrd;
+      G4double FragmentA;
+      G4double FragmentAthrd;
+
+
+};
+#endif
 
 
 
-  virtual G4double GetAlpha()
-  {
-    G4double C = 0.0;
-    G4double aZ = GetZ() + GetRestZ();
-    if (aZ <= 30) 
-      {
-        C = 0.10;
-      }
-    else if (aZ <= 50) 
-      {
-        C = 0.1 + -((aZ-50.)/20.)*0.02;
-      } 
-    else if (aZ < 70) 
-      {
-        C = 0.08 + -((aZ-70.)/20.)*0.02;
-      }
-    else 
-      {
-        C = 0.06;
-      }
-    return 1.0 + C*(4.0/3.0);
-  }
 
-  virtual G4double GetBeta()
-  {
-    return -GetCoulombBarrier();
-  }
 
-  virtual G4double FactorialFactor(const G4double N, const G4double P)
-  {
-    return 
-      (N-3.0)*(P-2.0)*(
-                       (((N-2.0)*(P-1.0))/2.0) *(
-                                                 (((N-1.0)*P)/3.0) 
-                                                 )
-                       );
-  }
-
-  virtual G4double CoalescenceFactor(const G4double A)
-  {
-    return 243.0/(A*A);
-  }    
-private:
-
-  G4He3CoulombBarrier theHe3CoulombBarrier;
-
-};
-
-#endif
  
-