ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4MuBremsstrahlungModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4MuBremsstrahlungModel.hh
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 //
32 // File name: G4MuBremsstrahlungModel
33 //
34 // Author: Vladimir Ivanchenko on base of Laszlo Urban code
35 //
36 // Creation date: 18.05.2002
37 //
38 // Modifications:
39 //
40 // 23-12-02 Change interface in order to move to cut per region (V.Ivanchenko)
41 // 27-01-03 Make models region aware (V.Ivanchenko)
42 // 13-02-03 Add name (V.Ivanchenko)
43 // 10-02-04 Add lowestKinEnergy (V.Ivanchenko)
44 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
45 // 13-02-06 Add ComputeCrossSectionPerAtom (mma)
46 // 11-10-07 Add ignoreCut flag (V.Ivanchenko)
47 // 28-02-08 Reorganized protected methods and members (V.Ivanchenko)
48 // 06-03-08 Remove obsolete methods and members (V.Ivanchenko)
49 // 31-05-13 Use element selectors instead of local data (V.Ivanchenko)
50 //
51 
52 //
53 // Class Description:
54 //
55 // Implementation of bremssrahlung by muons
56 
57 // -------------------------------------------------------------------
58 //
59 
60 #ifndef G4MuBremsstrahlungModel_h
61 #define G4MuBremsstrahlungModel_h 1
62 
63 #include <CLHEP/Units/PhysicalConstants.h>
64 
65 #include "G4VEmModel.hh"
66 #include "G4NistManager.hh"
67 
68 class G4Element;
69 class G4ParticleChangeForLoss;
70 
71 class G4MuBremsstrahlungModel : public G4VEmModel
72 {
73 
74 public:
75 
76  explicit G4MuBremsstrahlungModel(const G4ParticleDefinition* p = nullptr,
77  const G4String& nam = "MuBrem");
78 
79  ~G4MuBremsstrahlungModel() = default;
80 
81  virtual void Initialise(const G4ParticleDefinition*,
82  const G4DataVector&) override;
83 
84  virtual void InitialiseLocal(const G4ParticleDefinition*,
85  G4VEmModel* masterModel) override;
86 
87  virtual G4double MinEnergyCut(const G4ParticleDefinition*,
88  const G4MaterialCutsCouple*) override;
89 
90  virtual G4double ComputeCrossSectionPerAtom(
91  const G4ParticleDefinition*,
92  G4double kineticEnergy,
93  G4double Z, G4double A,
94  G4double cutEnergy,
95  G4double maxEnergy) override;
96 
97  virtual G4double ComputeDEDXPerVolume(const G4Material*,
98  const G4ParticleDefinition*,
99  G4double kineticEnergy,
100  G4double cutEnergy) override;
101 
102  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
103  const G4MaterialCutsCouple*,
104  const G4DynamicParticle*,
105  G4double tmin,
106  G4double maxEnergy) override;
107 
108  inline void SetLowestKineticEnergy(G4double e);
109 
110  virtual G4double MinPrimaryEnergy(const G4Material*,
111  const G4ParticleDefinition*,
112  G4double) override;
113 
114 protected:
115 
116  G4double ComputMuBremLoss(G4double Z, G4double tkin, G4double cut);
117 
118  G4double ComputeMicroscopicCrossSection(G4double tkin,
119  G4double Z,
120  G4double cut);
121 
122  virtual G4double ComputeDMicroscopicCrossSection(G4double tkin,
123  G4double Z,
124  G4double gammaEnergy);
125 
126  inline void SetParticle(const G4ParticleDefinition*);
127 
128 private:
129 
130  // hide assignment operator
131  G4MuBremsstrahlungModel &
132  operator=(const G4MuBremsstrahlungModel &right) = delete;
133  G4MuBremsstrahlungModel(const G4MuBremsstrahlungModel&) = delete;
134 
135 protected:
136 
137  const G4ParticleDefinition* particle;
138  G4NistManager* nist;
139  G4double mass;
140  G4double rmass;
141  G4double cc;
142  G4double coeff;
143  G4double sqrte;
144  G4double bh;
145  G4double bh1;
146  G4double btf;
147  G4double btf1;
148 
149  G4ParticleDefinition* theGamma;
150  G4ParticleChangeForLoss* fParticleChange;
151 
152  G4double lowestKinEnergy;
153  G4double minThreshold;
154 
155  static const G4double xgi[6];
156  static const G4double wgi[6];
157 
158  static G4double fDN[93];
159 };
160 
161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
162 
163 inline void G4MuBremsstrahlungModel::SetLowestKineticEnergy(G4double e)
164 {
165  lowestKinEnergy = e;
166 }
167 
168 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
169 
170 inline
171 void G4MuBremsstrahlungModel::SetParticle(const G4ParticleDefinition* p)
172 {
173  if(!particle) {
174  particle = p;
175  mass = particle->GetPDGMass();
176  rmass=mass/CLHEP::electron_mass_c2 ;
177  cc=CLHEP::classic_electr_radius/rmass ;
178  coeff= 16.*CLHEP::fine_structure_const*cc*cc/3. ;
179  }
180 }
181 
182 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
183 
184 #endif