Geant4  10.02.p01
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  G4IonCoulombScatteringModel(const G4String& nam = "IonCoulombScattering");
76 
78 
79  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
80 
82  const G4ParticleDefinition*,
83  G4double kinEnergy,
84  G4double Z,
85  G4double A,
86  G4double cut,
87  G4double emax);
88 
89  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
90  const G4MaterialCutsCouple*,
91  const G4DynamicParticle*,
92  G4double tmin,
93  G4double maxEnergy);
94 
95 
96 
97  inline void SetRecoilThreshold(G4double eth);
98 
99  inline void SetHeavyIonCorr(G4int b);
100 
101  inline G4int GetHeavyIonCorr();
102 
103  //protected:
104 
105 private:
106 
107  inline void DefineMaterial(const G4MaterialCutsCouple*);
108 
109  inline void SetupParticle(const G4ParticleDefinition*);
110 
111  // hide assignment operator
114 
115  //protected:
116 
121 
122  const std::vector<G4double>* pCuts;
127 
129 
132 
133  // projectile
137 
138 };
139 
140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
141 
142 inline void
144 {
145  if(cup != currentCouple) {
146  currentCouple = cup;
147  currentMaterial = cup->GetMaterial();
149  }
150 }
151 
152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
153 
154 inline
156 {
157  if(p != particle) {
158  particle = p;
159  mass = particle->GetPDGMass();
161  }
162 }
163 
164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
165 
167 {
168  recoilThreshold = eth;
169 }
170 
171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
172 
174 {
175  heavycorr = b;
176 }
177 
178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
179 
181 {
182  return heavycorr;
183 }
184 
185 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
186 
187 #endif
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax)
const std::vector< G4double > * pCuts
int G4int
Definition: G4Types.hh:78
G4ParticleChangeForGamma * fParticleChange
G4IonCoulombScatteringModel(const G4String &nam="IonCoulombScattering")
double A(double temperature)
void SetupParticle(const G4ParticleDefinition *)
const G4ParticleDefinition * theProton
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
static const G4double emax
G4double GetPDGMass() const
G4IonCoulombScatteringModel & operator=(const G4IonCoulombScatteringModel &right)
const G4ParticleDefinition * particle
void DefineMaterial(const G4MaterialCutsCouple *)
void SetupParticle(const G4ParticleDefinition *)
const G4MaterialCutsCouple * currentCouple
double G4double
Definition: G4Types.hh:76
const G4Material * GetMaterial() const