include/Geant4/G4StatMFMicroPartition.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 //
0028 // Hadronic Process: Nuclear De-excitations
0029 // by V. Lara
0030
0031 #ifndef G4StatMFMicroPartition_h
0032 #define G4StatMFMicroPartition_h 1
0033
0034 #include  <vector>
0035
0036 #include "globals.hh"
0037 #include "G4StatMFParameters.hh"
0038 #include "G4StatMFChannel.hh"
0039
0040 class G4StatMFMicroPartition {
0041
0042 public:
0043     // Constructor
0044     G4StatMFMicroPartition(G4int A, G4int Z) :
0045     theA(A), theZ(Z), _Probability(0.0), _Temperature(0.0),
0046     _Entropy(0.0) {};
0047
0048
0049     // Destructor
0050     ~G4StatMFMicroPartition() {};
0051
0052
0053 private:
0054     // Default constructor
0055     G4StatMFMicroPartition() {};
0056
0057     // Copy constructor
0058     G4StatMFMicroPartition(const G4StatMFMicroPartition & right);
0059
0060     // operators
0061     G4StatMFMicroPartition & operator=(const G4StatMFMicroPartition & right);
0062 public:
0063     G4bool operator==(const G4StatMFMicroPartition & right) const;
0064     G4bool operator!=(const G4StatMFMicroPartition & right) const;
0065
0066 public:
0067
0068     // Gives fragments charges
0069     G4StatMFChannel * ChooseZ(G4int A0, G4int Z0, G4double MeanT);
0070
0071     G4double GetProbability(void)
0072     { return _Probability; }
0073
0074     void SetPartitionFragment(G4int anA)
0075     {
0076         _thePartition.push_back(anA);
0077         CoulombFreeEnergy(anA);
0078     }
0079
0080     void Normalize(G4double Normalization)
0081     { _Probability /= Normalization; }
0082
0083     G4double CalcPartitionProbability(G4double U,
0084                       G4double FreeInternalE0,
0085                       G4double SCompound);
0086
0087     G4double GetTemperature(void)
0088     {
0089         return _Temperature;
0090     }
0091
0092     G4double GetEntropy(void)
0093     {
0094         return _Entropy;
0095     }
0096
0097 private:
0098
0099     void CoulombFreeEnergy(G4int anA);
0100
0101     G4double CalcPartitionTemperature(G4double U,
0102                       G4double FreeInternalE0);
0103
0104     G4double GetPartitionEnergy(G4double T);
0105
0106     G4double GetCoulombEnergy(void);
0107
0108     G4double GetDegeneracyFactor(G4int A);
0109
0110     G4double InvLevelDensity(G4double Af)
0111     {
0112         // Calculate Inverse Density Level
0113         // Epsilon0*(1 + 3 /(Af - 1))
0114         if (Af < 1.5) return 0.0;
0115         else return G4StatMFParameters::GetEpsilon0()*(1.0+3.0/(Af - 1.0));
0116     }
0117
0118 private:
0119
0120     // A and Z of initial nucleus
0121     G4int theA;
0122     G4int theZ;
0123
0124     // Partition probability
0125     G4double _Probability;
0126
0127     // Partition temperature
0128     G4double _Temperature;
0129
0130     // Partition entropy
0131     G4double _Entropy;
0132
0133     // The partition itself
0134     std::vector<G4int> _thePartition;
0135
0136     std::vector<G4double> _theCoulombFreeEnergy;
0137
0138 };
0139
0140 #endif