57 photonEvaporation(nullptr),lab4mom(0.,0.,0.,0.)
119 G4cout <<
"G4NeutronRadCapture::DoIt: Eini(MeV)="
122 <<
" Z= " << Z <<
" A= " << A <<
G4endl;
124 G4double e1 = (M - mass)*(M + mass)/(2*M);
138 else if (Z == 2 && A == 3) {theDef =
G4He3::He3();}
143 G4cout <<
"Gamma 4-mom: " << lv2 <<
" "
159 G4double ptot = std::sqrt((etot - mass)*(etot + mass));
167 G4cout <<
"G4NeutronRadCapture::ApplyYourself initial G4Fragmet:"
177 fv->push_back(aFragment);
178 size_t n = fv->size();
181 G4cout <<
"G4NeutronRadCapture: " << n <<
" final particle icID= " <<
icID <<
G4endl;
183 for(
size_t i=0; i<
n; ++i) {
195 else if (Z == 2 && A == 3) {theDef =
G4He3::He3();}
210 <<
" Ekin(MeV)= " << etot/
MeV
219 if(timeF < 0.0) { timeF = 0.0; }
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double GetCreationTime() const
const G4LorentzVector & GetMomentum() const
virtual ~G4NeutronRadCapture()
void SetMinEnergy(G4double anEnergy)
void SetMaxEnergy(const G4double anEnergy)
static G4ParticleTable * GetParticleTable()
static constexpr double MeV
G4DeexPrecoParameters * GetParameters()
G4IonTable * theTableOfIons
static G4NuclearLevelData * GetInstance()
G4int GetInternalConversionID() const
const G4String & GetParticleName() const
G4IonTable * GetIonTable() const
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
G4ThreeVector G4RandomDirection()
static constexpr double TeV
G4double GetMinExcitation() const
G4double GetPDGMass() const
G4double lowestEnergyLimit
G4FragmentVector * BreakUpFragment(G4Fragment *theNucleus)
static G4Deuteron * Deuteron()
G4double GetGlobalTime() const
virtual void InitialiseModel() final
void SetCreatorModelType(G4int idx)
double A(double temperature)
static G4double GetNuclearMass(const G4double A, const G4double Z)
static constexpr double eV
void set(double x, double y, double z, double t)
static G4Electron * Electron()
const G4LorentzVector & Get4Momentum() const
virtual void SetICM(G4bool)
G4double GetKineticEnergy() const
virtual void Initialise()
static G4Triton * Triton()
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus) final
G4VEvaporationChannel * photonEvaporation
G4GLOB_DLL std::ostream G4cout
Hep3Vector boostVector() const
const G4ParticleDefinition * GetParticleDefinition() const
void SetTime(G4double aT)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
G4double GetExcitationEnergy() const
static constexpr double keV
std::vector< G4Fragment * > G4FragmentVector
static constexpr double GeV
G4HadFinalState theParticleChange
HepLorentzVector & boost(double, double, double)
void SetStatusChange(G4HadFinalStateStatus aS)
const G4ParticleDefinition * electron