61 #ifdef G4MULTITHREADED
66 : fLevelManager(nullptr), fTransition(p), fPolarization(nullptr),
67 fVerbose(0), fPoints(0), vShellNumber(-1), fIndex(0),
69 fICM(true), fRDM(false), fSampleTime(true),
70 fCorrelatedGamma(false), isInitialised(false)
98 G4cout <<
"### G4PhotonEvaporation is initialized " <<
this <<
G4endl;
114 #ifdef G4MULTITHREADED
115 G4MUTEXLOCK(&G4PhotonEvaporation::PhotonEvaporationMutex);
119 static const G4float GRWfactor = 0.3f;
125 #ifdef G4MULTITHREADED
149 G4cout <<
"G4PhotonEvaporation::EmittedFragment: "
160 G4cout <<
"G4PhotonEvaporation::EmittedFragment: remove "
169 G4cout <<
"G4PhotonEvaporation::EmittedFragment: RDM= "
186 products->push_back(aNucleus);
195 G4cout <<
"G4PhotonEvaporation::BreakUpChain RDM= " <<
fRDM <<
" "
212 products->push_back(gamma);
214 G4cout <<
"G4PhotonEvaporation::BreakUpChain: "
243 G4cout <<
"G4PhotonEvaporation::GetEmissionProbability: Z="
255 static const G4float GREfactor = 5.0f;
290 G4double p0 =
G4Exp(-2.0*xsqr)*gammaR2*gammaE2/(egdp2*egdp2 + gammaR2);
296 gammaR2 = gammaE2*wres2;
297 egdp2 = gammaE2 - eres2;
299 *gammaR2*gammaE2/(egdp2*egdp2 + gammaR2);
351 G4bool isDiscrete =
false;
375 if(ntrans > 0) { isDiscrete =
true; }
385 <<
" Exc= " << eexc <<
" Emax= "
388 <<
" Ntr= " << ntrans <<
" discrete: " << isDiscrete
419 if(efinal >= eexc && 0 <
fIndex) {
433 G4cout <<
"Continues emission efinal(MeV)= " << efinal <<
G4endl;
438 G4cout <<
"Discrete emission from level Index= " <<
fIndex
454 G4cout <<
"Ntrans= " << ntrans <<
" idx= " << idx
455 <<
" ICM= " <<
fICM <<
" JP1= " << JP1 <<
G4endl;
459 if(
fICM && prob < 1.0) {
463 rndm = (rndm - prob)/(1.0 - prob);
486 isDiscrete, isGamma);
495 if(efinal == 0.0 &&
fIndex > 0) {
501 G4cout <<
"Final level E= " << efinal <<
" time= " << time
502 <<
" idxFinal= " <<
fIndex <<
" isDiscrete: " << isDiscrete
503 <<
" isGamma: " << isGamma <<
" multiP= " << multiP
505 <<
" JP1= " << JP1 <<
" JP2= " << JP2 <<
G4endl;
static G4NuclearPolarizationStore * GetInstance()
G4int SampleShell(size_t idx, G4double rndm) const
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double GetCreationTime() const
void SetNuclearPolarization(G4NuclearPolarization *)
virtual G4bool BreakUpChain(G4FragmentVector *theResult, G4Fragment *theNucleus) final
virtual void SetICM(G4bool)
virtual G4Fragment * EmittedFragment(G4Fragment *theNucleus) final
void SetTwoJMAX(G4int val)
void SetCreationTime(G4double time)
G4DeexPrecoParameters * GetParameters()
G4bool GetInternalConversionFlag() const
static G4float GRWidth[MAXGRDATA]
G4NuclearPolarization * GetNuclearPolarization() const
size_t SampleGammaTransition(G4double rndm) const
static G4NuclearLevelData * GetInstance()
static constexpr double hbarc
static G4float GREnergy[MAXGRDATA]
G4PhotonEvaporation(G4GammaTransition *ptr=nullptr)
static constexpr double pi2
void SetExcitationEnergy(G4double val)
static const G4double emax
G4double GetGroundStateMass() const
G4double GetMinExcitation() const
G4FragmentVector * BreakItUp(const G4Fragment &theNucleus)
virtual ~G4PhotonEvaporation()
G4float MultipolarityRatio(size_t idx) const
G4double G4Log(G4double x)
const G4NucLevel * GetLevel(size_t i) const
void SetPolarizationFlag(G4bool val)
static constexpr double neutron_mass_c2
G4double powZ(G4int Z, G4double y) const
G4int FloatingLevel(size_t i) const
G4double GetMaxLifeTime() const
size_t FinalExcitationIndex(size_t idx) const
#define G4MUTEX_INITIALIZER
static G4Pow * GetInstance()
G4double NearestLevelEnergy(G4double energy, size_t index=0) const
G4float GammaProbability(size_t idx) const
G4NuclearPolarizationStore * fNucPStore
void SetSpin(G4double value)
G4double ComputeGroundStateMass(G4int Z, G4int A) const
void SetGammaTransition(G4GammaTransition *)
static constexpr double MeV
double A(double temperature)
virtual void Initialise() final
static constexpr double eV
G4double LifeTime(size_t i) const
virtual G4double GetEmissionProbability(G4Fragment *theNucleus) final
G4NuclearLevelData * fNuclearLevelData
#define G4MUTEXUNLOCK(mutex)
G4Fragment * GenerateGamma(G4Fragment *nucleus)
G4double G4ParticleHPJENDLHEData::G4double result
void SetFloatingLevelNumber(G4int value)
virtual G4Fragment * SampleTransition(G4Fragment *nucleus, G4double newExcEnergy, G4double mpRatio, G4int JP1, G4int JP2, G4int MP, G4int shell, G4bool isDiscrete, G4bool isGamma)
G4double GetLevelDensity() const
void RemoveMe(G4NuclearPolarization *ptr)
G4int SpinTwo(size_t i) const
void SetVerbose(G4int val)
G4NuclearPolarization * FindOrBuild(G4int Z, G4int A, G4double Eexc)
virtual G4double GetFinalLevelEnergy(G4int Z, G4int A, G4double energy) final
#define G4MUTEXLOCK(mutex)
G4GammaTransition * fTransition
G4GLOB_DLL std::ostream G4cout
G4double fCummProbability[MAXDEPOINT]
const G4LevelManager * fLevelManager
virtual void RDMForced(G4bool)
size_t NumberOfTransitions() const
G4NuclearPolarization * fPolarization
G4double LevelEnergy(size_t i) const
G4double GetExcitationEnergy() const
std::vector< G4Fragment * > G4FragmentVector
void InitialiseLevelManager(G4int Z, G4int A)
G4bool CorrelatedGamma() const
G4int TransitionType(size_t idx) const
static constexpr double millibarn
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
size_t NearestLevelIndex(G4double energy, size_t index=0) const
virtual G4double GetUpperLevelEnergy(G4int Z, G4int A) final