62 for(
G4int i=0; i<6; ++i) {
89 G4cout <<
"#Unstable decay " <<
" Z= " <<
Z <<
" A= " <<
A
98 for(
G4int i=0; i<6; ++i) {
101 if(Zres >= 0 && Ares >= Zres && Ares >=
Afr[i]) {
103 for(
G4int j=0; j<6; ++j) {
104 if(Zres ==
Zfr[j] && Ares ==
Afr[j]) {
116 if(delm > 0.0) { isChannel =
true; }
122 if(isChannel) {
break; }
128 if(lman && e >= 0.0) {
140 if(e > 0.0) { mass1 += e; }
148 if(mass + dmlimit < massmin) {
return false; }
150 G4cout <<
"#Unstable decay correction: Z= " <<
Z <<
" A= " <<
A
152 <<
" deltaM(MeV)= " << mass - massmin
157 G4double mom = std::sqrt((e - mass)*(e + mass));
163 G4double e2 = 0.5*((mass - mass1)*(mass + mass1) + mass2*mass2)/mass;
175 results->push_back(frag);
static const G4double e2[44]
std::vector< G4Fragment * > G4FragmentVector
CLHEP::HepLorentzVector G4LorentzVector
G4ThreeVector G4RandomDirection()
G4GLOB_DLL std::ostream G4cout
Hep3Vector boostVector() const
HepLorentzVector & boost(double, double, double)
void set(double x, double y, double z, double t)
G4double GetExcitationEnergy() const
const G4LorentzVector & GetMomentum() const
void SetCreatorModelID(G4int value)
G4double GetCreationTime() const
void SetCreationTime(G4double time)
void SetMomentum(const G4LorentzVector &value)
void SetZandA_asInt(G4int Znew, G4int Anew)
const G4LevelManager * GetLevelManager(G4int Z, G4int A)
static G4NuclearLevelData * GetInstance()
static G4double GetNuclearMass(const G4double A, const G4double Z)
static G4int GetModelID(const G4int modelIndex)
G4NuclearLevelData * fLevelData
G4UnstableFragmentBreakUp()
static const G4int Afr[6]
virtual G4bool BreakUpChain(G4FragmentVector *, G4Fragment *) final
virtual G4double GetEmissionProbability(G4Fragment *fragment) final
virtual ~G4UnstableFragmentBreakUp()
static const G4int Zfr[6]
static constexpr double MeV
static constexpr double eV
T max(const T t1, const T t2)
brief Return the largest of the two arguments