Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4GNASHTransitions.cc
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: G4GNASHTransitions.cc 96603 2016-04-25 13:29:51Z gcosmo $
27 //
28 // 20.08.2010 V.Ivanchenko move constructor and destructor to the source
29 
30 #include "G4GNASHTransitions.hh"
31 #include "G4PhysicalConstants.hh"
32 #include "G4SystemOfUnits.hh"
33 #include "G4NuclearLevelData.hh"
34 #include "G4DeexPrecoParameters.hh"
35 #include "G4HadronicException.hh"
36 #include "Randomize.hh"
37 
39 {
41 }
42 
44 {}
45 
47 CalculateProbability(const G4Fragment & aFragment)
48 {
49  static const G4double k = 135.0 *CLHEP::MeV*CLHEP::MeV*CLHEP::MeV;
50  G4double E = aFragment.GetExcitationEnergy();
51  G4double P = aFragment.GetNumberOfParticles();
52  G4double H = aFragment.GetNumberOfHoles();
53  G4double N = P + H;
54  G4double A = aFragment.GetA_asInt();
55 
56  G4double theMatrixElement(k*N/(A*A*A*E));
57  G4double x = E/(N*CLHEP::MeV);
58  static const G4double xf = std::sqrt(2.0/7.0);
59  if ( x < 2.0) { x *= xf; }
60  else if ( x < 7.0) { x *= std::sqrt(x/7.0); }
61  else if ( x > 15.0){ x *= std::sqrt(15.0/x); }
62  theMatrixElement *= x;
63 
64  G4double gg = (6.0/pi2)*theParameters->GetLevelDensity()*A;
65 
66  G4double Epauli = ((P+1.0)*(P+1.0) + (H+1.0)*(H+1.0) + (P+1.0) - 3.0*(H-1.0))*0.25;
67 
68  G4double Probability = gg*gg*gg *(E-Epauli)*(E-Epauli);
69  Probability *= theMatrixElement/(2.0*(N+1.0)*CLHEP::h_Planck);
70 
71  return Probability;
72 }
73 
75 {
76  result.SetNumberOfParticles(result.GetNumberOfParticles()+1);
77  result.SetNumberOfHoles(result.GetNumberOfHoles()+1);
78  if (G4UniformRand()*result.GetA_asInt() <= G4double(result.GetZ_asInt()))
79  {
80  result.SetNumberOfCharged(result.GetNumberOfCharged()+1);
81  }
82 
83  if (result.GetNumberOfParticles() < result.GetNumberOfCharged())
84  {
86  }
87 }
G4double G4ParticleHPJENDLHEData::G4double result
static constexpr double h_Planck
virtual void PerformTransition(G4Fragment &aFragment)
void SetNumberOfHoles(G4int valueTot, G4int valueP=0)
Definition: G4Fragment.hh:375
virtual G4double CalculateProbability(const G4Fragment &aFragment)
G4int GetNumberOfParticles() const
Definition: G4Fragment.hh:345
tuple x
Definition: test.py:50
G4double GetLevelDensity() const
static double P[]
G4int GetNumberOfHoles() const
Definition: G4Fragment.hh:365
#define G4UniformRand()
Definition: Randomize.hh:97
double A(double temperature)
G4int GetA_asInt() const
Definition: G4Fragment.hh:266
static constexpr double MeV
G4DeexPrecoParameters * GetParameters()
void SetNumberOfParticles(G4int value)
Definition: G4Fragment.hh:384
G4int GetZ_asInt() const
Definition: G4Fragment.hh:271
**D E S C R I P T I O N
Definition: HEPEvtcom.cc:77
void SetNumberOfCharged(G4int value)
Definition: G4Fragment.hh:389
double G4double
Definition: G4Types.hh:76
G4int GetNumberOfCharged() const
Definition: G4Fragment.hh:350
static G4NuclearLevelData * GetInstance()
G4double GetExcitationEnergy() const
Definition: G4Fragment.hh:283
static constexpr double pi2
Definition: G4SIunits.hh:78