ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PenelopeBremsstrahlungModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4PenelopeBremsstrahlungModel.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 // Author: Luciano Pandola
28 //
29 // History:
30 // -----------
31 // 23 Aug 2010 L. Pandola 1st implementation.
32 // 24 May 2011 L. Pandola Renamed (make v2008 as default Penelope)
33 // 02 Oct 2013 L. Pandola Migrated to multi-thread
34 //
35 // -------------------------------------------------------------------
36 //
37 // Class description:
38 // Low Energy Electromagnetic Physics, e+ and e- bremsstrahlung
39 // with Penelope Model, version 2008
40 // -------------------------------------------------------------------
41 
42 #ifndef G4PENELOPEBREMSSTRAHLUNGMODEL_HH
43 #define G4PENELOPEBREMSSTRAHLUNGMODEL_HH 1
44 
45 #include "globals.hh"
46 #include "G4VEmModel.hh"
47 #include "G4DataVector.hh"
49 
50 
52 class G4PhysicsLogVector;
54 class G4DynamicParticle;
56 class G4Material;
61 
63 {
64 
65 public:
66 
68  const G4String& processName ="PenBrem");
69 
71 
72  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
73  virtual void InitialiseLocal(const G4ParticleDefinition*,
74  G4VEmModel*);
75 
76  //DUMMY METHOD
77  virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition* theParticle,
78  G4double kinEnergy,
79  G4double Z,
80  G4double A=0,
81  G4double cut=0,
83 
84 
86  const G4ParticleDefinition* theParticle,
87  G4double kineticEnergy,
88  G4double cutEnergy,
89  G4double maxEnergy = DBL_MAX);
90 
91 
92  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
93  const G4MaterialCutsCouple*,
94  const G4DynamicParticle*,
95  G4double tmin,
96  G4double maxEnergy);
97 
99  const G4ParticleDefinition*,
100  G4double kineticEnergy,
101  G4double cutEnergy);
102 
103  // Min cut in kinetic energy allowed by the model
105  const G4MaterialCutsCouple*);
106 
109 
110 protected:
113 
114 private:
115  void ClearTables();
116 
119 
121  const G4Material*,G4double cut);
122  void SetParticle(const G4ParticleDefinition*);
123 
124  //Intrinsic energy limits of the model: cannot be extended by the parent process
127 
129 
131 
133 
134  //
135  //Members to handle and store cross section tables
136  void BuildXSTable(const G4Material* material,G4double cut);
137 
138  //This is the main energy grid
140  size_t nBins;
141  //G4PenelopeCrossSection takes care of the logs
142  std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*> *XSTableElectron;
143  std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*> *XSTablePositron;
145 
146  //Helpers
149 
150  //Used only for G4EmCalculator and Unit Tests
152 
153 };
154 
155 #endif
156