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)) {
150 availableEnergy -= neutron_mass_c2 - nenergy;
157 residualMass*residualMass/availableEnergy);
181 G4int reentryCount = 0;
188 pDef = nucleons[index].GetDefinition();
197 G4double Enu = 0.5*(Ecms - neutron_mass_c2*neutron_mass_c2/Ecms);
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";
G4MuMinusCapturePrecompound(G4VPreCompoundModel *ptr=0)
static G4double GetNuclearMass(const G4double A, const G4double Z)
std::ostringstream G4ExceptionDescription
const G4ParticleDefinition * fNeutron
CLHEP::Hep3Vector G4ThreeVector
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)=0
const G4ParticleDefinition * fProton
void ModelDescription(std::ostream &outFile) const
G4ThreeVector G4RandomDirection()
const G4String & GetModelName() const
void SetStatusChange(G4HadFinalStateStatus aS)
G4VPreCompoundModel * fPreCompound
std::vector< G4ReactionProduct * > G4ReactionProductVector
const G4ParticleDefinition * GetDefinition() const
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
static G4NeutrinoMu * NeutrinoMu()
G4double GetBoundEnergy() const
void Init(G4int theA, G4int theZ)
G4double GetGlobalTime() const
static G4Triton * Triton()
static G4Proton * Proton()
void AddNewParticle(const G4ParticleDefinition *aParticle, G4ThreeVector &direction, G4double kinEnergy)
static G4Neutron * Neutron()
const std::vector< G4Nucleon > & GetNucleons()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4double GetKineticEnergy() const
G4HadronicInteraction * FindModel(const G4String &name)
G4double GetPDGMass() const
G4Fancy3DNucleus fNucleus
~G4MuMinusCapturePrecompound()
static G4HadronicInteractionRegistry * Instance()
static constexpr double GeV
G4ThreeVector GetMomentum() const
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
static G4MuonMinus * MuonMinus()
static constexpr double MeV
static constexpr double keV
G4int GetNumberOfSecondaries() const
CLHEP::HepLorentzVector G4LorentzVector