Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4EmCaptureCascade.hh
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27//-----------------------------------------------------------------------------
28//
29// GEANT4 Class header file
30//
31// File name: G4EmCaptureCascade
32//
33// Author: V.Ivanchenko (Vladimir.Ivantchenko@cern.ch)
34//
35// Creation date: 22 April 2012 on base of G4MuMinusCaptureCascade
36//
37// Class Description:
38//
39// Simulation of electromagnetic cascade from capture level to K-shell
40// of the mesonic atom
41//
42// Probabilities of gamma and Auger transitions from
43// N.C.Mukhopadhyay Phys. Rep. 30 (1977) 1.
44//
45//-----------------------------------------------------------------------------
46//
47// Modifications:
48//
49//-----------------------------------------------------------------------------
50
51#ifndef G4EmCaptureCascade_h
52#define G4EmCaptureCascade_h 1
53
54#include "globals.hh"
55#include "G4Nucleus.hh"
56#include "G4Track.hh"
57#include "G4HadProjectile.hh"
58#include "G4HadSecondary.hh"
59#include "G4HadFinalState.hh"
61#include "G4RandomDirection.hh"
63#include "G4DynamicParticle.hh"
64#include "G4ThreeVector.hh"
65
67{
68public:
69
70 explicit G4EmCaptureCascade();
71
72 virtual ~G4EmCaptureCascade();
73
74 virtual G4HadFinalState* ApplyYourself(const G4HadProjectile &aTrack,
75 G4Nucleus & targetNucleus );
76
77 virtual void ModelDescription(std::ostream& outFile) const;
78
79private:
80
81 inline void AddNewParticle(G4ParticleDefinition* aParticle,
82 G4double kinEnergy);
83
84 // hide assignment operator as private
85 G4EmCaptureCascade& operator=(const G4EmCaptureCascade &right) = delete;
86 G4EmCaptureCascade(const G4EmCaptureCascade& ) = delete;
87
88 G4HadFinalState result;
89 G4ParticleDefinition* theElectron;
90 G4ParticleDefinition* theGamma;
91 G4double fMuMass;
92 G4double fTime;
93 G4double fLevelEnergy[14];
94 G4double fKLevelEnergy[93];
95};
96
97//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
98
99inline void
100G4EmCaptureCascade::AddNewParticle(G4ParticleDefinition* aParticle,
101 G4double kinEnergy)
102{
103 G4DynamicParticle* dp = new G4DynamicParticle(aParticle,
105 kinEnergy);
106 G4HadSecondary hs(dp);
107 hs.SetTime(fTime);
108 result.AddSecondary(hs);
109}
110
111#endif
112
113
G4ThreeVector G4RandomDirection()
double G4double
Definition: G4Types.hh:83
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
virtual void ModelDescription(std::ostream &outFile) const
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)