Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4MuMinusCapturePrecompound Class Reference

#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 G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)=0
 
virtual G4double SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
 
virtual G4bool IsApplicable (const G4HadProjectile &, G4Nucleus &)
 
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)
 
const G4HadronicInteraction * GetMyPointer () const
 
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
 
G4bool operator== (const G4HadronicInteraction &right) const
 
G4bool operator!= (const G4HadronicInteraction &right) const
 
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels () const
 
virtual std::pair< G4double, G4double > GetEnergyMomentumCheckLevels () const
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
virtual void ModelDescription (std::ostream &outFile) const
 

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 65 of file G4MuMinusCapturePrecompound.hh.

Constructor & Destructor Documentation

◆ G4MuMinusCapturePrecompound()

G4MuMinusCapturePrecompound::G4MuMinusCapturePrecompound ( G4VPreCompoundModel *  ptr = 0)

Definition at line 64 of file G4MuMinusCapturePrecompound.cc.

66 : G4HadronicInteraction("muMinusNuclearCapture")
67{
68 fMuMass = G4MuonMinus::MuonMinus()->GetPDGMass();
69 fProton = G4Proton::Proton();
70 fNeutron = G4Neutron::Neutron();
71 fThreshold = 10*MeV;
72 fPreCompound = ptr;
73 if(!ptr) {
74 G4HadronicInteraction* p =
75 G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
76 ptr = static_cast<G4VPreCompoundModel*>(p);
77 fPreCompound = ptr;
78 if(!ptr) { fPreCompound = new G4PreCompoundModel(); }
79 }
80}
G4HadronicInteractionRegistry::FindModel
G4HadronicInteraction * FindModel(const G4String &name)
Definition: G4HadronicInteractionRegistry.cc:104
G4HadronicInteractionRegistry::Instance
static G4HadronicInteractionRegistry * Instance()
Definition: G4HadronicInteractionRegistry.cc:36
G4MuonMinus::MuonMinus
static G4MuonMinus * MuonMinus()
Definition: G4MuonMinus.cc:100
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93

◆ ~G4MuMinusCapturePrecompound()

G4MuMinusCapturePrecompound::~G4MuMinusCapturePrecompound ( )

Definition at line 84 of file G4MuMinusCapturePrecompound.cc.

85{
86 result.Clear();
87}

Member Function Documentation

◆ ApplyYourself()

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

Implements G4HadronicInteraction.

Definition at line 92 of file G4MuMinusCapturePrecompound.cc.

94{
95 result.Clear();
96 result.SetStatusChange(stopAndKill);
97 fTime = projectile.GetGlobalTime();
98 G4double time0 = fTime;
99
100 G4double muBindingEnergy = projectile.GetBoundEnergy();
101
102 G4int Z = targetNucleus.GetZ_asInt();
103 G4int A = targetNucleus.GetA_asInt();
104 G4double massA = G4NucleiProperties::GetNuclearMass(A, Z);
105
106 /*
107 G4cout << "G4MuMinusCapturePrecompound::ApplyYourself: Emu= "
108 << muBindingEnergy << G4endl;
109 */
110 // Energy on K-shell
111 G4double muEnergy = fMuMass + muBindingEnergy;
112 G4double muMom = std::sqrt(muBindingEnergy*(muBindingEnergy + 2.0*fMuMass));
113 G4double availableEnergy = massA + fMuMass - muBindingEnergy;
114 G4double residualMass = G4NucleiProperties::GetNuclearMass(A, Z - 1);
115
116 G4ThreeVector vmu = muMom*G4RandomDirection();
117 G4LorentzVector aMuMom(vmu, muEnergy);
118
119 // p or 3He as a target
120 // two body reaction mu- + A(Z,A) -> nuMu + A(Z-1,A)
121 if((1 == Z && 1 == A) || (2 == Z && 3 == A)) {
122
123 G4ParticleDefinition* pd = 0;
124 if(1 == Z) { pd = fNeutron; }
125 else { pd = G4Triton::Triton(); }
126
127 //
128 // Computation in assumption of CM reaction
129 //
130 G4double e = 0.5*(availableEnergy -
131 residualMass*residualMass/availableEnergy);
132
133 G4ThreeVector nudir = G4RandomDirection();
134 AddNewParticle(G4NeutrinoMu::NeutrinoMu(), nudir, e);
135 nudir *= -1.0;
136 AddNewParticle(pd, nudir, availableEnergy - e - residualMass);
137
138
139 } else {
140 // sample mu- + p -> nuMu + n reaction in CM of muonic atom
141
142 // muon
143//
144// NOTE by K.Genser and J.Yarba:
145// The code below isn't working because emu always turns smaller than fMuMass
146// For this reason the sqrt is producing a NaN
147//
148// G4double emu = (availableEnergy*availableEnergy - massA*massA
149// + fMuMass*fMuMass)/(2*availableEnergy);
150// G4ThreeVector mudir = G4RandomDirection();
151// G4LorentzVector momMuon(std::sqrt(emu*emu - fMuMass*fMuMass)*mudir, emu);
152
153 // nucleus
154 G4LorentzVector momInitial(0.0,0.0,0.0,availableEnergy);
155 G4LorentzVector momResidual, momNu;
156
157 // pick random proton inside nucleus
158 G4double eEx;
159 fNucleus.Init(A, Z);
160 const std::vector<G4Nucleon>& nucleons= fNucleus.GetNucleons();
161 G4ParticleDefinition* pDef;
162
163 G4int nneutrons = 1;
164 G4int reentryCount = 0;
165
166 do {
167 ++reentryCount;
168 G4int index = 0;
169 do {
170 index=G4int(A*G4UniformRand());
171 pDef = nucleons[index].GetDefinition();
172 } while(pDef != fProton);
173 G4LorentzVector momP = nucleons[index].Get4Momentum();
174
175 // Get CMS kinematics
176 G4LorentzVector theCMS = momP + aMuMom;
177 G4ThreeVector bst = theCMS.boostVector();
178
179 G4double Ecms = theCMS.mag();
180 G4double Enu = 0.5*(Ecms - neutron_mass_c2*neutron_mass_c2/Ecms);
181 eEx = 0.0;
182
183 if(Enu > 0.0) {
184 // make the nu, and transform to lab;
185 momNu.set(Enu*G4RandomDirection(), Enu);
186
187 // nu in lab.
188 momNu.boost(bst);
189 momResidual = momInitial - momNu;
190 eEx = momResidual.mag() - residualMass;
191
192 // release neutron
193
194 if(eEx > 0.0) {
195 G4double eth = residualMass - massA + fThreshold + 2*neutron_mass_c2;
196 if(Ecms - Enu > eth) {
197 theCMS -= momNu;
198 G4double ekin = theCMS.e() - eth;
199 G4ThreeVector dir = theCMS.vect().unit();
200 AddNewParticle(fNeutron, dir, ekin);
201 momResidual -=
202 result.GetSecondary(0)->GetParticle()->Get4Momentum();
203 --Z;
204 --A;
205 residualMass = G4NucleiProperties::GetNuclearMass(A, Z);
206 nneutrons = 0;
207 }
208 }
209 }
210 if(Enu <= 0.0 && eEx <= 0.0 && reentryCount > 100) {
211 G4ExceptionDescription ed;
212 ed << "Call for " << GetModelName() << G4endl;
213 ed << "Target Z= " << Z
214 << " A= " << A << G4endl;
215 ed << " ApplyYourself does not completed after 100 attempts" << G4endl;
216 G4Exception("G4MuMinusCapturePrecompound::AtRestDoIt", "had006",
217 FatalException, ed);
218 }
219 } while(eEx <= 0.0);
220
221 G4ThreeVector dir = momNu.vect().unit();
222 AddNewParticle(G4NeutrinoMu::NeutrinoMu(), dir, momNu.e());
223
224 G4Fragment initialState(A, Z, momResidual);
225 initialState.SetNumberOfExcitedParticle(nneutrons,0);
226 initialState.SetNumberOfHoles(1,1);
227
228 // decay time for pre-compound/de-excitation starts from zero
229 G4ReactionProductVector* rpv = fPreCompound->DeExcite(initialState);
230 size_t n = rpv->size();
231 for(size_t i=0; i<n; ++i) {
232 G4ReactionProduct* rp = (*rpv)[i];
233
234 // reaction time
235 fTime = time0 + rp->GetTOF();
236 G4ThreeVector direction = rp->GetMomentum().unit();
237 AddNewParticle(rp->GetDefinition(), direction, rp->GetKineticEnergy());
238 delete rp;
239 }
240 delete rpv;
241 }
242 if(verboseLevel > 1)
243 G4cout << "G4MuMinusCapturePrecompound::ApplyYourself: Nsec= "
244 << result.GetNumberOfSecondaries()
245 <<" E0(MeV)= " <<availableEnergy/MeV
246 <<" Mres(GeV)= " <<residualMass/GeV
247 <<G4endl;
248
249 return &result;
250}
FatalException
@ FatalException
Definition: G4ExceptionSeverity.hh:60
stopAndKill
@ stopAndKill
Definition: G4HadFinalState.hh:42
G4RandomDirection
G4ThreeVector G4RandomDirection()
Definition: G4RandomDirection.hh:55
G4ReactionProductVector
std::vector< G4ReactionProduct * > G4ReactionProductVector
Definition: G4ReactionProductVector.hh:43
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:53
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::HepLorentzVector::boostVector
Hep3Vector boostVector() const
Definition: LorentzVector.cc:177
CLHEP::HepLorentzVector::boost
HepLorentzVector & boost(double, double, double)
Definition: LorentzVector.cc:59
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
CLHEP::HepLorentzVector::set
void set(double x, double y, double z, double t)
G4DynamicParticle::Get4Momentum
G4LorentzVector Get4Momentum() const
G4Fancy3DNucleus::GetNucleons
const std::vector< G4Nucleon > & GetNucleons()
Definition: G4Fancy3DNucleus.cc:130
G4Fancy3DNucleus::Init
void Init(G4int theA, G4int theZ)
Definition: G4Fancy3DNucleus.cc:76
G4HadFinalState::SetStatusChange
void SetStatusChange(G4HadFinalStateStatus aS)
Definition: G4HadFinalState.cc:81
G4HadFinalState::GetNumberOfSecondaries
G4int GetNumberOfSecondaries() const
Definition: G4HadFinalState.cc:40
G4HadFinalState::GetSecondary
G4HadSecondary * GetSecondary(size_t i)
Definition: G4HadFinalState.cc:109
G4HadSecondary::GetParticle
G4DynamicParticle * GetParticle()
Definition: G4HadSecondary.hh:39
G4HadronicInteraction::GetModelName
const G4String & GetModelName() const
Definition: G4HadronicInteraction.hh:119
G4NeutrinoMu::NeutrinoMu
static G4NeutrinoMu * NeutrinoMu()
Definition: G4NeutrinoMu.cc:85
G4NucleiProperties::GetNuclearMass
static G4double GetNuclearMass(const G4double A, const G4double Z)
Definition: G4NucleiProperties.cc:53
G4Nucleus::GetA_asInt
G4int GetA_asInt() const
Definition: G4Nucleus.hh:109
G4Nucleus::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Nucleus.hh:115
G4ReactionProduct::GetKineticEnergy
G4double GetKineticEnergy() const
Definition: G4ReactionProduct.hh:135
G4ReactionProduct::GetMomentum
G4ThreeVector GetMomentum() const
Definition: G4ReactionProduct.hh:120
G4ReactionProduct::GetTOF
G4double GetTOF() const
Definition: G4ReactionProduct.hh:153
G4ReactionProduct::GetDefinition
G4ParticleDefinition * GetDefinition() const
Definition: G4ReactionProduct.hh:104
G4Triton::Triton
static G4Triton * Triton()
Definition: G4Triton.cc:95
G4VPreCompoundModel::DeExcite
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)=0
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
G4ExceptionDescription
std::ostringstream G4ExceptionDescription
Definition: globals.hh:76

◆ ModelDescription()

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

Reimplemented from G4HadronicInteraction.

Definition at line 254 of file G4MuMinusCapturePrecompound.cc.

255{
256 outFile << "Sampling of mu- capture by atomic nucleus from K-shell"
257 << " mesoatom orbit.\n"
258 << "Primary reaction mu- + p -> n + neutrino, neutron providing\n"
259 << " initial excitation of the target nucleus and PreCompound"
260 << " model samples final state\n";
261}
The documentation for this class was generated from the following files: