Changeset 1340 for trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4PreCompoundNucleon.hh

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4PreCompoundNucleon.hh (modified) (2 diffs)

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

@@ -24 +24 @@
 // ********************************************************************
 //
-//J. M. Quesada (August 2008).  
-//Based  on previous work by V. Lara
+// $Id: G4PreCompoundNucleon.hh,v 1.7 2010/08/28 15:16:55 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-03-ref-09 $
 //
+// J. M. Quesada (August 2008).  
+// Based  on previous work by V. Lara
+//
+// Modified:
+// 20.08.2010 V.Ivanchenko added int Z and A and cleanup; added 
+//                        G4ParticleDefinition to constructor
 
 
@@ -36 +42 @@
 class G4PreCompoundNucleon : public G4PreCompoundFragment
 {
+public:  
+
+  G4PreCompoundNucleon(const G4ParticleDefinition*,
+                       G4VCoulombBarrier * aCoulombBarrier);
+  
+  virtual ~G4PreCompoundNucleon();
+    
+protected:
+
+  virtual G4double 
+  ProbabilityDistributionFunction(G4double eKin,
+                                  const G4Fragment& aFragment);
+
+  virtual G4double CrossSection(G4double ekin) = 0; 
+
+  virtual G4double 
+  GetRj(G4int NumberParticles, G4int NumberCharged) = 0; 
+
+  virtual G4double GetAlpha() = 0;
+
+  virtual G4double GetBeta() = 0;
+
+  inline G4double GetOpt0(G4double ekin);
+
 private:
+
   // default constructor
-  G4PreCompoundNucleon() {};
+  G4PreCompoundNucleon();
+  // operators
+  G4PreCompoundNucleon(const G4PreCompoundNucleon &right);
+  const G4PreCompoundNucleon& 
+  operator= (const G4PreCompoundNucleon &right);
+  G4int operator==(const G4PreCompoundNucleon &right) const;
+  G4int operator!=(const G4PreCompoundNucleon &right) const;    
 
-public:
-  
-  // copy constructor
-  G4PreCompoundNucleon(const G4PreCompoundNucleon &right): 
-    G4PreCompoundFragment(right) {}
+  G4double fact;
+};
 
-  // constructor  
-  G4PreCompoundNucleon(const G4double anA, 
-                       const G4double aZ, 
-                       G4VCoulombBarrier* aCoulombBarrier,                     
-                       const G4String & aName) :
-    G4PreCompoundFragment(anA,aZ,aCoulombBarrier,aName) {} 
-
-
-  virtual ~G4PreCompoundNucleon() {}
-
-  // operators  
-  const G4PreCompoundNucleon & 
-  operator=(const G4PreCompoundNucleon &right) 
-  {
-    if (&right != this) this->G4PreCompoundFragment::operator=(right);
-    return *this;
-  }
-  
-  G4bool operator==(const G4PreCompoundNucleon &right) const 
-  { 
-    return G4PreCompoundFragment::operator==(right);
-  }
-    
-  G4bool operator!=(const G4PreCompoundNucleon &right) const 
-  { 
-    return G4PreCompoundFragment::operator!=(right);
-  }
-    
-  virtual G4double ProbabilityDistributionFunction(const G4double eKin,
-                                                   const G4Fragment& aFragment);
-  
-  private:
-
-  G4bool IsItPossible(const G4Fragment&) ;   
-    
- protected:
-
-  virtual G4double CrossSection(const G4double ekin)=0; 
-
-  virtual G4double GetRj(const G4int NumberParticles, const G4int NumberCharged) = 0; 
-
- };
+// *********************** OPT=0 : Dostrovski's cross section  ***************
+inline G4double G4PreCompoundNucleon::GetOpt0(G4double K)
+{
+  G4double r0 = theParameters->Getr0()*ResidualA13();
+  // cross section is now given in mb (r0 is in mm) for the sake of consistency
+  //with the rest of the options
+  return 1.e+25*CLHEP::pi*r0*r0*ResidualA13()*GetAlpha()*(1.+GetBeta()/K);
+}
 
 #endif