Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PreCompoundModel.hh
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4PreCompoundModel.hh 96527 2016-04-20 08:51:00Z gcosmo $
27 //
28 // by V. Lara
29 //
30 // Class Description
31 // Model implementation for pre-equilibrium decay models in geant4.
32 // To be used in your physics list, in case you neeed this kind of physics.
33 // Can be used as a stand-allone model, but also in conjunction with an intra-nuclear
34 // transport, or any of the string-parton models.
35 // Class Description - End
36 //
37 // Modified:
38 // 03.09.2008 J.M.Quesada added external choice of inverse
39 // cross section option.(default OPTxs=3)
40 // 06.09.2008 J.M.Quesada external choices have been added for:
41 // - superimposed Coulomb barrier (if useSICB=true, default false)
42 // - "never go back" hipothesis (if useNGB=true, default false)
43 // - soft cutoff from preeq. to equlibrium (if useSCO=true, default false)
44 // - CEM transition probabilities (if useCEMtr=true)
45 // 30.10.2009 J.M.Quesada CEM transition probabilities are set as default
46 // 20.08.2010 V.Ivanchenko Cleanup of the code - changed data members and inline methods
47 // 03.01.2012 V.Ivanchenko Added pointer to G4ExcitationHandler to the
48 // constructor
49 
50 #ifndef G4PreCompoundModel_h
51 #define G4PreCompoundModel_h 1
52 
53 #include "G4VPreCompoundModel.hh"
54 #include "G4Fragment.hh"
56 #include "G4ReactionProduct.hh"
57 #include "G4ExcitationHandler.hh"
58 
62 
64 {
65 public:
66 
67  explicit G4PreCompoundModel(G4ExcitationHandler* ptr = nullptr);
68 
69  virtual ~G4PreCompoundModel();
70 
71  virtual G4HadFinalState * ApplyYourself(const G4HadProjectile & thePrimary,
72  G4Nucleus & theNucleus) final;
73 
74  virtual G4ReactionProductVector* DeExcite(G4Fragment& aFragment) final;
75 
76  virtual void BuildPhysicsTable(const G4ParticleDefinition&) final;
77 
78  virtual void InitialiseModel() final;
79 
80  virtual void ModelDescription(std::ostream& outFile) const final;
81  virtual void DeExciteModelDescription(std::ostream& outFile) const final;
82 
83  //====== obsolete Set methods =======
84  void UseHETCEmission();
85  void UseDefaultEmission();
86  void UseGNASHTransition();
87  void UseDefaultTransition();
88 
89  //for cross section selection
90  void SetOPTxs(G4int opt);
91 
92  //for the rest of external choices
93  void UseSICB();
94  void UseNGB();
95  void UseSCO();
96  void UseCEMtr();
97  //======================================
98 
99 private:
100 
101  inline
102  void PerformEquilibriumEmission(const G4Fragment & aFragment,
103  G4ReactionProductVector * theResult) const;
104 
105  void PrintWarning(const G4String& mname);
106 
107  G4PreCompoundModel(const G4PreCompoundModel &) = delete;
108  const G4PreCompoundModel& operator=(const G4PreCompoundModel &right) = delete;
109  G4bool operator==(const G4PreCompoundModel &right) const = delete;
110  G4bool operator!=(const G4PreCompoundModel &right) const = delete;
111 
112  //==============
113  // Data Members
114  //==============
115 
116  G4PreCompoundEmission* theEmission;
117  G4VPreCompoundTransitions* theTransition;
118 
119  const G4ParticleDefinition* proton;
120  const G4ParticleDefinition* neutron;
121 
122  G4double fLevelDensity;
123  G4double fLimitEnergy;
124 
125  //for the rest of external choices
126  G4bool useSCO;
127  G4bool isInitialised;
128 
129  G4int minZ;
130  G4int minA;
131 
132  G4HadFinalState theResult;
133 
134 };
135 
136 inline void
137 G4PreCompoundModel::PerformEquilibriumEmission(const G4Fragment & aFragment,
138  G4ReactionProductVector * Result) const
139 {
140  G4ReactionProductVector* theEquilibriumResult =
141  GetExcitationHandler()->BreakItUp(aFragment);
142  Result->insert(Result->end(),theEquilibriumResult->begin(), theEquilibriumResult->end());
143  delete theEquilibriumResult;
144 }
145 
146 #endif
147 
virtual void ModelDescription(std::ostream &outFile) const final
virtual void DeExciteModelDescription(std::ostream &outFile) const final
int G4int
Definition: G4Types.hh:78
G4ReactionProductVector * BreakItUp(const G4Fragment &theInitialState)
virtual void BuildPhysicsTable(const G4ParticleDefinition &) final
std::vector< G4ReactionProduct * > G4ReactionProductVector
G4ExcitationHandler * GetExcitationHandler() const
bool G4bool
Definition: G4Types.hh:79
virtual void InitialiseModel() final
double G4double
Definition: G4Types.hh:76
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment) final
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &thePrimary, G4Nucleus &theNucleus) final
void SetOPTxs(G4int opt)
G4PreCompoundModel(G4ExcitationHandler *ptr=nullptr)