ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4DNARuddIonisationExtendedModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4DNARuddIonisationExtendedModel.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 #ifndef G4DNARuddIonisationExtendedModel_h
29 #define G4DNARuddIonisationExtendedModel_h 1
30 
31 #include "G4VEmModel.hh"
32 #include "G4ParticleChangeForGamma.hh"
33 #include "G4ProductionCutsTable.hh"
34 
35 #include "G4DNAGenericIonsManager.hh"
36 #include "G4DNACrossSectionDataSet.hh"
37 #include "G4Electron.hh"
38 #include "G4Proton.hh"
39 #include "G4LogLogInterpolation.hh"
40 
41 #include "G4DNAWaterIonisationStructure.hh"
42 #include "G4VAtomDeexcitation.hh"
43 #include "G4NistManager.hh"
44 
45 class G4DNARuddIonisationExtendedModel : public G4VEmModel
46 {
47 
48 public:
49 
50  G4DNARuddIonisationExtendedModel(const G4ParticleDefinition* p = 0,
51  const G4String& nam = "DNARuddIonisationExtendedModel");
52 
53  virtual ~G4DNARuddIonisationExtendedModel();
54 
55  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
56 
57  virtual G4double CrossSectionPerVolume( const G4Material* material,
58  const G4ParticleDefinition* p,
59  G4double ekin,
60  G4double emin,
61  G4double emax);
62 
63  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
64  const G4MaterialCutsCouple*,
65  const G4DynamicParticle*,
66  G4double tmin,
67  G4double maxEnergy);
68 
69  inline void SelectStationary(G4bool input);
70 
71 protected:
72 
73  G4ParticleChangeForGamma* fParticleChangeForGamma;
74 
75 private:
76 
77  G4bool statCode;
78 
79  // Water density table
80  const std::vector<G4double>* fpWaterDensity;
81 
82  //deexcitation manager to produce fluo photns and e-
83  G4VAtomDeexcitation* fAtomDeexcitation;
84 
85  std::map<G4String,G4double,std::less<G4String> > lowEnergyLimit;
86  std::map<G4String,G4double,std::less<G4String> > highEnergyLimit;
87 
88  // ZF: 26-10-2010
89  std::map<G4double,G4double> lowEnergyLimitForA, lowEnergyLimitOfModelForA, killBelowEnergyForA;
90 
91  G4bool isInitialised;
92  G4int verboseLevel;
93 
94  // Cross section
95 
96  typedef std::map<G4String,G4String,std::less<G4String> > MapFile;
97  MapFile tableFile;
98 
99  typedef std::map<G4String,G4DNACrossSectionDataSet*,std::less<G4String> > MapData;
100  MapData tableData;
101 
102  // Final state
103 
104  G4DNAWaterIonisationStructure waterStructure;
105 
106  G4double RandomizeEjectedElectronEnergy(G4ParticleDefinition* particleDefinition,
107  G4double incomingParticleEnergy,
108  G4int shell);
109 
110  // SI: Not necessary anymore since we now use interface to angle generator
111  /*
112  void RandomizeEjectedElectronDirection(G4ParticleDefinition* particleDefinition,
113  G4double incomingParticleEnergy,
114  G4double outgoingParticleEnergy,
115  G4double & cosTheta,
116  G4double & phi, G4int shell);
117  */
118 
119  G4double RejectionFunction(G4ParticleDefinition* particle,
120  G4double k,
121  G4double proposed_ws,
122  G4int ionizationLevelIndex);
123 
124  G4double ProposedSampledEnergy(G4ParticleDefinition* particle,
125  G4double k,
126  G4int ionizationLevelIndex);
127 
128  G4double CorrectionFactor(G4ParticleDefinition* particleDefinition, G4double k, G4int shell);
129 
130  G4double S_1s(G4double t,
131  G4double energyTransferred,
132  G4double slaterEffectiveChg,
133  G4double shellNumber);
134 
135  G4double S_2s(G4double t,
136  G4double energyTransferred,
137  G4double slaterEffectiveChg,
138  G4double shellNumber);
139 
140 
141  G4double S_2p(G4double t,
142  G4double energyTransferred,
143  G4double slaterEffectiveChg,
144  G4double shellNumber);
145 
146  G4double R(G4double t,
147  G4double energyTransferred,
148  G4double slaterEffectiveChg,
149  G4double shellNumber) ;
150 
151  G4double slaterEffectiveCharge[3];
152  G4double sCoefficient[3];
153 
154  // Partial cross section
155 
156  G4double PartialCrossSection(const G4Track& track);
157 
158  G4double Sum(G4double energy, const G4String& particle);
159 
160  G4int RandomSelect(G4double energy,const G4String& particle );
161 
162  //
163 
164  G4DNARuddIonisationExtendedModel & operator=(const G4DNARuddIonisationExtendedModel &right);
165  G4DNARuddIonisationExtendedModel(const G4DNARuddIonisationExtendedModel&);
166 
167 };
168 
169 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
170 
171 inline void G4DNARuddIonisationExtendedModel::SelectStationary (G4bool input)
172 {
173  statCode = input;
174 }
175 
176 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
177 
178 #endif