63 char* env = getenv(
"G4UsePhotonEvaporationOLD");
107 G4cout <<
"G4NeutronRadCapture::DoIt: Eini(MeV)="
110 <<
" Z= " << Z <<
" A= " << A <<
G4endl;
115 if(cost > 1.0) {cost = 1.0;}
116 else if(cost < -1.0) {cost = -1.0;}
117 G4double sint = std::sqrt((1. - cost)*(1.0 + cost));
134 else if (Z == 2 && A == 3) {theDef =
G4He3::He3();}
139 G4cout <<
"Gamma 4-mom: " << lv2 <<
" "
155 G4double ptot = std::sqrt((etot - mass)*(etot + mass));
157 lab4mom.set(v.x()*ptot,v.y()*ptot,v.z()*ptot,etot);
163 G4cout <<
"G4NeutronRadCapture::ApplyYourself initial G4Fragmet:"
173 fv->push_back(aFragment);
174 size_t n = fv->size();
177 G4cout <<
"G4NeutronRadCapture: " << n <<
" final particle" <<
G4endl;
179 for(
size_t i=0; i<
n; ++i) {
191 else if (Z == 2 && A == 3) {theDef =
G4He3::He3();}
206 <<
" Ekin(MeV)= " << etot/
MeV
215 if(timeF < 0.0) { timeF = 0.0; }
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
virtual ~G4NeutronRadCapture()
static G4double GetNuclearMass(const G4double A, const G4double Z)
virtual G4FragmentVector * BreakUpFragment(G4Fragment *theNucleus)
CLHEP::Hep3Vector G4ThreeVector
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
const G4ParticleDefinition * GetParticleDefinition() const
virtual void SetICM(G4bool)
const G4String & GetParticleName() const
G4VEvaporationChannel * photonEvaporation
void SetStatusChange(G4HadFinalStateStatus aS)
void SetMinEnergy(G4double anEnergy)
G4IonTable * GetIonTable() const
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
G4double GetCreationTime() const
const G4LorentzVector & GetMomentum() const
G4double GetKineticEnergy() const
static const double twopi
G4double GetGlobalTime() const
std::vector< G4Fragment * > G4FragmentVector
static G4Triton * Triton()
const G4LorentzVector & Get4Momentum() const
static G4Deuteron * Deuteron()
void SetTime(G4double aT)
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
void SetMaxEnergy(const G4double anEnergy)
G4HadFinalState theParticleChange
G4double lowestEnergyLimit
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
G4IonTable * theTableOfIons
G4double GetExcitationEnergy() const
CLHEP::HepLorentzVector G4LorentzVector