#include <G4MuMinusCapturePrecompound.hh>

Inheritance diagram for G4MuMinusCapturePrecompound:

Public Member Functions
	G4MuMinusCapturePrecompound (G4VPreCompoundModel *ptr=0)

	~G4MuMinusCapturePrecompound ()

G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)

void	ModelDescription (std::ostream &outFile) const

Public Member Functions inherited from G4HadronicInteraction
	G4HadronicInteraction (const G4String &modelName="HadronicModel")

virtual	~G4HadronicInteraction ()

virtual G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)

virtual G4bool	IsApplicable (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)

G4double	GetMinEnergy () const

G4double	GetMinEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMinEnergy (G4double anEnergy)

void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)

G4double	GetMaxEnergy () const

G4double	GetMaxEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMaxEnergy (const G4double anEnergy)

void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)

G4int	GetVerboseLevel () const

void	SetVerboseLevel (G4int value)

const G4String &	GetModelName () const

void	DeActivateFor (const G4Material *aMaterial)

void	ActivateFor (const G4Material *aMaterial)

void	DeActivateFor (const G4Element *anElement)

void	ActivateFor (const G4Element *anElement)

G4bool	IsBlocked (const G4Material *aMaterial) const

G4bool	IsBlocked (const G4Element *anElement) const

void	SetRecoilEnergyThreshold (G4double val)

G4double	GetRecoilEnergyThreshold () const

virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const

virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const

void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	InitialiseModel ()

	G4HadronicInteraction (const G4HadronicInteraction &right)=delete

const G4HadronicInteraction &	operator= (const G4HadronicInteraction &right)=delete

G4bool	operator== (const G4HadronicInteraction &right) const =delete

G4bool	operator!= (const G4HadronicInteraction &right) const =delete

Additional Inherited Members
Protected Member Functions inherited from G4HadronicInteraction
void	SetModelName (const G4String &nam)

G4bool	IsBlocked () const

void	Block ()

Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState	theParticleChange

G4int	verboseLevel

G4double	theMinEnergy

G4double	theMaxEnergy

G4bool	isBlocked

Detailed Description

Definition at line 64 of file G4MuMinusCapturePrecompound.hh.

Constructor & Destructor Documentation

◆ G4MuMinusCapturePrecompound()

G4MuMinusCapturePrecompound::G4MuMinusCapturePrecompound ( G4VPreCompoundModel * ptr = 0 )

Definition at line 63 of file G4MuMinusCapturePrecompound.cc.

  : G4HadronicInteraction("muMinusNuclearCapture")
{ 
  fMuMass = G4MuonMinus::MuonMinus()->GetPDGMass(); 
  fProton = G4Proton::Proton();
  fNeutron = G4Neutron::Neutron();
  fThreshold = 10*MeV;
  fTime = 0.0;
  fPreCompound = ptr;
  if(!ptr) { 
    G4HadronicInteraction* p =
      G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
    ptr = static_cast<G4VPreCompoundModel*>(p); 
    fPreCompound = ptr;
    if(!ptr) { fPreCompound = new G4PreCompoundModel(); }
  }
}

◆ ~G4MuMinusCapturePrecompound()

G4MuMinusCapturePrecompound::~G4MuMinusCapturePrecompound ( )

Definition at line 84 of file G4MuMinusCapturePrecompound.cc.

{
  result.Clear();
}

Member Function Documentation

◆ ApplyYourself()

G4HadFinalState * G4MuMinusCapturePrecompound::ApplyYourself	(	const G4HadProjectile &	aTrack,
		G4Nucleus &	targetNucleus )

virtual

Reimplemented from G4HadronicInteraction.

Definition at line 92 of file G4MuMinusCapturePrecompound.cc.

{
  result.Clear();
  result.SetStatusChange(stopAndKill);
  fTime = projectile.GetGlobalTime();
  G4double time0 = fTime;
 
  G4double muBindingEnergy = projectile.GetBoundEnergy();
 
  G4int Z = targetNucleus.GetZ_asInt(); 
  G4int A = targetNucleus.GetA_asInt();
  G4double massA = G4NucleiProperties::GetNuclearMass(A, Z);
 
  /*
  G4cout << "G4MuMinusCapturePrecompound::ApplyYourself: Emu= "
     << muBindingEnergy << G4endl;
  */
  // Energy on K-shell
  G4double muEnergy = fMuMass + muBindingEnergy;
  G4double muMom =std::sqrt(muBindingEnergy*(muBindingEnergy + 2.0*fMuMass));
  G4double availableEnergy = massA + fMuMass - muBindingEnergy;
  G4double residualMass = G4NucleiProperties::GetNuclearMass(A, Z - 1);
 
  G4ThreeVector vmu = muMom*G4RandomDirection();
  G4LorentzVector aMuMom(vmu, muEnergy);
 
  const G4double nenergy = keV;    
 
  // p or 3He as a target 
  // two body reaction mu- + A(Z,A) -> nuMu + A(Z-1,A)
  if((1 == Z && 1 == A) || (2 == Z && 3 == A)) {
 
    const G4ParticleDefinition* pd = 0;
    if(1 == Z) { pd = fNeutron; }
    else { pd = G4Triton::Triton(); }
 
    //
    //  Computation in assumption of CM reaction
    //  
    G4double e = 0.5*(availableEnergy - 
              residualMass*residualMass/availableEnergy);
 
    G4ThreeVector nudir = G4RandomDirection();
    AddNewParticle(G4NeutrinoMu::NeutrinoMu(), nudir, e);
    nudir *= -1.0;
    AddNewParticle(pd, nudir, availableEnergy - e - residualMass);
 
  // d or 4He as a target 
  // three body reaction mu- + A(Z,A) -> nuMu + n + A(Z-1,A)
  // extra neutron produced at rest
  } else if((1 == Z && 2 == A) || (2 == Z && 4 == A)) {
 
    const G4ParticleDefinition* pd = 0;
    if(1 == Z) { pd = fNeutron; }
    else { pd = G4Triton::Triton(); }
 
    availableEnergy -= neutron_mass_c2 - nenergy;
    residualMass = pd->GetPDGMass();
 
    //
    //  Computation in assumption of CM reaction
    //  
    G4double e = 0.5*(availableEnergy - 
              residualMass*residualMass/availableEnergy);
 
    G4ThreeVector nudir = G4RandomDirection();
    AddNewParticle(G4NeutrinoMu::NeutrinoMu(), nudir, e);
    nudir *= -1.0;
    AddNewParticle(pd, nudir, availableEnergy - e - residualMass);
 
    // extra low-energy neutron
    nudir = G4RandomDirection();
    AddNewParticle(fNeutron, nudir, nenergy);
 
  } else {
    // sample mu- + p -> nuMu + n reaction in CM of muonic atom
 
    // nucleus
    G4LorentzVector momInitial(0.0,0.0,0.0,availableEnergy);
    G4LorentzVector momResidual, momNu;
 
    // pick random proton inside nucleus 
    G4double eEx;
    fNucleus.Init(A, Z);
    const std::vector<G4Nucleon>& nucleons= fNucleus.GetNucleons();
    const G4ParticleDefinition* pDef;
 
    G4int reentryCount = 0;
  
    do {
      ++reentryCount;
      G4int index = 0;
      do {
    index=G4int(A*G4UniformRand());
    pDef = nucleons[index].GetDefinition();
      } while(pDef != fProton);
      G4LorentzVector momP = nucleons[index].Get4Momentum();
 
      // Get CMS kinematics
      G4LorentzVector theCMS = momP + aMuMom;
      G4ThreeVector bst = theCMS.boostVector();
 
      G4double Ecms = theCMS.mag();
      G4double Enu  = 0.5*(Ecms - neutron_mass_c2*neutron_mass_c2/Ecms);
      eEx = 0.0;
 
      if(Enu > 0.0) {
    // make the nu, and transform to lab;
    momNu.set(Enu*G4RandomDirection(), Enu);
 
    // nu in lab.
    momNu.boost(bst);
    momResidual = momInitial - momNu;
    eEx = momResidual.mag() - residualMass;
        if(eEx < 0.0 && eEx + nenergy >= 0.0) {
          momResidual.set(0.0, 0.0, 0.0, residualMass);
          eEx = 0.0;
    }
      }
      // in the case of many iterations stop the loop
      // with zero excitation energy
      if(reentryCount > 100 && eEx < 0.0) {
    G4ExceptionDescription ed;
    ed << "Call for " << GetModelName() << G4endl;
    ed << "Target  Z= " << Z  
       << "  A= " << A << "  Eex(MeV)= " << eEx/MeV << G4endl;
    ed << " ApplyYourself does not completed after 100 attempts -"
       << " excitation energy is set to zero";
    G4Exception("G4MuMinusCapturePrecompound::ApplyYourself", "had006", 
            JustWarning, ed);
    momResidual.set(0.0, 0.0, 0.0, residualMass);
    eEx = 0.0;
      }
      // Loop checking, 06-Aug-2015, Vladimir Ivanchenko
    } while(eEx <= 0.0);
 
    G4ThreeVector dir = momNu.vect().unit();
    AddNewParticle(G4NeutrinoMu::NeutrinoMu(), dir, momNu.e());
 
    G4Fragment initialState(A, Z-1, momResidual);
    initialState.SetNumberOfExcitedParticle(2,0);
    initialState.SetNumberOfHoles(1,1);
 
    // decay time for pre-compound/de-excitation starts from zero
    G4ReactionProductVector* rpv = fPreCompound->DeExcite(initialState);
    size_t n = rpv->size();
    for(size_t i=0; i<n; ++i) {
      G4ReactionProduct* rp = (*rpv)[i];
 
      // reaction time
      fTime = time0 + rp->GetTOF();
      G4ThreeVector direction = rp->GetMomentum().unit();
      AddNewParticle(rp->GetDefinition(), direction, rp->GetKineticEnergy());
      delete rp;
    }
    delete rpv;
  } 
  if(verboseLevel > 1)
    G4cout << "G4MuMinusCapturePrecompound::ApplyYourself:  Nsec= " 
       << result.GetNumberOfSecondaries() 
       <<" E0(MeV)= " <<availableEnergy/MeV
       <<" Mres(GeV)= " <<residualMass/GeV
       <<G4endl;
 
  return &result;
}

◆ ModelDescription()

void G4MuMinusCapturePrecompound::ModelDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4HadronicInteraction.

Definition at line 262 of file G4MuMinusCapturePrecompound.cc.

{
  outFile << "Sampling of mu- capture by atomic nucleus from K-shell"
      << " mesoatom orbit.\n"
      << "Primary reaction mu- + p -> n + neutrino, neutron providing\n"
      << "  initial excitation of the target nucleus and PreCompound"
      << " model samples final state\n";
}

The documentation for this class was generated from the following files:

Public Member Functions