Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4GEMChannelVI.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: G4GEMChannelVI.hh 98577 2016-07-25 13:05:12Z vnivanch $
27 //
28 // GEM de-excitation model
29 // by V. Ivanchenko (July 2016)
30 //
31 
32 #ifndef G4GEMChannelVI_h
33 #define G4GEMChannelVI_h 1
34 
35 #include "G4VEvaporationChannel.hh"
37 #include "G4NucleiProperties.hh"
38 #include "G4VCoulombBarrier.hh"
39 #include "G4Exp.hh"
40 
41 class G4Pow;
43 class G4VCoulombBarrier;
44 class G4LevelManager;
45 class G4NuclearLevelData;
46 
47 const G4int NPOINTSGEM = 10;
48 
50 {
51 public:
52 
53  explicit G4GEMChannelVI(G4int theA, G4int theZ);
54 
55  virtual ~G4GEMChannelVI();
56 
57  virtual void Initialise() final;
58 
59  virtual G4double GetEmissionProbability(G4Fragment* theNucleus) final;
60 
61  virtual G4Fragment* EmittedFragment(G4Fragment* theNucleus) final;
62 
63  virtual void Dump() const;
64 
65 private:
66 
67  G4double IntegratedProbability(G4double exc);
68 
69  G4double ProbabilityDistributionFunction(G4double exc, G4double resExc);
70 
71  G4double FindLevel(const G4LevelManager*,
72  G4double exc, G4double exclim);
73 
74  inline G4double I0(G4double t);
75  inline G4double I1(G4double t, G4double tx);
76  inline G4double I2(G4double s0, G4double sx);
77  G4double I3(G4double s0, G4double sx);
78 
79  G4GEMChannelVI(const G4GEMChannelVI & right) = delete;
80  const G4GEMChannelVI & operator=(const G4GEMChannelVI & right) = delete;
81  G4bool operator==(const G4GEMChannelVI & right) const = delete;
82  G4bool operator!=(const G4GEMChannelVI & right) const = delete;
83 
84  G4Pow* fG4pow;
85 
86  const G4VCoulombBarrier* cBarrier;
87 
88  const G4PairingCorrection* pairingCorrection;
89 
90  const G4LevelManager* levelManager;
91 
92  G4NuclearLevelData* nData;
93 
94  G4int A;
95  G4int Z;
96  G4int resA;
97  G4int resZ;
98  G4int fragA;
99  G4int fragZ;
100  G4int nWarn;
101 
102  G4double massGround;
103  G4double maxLevelE;
104  G4double Z13;
105  G4double A13;
106 
107  G4double massFrag;
108  G4double eCBarrier;
109  G4double resMassGround;
110  G4double maxKinEnergy;
111  G4double resZ13;
112  G4double resA13;
113  G4double delta0;
114  G4double delta1;
115 
116  G4double alphaP;
117  G4double betaP;
118 
119  G4double maxExc;
120  G4double maxProb;
121  G4double coeff;
122  G4double levelDensity;
123 
124  static const G4double ws[NPOINTSGEM];
125  static const G4double xs[NPOINTSGEM];
126 
127 };
128 
129 inline G4double G4GEMChannelVI::I0(G4double t)
130 {
131  return G4Exp(t) - 1.0;
132 }
133 
134 inline G4double G4GEMChannelVI::I1(G4double t, G4double tx)
135 {
136  return (t - tx + 1.0)*G4Exp(tx) - t - 1.0;
137 }
138 
139 
140 inline G4double G4GEMChannelVI::I2(G4double s0, G4double sx)
141 {
142  G4double S = 1.0/std::sqrt(s0);
143  G4double Sx = 1.0/std::sqrt(sx);
144 
145  G4double p1 = S*S*S*( 1.0 + S*S*( 1.5 + 3.75*S*S) );
146  G4double p2 = Sx*Sx*Sx*( 1.0 + Sx*Sx*( 1.5 + 3.75*Sx*Sx) )*G4Exp(sx-s0);
147 
148  return p1-p2;
149 }
150 
151 #endif
double S(double temp)
Definition: G4Pow.hh:56
G4GEMChannelVI(G4int theA, G4int theZ)
virtual G4Fragment * EmittedFragment(G4Fragment *theNucleus) final
virtual void Initialise() final
int G4int
Definition: G4Types.hh:78
bool G4bool
Definition: G4Types.hh:79
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:183
virtual ~G4GEMChannelVI()
virtual void Dump() const
double G4double
Definition: G4Types.hh:76
const G4int NPOINTSGEM
virtual G4double GetEmissionProbability(G4Fragment *theNucleus) final