Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4IonCoulombScatteringModel.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/electromagnetic/standard/include/G4IonCoulombScatteringModel.hh (Version 11.3.0) and /processes/electromagnetic/standard/include/G4IonCoulombScatteringModel.hh (Version 10.2.p2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //  G4IonCoulombScatteringModel.hh                 26 //  G4IonCoulombScatteringModel.hh
 27 // -------------------------------------------     27 // -------------------------------------------------------------------
 28 //                                                 28 //
 29 // GEANT4 Class header file                        29 // GEANT4 Class header file
 30 //                                                 30 //
 31 // File name:    G4IonCoulombScatteringModel       31 // File name:    G4IonCoulombScatteringModel 
 32 //                                                 32 //
 33 // Author:  Cristina Consolandi                    33 // Author:  Cristina Consolandi
 34 //                                                 34 //
 35 // Creation date: 05.10.2010 from G4eCoulombSc     35 // Creation date: 05.10.2010 from G4eCoulombScatteringModel 
 36 //        & G4CoulombScatteringModel               36 //        & G4CoulombScatteringModel
 37 //                                                 37 //
 38 //                                                 38 //
 39 // Class Description:                              39 // Class Description:
 40 //  Single Scattering Model for                    40 //  Single Scattering Model for
 41 //      for protons, alpha and heavy Ions          41 //      for protons, alpha and heavy Ions
 42 //                                                 42 //
 43 // Reference:                                      43 // Reference:   
 44 //  M.J. Boschini et al. "Nuclear and Non-Ioni     44 //  M.J. Boschini et al. "Nuclear and Non-Ionizing Energy-Loss 
 45 //  for Coulomb ScatteredParticles from Low En     45 //  for Coulomb ScatteredParticles from Low Energy up to Relativistic 
 46 //  Regime in Space Radiation Environment"         46 //  Regime in Space Radiation Environment"
 47 //  Accepted for publication in the Proceeding     47 //  Accepted for publication in the Proceedings of  the  ICATPP Conference
 48 //  on Cosmic Rays for Particle and Astroparti     48 //  on Cosmic Rays for Particle and Astroparticle Physics, Villa  Olmo, 7-8
 49 //  October,  2010, to be published by World S     49 //  October,  2010, to be published by World Scientific (Singapore).
 50 //                                                 50 //
 51 //      Available for downloading at:              51 //      Available for downloading at:
 52 //      http://arxiv.org/abs/1011.4822             52 //      http://arxiv.org/abs/1011.4822
 53 //                                                 53 //
 54 // -------------------------------------------     54 // -------------------------------------------------------------------
 55 //                                                 55 //
 56                                                    56 
 57 #ifndef G4IonCoulombScatteringModel_h              57 #ifndef G4IonCoulombScatteringModel_h
 58 #define G4IonCoulombScatteringModel_h 1            58 #define G4IonCoulombScatteringModel_h 1
 59                                                    59 
 60 #include "G4VEmModel.hh"                           60 #include "G4VEmModel.hh"
 61 #include "globals.hh"                              61 #include "globals.hh"
 62 #include "G4NistManager.hh"                        62 #include "G4NistManager.hh"
 63 #include "G4IonCoulombCrossSection.hh"             63 #include "G4IonCoulombCrossSection.hh"
 64                                                    64 
 65 #include <vector>                                  65 #include <vector>
 66                                                    66 
 67 class G4ParticleChangeForGamma;                    67 class G4ParticleChangeForGamma;
 68 class G4ParticleDefinition;                        68 class G4ParticleDefinition;
 69 class G4IonTable;                                  69 class G4IonTable;
 70                                                    70 
 71 class G4IonCoulombScatteringModel : public G4V     71 class G4IonCoulombScatteringModel : public G4VEmModel
 72 {                                                  72 {
 73 public:                                            73 public:
 74                                                    74 
 75   explicit G4IonCoulombScatteringModel(const G <<  75   G4IonCoulombScatteringModel(const G4String& nam = "IonCoulombScattering");
 76                "IonCoulombScattering");        << 
 77                                                    76  
 78   ~G4IonCoulombScatteringModel() override;     <<  77   virtual ~G4IonCoulombScatteringModel();
 79                                                    78 
 80   void Initialise(const G4ParticleDefinition*, <<  79   virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
 81                                                    80  
 82   G4double ComputeCrossSectionPerAtom(         <<  81   virtual G4double ComputeCrossSectionPerAtom(
 83                                 const G4Partic     82                                 const G4ParticleDefinition*,
 84         G4double kinEnergy,                        83         G4double kinEnergy, 
 85         G4double Z,                                84         G4double Z, 
 86         G4double A,                                85         G4double A, 
 87         G4double cut,                              86         G4double cut,
 88         G4double emax) final;                  <<  87         G4double emax);
 89                                                    88 
 90   void SampleSecondaries(std::vector<G4Dynamic <<  89   virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
 91        const G4MaterialCutsCouple*,            <<  90          const G4MaterialCutsCouple*,
 92        const G4DynamicParticle*,               <<  91          const G4DynamicParticle*,
 93        G4double tmin,                          <<  92          G4double tmin,
 94        G4double maxEnergy) final;              <<  93          G4double maxEnergy);
 95                                                    94 
 96                                                    95 
 97                                                    96     
 98   inline void  SetRecoilThreshold(G4double eth     97   inline void  SetRecoilThreshold(G4double eth);
 99                                                    98 
100   inline void  SetHeavyIonCorr(G4int b);           99   inline void  SetHeavyIonCorr(G4int b);
101                                                   100 
102   inline G4int GetHeavyIonCorr();                 101   inline G4int GetHeavyIonCorr();
103                                                   102 
104   // hide assignment operator                  << 103   //protected: 
105   G4IonCoulombScatteringModel & operator=      << 
106   (const G4IonCoulombScatteringModel &right) = << 
107   G4IonCoulombScatteringModel(const  G4IonCoul << 
108                                                   104    
109 private:                                          105 private:
110                                                   106 
111   inline void DefineMaterial(const G4MaterialC    107   inline void DefineMaterial(const G4MaterialCutsCouple*);
112                                                   108   
113   inline void SetupParticle(const G4ParticleDe    109   inline void SetupParticle(const G4ParticleDefinition*);
                                                   >> 110 
                                                   >> 111   // hide assignment operator
                                                   >> 112   G4IonCoulombScatteringModel & operator=(const G4IonCoulombScatteringModel &right);
                                                   >> 113   G4IonCoulombScatteringModel(const  G4IonCoulombScatteringModel&);
                                                   >> 114 
                                                   >> 115   //protected:
114                                                   116 
115   G4IonTable*               theIonTable;          117   G4IonTable*               theIonTable;
116   G4ParticleChangeForGamma* fParticleChange;      118   G4ParticleChangeForGamma* fParticleChange; 
117   G4NistManager*            fNistManager;         119   G4NistManager*            fNistManager;
118   G4IonCoulombCrossSection* ioncross;             120   G4IonCoulombCrossSection* ioncross;   
119                                                   121 
120   const std::vector<G4double>* pCuts;             122   const std::vector<G4double>* pCuts;
121   const G4MaterialCutsCouple* currentCouple;      123   const G4MaterialCutsCouple* currentCouple; 
122   const G4Material*           currentMaterial;    124   const G4Material*           currentMaterial;
123   const G4Element*            currentElement;     125   const G4Element*            currentElement;
124   G4int                       currentMaterialI    126   G4int                       currentMaterialIndex;
125                                                   127   
126   G4int           heavycorr;                      128   G4int           heavycorr;
127                                                   129 
128   G4double                  cosThetaMin;          130   G4double                  cosThetaMin;
129   G4double                  recoilThreshold;      131   G4double                  recoilThreshold;
130                                                   132         
131   // projectile                                   133   // projectile
132   const G4ParticleDefinition* particle;           134   const G4ParticleDefinition* particle;   
133   const G4ParticleDefinition* theProton;          135   const G4ParticleDefinition* theProton;  
134   G4double                  mass;                 136   G4double                  mass;   
135                                                   137 
136 };                                                138 };
137                                                   139 
138 //....oooOO0OOooo........oooOO0OOooo........oo    140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
139                                                   141 
140 inline void                                       142 inline void 
141 G4IonCoulombScatteringModel::DefineMaterial(co    143 G4IonCoulombScatteringModel::DefineMaterial(const G4MaterialCutsCouple* cup)
142 {                                                 144 { 
143   if(cup != currentCouple) {                      145   if(cup != currentCouple) {
144     currentCouple = cup;                          146     currentCouple = cup;
145     currentMaterial = cup->GetMaterial();         147     currentMaterial = cup->GetMaterial();
146     currentMaterialIndex = currentCouple->GetI    148     currentMaterialIndex = currentCouple->GetIndex();
147   }                                               149   }
148 }                                                 150 }
149                                                   151 
150 //....oooOO0OOooo........oooOO0OOooo........oo    152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
151                                                   153 
152 inline                                            154 inline
153 void G4IonCoulombScatteringModel::SetupParticl    155 void G4IonCoulombScatteringModel::SetupParticle(const G4ParticleDefinition* p)
154 {                                                 156 {
155   if(p != particle) {                             157   if(p != particle) {
156     particle = p;                                 158     particle = p;
157     mass = particle->GetPDGMass();                159     mass = particle->GetPDGMass();
158     ioncross->SetupParticle(p);                   160     ioncross->SetupParticle(p);
159   }                                               161   }
160 }                                                 162 }
161                                                   163 
162 //....oooOO0OOooo........oooOO0OOooo........oo    164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
163                                                   165 
164 inline void G4IonCoulombScatteringModel::SetRe    166 inline void G4IonCoulombScatteringModel::SetRecoilThreshold(G4double eth)
165 {                                                 167 {
166   recoilThreshold = eth;                          168   recoilThreshold = eth;
167 }                                                 169 }
168                                                   170 
169 //....oooOO0OOooo........oooOO0OOooo........oo    171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
170                                                   172 
171 inline void G4IonCoulombScatteringModel::SetHe    173 inline void G4IonCoulombScatteringModel::SetHeavyIonCorr(G4int b) 
172 {                                                 174 {
173   heavycorr = b;                                  175   heavycorr = b; 
174 }                                                 176 }
175                                                   177 
176 //....oooOO0OOooo........oooOO0OOooo........oo    178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
177                                                   179 
178 inline G4int G4IonCoulombScatteringModel::GetH    180 inline G4int G4IonCoulombScatteringModel::GetHeavyIonCorr() 
179 {                                                 181 {
180   return heavycorr;                               182   return heavycorr; 
181 }                                                 183 }
182                                                   184 
183 //....oooOO0OOooo........oooOO0OOooo........oo    185 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
184                                                   186 
185 #endif                                            187 #endif
186                                                   188