ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4DNADingfelderChargeDecreaseModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4DNADingfelderChargeDecreaseModel.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 G4DNADingfelderChargeDecreaseModel_h
29 #define G4DNADingfelderChargeDecreaseModel_h 1
30 
31 #include "G4VEmModel.hh"
33 #include "G4ProductionCutsTable.hh"
34 
35 #include "G4Proton.hh"
37 #include "G4NistManager.hh"
38 
40 {
41 
42 public:
43 
45  const G4String& nam = "DNADingfelderChargeDecreaseModel");
46 
48 
49  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
50 
52  const G4ParticleDefinition* p,
53  G4double ekin,
54  G4double emin,
55  G4double emax);
56 
57  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
58  const G4MaterialCutsCouple*,
59  const G4DynamicParticle*,
60  G4double tmin,
61  G4double maxEnergy);
62 
63  inline void SelectStationary(G4bool input);
64 
65 protected:
66 
68 
69 private:
70 
72 
73  // Water density table
74  const std::vector<G4double>* fpMolWaterDensity;
75 
76  std::map<G4String,G4double,std::less<G4String> > lowEnergyLimit;
77  std::map<G4String,G4double,std::less<G4String> > highEnergyLimit;
78 
81 
82  // Partial cross section
83 
85 
87 
89 
90  G4int numberOfPartialCrossSections[3]; // 3 is the particle type index
91 
92  G4double f0[2][3];
93  G4double a0[2][3];
94  G4double a1[2][3];
95  G4double b0[2][3];
96  G4double b1[2][3];
97  G4double c0[2][3];
98  G4double d0[2][3];
99  G4double x0[2][3];
100  G4double x1[2][3];
101 
102  // Final state
103 
104  G4int NumberOfFinalStates(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
105 
106  G4ParticleDefinition* OutgoingParticleDefinition(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
107 
108  G4double WaterBindingEnergyConstant(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
109 
110  G4double OutgoingParticleBindingEnergyConstant(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
111 
112  //
113 
116 
117 };
118 
119 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
120 
122 {
123  statCode = input;
124 }
125 
126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
127 
128 #endif