113 G4double muMom = std::sqrt(muBindingEnergy*(muBindingEnergy + 2.0*
fMuMass));
122 if((1 == Z && 1 == A) || (2 == Z && 3 == A)) {
132 residualMass*residualMass/availableEnergy);
153 G4int reentryCount = 0;
160 pDef = nucleons[index].GetDefinition();
169 G4double Enu = 0.5*(Ecms - neutron_mass_c2*neutron_mass_c2/Ecms);
178 momResidual = momInitial - momNu;
179 eEx = momResidual.mag() - residualMass;
181 if(reentryCount > 1000) {
184 ed <<
"Target Z= " << Z
185 <<
" A= " << A <<
" Eex(MeV)= " << eEx/
MeV <<
G4endl;
186 ed <<
" ApplyYourself does not completed after 1000 attempts";
187 G4Exception(
"G4MuMinusCapturePrecompound::ApplyYourself",
"had006",
196 initialState.SetNumberOfExcitedParticle(2,0);
197 initialState.SetNumberOfHoles(1,1);
201 size_t n = rpv->size();
202 for(
size_t i=0; i<
n; ++i) {
214 G4cout <<
"G4MuMinusCapturePrecompound::ApplyYourself: Nsec= "
216 <<
" E0(MeV)= " <<availableEnergy/
MeV
217 <<
" Mres(GeV)= " <<residualMass/
GeV
227 outFile <<
"Sampling of mu- capture by atomic nucleus from K-shell"
228 <<
" mesoatom orbit.\n"
229 <<
"Primary reaction mu- + p -> n + neutrino, neutron providing\n"
230 <<
" initial excitation of the target nucleus and PreCompound"
231 <<
" model samples final state\n";
G4MuMinusCapturePrecompound(G4VPreCompoundModel *ptr=0)
static G4double GetNuclearMass(const G4double A, const G4double Z)
std::ostringstream G4ExceptionDescription
CLHEP::Hep3Vector G4ThreeVector
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)=0
void ModelDescription(std::ostream &outFile) const
void AddNewParticle(G4ParticleDefinition *aParticle, G4ThreeVector &direction, G4double kinEnergy)
G4ThreeVector G4RandomDirection()
const G4String & GetModelName() const
G4ParticleDefinition * fNeutron
G4ParticleDefinition * GetDefinition() const
void SetStatusChange(G4HadFinalStateStatus aS)
G4VPreCompoundModel * fPreCompound
std::vector< G4ReactionProduct * > G4ReactionProductVector
G4GLOB_DLL std::ostream G4cout
static G4NeutrinoMu * NeutrinoMu()
G4double GetBoundEnergy() const
void Init(G4int theA, G4int theZ)
G4double GetGlobalTime() const
static G4Triton * Triton()
static G4Proton * Proton()
static G4Neutron * Neutron()
static const G4double A[nN]
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()
G4ThreeVector GetMomentum() const
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
static G4MuonMinus * MuonMinus()
G4int GetNumberOfSecondaries() const
CLHEP::HepLorentzVector G4LorentzVector
G4ParticleDefinition * fProton