113 G4double muMom =std::sqrt(muBindingEnergy*(muBindingEnergy + 2.0*
fMuMass));
124 if((1 == Z && 1 == A) || (2 == Z && 3 == A)) {
134 residualMass*residualMass/availableEnergy);
144 }
else if((1 == Z && 2 == A) || (2 == Z && 4 == A)) {
157 residualMass*residualMass/availableEnergy);
181 G4int reentryCount = 0;
188 pDef = nucleons[index].GetDefinition();
206 momResidual = momInitial - momNu;
207 eEx = momResidual.
mag() - residualMass;
208 if(eEx < 0.0 && eEx + nenergy >= 0.0) {
209 momResidual.
set(0.0, 0.0, 0.0, residualMass);
215 if(reentryCount > 100 && eEx < 0.0) {
218 ed <<
"Target Z= " << Z
219 <<
" A= " << A <<
" Eex(MeV)= " << eEx/
MeV <<
G4endl;
220 ed <<
" ApplyYourself does not completed after 100 attempts -"
221 <<
" excitation energy is set to zero";
222 G4Exception(
"G4MuMinusCapturePrecompound::ApplyYourself",
"had006",
224 momResidual.
set(0.0, 0.0, 0.0, residualMass);
234 initialState.SetNumberOfExcitedParticle(2,0);
235 initialState.SetNumberOfHoles(1,1);
239 size_t n = rpv->size();
240 for(
size_t i=0; i<
n; ++i) {
252 G4cout <<
"G4MuMinusCapturePrecompound::ApplyYourself: Nsec= "
254 <<
" E0(MeV)= " <<availableEnergy/
MeV
255 <<
" Mres(GeV)= " <<residualMass/
GeV
265 outFile <<
"Sampling of mu- capture by atomic nucleus from K-shell"
266 <<
" mesoatom orbit.\n"
267 <<
"Primary reaction mu- + p -> n + neutrino, neutron providing\n"
268 <<
" initial excitation of the target nucleus and PreCompound"
269 <<
" model samples final state\n";
G4double GetBoundEnergy() const
void AddNewParticle(const G4ParticleDefinition *aParticle, G4ThreeVector &direction, G4double kinEnergy)
std::ostringstream G4ExceptionDescription
static constexpr double MeV
G4Fancy3DNucleus fNucleus
static constexpr double keV
const G4ParticleDefinition * fProton
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)=0
G4MuMinusCapturePrecompound(G4VPreCompoundModel *ptr=0)
const G4ParticleDefinition * GetDefinition() const
static G4MuonMinus * MuonMinus()
G4ThreeVector G4RandomDirection()
static G4Proton * Proton()
G4double GetPDGMass() const
~G4MuMinusCapturePrecompound()
static constexpr double neutron_mass_c2
G4double GetKineticEnergy() const
G4double GetGlobalTime() const
static G4NeutrinoMu * NeutrinoMu()
static G4HadronicInteractionRegistry * Instance()
void Init(G4int theA, G4int theZ)
double A(double temperature)
void ModelDescription(std::ostream &outFile) const
static G4double GetNuclearMass(const G4double A, const G4double Z)
const std::vector< G4Nucleon > & GetNucleons()
G4ThreeVector GetMomentum() const
void set(double x, double y, double z, double t)
std::vector< G4ReactionProduct * > G4ReactionProductVector
G4int GetNumberOfSecondaries() const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
static G4Triton * Triton()
static G4Neutron * Neutron()
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
const G4ParticleDefinition * fNeutron
G4GLOB_DLL std::ostream G4cout
Hep3Vector boostVector() const
G4HadronicInteraction * FindModel(const G4String &name)
static constexpr double GeV
HepLorentzVector & boost(double, double, double)
void SetStatusChange(G4HadFinalStateStatus aS)
const G4String & GetModelName() const
G4VPreCompoundModel * fPreCompound