ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PreCompoundModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4PreCompoundModel.hh
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 // by V. Lara
27 //
28 // Class Description
29 // Model implementation for pre-equilibrium decay models in geant4.
30 // To be used in your physics list, in case you neeed this kind of physics.
31 // Can be used as a stand-allone model, but also in conjunction with an intra-nuclear
32 // transport, or any of the string-parton models.
33 // Class Description - End
34 //
35 // Modified:
36 // 03.09.2008 J.M.Quesada added external choice of inverse
37 // cross section option.(default OPTxs=3)
38 // 06.09.2008 J.M.Quesada external choices have been added for:
39 // - superimposed Coulomb barrier (if useSICB=true, default false)
40 // - "never go back" hipothesis (if useNGB=true, default false)
41 // - soft cutoff from preeq. to equlibrium (if useSCO=true, default false)
42 // - CEM transition probabilities (if useCEMtr=true)
43 // 30.10.2009 J.M.Quesada CEM transition probabilities are set as default
44 // 20.08.2010 V.Ivanchenko Cleanup of the code
45 // 03.01.2012 V.Ivanchenko Added pointer to G4ExcitationHandler to the
46 // constructor
47 
48 #ifndef G4PreCompoundModel_h
49 #define G4PreCompoundModel_h 1
50 
51 #include "G4VPreCompoundModel.hh"
52 #include "G4Fragment.hh"
54 #include "G4ReactionProduct.hh"
55 #include "G4ExcitationHandler.hh"
56 
59 class G4NuclearLevelData;
61 
63 {
64 public:
65 
66  explicit G4PreCompoundModel(G4ExcitationHandler* ptr = nullptr);
67 
68  virtual ~G4PreCompoundModel();
69 
70  virtual G4HadFinalState * ApplyYourself(const G4HadProjectile & thePrimary,
71  G4Nucleus & theNucleus) final;
72 
73  virtual G4ReactionProductVector* DeExcite(G4Fragment& aFragment) final;
74 
75  virtual void BuildPhysicsTable(const G4ParticleDefinition&) final;
76 
77  virtual void InitialiseModel() final;
78 
79  virtual void ModelDescription(std::ostream& outFile) const final;
80  virtual void DeExciteModelDescription(std::ostream& outFile) const final;
81 
82  //====== obsolete Set methods =======
83  void UseHETCEmission();
84  void UseDefaultEmission();
85  void UseGNASHTransition();
86  void UseDefaultTransition();
87 
88  //for cross section selection
89  void SetOPTxs(G4int opt);
90 
91  //for the rest of external choices
92  void UseSICB();
93  void UseNGB();
94  void UseSCO();
95  void UseCEMtr();
96  //======================================
97 
98 private:
99 
100  inline
101  void PerformEquilibriumEmission(const G4Fragment & aFragment,
103 
104  void PrintWarning(const G4String& mname);
105 
106  G4PreCompoundModel(const G4PreCompoundModel &) = delete;
107  const G4PreCompoundModel& operator=(const G4PreCompoundModel &right) = delete;
108  G4bool operator==(const G4PreCompoundModel &right) const = delete;
109  G4bool operator!=(const G4PreCompoundModel &right) const = delete;
110 
111  //==============
112  // Data Members
113  //==============
114 
118 
121 
124 
125  //for the rest of external choices
129 
132 
134 
135 };
136 
137 inline void
139  G4ReactionProductVector * Result) const
140 {
141  G4ReactionProductVector* theEquilibriumResult =
142  GetExcitationHandler()->BreakItUp(aFragment);
143  Result->insert(Result->end(),theEquilibriumResult->begin(), theEquilibriumResult->end());
144  delete theEquilibriumResult;
145 }
146 
147 #endif
148