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

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

@@ -24 +24 @@
 // ********************************************************************
 //
-//
-// $Id: G4PreCompoundModel.hh,v 1.7 2009/11/19 10:19:31 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4PreCompoundModel.hh,v 1.9 2010/08/20 07:41:48 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-03-ref-09 $
 //
 // by V. Lara
-
+//
 // Class Description
 // Model implementation for pre-equilibrium decay models in geant4. 
@@ -37 +36 @@
 // Class Description - End
 //
-// Modif (03 September 2008) by J. M. Quesada for external choice of inverse 
-// cross section option.(default OPTxs=3)
-// JMQ (06 September 2008) Also external choices have been added for:
+// Modified:
+// 03.09.2008 J.M.Quesada added external choice of inverse 
+//            cross section option.(default OPTxs=3)
+// 06.09.2008 J.M.Quesada external choices have been added for:
 //                - superimposed Coulomb barrier (if useSICB=true, default false) 
 //                - "never go back"  hipothesis (if useNGB=true, default false) 
 //                - soft cutoff from preeq. to equlibrium (if useSCO=true, default false)
 //                - CEM transition probabilities (if useCEMtr=true)
-// J. M. Quesada (30.10.09) : CEM transition probabilities are set as default
+// 30.10.2009 J.M.Quesada CEM transition probabilities are set as default
+// 20.08.2010 V.Ivanchenko Cleanup of the code - changed data members and inline methods
 
 #ifndef G4PreCompoundModel_h
@@ -50 +51 @@
 
 #include "G4VPreCompoundModel.hh"
-#include "G4LorentzVector.hh"
-
-
-#include "G4NucleiProperties.hh"
-#include "G4PreCompoundParameters.hh"
-#include "G4ExcitationHandler.hh"
 #include "G4Fragment.hh"
-#include "Randomize.hh"
-
-//#include "G4PreCompoundEmission.hh"
-
-#include "G4DynamicParticle.hh"
 #include "G4ReactionProductVector.hh"
 #include "G4ReactionProduct.hh"
-#include "G4ParticleTypes.hh"
-#include "G4ParticleTable.hh"
+#include "G4ExcitationHandler.hh"
 
-//#define debug
-//#define verbose
+class G4PreCompoundParameters;
+class G4PreCompoundEmission;
+class G4VPreCompoundTransitions;
+class G4ParticleDefinition;
 
 class G4PreCompoundModel : public G4VPreCompoundModel
-{
- 
+{ 
 public:
 
@@ -79 +69 @@
   virtual ~G4PreCompoundModel();
 
-private:
+  virtual G4HadFinalState * ApplyYourself(const G4HadProjectile & thePrimary, 
+                                          G4Nucleus & theNucleus);
 
-  G4PreCompoundModel();
+  virtual G4ReactionProductVector* DeExcite(G4Fragment& aFragment);
 
-  G4PreCompoundModel(const G4PreCompoundModel &);
+  void UseHETCEmission();
+  void UseDefaultEmission();
+  void UseGNASHTransition();
+  void UseDefaultTransition();
 
-  const G4PreCompoundModel& operator=(const G4PreCompoundModel &right);
-  G4bool operator==(const G4PreCompoundModel &right) const;
-  G4bool operator!=(const G4PreCompoundModel &right) const;
+  //for cross section selection
+  void SetOPTxs(G4int opt);
 
-public:
-
-  G4HadFinalState * ApplyYourself(const G4HadProjectile & thePrimary, G4Nucleus & theNucleus);
-
-  G4ReactionProductVector* DeExcite(const G4Fragment& aFragment) const;
+  //for the rest of external choices
+  void UseSICB();
+  void UseNGB();
+  void UseSCO();
+  void UseCEMtr();
 
 #ifdef PRECOMPOUND_TEST
@@ -102 +95 @@
 #endif
 
-  inline void UseHETCEmission() { useHETCEmission = true; }
-  inline void UseDefaultEmission() { useHETCEmission = false; }
-  inline void UseGNASHTransition() { useGNASHTransition = true; }
-  inline void UseDefaultTransition() { useGNASHTransition = false; }
-
- //for cross section selection
-  inline void SetOPTxs(G4int opt) { OPTxs = opt; }
-//for the rest of external choices
-  inline void UseSICB() { useSICB = true; }
-  inline void UseNGB()  { useNGB = true; }
-  inline void UseSCO()  { useSCO = true; }
-  inline void UseCEMtr() { useCEMtr = true; }
 private:  
 
+  inline 
   void PerformEquilibriumEmission(const G4Fragment & aFragment, 
                                   G4ReactionProductVector * theResult) const;
 
-private:
+  G4PreCompoundModel();
+  G4PreCompoundModel(const G4PreCompoundModel &);
+  const G4PreCompoundModel& operator=(const G4PreCompoundModel &right);
+  G4bool operator==(const G4PreCompoundModel &right) const;
+  G4bool operator!=(const G4PreCompoundModel &right) const;
 
 #ifdef debug                              
@@ -131 +117 @@
   //==============
 
+  G4PreCompoundParameters*   theParameters;
+  G4PreCompoundEmission*     theEmission;
+  G4VPreCompoundTransitions* theTransition;
 
+  const G4ParticleDefinition* proton;
+  const G4ParticleDefinition* neutron;
 
-  G4bool           useHETCEmission;
-  G4bool           useGNASHTransition;
+  G4bool useHETCEmission;
+  G4bool useGNASHTransition;
 
-//for cross section options
+  //for cross section options
   G4int OPTxs;
-//for the rest of external choices
+
+  //for the rest of external choices
   G4bool useSICB;
   G4bool useNGB;
@@ -144 +136 @@
   G4bool useCEMtr;
 
-
-    G4HadFinalState theResult;
+  G4HadFinalState theResult;
 
 #ifdef PRECOMPOUND_TEST
@@ -151 +142 @@
   static std::vector<G4String*> theCreatorModels;
 #endif
+};
 
-};
+inline void 
+G4PreCompoundModel::PerformEquilibriumEmission(const G4Fragment & aFragment,
+                                               G4ReactionProductVector * Result) const 
+{
+  G4ReactionProductVector* theEquilibriumResult = 
+    GetExcitationHandler()->BreakItUp(aFragment);
+  Result->insert(Result->end(),theEquilibriumResult->begin(), theEquilibriumResult->end());
+  delete theEquilibriumResult;
+}
+
 #endif