ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
G4EMDissociationSpectrum.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4EMDissociationSpectrum.cc
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
// * *
21
// * Parts of this code which have been developed by QinetiQ Ltd *
22
// * under contract to the European Space Agency (ESA) are the *
23
// * intellectual property of ESA. Rights to use, copy, modify and *
24
// * redistribute this software for general public use are granted *
25
// * in compliance with any licensing, distribution and development *
26
// * policy adopted by the Geant4 Collaboration. This code has been *
27
// * written by QinetiQ Ltd for the European Space Agency, under ESA *
28
// * contract 17191/03/NL/LvH (Aurora Programme). *
29
// * *
30
// * By using, copying, modifying or distributing the software (or *
31
// * any work based on the software) you agree to acknowledge its *
32
// * use in resulting scientific publications, and indicate your *
33
// * acceptance of all terms of the Geant4 Software license. *
34
// ********************************************************************
35
//
36
// %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
37
//
38
// MODULE: G4EMDissociationSpectrum.cc
39
//
40
// Version: B.1
41
// Date: 15/04/04
42
// Author: P R Truscott
43
// Organisation: QinetiQ Ltd, UK
44
// Customer: ESA/ESTEC, NOORDWIJK
45
// Contract: 17191/03/NL/LvH
46
//
47
// %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
48
//
49
// CHANGE HISTORY
50
// --------------
51
//
52
// 6 October 2003, P R Truscott, QinetiQ Ltd, UK
53
// Created.
54
//
55
// 15 March 2004, P R Truscott, QinetiQ Ltd, UK
56
// Beta release
57
//
58
// %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
60
//
61
#include "
G4EMDissociationSpectrum.hh
"
62
#include "
globals.hh
"
63
#include "
G4PhysicalConstants.hh
"
64
#include "
G4SystemOfUnits.hh
"
65
#include "
G4Pow.hh
"
66
68
//
69
G4EMDissociationSpectrum::G4EMDissociationSpectrum
()
70
{
71
bessel
=
new
G4Bessel
();
72
}
74
//
75
G4EMDissociationSpectrum::~G4EMDissociationSpectrum
()
76
{
77
delete
bessel
;
78
}
80
//
81
G4double
G4EMDissociationSpectrum::GetGeneralE1Spectrum
82
(
G4double
Eg,
G4double
b
,
G4double
bmin)
83
{
84
G4double
b2 = b*
b
;
85
G4double
gg = 1.0/std::sqrt(1.0-b2);
86
G4double
xi = Eg * bmin / gg / b /
hbarc
;
87
G4double
K0 = bessel->K0(xi);
88
G4double
K1 = bessel->K1(xi);
89
G4double
n
= 2.0 *
fine_structure_const
/
pi
/ b2 / Eg *
90
(xi*K0*K1 - xi*xi*b2/2.0*(K1*K1-K0*K0));
91
return
n
;
92
}
94
//
95
G4double
G4EMDissociationSpectrum::GetGeneralE2Spectrum
96
(
G4double
Eg,
G4double
b,
G4double
bmin)
97
{
98
G4double
b2 = b *
b
;
99
G4double
b4 = b2 * b2;
100
G4double
gg = 1.0/std::sqrt(1.0-b2);
101
G4double
xi = Eg * bmin / gg / b /
hbarc
;
102
G4double
K0 = bessel->K0(xi);
103
G4double
K1 = bessel->K1(xi);
104
G4double
n
= 2.0 *
fine_structure_const
/
pi
/ b4 / Eg *
105
(2.0*(1.0-b2)*K1*K1 + xi*
G4Pow::GetInstance
()->
powA
((2.0-b2),2.0)*K0*K1 -
106
xi*xi*b4/2.0*(K1*K1-K0*K0));
107
return
n
;
108
}
110
//
111
G4double
G4EMDissociationSpectrum::GetClosestApproach
112
(
const
G4double
AP,
const
G4double
ZP,
G4double
AT,
G4double
ZT,
G4double
b)
113
{
114
G4double
bsq = b *
b
;
115
G4double
gg = 1.0/std::sqrt(1-bsq);
116
G4double
AProot3 =
G4Pow::GetInstance
()->
powA
(AP,1.0/3.0);
117
G4double
ATroot3 =
G4Pow::GetInstance
()->
powA
(AT,1.0/3.0);
118
G4double
bc = 1.34 *
fermi
* (AProot3+ATroot3 - 0.75 *(1.0/AProot3+1.0/ATroot3));
119
// G4double a0 = ZP * ZT * classic_electr_radius/bsq;
120
G4double
a0
= ZP * ZT *
elm_coupling
/ (AT*AP*
amu_c2
/(AT+AP)) / bsq;
121
G4double
bmin = 1.25 * bc +
halfpi
*a0/gg;
122
return
bmin;
123
}
125
//
geant4
tree
geant4-10.6-release
source
processes
hadronic
cross_sections
src
G4EMDissociationSpectrum.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:37
using
1.8.2 with
ECCE GitHub integration