Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4HETCProton.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//
27// by V. Lara
28//
29// Modified:
30// 23.08.2010 V.Ivanchenko general cleanup, move constructor and destructor
31// the source, use G4Pow
32
33#include "G4HETCProton.hh"
34#include "G4Proton.hh"
35#include "G4CoulombBarrier.hh"
36
40
42{
43 G4double C = 0.0;
44 if (theResZ <= 70)
45 {
46 C = 0.10;
47 }
48 else
49 {
50 C = ((((0.15417e-06*theResZ) - 0.29875e-04)*theResZ
51 + 0.21071e-02)*theResZ - 0.66612e-01)*theResZ + 0.98375;
52 }
53 return 1.0 + C;
54}
55
57{
58 // 2s+1
59 return 2.0;
60}
61
62G4double G4HETCProton::K(const G4Fragment& aFragment) const
63{
64 // Number of protons in emitted fragment
65 G4int Pa = theZ;
66
68
69 G4int P = aFragment.GetNumberOfParticles();
70 G4int H = aFragment.GetNumberOfHoles();
71
72 G4double result = 0.0;
73 if (P > 0)
74 {
75 result = (H*r + Pa)/(P*r);
76 }
77
78 return std::max(0.0,result);
79}
G4double C(G4double temp)
double G4double
Definition G4Types.hh:83
int G4int
Definition G4Types.hh:85
G4int GetNumberOfParticles() const
G4int GetNumberOfHoles() const
G4double GetSpinFactor() const override
G4double GetAlpha() const override
G4double K(const G4Fragment &aFragment) const override