Changeset 1196 for trunk/source/processes/hadronic/models/management/include/G4HadronicInteraction.hh

Timestamp:

Nov 25, 2009, 5:13:58 PM (15 years ago)

Author:

garnier

Message:

update CVS release candidate geant4.9.3.01

File:

• trunk/source/processes/hadronic/models/management/include/G4HadronicInteraction.hh (modified) (9 diffs)

trunk/source/processes/hadronic/models/management/include/G4HadronicInteraction.hh

@@ -25 +25 @@
 //
 //
-// $Id: G4HadronicInteraction.hh,v 1.9 2009/01/24 11:56:27 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03-beta-cand-01 $
+// $Id: G4HadronicInteraction.hh,v 1.13 2009/10/02 17:18:33 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-03-cand-01 $
 //
 // Hadronic Interaction  abstract base class
@@ -41 +41 @@
 // 23-Jan-2009 V.Ivanchenko move constructor and destructor to the body
 //                          and reorder methods in the header 
+// 29-Jun-2009 V.Ivanchenko add SampleInvariantT method
+// 29-Aug-2009 V.Ivanchenko moved G4ReactionDynamics to G4InelasticInteraction,
+//                          add const pointers, and added recoilEnergyThreshold
+//                          member and accesors
 
 // Class Description
@@ -52 +56 @@
  
 #include "G4HadFinalState.hh"
-#include "G4ReactionDynamics.hh"
 #include "G4Material.hh"
 #include "G4Nucleus.hh"
 #include "G4Track.hh"
 #include "G4HadProjectile.hh"
+#include "G4ReactionDynamics.hh"
 
 class G4HadronicInteraction
 {
 public: // With description
-
     
   G4HadronicInteraction(const G4String& modelName = "HadronicModel");
@@ -69 +72 @@
   virtual G4HadFinalState *ApplyYourself(const G4HadProjectile &aTrack, 
                                          G4Nucleus & targetNucleus ) = 0;
-  // This is the interface to implement for final state production code.
+  // The interface to implement for final state production code.
+
+  virtual G4double SampleInvariantT(const G4ParticleDefinition* p, 
+                                    G4double plab,
+                                    G4int Z, G4int A);
+  // The interface to implement sampling of scattering or change exchange
    
   virtual G4bool IsApplicable(const G4HadProjectile &/*aTrack*/, 
                               G4Nucleus & /*targetNucleus*/)
-  {  return true;}
+  { return true;}
  
   inline G4double GetMinEnergy() const
   { return theMinEnergy; }
     
-  virtual G4double GetMinEnergy( const G4Material *aMaterial,
-                                 const G4Element *anElement ) const;
+  G4double GetMinEnergy( const G4Material *aMaterial,
+                         const G4Element *anElement ) const;
    
-  inline void SetMinEnergy( const G4double anEnergy )
+  inline void SetMinEnergy( G4double anEnergy )
   { theMinEnergy = anEnergy; }
     
-  virtual void SetMinEnergy( G4double anEnergy, G4Element *anElement );
-    
-  virtual void SetMinEnergy( G4double anEnergy, G4Material *aMaterial );
+  void SetMinEnergy( G4double anEnergy, const G4Element *anElement );
+    
+  void SetMinEnergy( G4double anEnergy, const G4Material *aMaterial );
     
   inline G4double GetMaxEnergy() const
   { return theMaxEnergy; }
     
-  virtual G4double GetMaxEnergy( const G4Material *aMaterial,
-                                 const G4Element *anElement ) const;
+  G4double GetMaxEnergy( const G4Material *aMaterial,
+                         const G4Element *anElement ) const;
     
   inline void SetMaxEnergy( const G4double anEnergy )
   { theMaxEnergy = anEnergy; }
     
-  virtual void SetMaxEnergy( G4double anEnergy, G4Element *anElement );
-    
-  virtual void SetMaxEnergy( G4double anEnergy, G4Material *aMaterial );
+  void SetMaxEnergy( G4double anEnergy, const G4Element *anElement );
+    
+  void SetMaxEnergy( G4double anEnergy, const G4Material *aMaterial );
   
   inline const G4HadronicInteraction *GetMyPointer() const
@@ -113 +121 @@
   { return theModelName; }
     
-public: // Without description
-
-  virtual void DeActivateFor( G4Material *aMaterial );
-    
-  virtual void ActivateFor( G4Material *aMaterial ) 
+  void DeActivateFor( const G4Material *aMaterial );
+    
+  inline void ActivateFor( const G4Material *aMaterial ) 
   { 
     Block(); 
@@ -124 +130 @@
   }
 
-  virtual void DeActivateFor( G4Element *anElement ); 
-  virtual void ActivateFor( G4Element *anElement )
+  void DeActivateFor( const G4Element *anElement ); 
+  inline void ActivateFor( const G4Element *anElement )
   { 
     Block(); 
@@ -132 +138 @@
   }
 
-  virtual G4bool IsBlocked( const G4Material *aMaterial ) const;
-
-  virtual G4bool IsBlocked( const G4Element *anElement) const;
+  G4bool IsBlocked( const G4Material *aMaterial ) const;
+
+  G4bool IsBlocked( const G4Element *anElement) const;
+
+  inline void SetRecoilEnergyThreshold(G4double val) 
+  { recoilEnergyThreshold = val; }
+
+  G4double GetRecoilEnergyThreshold() const 
+  { return recoilEnergyThreshold;}
 
   inline G4bool operator==(const G4HadronicInteraction &right ) const
@@ -149 +161 @@
 
 protected:
+
+  inline void SetModelName(const G4String& nam) 
+  { theModelName = nam; }
 
   inline G4bool IsBlocked() const { return isBlocked;}
@@ -164 +179 @@
   // 2: more
   // (instead of verboseLevel as found in G4VProcess)
-    
-  G4ReactionDynamics theReactionDynamics;
-    
-  // these two have global validity
-  // units are assumed to be MeV
-    
+
+  // these two have global validity energy range    
   G4double theMinEnergy;
   G4double theMaxEnergy;
-    
+
   G4bool isBlocked;
 
+private:       
+    
+  G4double recoilEnergyThreshold;
+
   G4String theModelName;
     
-private:
-        
-   std::vector<std::pair<G4double, G4Material *> > theMinEnergyList;
-   std::vector<std::pair<G4double, G4Material *> > theMaxEnergyList;
-   std::vector<std::pair<G4double, G4Element *> > theMinEnergyListElements;
-   std::vector<std::pair<G4double, G4Element *> > theMaxEnergyListElements;
-   std::vector<G4Material *> theBlockedList;
-   std::vector<G4Element *> theBlockedListElements;
+  std::vector<std::pair<G4double, const G4Material *> > theMinEnergyList;
+  std::vector<std::pair<G4double, const G4Material *> > theMaxEnergyList;
+  std::vector<std::pair<G4double, const G4Element *> > theMinEnergyListElements;
+  std::vector<std::pair<G4double, const G4Element *> > theMaxEnergyListElements;
+  std::vector<const G4Material *> theBlockedList;
+  std::vector<const G4Element *> theBlockedListElements;
 };