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
G4HadronPhysicsFTFP_BERT_TRV.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4HadronPhysicsFTFP_BERT_TRV.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
// * 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
// ClassName:
30
//
31
// Author: 2007 Gunter Folger
32
// created from G4HadronPhysicsFTFP
33
//
34
// Modified:
35
//
36
//----------------------------------------------------------------------------
37
//
38
#include <iomanip>
39
40
#include "
G4HadronPhysicsFTFP_BERT_TRV.hh
"
41
#include "
G4PiKBuilder.hh
"
42
#include "
G4FTFPPiKBuilder.hh
"
43
#include "
G4BertiniPiKBuilder.hh
"
44
#include "
G4CrossSectionInelastic.hh
"
45
#include "
G4HadronCaptureProcess.hh
"
46
#include "
G4NeutronRadCapture.hh
"
47
#include "
G4NeutronInelasticXS.hh
"
48
#include "
G4NeutronCaptureXS.hh
"
49
50
#include "
globals.hh
"
51
#include "
G4ios.hh
"
52
#include "
G4SystemOfUnits.hh
"
53
54
#include "
G4PhysListUtil.hh
"
55
#include "
G4Threading.hh
"
56
57
#include "
G4HadronicParameters.hh
"
58
59
// factory
60
#include "
G4PhysicsConstructorFactory.hh
"
61
//
62
G4_DECLARE_PHYSCONSTR_FACTORY
(
G4HadronPhysicsFTFP_BERT_TRV
);
63
64
G4HadronPhysicsFTFP_BERT_TRV::G4HadronPhysicsFTFP_BERT_TRV
(
G4int
) :
65
G4HadronPhysicsFTFP_BERT_TRV
(
"hInelastic FTFP_BERT_TRV"
,
false
)
66
{}
67
68
G4HadronPhysicsFTFP_BERT_TRV::G4HadronPhysicsFTFP_BERT_TRV
(
const
G4String
&
name
,
G4bool
quasiElastic)
69
:
G4HadronPhysicsFTFP_BERT
(name,quasiElastic)
70
{
71
minFTFP_pion
=
G4HadronicParameters::Instance
()->
GetMinEnergyTransitionFTF_Cascade
();
72
maxBERT_pion
=
G4HadronicParameters::Instance
()->
GetMaxEnergyTransitionFTF_Cascade
();
73
minFTFP_kaon
=
G4HadronicParameters::Instance
()->
GetMinEnergyTransitionFTF_Cascade
();
74
maxBERT_kaon
=
G4HadronicParameters::Instance
()->
GetMaxEnergyTransitionFTF_Cascade
();
75
minFTFP_proton
=
G4HadronicParameters::Instance
()->
GetMinEnergyTransitionFTF_Cascade
();
76
maxBERT_proton
=
G4HadronicParameters::Instance
()->
GetMaxEnergyTransitionFTF_Cascade
();
77
minFTFP_neutron
=
G4HadronicParameters::Instance
()->
GetMinEnergyTransitionFTF_Cascade
();
78
maxBERT_neutron
=
G4HadronicParameters::Instance
()->
GetMaxEnergyTransitionFTF_Cascade
();
79
QuasiElastic
=
false
;
80
}
81
82
/*
83
void G4HadronPhysicsFTFP_BERT_TRV::DumpBanner()
84
{
85
G4cout << " Revised FTFTP_BERT_TRV - new threshold between BERT and FTFP "
86
<< " is over the interval " << minFTFP_pion/GeV << " to " << maxBERT_pion/GeV
87
<< " GeV. " << G4endl;
88
G4cout << " -- quasiElastic was asked to be " << QuasiElastic
89
<< " and it is reset to " << false << G4endl;
90
}
91
*/
92
93
void
G4HadronPhysicsFTFP_BERT_TRV::Pion
()
94
{
95
auto
pik =
new
G4PiKBuilder
;
96
AddBuilder
(pik);
97
auto
ftfppik =
new
G4FTFPPiKBuilder
(
QuasiElastic
);
98
AddBuilder
(ftfppik);
99
ftfppik->SetMinEnergy(
minFTFP_pion
);
100
pik->RegisterMe(ftfppik);
101
auto
bertpik =
new
G4BertiniPiKBuilder
();
102
AddBuilder
(bertpik);
103
bertpik->SetMaxEnergy(
maxBERT_pion
);
104
pik->RegisterMe(bertpik);
105
pik->Build();
106
}
107
108
void
G4HadronPhysicsFTFP_BERT_TRV::Kaon
() {
109
//Use combined with pions
110
}
111
112
113
void
G4HadronPhysicsFTFP_BERT_TRV::ExtraConfiguration
()
114
{
115
//Modify Neutrons
116
const
G4ParticleDefinition
*
neutron
=
G4Neutron::Neutron
();
117
G4HadronicProcess
* inel =
G4PhysListUtil::FindInelasticProcess
(neutron);
118
if
(inel) { inel->
AddDataSet
(
new
G4NeutronInelasticXS
()); }
119
G4HadronicProcess
* capture =
G4PhysListUtil::FindCaptureProcess
(neutron);
120
if
(capture) {
121
capture->
RegisterMe
(
new
G4NeutronRadCapture
());
122
}
123
}
124
geant4
tree
geant4-10.6-release
source
physics_lists
constructors
hadron_inelastic
src
G4HadronPhysicsFTFP_BERT_TRV.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:27
using
1.8.2 with
ECCE GitHub integration