Geant4_10
G4StatMFMicroPartition.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 //
27 // $Id: G4StatMFMicroPartition.hh 67983 2013-03-13 10:42:03Z gcosmo $
28 //
29 // Hadronic Process: Nuclear De-excitations
30 // by V. Lara
31 
32 #ifndef G4StatMFMicroPartition_h
33 #define G4StatMFMicroPartition_h 1
34 
35 #include <vector>
36 
37 #include "globals.hh"
38 #include "G4StatMFParameters.hh"
39 #include "G4StatMFChannel.hh"
40 
42 
43 public:
44  // Constructor
46  theA(A), theZ(Z), _Probability(0.0), _Temperature(0.0),
47  _Entropy(0.0) {};
48 
49 
50  // Destructor
52 
53 
54 private:
55  // Default constructor
57 
58  // Copy constructor
60 
61  // operators
63 public:
66 
67 public:
68 
69  // Gives fragments charges
70  G4StatMFChannel * ChooseZ(G4int A0, G4int Z0, G4double MeanT);
71 
73  { return _Probability; }
74 
76  {
77  _thePartition.push_back(anA);
78  CoulombFreeEnergy(anA);
79  }
80 
81  void Normalize(G4double Normalization)
82  { _Probability /= Normalization; }
83 
85  G4double FreeInternalE0,
86  G4double SCompound);
87 
89  {
90  return _Temperature;
91  }
92 
94  {
95  return _Entropy;
96  }
97 
98 private:
99 
100  void CoulombFreeEnergy(G4int anA);
101 
102  G4double CalcPartitionTemperature(G4double U,
103  G4double FreeInternalE0);
104 
105  G4double GetPartitionEnergy(G4double T);
106 
107  G4double GetCoulombEnergy(void);
108 
109  G4double GetDegeneracyFactor(G4int A);
110 
111  G4double InvLevelDensity(G4double Af)
112  {
113  // Calculate Inverse Density Level
114  // Epsilon0*(1 + 3 /(Af - 1))
115  if (Af < 1.5) return 0.0;
116  else return G4StatMFParameters::GetEpsilon0()*(1.0+3.0/(Af - 1.0));
117  }
118 
119 private:
120 
121  // A and Z of initial nucleus
122  G4int theA;
123  G4int theZ;
124 
125  // Partition probability
126  G4double _Probability;
127 
128  // Partition temperature
129  G4double _Temperature;
130 
131  // Partition entropy
132  G4double _Entropy;
133 
134  // The partition itself
135  std::vector<G4int> _thePartition;
136 
137  std::vector<G4double> _theCoulombFreeEnergy;
138 
139 };
140 
141 #endif
G4bool operator!=(const G4StatMFMicroPartition &right) const
G4bool operator==(const G4StatMFMicroPartition &right) const
int G4int
Definition: G4Types.hh:78
G4StatMFChannel * ChooseZ(G4int A0, G4int Z0, G4double MeanT)
void Normalize(G4double Normalization)
G4StatMFMicroPartition(G4int A, G4int Z)
Float_t Z
Definition: plot.C:39
bool G4bool
Definition: G4Types.hh:79
G4double CalcPartitionProbability(G4double U, G4double FreeInternalE0, G4double SCompound)
static G4double GetEpsilon0()
double G4double
Definition: G4Types.hh:76
void SetPartitionFragment(G4int anA)