Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4PreCompoundTriton.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// -------------------------------------------------------------------
28//
29// GEANT4 Class file
30//
31//
32// File name: G4PreCompoundTriton
33//
34// Author: V.Lara
35//
36// Modified:
37// 21.08.2008 J. M. Quesada add choice of options
38// 20.08.2010 V.Ivanchenko added G4Pow and G4PreCompoundParameters pointers
39// use int Z and A and cleanup
40// 05.07.2013 J.M. Quesada FactorialFactor fixed
41//
42
44#include "G4Triton.hh"
45#include "G4CoulombBarrier.hh"
46
50
52{
53 return static_cast<G4double>(((N-3)*(P-2)*(N-2))*((P-1)*(N-1)*P))/12.0;
54}
55
57{
58 return 243.0/static_cast<G4double>(A*A);
59}
60
62{
63 G4double rj = 0.0;
64 if(nCharged >= 1 && (nParticles-nCharged) >= 2) {
65 G4double denominator = (G4double)(nParticles*(nParticles-1)*(nParticles-2));
66 rj = (3*nCharged*(nParticles-nCharged)*(nParticles-nCharged-1))
67 /denominator;
68 }
69 return rj;
70}
71
73{
74 G4double C = 0.0;
75 if (theFragZ <= 70)
76 {
77 C = 0.10;
78 }
79 else
80 {
81 C = ((((0.15417e-06*theFragZ) - 0.29875e-04)*theFragZ
82 + 0.21071e-02)*theFragZ - 0.66612e-01)*theFragZ + 0.98375;
83 }
84
85 return 1.0 + C/3.0;
86}
87
G4double C(G4double temp)
double G4double
Definition G4Types.hh:83
int G4int
Definition G4Types.hh:85
const G4double A[17]
G4double CoalescenceFactor(G4int A) const override
G4double FactorialFactor(G4int N, G4int P) const override
G4double GetAlpha() const override
G4double GetRj(G4int NumberParticles, G4int NumberCharged) const override
#define N
Definition crc32.c:57