Geant4-11
Public Member Functions | Private Attributes
G4teoCrossSection Class Reference

#include <G4teoCrossSection.hh>

Inheritance diagram for G4teoCrossSection:
G4VhShellCrossSection

Public Member Functions

G4double CrossSection (G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass, const G4Material *mat) override
 
 G4teoCrossSection (const G4String &name)
 
 G4teoCrossSection (const G4teoCrossSection &)=delete
 
std::vector< G4double > GetCrossSection (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy=0, const G4Material *mat=nullptr) override
 
const G4String & GetName () const
 
G4teoCrossSection & operator= (const G4teoCrossSection &right)=delete
 
std::vector< G4double > Probabilities (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy=0, const G4Material *mat=nullptr) override
 
G4int SelectRandomShell (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
 
void SetTotalCS (G4double) override
 
virtual ~G4teoCrossSection ()
 

Private Attributes

G4VecpssrKModel * ecpssrShellK
 
G4VecpssrLiModel * ecpssrShellLi
 
G4VecpssrMiModel * ecpssrShellMi
 
G4String name
 
G4double totalCS
 

Detailed Description

Definition at line 54 of file G4teoCrossSection.hh.

Constructor & Destructor Documentation

◆ G4teoCrossSection() [1/2]

G4teoCrossSection::G4teoCrossSection ( const G4String &  name)
explicit

Definition at line 54 of file G4teoCrossSection.cc.

55 :G4VhShellCrossSection(nam),totalCS(0.0)
56{
57 ecpssrShellMi = nullptr;
58 if (nam == "ECPSSR_Analytical")
59 {
60 ecpssrShellK = new G4ecpssrBaseKxsModel();
61 ecpssrShellLi = new G4ecpssrBaseLixsModel();
62 }
63 else if (nam == "ECPSSR_FormFactor")
64 {
65 ecpssrShellK = new G4ecpssrFormFactorKxsModel();
66 ecpssrShellLi = new G4ecpssrFormFactorLixsModel();
67 ecpssrShellMi = new G4ecpssrFormFactorMixsModel();
68 }
69 else if (nam == "ECPSSR_ANSTO")
70 {
71 ecpssrShellK = new G4ANSTOecpssrKxsModel();
72 ecpssrShellLi = new G4ANSTOecpssrLixsModel();
73 ecpssrShellMi = new G4ANSTOecpssrMixsModel();
74 }
75 else
76 {
77 G4cout << "G4teoCrossSection::G4teoCrossSection: ERROR "
78 << " in cross section name ECPSSR_Analytical is used"
79 << G4endl;
80 ecpssrShellK = new G4ecpssrBaseKxsModel();
81 ecpssrShellLi = new G4ecpssrBaseLixsModel();
82 }
83}
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4VhShellCrossSection::G4VhShellCrossSection
G4VhShellCrossSection(const G4String &xname="")
Definition: G4VhShellCrossSection.cc:48
G4teoCrossSection::ecpssrShellLi
G4VecpssrLiModel * ecpssrShellLi
Definition: G4teoCrossSection.hh:84
G4teoCrossSection::ecpssrShellMi
G4VecpssrMiModel * ecpssrShellMi
Definition: G4teoCrossSection.hh:85
G4teoCrossSection::ecpssrShellK
G4VecpssrKModel * ecpssrShellK
Definition: G4teoCrossSection.hh:83

References ecpssrShellK, ecpssrShellLi, ecpssrShellMi, G4cout, and G4endl.

◆ ~G4teoCrossSection()

G4teoCrossSection::~G4teoCrossSection ( )
virtual

Definition at line 87 of file G4teoCrossSection.cc.

88{
89 delete ecpssrShellK;
90 delete ecpssrShellLi;
91 delete ecpssrShellMi;
92}

References ecpssrShellK, ecpssrShellLi, and ecpssrShellMi.

◆ G4teoCrossSection() [2/2]

G4teoCrossSection::G4teoCrossSection ( const G4teoCrossSection &  )
delete

Member Function Documentation

◆ CrossSection()

G4double G4teoCrossSection::CrossSection ( G4int  Z,
G4AtomicShellEnumerator  shell,
G4double  incidentEnergy,
G4double  mass,
const G4Material *  mat 
)
overridevirtual

Implements G4VhShellCrossSection.

Definition at line 122 of file G4teoCrossSection.cc.

126{
127 G4double res = 0.0;
128 if(shell > 3 && !ecpssrShellMi) {
129 return res;
130 }
131 else if(shell > 8) {
132 return res;
133 }
134 else if(fKShell == shell)
135 {
136 res = ecpssrShellK->CalculateCrossSection(Z, mass, incidentEnergy);
137 }
138 else if(fL1Shell == shell)
139 {
140 res = ecpssrShellLi->CalculateL1CrossSection(Z, mass, incidentEnergy);
141 }
142 else if(fL2Shell == shell)
143 {
144 res = ecpssrShellLi->CalculateL2CrossSection(Z, mass, incidentEnergy);
145 }
146 else if(fL3Shell == shell)
147 {
148 res = ecpssrShellLi->CalculateL3CrossSection(Z, mass, incidentEnergy);
149 }
150 else if(fM1Shell == shell)
151 {
152 res = ecpssrShellMi->CalculateM1CrossSection(Z, mass, incidentEnergy);
153 }
154 else if(fM2Shell == shell)
155 {
156 res = ecpssrShellMi->CalculateM2CrossSection(Z, mass, incidentEnergy);
157 }
158 else if(fM3Shell == shell)
159 {
160 res = ecpssrShellMi->CalculateM3CrossSection(Z, mass, incidentEnergy);
161 }
162 else if(fM4Shell == shell)
163 {
164 res = ecpssrShellMi->CalculateM4CrossSection(Z, mass, incidentEnergy);
165 }
166 else if(fM5Shell == shell)
167 {
168 res = ecpssrShellMi->CalculateM5CrossSection(Z, mass, incidentEnergy);
169 }
170 return res;
171}
G4double
double G4double
Definition: G4Types.hh:83
Z
const G4int Z[17]
Definition: G4WaterStopping.cc:51
G4VecpssrKModel::CalculateCrossSection
virtual G4double CalculateCrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrLiModel::CalculateL2CrossSection
virtual G4double CalculateL2CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrLiModel::CalculateL3CrossSection
virtual G4double CalculateL3CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrLiModel::CalculateL1CrossSection
virtual G4double CalculateL1CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrMiModel::CalculateM5CrossSection
virtual G4double CalculateM5CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrMiModel::CalculateM1CrossSection
virtual G4double CalculateM1CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrMiModel::CalculateM2CrossSection
virtual G4double CalculateM2CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrMiModel::CalculateM4CrossSection
virtual G4double CalculateM4CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0
G4VecpssrMiModel::CalculateM3CrossSection
virtual G4double CalculateM3CrossSection(G4int zTarget, G4double massIncident, G4double energyIncident)=0

References G4VecpssrKModel::CalculateCrossSection(), G4VecpssrLiModel::CalculateL1CrossSection(), G4VecpssrLiModel::CalculateL2CrossSection(), G4VecpssrLiModel::CalculateL3CrossSection(), G4VecpssrMiModel::CalculateM1CrossSection(), G4VecpssrMiModel::CalculateM2CrossSection(), G4VecpssrMiModel::CalculateM3CrossSection(), G4VecpssrMiModel::CalculateM4CrossSection(), G4VecpssrMiModel::CalculateM5CrossSection(), ecpssrShellK, ecpssrShellLi, ecpssrShellMi, fKShell, fL1Shell, fL2Shell, fL3Shell, fM1Shell, fM2Shell, fM3Shell, fM4Shell, fM5Shell, and Z.

◆ GetCrossSection()

std::vector< G4double > G4teoCrossSection::GetCrossSection ( G4int  Z,
G4double  incidentEnergy,
G4double  mass,
G4double  deltaEnergy = 0,
const G4Material *  mat = nullptr 
)
overridevirtual

Implements G4VhShellCrossSection.

Definition at line 96 of file G4teoCrossSection.cc.

101{
102 std::vector<G4double> crossSections;
103
104 crossSections.push_back( ecpssrShellK->CalculateCrossSection(Z, mass, incidentEnergy) );
105
106 crossSections.push_back( ecpssrShellLi->CalculateL1CrossSection(Z, mass, incidentEnergy) );
107 crossSections.push_back( ecpssrShellLi->CalculateL2CrossSection(Z, mass, incidentEnergy) );
108 crossSections.push_back( ecpssrShellLi->CalculateL3CrossSection(Z, mass, incidentEnergy) );
109
110 if (ecpssrShellMi) {
111 crossSections.push_back( ecpssrShellMi->CalculateM1CrossSection(Z, mass, incidentEnergy) );
112 crossSections.push_back( ecpssrShellMi->CalculateM2CrossSection(Z, mass, incidentEnergy) );
113 crossSections.push_back( ecpssrShellMi->CalculateM3CrossSection(Z, mass, incidentEnergy) );
114 crossSections.push_back( ecpssrShellMi->CalculateM4CrossSection(Z, mass, incidentEnergy) );
115 crossSections.push_back( ecpssrShellMi->CalculateM5CrossSection(Z, mass, incidentEnergy) );
116 }
117 return crossSections;
118}

References G4VecpssrKModel::CalculateCrossSection(), G4VecpssrLiModel::CalculateL1CrossSection(), G4VecpssrLiModel::CalculateL2CrossSection(), G4VecpssrLiModel::CalculateL3CrossSection(), G4VecpssrMiModel::CalculateM1CrossSection(), G4VecpssrMiModel::CalculateM2CrossSection(), G4VecpssrMiModel::CalculateM3CrossSection(), G4VecpssrMiModel::CalculateM4CrossSection(), G4VecpssrMiModel::CalculateM5CrossSection(), ecpssrShellK, ecpssrShellLi, ecpssrShellMi, and Z.

Referenced by Probabilities().

◆ GetName()

const G4String & G4VhShellCrossSection::GetName ( ) const
inlineinherited

Definition at line 102 of file G4VhShellCrossSection.hh.

103{
104 return name;
105}

References G4VhShellCrossSection::name.

Referenced by G4UAtomicDeexcitation::InitialiseForNewRun().

◆ operator=()

G4teoCrossSection & G4teoCrossSection::operator= ( const G4teoCrossSection &  right)
delete

◆ Probabilities()

std::vector< G4double > G4teoCrossSection::Probabilities ( G4int  Z,
G4double  incidentEnergy,
G4double  mass,
G4double  deltaEnergy = 0,
const G4Material *  mat = nullptr 
)
overridevirtual

Implements G4VhShellCrossSection.

Definition at line 175 of file G4teoCrossSection.cc.

180{
181 std::vector<G4double> crossSections =
182 GetCrossSection(Z, incidentEnergy, mass, deltaEnergy);
183
184 for (size_t i=0; i<crossSections.size(); ++i ) {
185 if (totalCS) {
186 crossSections[i] = crossSections[i]/totalCS;
187 }
188 }
189 return crossSections;
190}
G4teoCrossSection::GetCrossSection
std::vector< G4double > GetCrossSection(G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy=0, const G4Material *mat=nullptr) override
Definition: G4teoCrossSection.cc:96

References GetCrossSection(), totalCS, and Z.

◆ SelectRandomShell()

G4int G4VhShellCrossSection::SelectRandomShell ( G4int  Z,
G4double  incidentEnergy,
G4double  mass,
G4double  deltaEnergy,
const G4Material *  mat 
)
inherited

Definition at line 63 of file G4VhShellCrossSection.cc.

70{
71 std::vector<G4double> p =
72 Probabilities(Z,incidentEnergy,mass,deltaEnergy,mat);
73 G4int shell = -1;
74 size_t nShells = p.size();
75 G4double q = G4UniformRand();
76 for (size_t i=0; i<nShells; i++) {
77
78 if (p[i] >= q) {
79 shell = i;
80 break;
81 }
82 q -= p[i];
83 }
84 return shell;
85}
G4int
int G4int
Definition: G4Types.hh:85
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
G4VhShellCrossSection::Probabilities
virtual std::vector< G4double > Probabilities(G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)=0

References G4UniformRand, G4VhShellCrossSection::Probabilities(), and Z.

◆ SetTotalCS()

void G4teoCrossSection::SetTotalCS ( G4double  val)
overridevirtual

Reimplemented from G4VhShellCrossSection.

Definition at line 194 of file G4teoCrossSection.cc.

194 {
195 totalCS = val;
196}

References totalCS.

Field Documentation

◆ ecpssrShellK

G4VecpssrKModel* G4teoCrossSection::ecpssrShellK
private

Definition at line 83 of file G4teoCrossSection.hh.

Referenced by CrossSection(), G4teoCrossSection(), GetCrossSection(), and ~G4teoCrossSection().

◆ ecpssrShellLi

G4VecpssrLiModel* G4teoCrossSection::ecpssrShellLi
private

Definition at line 84 of file G4teoCrossSection.hh.

Referenced by CrossSection(), G4teoCrossSection(), GetCrossSection(), and ~G4teoCrossSection().

◆ ecpssrShellMi

G4VecpssrMiModel* G4teoCrossSection::ecpssrShellMi
private

Definition at line 85 of file G4teoCrossSection.hh.

Referenced by CrossSection(), G4teoCrossSection(), GetCrossSection(), and ~G4teoCrossSection().

◆ name

G4String G4VhShellCrossSection::name
privateinherited

Definition at line 98 of file G4VhShellCrossSection.hh.

Referenced by source.g4viscp.G4Scene::create_scene(), G4VhShellCrossSection::GetName(), mcscore.MCParticle::printout(), and source.g4viscp.G4Scene::update_scene().

◆ totalCS

G4double G4teoCrossSection::totalCS
private

Definition at line 87 of file G4teoCrossSection.hh.

Referenced by Probabilities(), and SetTotalCS().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3