Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4IonCoulombScatteringModel.hh
Go to the documentation of this file.
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 // G4IonCoulombScatteringModel.hh
27 // -------------------------------------------------------------------
28 //
29 // GEANT4 Class header file
30 //
31 // File name: G4IonCoulombScatteringModel
32 //
33 // Author: Cristina Consolandi
34 //
35 // Creation date: 05.10.2010 from G4eCoulombScatteringModel
36 // & G4CoulombScatteringModel
37 //
38 //
39 // Class Description:
40 // Single Scattering Model for
41 // for protons, alpha and heavy Ions
42 //
43 // Reference:
44 // M.J. Boschini et al. "Nuclear and Non-Ionizing Energy-Loss
45 // for Coulomb ScatteredParticles from Low Energy up to Relativistic
46 // Regime in Space Radiation Environment"
47 // Accepted for publication in the Proceedings of the ICATPP Conference
48 // on Cosmic Rays for Particle and Astroparticle Physics, Villa Olmo, 7-8
49 // October, 2010, to be published by World Scientific (Singapore).
50 //
51 // Available for downloading at:
52 // http://arxiv.org/abs/1011.4822
53 //
54 // -------------------------------------------------------------------
55 //
56 
57 #ifndef G4IonCoulombScatteringModel_h
58 #define G4IonCoulombScatteringModel_h 1
59 
60 #include "G4VEmModel.hh"
61 #include "globals.hh"
62 #include "G4NistManager.hh"
64 
65 #include <vector>
66 
69 class G4IonTable;
70 
72 {
73 public:
74 
75  explicit G4IonCoulombScatteringModel(const G4String& nam =
76  "IonCoulombScattering");
77 
79 
80  virtual void Initialise(const G4ParticleDefinition*,
81  const G4DataVector&) final;
82 
84  const G4ParticleDefinition*,
85  G4double kinEnergy,
86  G4double Z,
87  G4double A,
88  G4double cut,
89  G4double emax) final;
90 
91  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
92  const G4MaterialCutsCouple*,
93  const G4DynamicParticle*,
94  G4double tmin,
95  G4double maxEnergy) final;
96 
97 
98 
99  inline void SetRecoilThreshold(G4double eth);
100 
101  inline void SetHeavyIonCorr(G4int b);
102 
103  inline G4int GetHeavyIonCorr();
104 
105  //protected:
106 
107 private:
108 
109  inline void DefineMaterial(const G4MaterialCutsCouple*);
110 
111  inline void SetupParticle(const G4ParticleDefinition*);
112 
113  // hide assignment operator
114  G4IonCoulombScatteringModel & operator=
115  (const G4IonCoulombScatteringModel &right) = delete;
117 
118  //protected:
119 
120  G4IonTable* theIonTable;
121  G4ParticleChangeForGamma* fParticleChange;
122  G4NistManager* fNistManager;
123  G4IonCoulombCrossSection* ioncross;
124 
125  const std::vector<G4double>* pCuts;
126  const G4MaterialCutsCouple* currentCouple;
127  const G4Material* currentMaterial;
128  const G4Element* currentElement;
129  G4int currentMaterialIndex;
130 
131  G4int heavycorr;
132 
133  G4double cosThetaMin;
134  G4double recoilThreshold;
135 
136  // projectile
137  const G4ParticleDefinition* particle;
138  const G4ParticleDefinition* theProton;
139  G4double mass;
140 
141 };
142 
143 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
144 
145 inline void
146 G4IonCoulombScatteringModel::DefineMaterial(const G4MaterialCutsCouple* cup)
147 {
148  if(cup != currentCouple) {
149  currentCouple = cup;
150  currentMaterial = cup->GetMaterial();
151  currentMaterialIndex = currentCouple->GetIndex();
152  }
153 }
154 
155 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
156 
157 inline
158 void G4IonCoulombScatteringModel::SetupParticle(const G4ParticleDefinition* p)
159 {
160  if(p != particle) {
161  particle = p;
162  mass = particle->GetPDGMass();
163  ioncross->SetupParticle(p);
164  }
165 }
166 
167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
168 
170 {
171  recoilThreshold = eth;
172 }
173 
174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
175 
177 {
178  heavycorr = b;
179 }
180 
181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
182 
184 {
185  return heavycorr;
186 }
187 
188 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
189 
190 #endif
const char * p
Definition: xmltok.h:285
int G4int
Definition: G4Types.hh:78
G4IonCoulombScatteringModel(const G4String &nam="IonCoulombScattering")
tuple b
Definition: test.py:12
double A(double temperature)
void SetupParticle(const G4ParticleDefinition *)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax) final
static const G4double emax
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) final
G4double GetPDGMass() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) final
double G4double
Definition: G4Types.hh:76
const G4Material * GetMaterial() const