56 lowestEnergyLimit = 0.1*eV;
65 delete photonEvaporation;
90 if(M - mass <= lowestEnergyLimit) {
95 G4cout <<
"G4NeutronRadCapture::DoIt: Eini(MeV)="
98 <<
" Z= " << Z <<
" A= " << A <<
G4endl;
100 G4double e1 = (M - mass)*(M + mass)/(2*M);
102 if(cost > 1.0) {cost = 1.0;}
103 else if(cost < -1.0) {cost = -1.0;}
104 G4double sint = std::sqrt((1. - cost)*(1.0 + cost));
106 G4LorentzVector lv2(e1*sint*std::cos(phi),e1*sint*std::sin(phi),e1*cost,e1);
114 else if (Z == 2 && A == 3) {theDef =
G4He3::He3();}
123 G4cout <<
"Gamma 4-mom: " << lv2 <<
" "
136 G4cout <<
"G4NeutronRadCapture::ApplyYourself initial G4Fragmet:" <<
G4endl;
145 fv->push_back(aFragment);
146 size_t n = fv->size();
149 G4cout <<
"G4NeutronRadCapture: " << n <<
" final particle" <<
G4endl;
151 for(
size_t i=0; i<n; ++i) {
166 <<
" " << lvres <<
G4endl;
std::vector< G4Fragment * > G4FragmentVector
G4DLLIMPORT std::ostream G4cout
Hep3Vector boostVector() const
HepLorentzVector & boost(double, double, double)
static G4Deuteron * Deuteron()
const G4LorentzVector & GetMomentum() const
G4ParticleDefinition * GetParticleDefinition() const
void SetStatusChange(G4HadFinalStateStatus aS)
void AddSecondary(G4DynamicParticle *aP)
G4double GetKineticEnergy() const
const G4LorentzVector & Get4Momentum() const
G4HadFinalState theParticleChange
void SetMinEnergy(G4double anEnergy)
void SetMaxEnergy(const G4double anEnergy)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int J=0)
virtual ~G4NeutronRadCapture()
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
static G4double GetNuclearMass(const G4double A, const G4double Z)
const G4String & GetParticleName() const
static G4ParticleTable * GetParticleTable()
G4IonTable * GetIonTable()
virtual G4FragmentVector * BreakUpFragment(G4Fragment *theNucleus)
static G4Triton * Triton()