ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PenelopeRayleighModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4PenelopeRayleighModel.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 // 03 Dec 2009 L. Pandola 1st implementation.
32 // 25 May 2011 L. Pandola Renamed (make v2008 as default Penelope)
33 // 27 Sep 2013 L. Pandola Migration to MT paradigm
34 //
35 // -------------------------------------------------------------------
36 //
37 // Class description:
38 // Low Energy Electromagnetic Physics, Rayleigh Scattering
39 // with the model from Penelope, version 2008
40 // -------------------------------------------------------------------
41 
42 #ifndef G4PENELOPERAYLEIGHMODEL_HH
43 #define G4PENELOPERAYLEIGHMODEL_HH 1
44 
45 #include "globals.hh"
46 #include "G4VEmModel.hh"
47 #include "G4DataVector.hh"
49 
51 class G4DynamicParticle;
53 class G4Material;
56 
58 {
59 
60 public:
62  const G4String& processName ="PenRayleigh");
63 
64  virtual ~G4PenelopeRayleighModel();
65 
66  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
67  virtual void InitialiseLocal(const G4ParticleDefinition*,
68  G4VEmModel *masterModel);
69 
71  G4double kinEnergy,
72  G4double Z,
73  G4double A=0,
74  G4double cut=0,
76 
77  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
78  const G4MaterialCutsCouple*,
79  const G4DynamicParticle*,
80  G4double tmin,
81  G4double maxEnergy);
82 
85 
86  //Testing purposes
87  void DumpFormFactorTable(const G4Material*);
88 
89 protected:
92 
93 private:
96 
97  void SetParticle(const G4ParticleDefinition*);
98 
99  //Intrinsic energy limits of the model: cannot be extended by
100  //the parent process
103 
106 
107  //Internal tables and manager methods
108  std::map<G4int,G4PhysicsFreeVector*> *logAtomicCrossSection;
109  std::map<G4int,G4PhysicsFreeVector*> *atomicFormFactor;
110 
111 
112  G4DataVector logQSquareGrid; //log(Q^2) grid for interpolation
113  std::map<const G4Material*,G4PhysicsFreeVector*> *logFormFactorTable; //log(Q^2) vs. log(F^2)
114 
115  G4DataVector logEnergyGridPMax; //energy grid for PMac (and originally for the x-section)
116  std::map<const G4Material*,G4PhysicsFreeVector*> *pMaxTable; //E vs. Pmax
117 
118  std::map<const G4Material*,G4PenelopeSamplingData*> *samplingTable;
119 
120  //Helper methods
121  void ReadDataFile(G4int);
122  void ClearTables();
123  void BuildFormFactorTable(const G4Material*);
124  void GetPMaxTable(const G4Material*);
125 
126  G4double GetFSquared(const G4Material*,const G4double);
128 
129  //Used only for G4EmCalculator and Unit Tests
131 
132 };
133 
134 #endif
135