Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4ionIonisation.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/G4ionIonisation.hh (Version 11.3.0) and /processes/electromagnetic/standard/include/G4ionIonisation.hh (Version 9.1.p1)


  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 // $Id: G4ionIonisation.hh,v 1.50 2007/11/09 11:45:45 vnivanch Exp $
                                                   >>  27 // GEANT4 tag $Name: geant4-09-01-patch-01 $
 26 //                                                 28 //
 27 // -------------------------------------------     29 // -------------------------------------------------------------------
 28 //                                                 30 //
 29 // GEANT4 Class header file                        31 // GEANT4 Class header file
 30 //                                                 32 //
 31 //                                                 33 //
 32 // File name:     G4ionIonisation                  34 // File name:     G4ionIonisation
 33 //                                                 35 //
 34 // Author:        Vladimir Ivanchenko              36 // Author:        Vladimir Ivanchenko
 35 //                                                 37 //
 36 // Creation date: 07.05.2002                       38 // Creation date: 07.05.2002
 37 //                                                 39 //
 38 // Modifications:                                  40 // Modifications:
 39 //                                                 41 //
 40 // 26-12-02 Secondary production moved to deri     42 // 26-12-02 Secondary production moved to derived classes (VI)
 41 // 24-01-03 Make models region aware (V.Ivanch     43 // 24-01-03 Make models region aware (V.Ivanchenko)
 42 // 05-02-03 Fix compilation warnings (V.Ivanch     44 // 05-02-03 Fix compilation warnings (V.Ivanchenko)
 43 // 13-02-03 SubCutoff regime is assigned to a      45 // 13-02-03 SubCutoff regime is assigned to a region (V.Ivanchenko)
 44 // 15-02-03 Add control on delta pointer (V.Iv     46 // 15-02-03 Add control on delta pointer (V.Ivanchenko)
 45 // 23-05-03 Add fluctuation model as a member      47 // 23-05-03 Add fluctuation model as a member function (V.Ivanchenko)
 46 // 03-08-03 Add effective charge and saturatio     48 // 03-08-03 Add effective charge and saturation of tmax (V.Ivanchenko)
 47 // 12-11-03 Fix problem of negative effective      49 // 12-11-03 Fix problem of negative effective charge (V.Ivanchenko)
 48 // 12-11-03 G4EnergyLossSTD -> G4EnergyLossPro     50 // 12-11-03 G4EnergyLossSTD -> G4EnergyLossProcess (V.Ivanchenko)
 49 // 21-01-04 Migrade to G4ParticleChangeForLoss     51 // 21-01-04 Migrade to G4ParticleChangeForLoss (V.Ivanchenko)
 50 // 08-11-04 Migration to new interface of Stor     52 // 08-11-04 Migration to new interface of Store/Retrieve tables (V.Ivantchenko)
 51 // 08-04-05 Major optimisation of internal int     53 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
 52 // 11-04-05 Move MaxSecondary energy to model      54 // 11-04-05 Move MaxSecondary energy to model (V.Ivanchneko)
 53 // 11-04-04 Move MaxSecondaryEnergy to models      55 // 11-04-04 Move MaxSecondaryEnergy to models (V.Ivanchenko)
 54 // 10-05-06 Add a possibility to download user     56 // 10-05-06 Add a possibility to download user data (V.Ivantchenko)
 55 // 22-07-06 Remove obsolete method (V.Ivantche     57 // 22-07-06 Remove obsolete method (V.Ivantchenko)
 56 // 07-11-07 Moved CorrectionsAlongStep to cc (     58 // 07-11-07 Moved CorrectionsAlongStep to cc (V.Ivantchenko)
 57 // 12-09-08 Removed InitialiseMassCharge and C << 
 58 //                                                 59 //
 59 // Class Description:                              60 // Class Description:
 60 //                                                 61 //
 61 // This class manages the ionisation process f     62 // This class manages the ionisation process for ions. Effective charge,
 62 // nuclear stopping power, energy loss correct     63 // nuclear stopping power, energy loss corrections are taken into account.
 63 // It inherites from G4VEnergyLossLoss.            64 // It inherites from G4VEnergyLossLoss.
 64 //                                                 65 //
 65                                                    66 
 66 // -------------------------------------------     67 // -------------------------------------------------------------------
 67 //                                                 68 //
 68                                                    69 
 69 #ifndef G4ionIonisation_h                          70 #ifndef G4ionIonisation_h
 70 #define G4ionIonisation_h 1                        71 #define G4ionIonisation_h 1
 71                                                    72 
 72 #include "G4VEnergyLossProcess.hh"                 73 #include "G4VEnergyLossProcess.hh"
                                                   >>  74 #include "G4ionEffectiveCharge.hh"
                                                   >>  75 #include "G4VEmModel.hh"
                                                   >>  76 #include "G4EmCorrections.hh"
 73                                                    77 
 74 class G4Material;                                  78 class G4Material;
                                                   >>  79 class G4PhysicsVector;
                                                   >>  80 class G4BraggIonModel;
 75                                                    81 
 76 class G4ionIonisation : public G4VEnergyLossPr     82 class G4ionIonisation : public G4VEnergyLossProcess
 77 {                                                  83 {
 78 public:                                            84 public:
 79                                                    85 
 80   explicit G4ionIonisation(const G4String& nam <<  86   G4ionIonisation(const G4String& name = "ionIoni");
 81                                                    87 
 82   ~G4ionIonisation() override = default;       <<  88   virtual ~G4ionIonisation();
 83                                                    89 
 84   G4bool IsApplicable(const G4ParticleDefiniti <<  90   inline G4bool IsApplicable(const G4ParticleDefinition& p);
                                                   >>  91 
                                                   >>  92   // Print out of the class parameters
                                                   >>  93   virtual void PrintInfo();
 85                                                    94 
 86   void AddStoppingData(G4int Z, G4int A, const     95   void AddStoppingData(G4int Z, G4int A, const G4String& materialName,
 87            G4PhysicsVector* dVector);          <<  96            G4PhysicsVector& dVector);
 88                                                    97 
 89   void ActivateStoppingData(G4bool);               98   void ActivateStoppingData(G4bool);
 90                                                    99 
 91   // print documentation in html format        << 100   void ActivateNuclearStopping(G4bool);
 92   void ProcessDescription(std::ostream&) const << 
 93                                                << 
 94   // hide assignment operator                  << 
 95   G4ionIonisation & operator=(const G4ionIonis << 
 96   G4ionIonisation(const G4ionIonisation&) = de << 
 97                                                   101 
 98 protected:                                        102 protected:
 99                                                   103 
100   void InitialiseEnergyLossProcess(const G4Par << 104   virtual void InitialiseEnergyLossProcess(const G4ParticleDefinition*,
101            const G4ParticleDefinition*) overri << 105              const G4ParticleDefinition*);
                                                   >> 106 
                                                   >> 107   virtual void CorrectionsAlongStep(const G4MaterialCutsCouple*,
                                                   >> 108             const G4DynamicParticle*,
                                                   >> 109             G4double& eloss,
                                                   >> 110             G4double& length);
                                                   >> 111 
                                                   >> 112   inline void InitialiseMassCharge(const G4Track&);
102                                                   113 
103   G4double MinPrimaryEnergy(const G4ParticleDe << 114   inline G4double MinPrimaryEnergy(const G4ParticleDefinition* p,
104           const G4Material*, G4double cut) fin << 115            const G4Material*, G4double cut);
105                                                   116 
106   inline G4double BetheBlochEnergyThreshold();    117   inline G4double BetheBlochEnergyThreshold();
107                                                   118 
108 private:                                       << 119   inline G4bool NuclearStoppingFlag();
                                                   >> 120 
                                                   >> 121   // protected pointers 
                                                   >> 122   G4ionEffectiveCharge*       effCharge;
                                                   >> 123   G4EmCorrections*            corr;
109                                                   124 
                                                   >> 125 private:
110                                                   126 
111   const G4ParticleDefinition* theParticle = nu << 127   // hide assignment operator
112   G4double eth;                                << 128   G4ionIonisation & operator=(const G4ionIonisation &right);
                                                   >> 129   G4ionIonisation(const G4ionIonisation&);
113                                                   130 
114   G4bool isInitialised = false;                << 131   // cash
                                                   >> 132   const G4Material*           curMaterial;
                                                   >> 133   const G4ParticleDefinition* curParticle;
                                                   >> 134   const G4ParticleDefinition* theParticle;
                                                   >> 135   const G4ParticleDefinition* theBaseParticle;
                                                   >> 136 
                                                   >> 137   G4double                    preKinEnergy;
                                                   >> 138 
                                                   >> 139   G4double                    eth;
                                                   >> 140   G4double                    baseMass;
                                                   >> 141   G4double                    massRatio;
                                                   >> 142   G4double                    massFactor;
                                                   >> 143   G4double                    charge2;
                                                   >> 144 
                                                   >> 145   G4bool                      isInitialised;
                                                   >> 146   G4bool                      stopDataActive;
                                                   >> 147   G4bool                      nuclearStopping;
115 };                                                148 };
116                                                   149 
117 //....oooOO0OOooo........oooOO0OOooo........oo    150 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 152 
                                                   >> 153 inline G4bool G4ionIonisation::IsApplicable(const G4ParticleDefinition& p)
                                                   >> 154 {
                                                   >> 155   return (p.GetPDGCharge() != 0.0 && !p.IsShortLived() &&
                                                   >> 156           p.GetParticleType() == "nucleus");
                                                   >> 157 }
                                                   >> 158 
                                                   >> 159 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 160 
                                                   >> 161 inline G4double G4ionIonisation::MinPrimaryEnergy(
                                                   >> 162           const G4ParticleDefinition*, const G4Material*, G4double cut)
                                                   >> 163 {
                                                   >> 164   G4double x = 0.5*cut/electron_mass_c2;
                                                   >> 165   G4double g = std::sqrt(1. + x);
                                                   >> 166   return proton_mass_c2*(g - 1.0);
                                                   >> 167 }
                                                   >> 168 
                                                   >> 169 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 170 
                                                   >> 171 inline void G4ionIonisation::InitialiseMassCharge(const G4Track& track)
                                                   >> 172 {
                                                   >> 173   preKinEnergy = track.GetKineticEnergy();
                                                   >> 174   massRatio    = baseMass/track.GetDynamicParticle()->GetMass();
                                                   >> 175   charge2      = effCharge->EffectiveChargeSquareRatio(track.GetDefinition(),
                                                   >> 176                    track.GetMaterial(),
                                                   >> 177                    preKinEnergy);
                                                   >> 178   SetDynamicMassCharge(massRatio, charge2);
                                                   >> 179 }
                                                   >> 180 
                                                   >> 181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 182 
                                                   >> 183 inline void G4ionIonisation::ActivateStoppingData(G4bool val)
                                                   >> 184 {
                                                   >> 185   stopDataActive = val;
                                                   >> 186 }
                                                   >> 187 
                                                   >> 188 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 189 
                                                   >> 190 inline void G4ionIonisation::ActivateNuclearStopping(G4bool val)
                                                   >> 191 {
                                                   >> 192   nuclearStopping = val;
                                                   >> 193 }
                                                   >> 194 
                                                   >> 195 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
118                                                   196 
119 inline G4double G4ionIonisation::BetheBlochEne    197 inline G4double G4ionIonisation::BetheBlochEnergyThreshold()
120 {                                                 198 {
121   return eth;                                     199   return eth;
                                                   >> 200 }
                                                   >> 201 
                                                   >> 202 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 203 
                                                   >> 204 inline G4bool G4ionIonisation::NuclearStoppingFlag()
                                                   >> 205 {
                                                   >> 206   return nuclearStopping;
122 }                                                 207 }
123                                                   208 
124 //....oooOO0OOooo........oooOO0OOooo........oo    209 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
125                                                   210 
126 #endif                                            211 #endif
127                                                   212