Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
G4BertiniEvaporationChannel Class Referenceabstract

#include <G4BertiniEvaporationChannel.hh>

+ Inheritance diagram for G4BertiniEvaporationChannel:

Public Member Functions

 G4BertiniEvaporationChannel ()
 
virtual ~G4BertiniEvaporationChannel ()
 
void setVerboseLevel (G4int verbose)
 
void setNucleusA (G4int inputA)
 
void setNucleusZ (G4int inputZ)
 
G4int getParticleA ()
 
G4int getParticleZ ()
 
void setExcitationEnergy (G4double inputE)
 
void setQ (G4double inputQ)
 
void setPairingCorrection (G4int isCorrection)
 
G4double getLevelDensityParameter ()
 
G4String getName ()
 
virtual G4double getProbability ()
 
virtual void setProbability (G4double newProb)
 
virtual void calculateProbability ()=0
 
virtual G4double qmFactor ()
 
virtual G4double getQ ()
 
virtual G4double getCoulomb ()
 
virtual G4double getThresh ()
 
virtual G4DynamicParticle * emit ()=0
 

Protected Member Functions

G4double Q (G4double a, G4double z)
 
G4double pairingEnergyProtons (G4int A)
 
G4double pairingEnergyNeutrons (G4int N)
 
G4double cameron (G4double a, G4double z)
 
G4double cameronShellCorrectionP (G4int Z)
 
G4double cameronShellCorrectionN (G4int N)
 
void isotropicCosines (G4double &, G4double &, G4double &)
 

Protected Attributes

G4String name
 
G4int verboseLevel
 
G4int nucleusA
 
G4int nucleusZ
 
G4int particleA
 
G4int particleZ
 
G4double exmass
 
G4double emissionProbability
 
G4double rho
 
G4double correction
 
G4double excitationEnergy
 
G4int spin
 

Detailed Description

Definition at line 38 of file G4BertiniEvaporationChannel.hh.

Constructor & Destructor Documentation

◆ G4BertiniEvaporationChannel()

G4BertiniEvaporationChannel::G4BertiniEvaporationChannel ( )

Definition at line 46 of file G4BertiniEvaporationChannel.cc.

47{
48 verboseLevel = 0;
49}

◆ ~G4BertiniEvaporationChannel()

G4BertiniEvaporationChannel::~G4BertiniEvaporationChannel ( )
virtual

Definition at line 52 of file G4BertiniEvaporationChannel.cc.

53{
54}

Member Function Documentation

◆ calculateProbability()

virtual void G4BertiniEvaporationChannel::calculateProbability ( )
pure virtual

Implemented in G4BEChargedChannel, and G4BENeutronChannel.

◆ cameron()

G4double G4BertiniEvaporationChannel::cameron ( G4double  a,
G4double  z 
)
protected

◆ cameronShellCorrectionN()

G4double G4BertiniEvaporationChannel::cameronShellCorrectionN ( G4int  N)
protected

Definition at line 353 of file G4BertiniEvaporationChannel.cc.

354{
355 // Gives the binding energy correction depending in N
356 // due to shell correction and pairing energies in MeV
357 //
358 // see Cameron, Canadian Journal of Physics,
359 // vol. 35, 1957, p.1022
360 G4double table[200] = {
361 -0.83199997E+01, -0.15899999E+02, -0.11509999E+02, -0.14309999E+02, -0.11570000E+02,
362 -0.15899999E+02, -0.13909999E+02, -0.16029999E+02, -0.12129999E+02, -0.13869999E+02,
363 -0.12249998E+02, -0.14399999E+02, -0.13069999E+02, -0.15799998E+02, -0.13809999E+02,
364 -0.14980000E+02, -0.12629999E+02, -0.13759999E+02, -0.11369999E+02, -0.12379998E+02,
365 -0.92299995E+01, -0.96499996E+01, -0.76399994E+01, -0.91699991E+01, -0.80499992E+01,
366 -0.97199993E+01, -0.88699989E+01, -0.10759999E+02, -0.86399994E+01, -0.88899994E+01,
367 -0.65999994E+01, -0.71299992E+01, -0.47699995E+01, -0.53299990E+01, -0.30599995E+01,
368 -0.37899990E+01, -0.17199993E+01, -0.27899990E+01, -0.92999995E+00, -0.21899996E+01,
369 -0.51999998E+00, -0.18999996E+01, -0.44999999E+00, -0.21999998E+01, -0.12199993E+01,
370 -0.30699997E+01, -0.24199991E+01, -0.43699999E+01, -0.39399996E+01, -0.60799999E+01,
371 -0.44899998E+01, -0.45000000E+01, -0.31399994E+01, -0.29299994E+01, -0.10399990E+01,
372 -0.13599997E+01, 0.69000000E+00, 0.20999998E+00, 0.21099997E+01, 0.13299990E+01,
373 0.32900000E+01, 0.24599991E+01, 0.42999992E+01, 0.33199997E+01, 0.47899990E+01,
374 0.36199999E+01, 0.49699993E+01, 0.36399994E+01, 0.46299992E+01, 0.30699997E+01,
375 0.40599995E+01, 0.24899998E+01, 0.32999992E+01, 0.14599991E+01, 0.20599995E+01,
376 0.50999999E+00, 0.73999995E+00, -0.11799994E+01, -0.12599993E+01, -0.35399990E+01,
377 -0.39699993E+01, -0.52599993E+01, -0.41799994E+01, -0.37099991E+01, -0.20999994E+01,
378 -0.16999998E+01, -0.79999983E-01, -0.17999995E+00, 0.94000000E+00, 0.26999998E+00,
379 0.11299992E+01, 0.79999983E-01, 0.90999997E+00, -0.30999994E+00, 0.48999995E+00,
380 -0.77999997E+00, 0.79999983E-01, -0.11499996E+01, -0.22999996E+00, -0.14099998E+01,
381 -0.41999996E+00, -0.15499992E+01, -0.54999995E+00, -0.16599998E+01, -0.65999997E+00,
382 -0.17299995E+01, -0.75000000E+00, -0.17399998E+01, -0.77999997E+00, -0.16899996E+01,
383 -0.77999997E+00, -0.15999994E+01, -0.75000000E+00, -0.14599991E+01, -0.66999996E+00,
384 -0.12599993E+01, -0.50999999E+00, -0.10399990E+01, -0.52999997E+00, -0.18399992E+01,
385 -0.24199991E+01, -0.45199995E+01, -0.47599993E+01, -0.63299990E+01, -0.67599993E+01,
386 -0.78099995E+01, -0.57999992E+01, -0.53699999E+01, -0.36299992E+01, -0.33499994E+01,
387 -0.17500000E+01, -0.18799992E+01, -0.60999995E+00, -0.89999998E+00, 0.89999974E-01,
388 -0.31999993E+00, 0.54999995E+00, -0.13000000E+00, 0.69999999E+00, -0.59999999E-01,
389 0.48999995E+00, -0.19999999E+00, 0.39999998E+00, -0.21999997E+00, 0.35999995E+00,
390 -0.89999974E-01, 0.57999998E+00, 0.11999995E+00, 0.75000000E+00, 0.14999998E+00,
391 0.69999999E+00, 0.16999996E+00, 0.11099997E+01, 0.88999999E+00, 0.18499994E+01,
392 0.16199999E+01, 0.25399990E+01, 0.22899990E+01, 0.31999998E+01, 0.29099998E+01,
393 0.38399992E+01, 0.35299997E+01, 0.44799995E+01, 0.41499996E+01, 0.51199999E+01,
394 0.47799997E+01, 0.57500000E+01, 0.53899994E+01, 0.63099995E+01, 0.59099998E+01,
395 0.68699999E+01, 0.63299990E+01, 0.71299992E+01, 0.66099997E+01, 0.72999992E+01,
396 0.63099995E+01, 0.62699995E+01, 0.48299999E+01, 0.44899998E+01, 0.28499994E+01,
397 0.23199997E+01, 0.57999998E+00, -0.10999995E+00, -0.97999996E+00, 0.80999994E+00,
398 0.17699995E+01, 0.33699999E+01, 0.41299992E+01, 0.55999994E+01, 0.61499996E+01,
399 0.72899990E+01, 0.73499994E+01, 0.79499998E+01, 0.76699991E+01, 0.81599998E+01,
400 0.78299990E+01, 0.83099995E+01, 0.80099993E+01, 0.85299997E+01, 0.82699995E+01
401 };
402 if ( N > 130 ) throw G4HadronicException(__FILE__, __LINE__, " G4BertiniEvaporationChannel: shell correction for protons called with too large N " );
403 return table[ N-1 ]*MeV;
404}
G4double
double G4double
Definition: G4Types.hh:64

◆ cameronShellCorrectionP()

G4double G4BertiniEvaporationChannel::cameronShellCorrectionP ( G4int  Z)
protected

Definition at line 313 of file G4BertiniEvaporationChannel.cc.

314{
315 // Gives the binding energy correction depending in Z
316 // due to shell correction and pairing energies in MeV
317 //
318 // see Cameron, Canadian Journal of Physics,
319 // vol. 35, 1957, p.1022
320 G4double table[130] = {
321 0.26169998E+02, 0.19250000E+02, 0.24209991E+02, 0.20919998E+02, 0.23149994E+02,
322 0.18009995E+02, 0.19549988E+02, 0.16939987E+02, 0.19729996E+02, 0.17069992E+02,
323 0.18209991E+02, 0.14990000E+02, 0.16009995E+02, 0.12040000E+02, 0.13270000E+02,
324 0.11089998E+02, 0.12169999E+02, 0.10259998E+02, 0.11040000E+02, 0.84099998E+01,
325 0.97899990E+01, 0.73599997E+01, 0.81499996E+01, 0.56299992E+01, 0.58799992E+01,
326 0.31699991E+01, 0.33199997E+01, 0.81999993E+00, 0.18299999E+01, 0.96999997E+00,
327 0.23299999E+01, 0.12699995E+01, 0.29199991E+01, 0.16099997E+01, 0.29099998E+01,
328 0.13499994E+01, 0.23999996E+01, 0.88999999E+00, 0.17399998E+01, 0.35999995E+00,
329 0.94999999E+00, -0.64999998E+00, -0.39999995E-01, -0.17299995E+01, -0.95999998E+00,
330 -0.28699999E+01, -0.20499992E+01, -0.40499992E+01, -0.33999996E+01, -0.57199993E+01,
331 -0.37499990E+01, -0.41299992E+01, -0.24199991E+01, -0.28499994E+01, -0.10099993E+01,
332 -0.13299990E+01, 0.53999996E+00, -0.20000000E-01, 0.17399998E+01, 0.75000000E+00,
333 0.22399998E+01, 0.99999988E+00, 0.19799995E+01, 0.78999996E+00, 0.15399990E+01,
334 0.38999999E+00, 0.10799999E+01, 0.00000000E+00, 0.77999997E+00, -0.34999996E+00,
335 0.57999998E+00, -0.54999995E+00, 0.58999997E+00, -0.60999995E+00, 0.58999997E+00,
336 -0.34999996E+00, 0.31999993E+00, -0.95999998E+00, -0.51999998E+00, -0.20799999E+01,
337 -0.24599991E+01, -0.36399994E+01, -0.15499992E+01, -0.95999998E+00, 0.96999997E+00,
338 0.88000000E+00, 0.23699999E+01, 0.17500000E+01, 0.27199993E+01, 0.18999996E+01,
339 0.25499992E+01, 0.14599991E+01, 0.19299994E+01, 0.85999995E+00, 0.11699991E+01,
340 0.79999983E-01, 0.38999999E+00, -0.75999999E+00, -0.38999999E+00, -0.15099993E+01,
341 -0.11699991E+01, -0.23599997E+01, -0.19499998E+01, -0.30599995E+01, -0.26199999E+01,
342 -0.35499992E+01, -0.29499998E+01, -0.37499990E+01, -0.30699997E+01, -0.37899990E+01,
343 -0.30599995E+01, -0.37699995E+01, -0.30499992E+01, -0.37799997E+01, -0.31199999E+01,
344 -0.38999996E+01, -0.33499994E+01, -0.42399998E+01, -0.38599997E+01, -0.49199991E+01,
345 -0.50599995E+01, -0.67699995E+01, -0.74099998E+01, -0.91799994E+01, -0.10160000E+02,
346 -0.11120000E+02, -0.97599993E+01, -0.92299995E+01, -0.79599991E+01, -0.76499996E+01,
347 };
348 if ( Z > 130 ) throw G4HadronicException(__FILE__, __LINE__, " G4BertiniEvaporationChannel: shell correction for protons called with too large Z " );
349 return table[ Z-1 ]*MeV;
350}

◆ emit()

virtual G4DynamicParticle * G4BertiniEvaporationChannel::emit ( )
pure virtual

Implemented in G4BEDeuteronChannel, G4BEHe3Channel, G4BEHe4Channel, G4BENeutronChannel, G4BEProtonChannel, G4BETritonChannel, and G4BEChargedChannel.

Referenced by G4BertiniEvaporation::BreakItUp().

◆ getCoulomb()

G4double G4BertiniEvaporationChannel::getCoulomb ( )
virtual

Definition at line 159 of file G4BertiniEvaporationChannel.cc.

160{
161 G4int residualZ = nucleusZ - particleZ;
162 G4int residualA = nucleusA - particleA;
163 const G4double factor = 0.84696; // = e / ( 4 pi epsilon_0 r0 ) * 10^-6, r0=1.7E-15
164 // In HETC88 this factor was 0.88235, perhaps due to different r0
165
166 G4double coulombBarrier = factor * particleZ * qmFactor() * residualZ /
167 ( std::pow( G4double(residualA), 0.33333333 ) + rho ) * MeV;
168
169 if ( verboseLevel >= 10 )
170 G4cout << " G4BertiniEvaporationChannel::getThresh() " << G4endl
171 << " residualZ " << residualZ << G4endl
172 << " residualA " << residualA << G4endl
173 << " qmFactor " << qmFactor() << G4endl
174 << " Q " << getQ() << G4endl
175 << " rho " << rho << G4endl
176 << " part Z " << particleZ << G4endl
177 << " (correction) " << correction << G4endl;
178
179 return coulombBarrier;
180}
G4int
int G4int
Definition: G4Types.hh:66
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout

Referenced by getThresh(), and G4BEChargedChannel::sampleKineticEnergy().

◆ getLevelDensityParameter()

G4double G4BertiniEvaporationChannel::getLevelDensityParameter ( )

Definition at line 192 of file G4BertiniEvaporationChannel.cc.

193{
194 G4int residualZ = nucleusZ - particleZ;
195 G4int residualA = nucleusA - particleA;
196 G4double b0 = 8;
197 G4double y0 = 1.5;
198
199 G4double temp = ( residualA - 2.0 * residualZ ) / residualA;
200 G4double smallA = residualA * ( 1.0 + y0 * std::pow( temp, 2. ) ) / b0 / MeV;
201
202 // In HETC98 b0 = b0(E).
203
204 return smallA;
205}

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), G4BEChargedChannel::sampleKineticEnergy(), and G4BENeutronChannel::sampleKineticEnergy().

◆ getName()

G4String G4BertiniEvaporationChannel::getName ( )

Definition at line 116 of file G4BertiniEvaporationChannel.cc.

117{
118 return name;
119}

Referenced by G4BertiniEvaporation::BreakItUp(), and G4BEChargedChannel::calculateProbability().

◆ getParticleA()

G4int G4BertiniEvaporationChannel::getParticleA ( )

Definition at line 75 of file G4BertiniEvaporationChannel.cc.

76{
77 return particleA;
78}

Referenced by G4BertiniEvaporation::BreakItUp().

◆ getParticleZ()

G4int G4BertiniEvaporationChannel::getParticleZ ( )

Definition at line 81 of file G4BertiniEvaporationChannel.cc.

82{
83 return particleZ;
84}

Referenced by G4BertiniEvaporation::BreakItUp().

◆ getProbability()

G4double G4BertiniEvaporationChannel::getProbability ( )
virtual

Definition at line 106 of file G4BertiniEvaporationChannel.cc.

107{
108 if ( verboseLevel >= 4 )
109 G4cout << "G4BerEvaporationChannel : probability of particle " << name
110 << " is " << emissionProbability << G4endl;
111
112 return emissionProbability;
113}

◆ getQ()

G4double G4BertiniEvaporationChannel::getQ ( )
virtual

Definition at line 129 of file G4BertiniEvaporationChannel.cc.

130{
131 // The reaction Q is calculated by using the excess masses of the
132 // secondary particle and the initial and residual nuclei
133 G4int residualZ = nucleusZ - particleZ;
134 G4int residualA = nucleusA - particleA;
135// G4int residualN = residualA - residualZ;
136// G4int nucleusN = nucleusA - nucleusZ;
137
138 // All the following are by default in MeV
139 const G4double e1 = G4NucleiProperties::GetMassExcess( residualA, residualZ );
140 const G4double e2 = G4NucleiProperties::GetMassExcess( nucleusA, nucleusZ );
141
142 return e1 - e2 + G4NucleiProperties::GetMassExcess( particleA, particleZ );
143 // In HETC88 : Mass excesses were calculated with the Cameron mass excess formula,
144 // see Cameron, Canadian Journal of Physics,
145 // vol. 35, 1957, p.1022
146 // Also mass excesses of particles were
147 // 8.3675, 7.5851, 13.7279, 15.8381, 15.8195, 3.6092 for
148 // n, p, D, T, 3He, 4He
149}
G4NucleiProperties::GetMassExcess
static G4double GetMassExcess(const G4int A, const G4int Z)
Definition: G4NucleiProperties.cc:164

Referenced by G4BertiniEvaporation::BreakItUp(), getCoulomb(), and getThresh().

◆ getThresh()

G4double G4BertiniEvaporationChannel::getThresh ( )
virtual

Definition at line 152 of file G4BertiniEvaporationChannel.cc.

153{
154 G4double threshold = getQ() + getCoulomb();
155 return threshold;
156}

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), G4BEChargedChannel::sampleKineticEnergy(), and G4BENeutronChannel::sampleKineticEnergy().

◆ isotropicCosines()

void G4BertiniEvaporationChannel::isotropicCosines ( G4double &  u,
G4double &  v,
G4double &  w 
)
protected

Definition at line 208 of file G4BertiniEvaporationChannel.cc.

209{
210 // Samples isotropic random direction cosines.
211 G4double CosTheta = 1.0 - 2.0 * G4UniformRand();
212 G4double SinTheta = std::sqrt( 1.0 - CosTheta * CosTheta );
213 G4double Phi = twopi * G4UniformRand();
214
215 u = std::cos( Phi ) * SinTheta;
216 v = std::cos( Phi ) * CosTheta,
217 w = std::sin( Phi );
218
219 return;
220}
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:53

Referenced by G4BEDeuteronChannel::emit(), G4BEHe3Channel::emit(), G4BEHe4Channel::emit(), G4BENeutronChannel::emit(), G4BEProtonChannel::emit(), and G4BETritonChannel::emit().

◆ pairingEnergyNeutrons()

G4double G4BertiniEvaporationChannel::pairingEnergyNeutrons ( G4int  N)
protected

Definition at line 261 of file G4BertiniEvaporationChannel.cc.

262{
263// Pairing energy for neutrons P(N), see
264// Cameron, Nuclear Level Spacings, p. 1040
265// Canadian Journal of Physics, vol 36, 1958
266 G4double table[200] = {
267 0.00000000E+00, 0.59799995E+01, 0.00000000E+00, 0.27699995E+01, 0.00000000E+00,
268 0.31599998E+01, 0.00000000E+00, 0.30099993E+01, 0.00000000E+00, 0.16799994E+01,
269 0.00000000E+00, 0.17299995E+01, 0.00000000E+00, 0.21699991E+01, 0.00000000E+00,
270 0.17399998E+01, 0.00000000E+00, 0.17500000E+01, 0.00000000E+00, 0.17199993E+01,
271 0.00000000E+00, 0.16299992E+01, 0.00000000E+00, 0.14099998E+01, 0.00000000E+00,
272 0.12899990E+01, 0.00000000E+00, 0.14699993E+01, 0.00000000E+00, 0.13199997E+01,
273 0.00000000E+00, 0.14599991E+01, 0.00000000E+00, 0.14399996E+01, 0.00000000E+00,
274 0.14599991E+01, 0.00000000E+00, 0.15199995E+01, 0.00000000E+00, 0.15099993E+01,
275 0.00000000E+00, 0.14699993E+01, 0.00000000E+00, 0.14499998E+01, 0.00000000E+00,
276 0.12799997E+01, 0.00000000E+00, 0.12299995E+01, 0.00000000E+00, 0.12699995E+01,
277 0.00000000E+00, 0.61999995E+00, 0.00000000E+00, 0.75999999E+00, 0.00000000E+00,
278 0.12299995E+01, 0.00000000E+00, 0.12199993E+01, 0.00000000E+00, 0.13999996E+01,
279 0.00000000E+00, 0.13599997E+01, 0.00000000E+00, 0.12999992E+01, 0.00000000E+00,
280 0.12899990E+01, 0.00000000E+00, 0.12399998E+01, 0.00000000E+00, 0.12799997E+01,
281 0.00000000E+00, 0.12399998E+01, 0.00000000E+00, 0.11999998E+01, 0.00000000E+00,
282 0.94000000E+00, 0.00000000E+00, 0.99999988E+00, 0.00000000E+00, 0.10499992E+01,
283 0.00000000E+00, 0.53999996E+00, 0.00000000E+00, 0.59999996E+00, 0.00000000E+00,
284 0.75000000E+00, 0.00000000E+00, 0.75000000E+00, 0.00000000E+00, 0.84999996E+00,
285 0.00000000E+00, 0.96999997E+00, 0.00000000E+00, 0.10199995E+01, 0.00000000E+00,
286 0.10499992E+01, 0.00000000E+00, 0.10599995E+01, 0.00000000E+00, 0.10699997E+01,
287 0.00000000E+00, 0.10599995E+01, 0.00000000E+00, 0.10499992E+01, 0.00000000E+00,
288 0.10199995E+01, 0.00000000E+00, 0.96999997E+00, 0.00000000E+00, 0.90999997E+00,
289 0.00000000E+00, 0.82999998E+00, 0.00000000E+00, 0.73999995E+00, 0.00000000E+00,
290 0.65999997E+00, 0.00000000E+00, 0.60999995E+00, 0.00000000E+00, 0.60999995E+00,
291 0.00000000E+00, 0.89999998E+00, 0.00000000E+00, 0.51999998E+00, 0.00000000E+00,
292 0.80999994E+00, 0.00000000E+00, 0.67999995E+00, 0.00000000E+00, 0.71999997E+00,
293 0.00000000E+00, 0.76999998E+00, 0.00000000E+00, 0.67999995E+00, 0.00000000E+00,
294 0.66999996E+00, 0.00000000E+00, 0.79999995E+00, 0.00000000E+00, 0.67999995E+00,
295 0.00000000E+00, 0.63999999E+00, 0.00000000E+00, 0.57999998E+00, 0.00000000E+00,
296 0.54999995E+00, 0.00000000E+00, 0.56999993E+00, 0.00000000E+00, 0.56999993E+00,
297 0.00000000E+00, 0.54999995E+00, 0.00000000E+00, 0.59999996E+00, 0.00000000E+00,
298 0.57999998E+00, 0.00000000E+00, 0.57999998E+00, 0.00000000E+00, 0.60999995E+00,
299 0.00000000E+00, 0.63000000E+00, 0.00000000E+00, 0.64999998E+00, 0.00000000E+00,
300 0.65999997E+00, 0.00000000E+00, 0.64999998E+00, 0.00000000E+00, 0.64999998E+00,
301 0.00000000E+00, 0.63999999E+00, 0.00000000E+00, 0.63999999E+00, 0.00000000E+00,
302 0.63000000E+00, 0.00000000E+00, 0.60999995E+00, 0.00000000E+00, 0.58999997E+00,
303 0.00000000E+00, 0.54999995E+00, 0.00000000E+00, 0.38999999E+00, 0.00000000E+00,
304 0.35999995E+00, 0.00000000E+00, 0.39999998E+00, 0.00000000E+00, 0.39999998E+00,
305 0.00000000E+00, 0.39999998E+00, 0.00000000E+00, 0.39999998E+00, 0.00000000E+00,
306 0.39999998E+00, 0.00000000E+00, 0.39999998E+00, 0.00000000E+00, 0.39999998E+00
307 };
308 if ( N > 200 ) throw G4HadronicException(__FILE__, __LINE__, " G4BertiniEvaporationChannel: pairing energy for neutrons called with too large Z " );
309 return table[ N-1 ]*MeV;
310}

Referenced by setPairingCorrection().

◆ pairingEnergyProtons()

G4double G4BertiniEvaporationChannel::pairingEnergyProtons ( G4int  A)
protected

Definition at line 223 of file G4BertiniEvaporationChannel.cc.

224{
225 // Pairing energy for protons P(Z), see
226 // Cameron, Nuclear Level Spacings, p. 1040
227 // Canadian Journal of Physics, vol 36, 1958
228 G4double table [130] = {
229 0.00000000E+00, 0.54399996E+01, 0.00000000E+00, 0.27599993E+01, 0.00000000E+00,
230 0.33399992E+01, 0.00000000E+00, 0.26999998E+01, 0.00000000E+00, 0.18999996E+01,
231 0.00000000E+00, 0.21199999E+01, 0.00000000E+00, 0.21299992E+01, 0.00000000E+00,
232 0.15399990E+01, 0.00000000E+00, 0.14199991E+01, 0.00000000E+00, 0.15099993E+01,
233 0.00000000E+00, 0.17299995E+01, 0.00000000E+00, 0.14399996E+01, 0.00000000E+00,
234 0.14499998E+01, 0.00000000E+00, 0.13699999E+01, 0.00000000E+00, 0.10899992E+01,
235 0.00000000E+00, 0.13599997E+01, 0.00000000E+00, 0.14199991E+01, 0.00000000E+00,
236 0.13299990E+01, 0.00000000E+00, 0.11999998E+01, 0.00000000E+00, 0.99999988E+00,
237 0.00000000E+00, 0.11599998E+01, 0.00000000E+00, 0.12799997E+01, 0.00000000E+00,
238 0.13799992E+01, 0.00000000E+00, 0.13799992E+01, 0.00000000E+00, 0.13199997E+01,
239 0.00000000E+00, 0.10399990E+01, 0.00000000E+00, 0.11099997E+01, 0.00000000E+00,
240 0.11299992E+01, 0.00000000E+00, 0.12099991E+01, 0.00000000E+00, 0.14299994E+01,
241 0.00000000E+00, 0.11499996E+01, 0.00000000E+00, 0.98999995E+00, 0.00000000E+00,
242 0.90999997E+00, 0.00000000E+00, 0.91999996E+00, 0.00000000E+00, 0.99999988E+00,
243 0.00000000E+00, 0.11099997E+01, 0.00000000E+00, 0.12299995E+01, 0.00000000E+00,
244 0.84999996E+00, 0.00000000E+00, 0.97999996E+00, 0.00000000E+00, 0.71999997E+00,
245 0.00000000E+00, 0.79999995E+00, 0.00000000E+00, 0.76999998E+00, 0.00000000E+00,
246 0.88999999E+00, 0.00000000E+00, 0.91999996E+00, 0.00000000E+00, 0.79999995E+00,
247 0.00000000E+00, 0.80999994E+00, 0.00000000E+00, 0.69000000E+00, 0.00000000E+00,
248 0.69999999E+00, 0.00000000E+00, 0.75999999E+00, 0.00000000E+00, 0.72999996E+00,
249 0.00000000E+00, 0.79999995E+00, 0.00000000E+00, 0.73999995E+00, 0.00000000E+00,
250 0.72999996E+00, 0.00000000E+00, 0.71999997E+00, 0.00000000E+00, 0.71999997E+00,
251 0.00000000E+00, 0.71999997E+00, 0.00000000E+00, 0.70999998E+00, 0.00000000E+00,
252 0.69000000E+00, 0.00000000E+00, 0.67999995E+00, 0.00000000E+00, 0.65999997E+00,
253 0.00000000E+00, 0.60999995E+00, 0.00000000E+00, 0.41999996E+00, 0.00000000E+00,
254 0.35999995E+00, 0.00000000E+00, 0.40999997E+00, 0.00000000E+00, 0.48999995E+00
255 };
256 if ( Z>130 ) throw G4HadronicException(__FILE__, __LINE__, " G4BertiniEvaporationChannel: pairing energy for protons called with too large Z " );
257 return table[ Z-1 ]*MeV;
258}

Referenced by setPairingCorrection().

◆ Q()

G4double G4BertiniEvaporationChannel::Q ( G4double  a,
G4double  z 
)
protected

◆ qmFactor()

G4double G4BertiniEvaporationChannel::qmFactor ( )
virtual

Reimplemented in G4BEDeuteronChannel, G4BEHe3Channel, G4BEHe4Channel, G4BEProtonChannel, and G4BETritonChannel.

Definition at line 183 of file G4BertiniEvaporationChannel.cc.

184{
185 // Coefficient k_p for empirical cross section
186 // formula presented in Dostrovsky, Phys. Rev.,
187 // vol. 116, 1959
188 return 0;
189}

Referenced by G4BEChargedChannel::calculateProbability(), and getCoulomb().

◆ setExcitationEnergy()

void G4BertiniEvaporationChannel::setExcitationEnergy ( G4double  inputE)

Definition at line 87 of file G4BertiniEvaporationChannel.cc.

88{
89 excitationEnergy = energy;
90}

◆ setNucleusA()

void G4BertiniEvaporationChannel::setNucleusA ( G4int  inputA)

Definition at line 63 of file G4BertiniEvaporationChannel.cc.

64{
65 nucleusA = a;
66}

◆ setNucleusZ()

void G4BertiniEvaporationChannel::setNucleusZ ( G4int  inputZ)

Definition at line 69 of file G4BertiniEvaporationChannel.cc.

70{
71 nucleusZ = z;
72}

◆ setPairingCorrection()

void G4BertiniEvaporationChannel::setPairingCorrection ( G4int  isCorrection)

Definition at line 93 of file G4BertiniEvaporationChannel.cc.

94{
95 G4int residualZ = nucleusZ - particleZ;
96 G4int residualA = nucleusA - particleA;
97
98 // Note: pairing energy is zero for odd nucleon number.
99
100 if ( isCorrection ) correction = pairingEnergyProtons( residualZ )
101 + pairingEnergyNeutrons( residualA - residualZ );
102 else correction = 0;
103}

◆ setProbability()

void G4BertiniEvaporationChannel::setProbability ( G4double  newProb)
virtual

Definition at line 122 of file G4BertiniEvaporationChannel.cc.

123{
124 emissionProbability = newProb;
125 return;
126}

Referenced by G4BertiniEvaporation::BreakItUp().

◆ setQ()

void G4BertiniEvaporationChannel::setQ ( G4double  inputQ)

◆ setVerboseLevel()

void G4BertiniEvaporationChannel::setVerboseLevel ( G4int  verbose)

Definition at line 57 of file G4BertiniEvaporationChannel.cc.

58{
59 verboseLevel = level;
60}

Member Data Documentation

◆ correction

G4double G4BertiniEvaporationChannel::correction
protected

Definition at line 75 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), getCoulomb(), G4BEChargedChannel::sampleKineticEnergy(), G4BENeutronChannel::sampleKineticEnergy(), and setPairingCorrection().

◆ emissionProbability

G4double G4BertiniEvaporationChannel::emissionProbability
protected

Definition at line 73 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), getProbability(), and setProbability().

◆ excitationEnergy

G4double G4BertiniEvaporationChannel::excitationEnergy
protected

Definition at line 76 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), G4BEChargedChannel::sampleKineticEnergy(), G4BENeutronChannel::sampleKineticEnergy(), and setExcitationEnergy().

◆ exmass

G4double G4BertiniEvaporationChannel::exmass
protected

Definition at line 72 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEDeuteronChannel::G4BEDeuteronChannel().

◆ name

G4String G4BertiniEvaporationChannel::name
protected

Definition at line 66 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEDeuteronChannel::G4BEDeuteronChannel(), G4BEHe3Channel::G4BEHe3Channel(), G4BEHe4Channel::G4BEHe4Channel(), G4BENeutronChannel::G4BENeutronChannel(), G4BEProtonChannel::G4BEProtonChannel(), G4BETritonChannel::G4BETritonChannel(), getName(), and getProbability().

◆ nucleusA

G4int G4BertiniEvaporationChannel::nucleusA
protected

Definition at line 68 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), getCoulomb(), getLevelDensityParameter(), getQ(), setNucleusA(), and setPairingCorrection().

◆ nucleusZ

G4int G4BertiniEvaporationChannel::nucleusZ
protected

Definition at line 69 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), G4BEChargedChannel::coulombFactorForProton(), getCoulomb(), getLevelDensityParameter(), getQ(), G4BEChargedChannel::qmFactorForAlpha(), G4BEChargedChannel::qmFactorForProton(), setNucleusZ(), and setPairingCorrection().

◆ particleA

G4int G4BertiniEvaporationChannel::particleA
protected

Definition at line 70 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), G4BEDeuteronChannel::G4BEDeuteronChannel(), G4BEHe3Channel::G4BEHe3Channel(), G4BEHe4Channel::G4BEHe4Channel(), G4BENeutronChannel::G4BENeutronChannel(), G4BEProtonChannel::G4BEProtonChannel(), G4BETritonChannel::G4BETritonChannel(), getCoulomb(), getLevelDensityParameter(), getParticleA(), getQ(), and setPairingCorrection().

◆ particleZ

G4int G4BertiniEvaporationChannel::particleZ
protected

Definition at line 71 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), G4BEChargedChannel::coulombFactorForProton(), G4BEDeuteronChannel::G4BEDeuteronChannel(), G4BEHe3Channel::G4BEHe3Channel(), G4BEHe4Channel::G4BEHe4Channel(), G4BENeutronChannel::G4BENeutronChannel(), G4BEProtonChannel::G4BEProtonChannel(), G4BETritonChannel::G4BETritonChannel(), getCoulomb(), getLevelDensityParameter(), getParticleZ(), getQ(), G4BEChargedChannel::qmFactorForAlpha(), G4BEChargedChannel::qmFactorForProton(), and setPairingCorrection().

◆ rho

G4double G4BertiniEvaporationChannel::rho
protected

Definition at line 74 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEDeuteronChannel::G4BEDeuteronChannel(), G4BEHe3Channel::G4BEHe3Channel(), G4BEHe4Channel::G4BEHe4Channel(), G4BENeutronChannel::G4BENeutronChannel(), G4BEProtonChannel::G4BEProtonChannel(), G4BETritonChannel::G4BETritonChannel(), and getCoulomb().

◆ spin

G4int G4BertiniEvaporationChannel::spin
protected

Definition at line 78 of file G4BertiniEvaporationChannel.hh.

◆ verboseLevel

G4int G4BertiniEvaporationChannel::verboseLevel
protected

Definition at line 67 of file G4BertiniEvaporationChannel.hh.

Referenced by G4BEChargedChannel::calculateProbability(), G4BENeutronChannel::calculateProbability(), G4BEChargedChannel::G4BEChargedChannel(), G4BEDeuteronChannel::G4BEDeuteronChannel(), G4BEHe3Channel::G4BEHe3Channel(), G4BEHe4Channel::G4BEHe4Channel(), G4BENeutronChannel::G4BENeutronChannel(), G4BEProtonChannel::G4BEProtonChannel(), G4BertiniEvaporationChannel(), G4BETritonChannel::G4BETritonChannel(), getCoulomb(), getProbability(), G4BEChargedChannel::sampleKineticEnergy(), G4BENeutronChannel::sampleKineticEnergy(), and setVerboseLevel().

The documentation for this class was generated from the following files: