Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/utils/include/G4ElectronIonPair.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/utils/include/G4ElectronIonPair.hh (Version 11.3.0) and /processes/electromagnetic/utils/include/G4ElectronIonPair.hh (Version 9.4.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 // $Id: G4ElectronIonPair.hh,v 1.5 2010-10-25 17:23:01 vnivanch Exp $
                                                   >>  27 // GEANT4 tag $Name: geant4-09-04-patch-02 $
                                                   >>  28 //
                                                   >>  29 //
 26 #ifndef G4ElectronIonPair_h                        30 #ifndef G4ElectronIonPair_h
 27 #define G4ElectronIonPair_h 1                      31 #define G4ElectronIonPair_h 1
 28                                                    32 
 29 // -------------------------------------------     33 // -------------------------------------------------------------
 30 //                                                 34 //
 31 // GEANT4 Class header file                        35 // GEANT4 Class header file
 32 //                                                 36 //
 33 //                                                 37 //
 34 // File name:     G4ElectronIonPair                38 // File name:     G4ElectronIonPair
 35 //                                                 39 //
 36 // Author:        Vladimir Ivanchenko              40 // Author:        Vladimir Ivanchenko
 37 //                                                 41 //
 38 // Creation date: 08.07.2008                       42 // Creation date: 08.07.2008
 39 //                                                 43 //
 40 // Modifications:                                  44 // Modifications:
 41 //                                                 45 //
 42 //                                                 46 //
 43 // Class Description:                              47 // Class Description:
 44 //   Compution on number of electon-ion or ele     48 //   Compution on number of electon-ion or electorn-hole pairs
 45 //   at the step of a particle and sampling io     49 //   at the step of a particle and sampling ionisation points 
 46 //   in space                                      50 //   in space
 47 //                                                 51 //
 48 // Based on ICRU Report 31, 1979                   52 // Based on ICRU Report 31, 1979
 49 // "Average Energy Required to Produce an Ion  <<  53 // "Average Energy Required to Produce an Ion Pair"  
 50 //                                             << 
 51 // 06.04.2010 V. Grichine, substitute Gauss by << 
 52 //                         distribution at fix << 
 53 //                                                 54 // 
 54 // -------------------------------------------     55 // -------------------------------------------------------------
 55                                                    56 
 56 //....oooOO0OOooo........oooOO0OOooo........oo     57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 57 //....oooOO0OOooo........oooOO0OOooo........oo     58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 58                                                    59 
 59 #include "globals.hh"                              60 #include "globals.hh"
 60 #include "G4Step.hh"                               61 #include "G4Step.hh"
 61 #include "G4ParticleDefinition.hh"                 62 #include "G4ParticleDefinition.hh"
 62 #include "G4ThreeVector.hh"                        63 #include "G4ThreeVector.hh"
 63 #include "G4VProcess.hh"                           64 #include "G4VProcess.hh"
 64 #include "Randomize.hh"                        << 
 65 #include <vector>                                  65 #include <vector>
 66                                                    66 
 67 //....oooOO0OOooo........oooOO0OOooo........oo     67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 68                                                    68 
 69 class G4Material;                                  69 class G4Material;
 70                                                    70 
 71 class G4ElectronIonPair                            71 class G4ElectronIonPair
 72 {                                                  72 {
 73 public:                                            73 public: 
 74                                                    74 
 75   explicit G4ElectronIonPair(G4int verb);      <<  75   G4ElectronIonPair();
 76                                                    76 
 77   virtual ~G4ElectronIonPair();                    77   virtual ~G4ElectronIonPair();
 78                                                    78 
 79   // compute mean number of ionisation points      79   // compute mean number of ionisation points at a step
 80   G4double MeanNumberOfIonsAlongStep(const G4P     80   G4double MeanNumberOfIonsAlongStep(const G4ParticleDefinition*, 
 81              const G4Material*,                    81              const G4Material*,
 82              G4double edepTotal,                   82              G4double edepTotal,
 83              G4double edepNIEL = 0.0);             83              G4double edepNIEL = 0.0);
 84                                                    84 
 85   inline G4double MeanNumberOfIonsAlongStep(co     85   inline G4double MeanNumberOfIonsAlongStep(const G4Step*); 
 86                                                    86 
 87   inline G4int SampleNumberOfIonsAlongStep(con     87   inline G4int SampleNumberOfIonsAlongStep(const G4Step*); 
 88                                                    88 
 89   // returns pointer to the new vector of posi     89   // returns pointer to the new vector of positions of
 90   // ionisation points in the World coordinate     90   // ionisation points in the World coordinate system 
 91   std::vector<G4ThreeVector>* SampleIonsAlongS     91   std::vector<G4ThreeVector>* SampleIonsAlongStep(const G4Step*);
 92                                                    92 
 93   // compute number of holes in the atom after     93   // compute number of holes in the atom after PostStep interaction
 94   G4int ResidualeChargePostStep(const G4Partic     94   G4int ResidualeChargePostStep(const G4ParticleDefinition*,
 95         const G4TrackVector* secondary = nullp <<  95         const G4TrackVector* secondary = 0,
 96         G4int processSubType = -1) const;      <<  96         G4int processSubType = -1);
 97                                                    97 
 98   inline G4int ResidualeChargePostStep(const G <<  98   inline G4int ResidualeChargePostStep(const G4Step*);
 99                                                    99 
100   // find mean energies per ionisation            100   // find mean energies per ionisation 
101   G4double FindG4MeanEnergyPerIonPair(const G4 << 101   G4double FindG4MeanEnergyPerIonPair(const G4Material*);
102                                                   102 
103   // dump mean energies per ionisation used in    103   // dump mean energies per ionisation used in run time
104   void DumpMeanEnergyPerIonPair() const;       << 104   void DumpMeanEnergyPerIonPair();
105                                                   105 
106   // dump G4 list                                 106   // dump G4 list
107   void DumpG4MeanEnergyPerIonPair() const;     << 107   void DumpG4MeanEnergyPerIonPair();
108                                                   108 
109   inline void SetVerbose(G4int);                  109   inline void SetVerbose(G4int);
110                                                   110 
111   // hide assignment operator                  << 
112   G4ElectronIonPair & operator=(const G4Electr << 
113   G4ElectronIonPair(const G4ElectronIonPair&)  << 
114                                                << 
115 private:                                          111 private:
116                                                   112 
117   void Initialise();                              113   void Initialise();
118                                                   114 
119   G4double FindMeanEnergyPerIonPair(const G4Ma << 115   G4double FindMeanEnergyPerIonPair(const G4Material*);
                                                   >> 116 
                                                   >> 117   // hide assignment operator
                                                   >> 118   G4ElectronIonPair & operator=(const G4ElectronIonPair &right);
                                                   >> 119   G4ElectronIonPair(const G4ElectronIonPair&);
120                                                   120 
121   // cache                                     << 121   // cash
122   const G4Material*  curMaterial;                 122   const G4Material*  curMaterial;
123   G4double           curMeanEnergy;               123   G4double           curMeanEnergy;
124                                                   124 
125   G4double invFanoFactor;                      << 125   G4double FanoFactor;
126                                                   126   
127   G4int    verbose;                               127   G4int    verbose;             
128   G4int    nMaterials;                            128   G4int    nMaterials;
129                                                   129 
130   // list of G4 NIST materials with mean energ    130   // list of G4 NIST materials with mean energy per ion defined 
131   std::vector<G4double>  g4MatData;               131   std::vector<G4double>  g4MatData;
132   std::vector<G4String>  g4MatNames;              132   std::vector<G4String>  g4MatNames;
133 };                                                133 };
134                                                   134 
135 inline G4double                                   135 inline G4double 
136 G4ElectronIonPair::MeanNumberOfIonsAlongStep(c << 136 G4ElectronIonPair::MeanNumberOfIonsAlongStep(const G4Step* step)
137 {                                                 137 {
138   return MeanNumberOfIonsAlongStep(step->GetTr    138   return MeanNumberOfIonsAlongStep(step->GetTrack()->GetParticleDefinition(),
139            step->GetPreStepPoint()->GetMateria    139            step->GetPreStepPoint()->GetMaterial(),
140            step->GetTotalEnergyDeposit(),         140            step->GetTotalEnergyDeposit(),
141            step->GetNonIonizingEnergyDeposit()    141            step->GetNonIonizingEnergyDeposit());
142 }                                                 142 }
143                                                   143 
144 inline G4int                                   << 144 inline 
145 G4ElectronIonPair::SampleNumberOfIonsAlongStep << 145 G4int G4ElectronIonPair::SampleNumberOfIonsAlongStep(const G4Step* step)
146 {                                                 146 {
147   // use gamma distribution with mean value n= << 
148   // dispersion D=meanion/invFanoFactor        << 
149   G4double meanion = MeanNumberOfIonsAlongStep    147   G4double meanion = MeanNumberOfIonsAlongStep(step);
150   return G4lrint(G4RandGamma::shoot(meanion*in << 148   G4double sig = FanoFactor*std::sqrt(meanion);
                                                   >> 149   G4int nion = G4int(G4RandGauss::shoot(meanion,sig) + 0.5);
                                                   >> 150   return nion;
151 }                                                 151 } 
152                                                   152 
153 inline G4int                                   << 153 inline 
154 G4ElectronIonPair::ResidualeChargePostStep(con << 154 G4int G4ElectronIonPair::ResidualeChargePostStep(const G4Step* step)
155 {                                                 155 {
156   G4int subtype = -1;                             156   G4int subtype = -1;
157   const G4VProcess* proc = step->GetPostStepPo    157   const G4VProcess* proc = step->GetPostStepPoint()->GetProcessDefinedStep();
158   if(proc) { subtype = proc->GetProcessSubType    158   if(proc) { subtype = proc->GetProcessSubType(); }
159   return ResidualeChargePostStep(step->GetTrac    159   return ResidualeChargePostStep(step->GetTrack()->GetParticleDefinition(),
160          step->GetSecondary(),                    160          step->GetSecondary(),
161          subtype);                                161          subtype);
162 }                                                 162 }
163                                                   163 
164 inline void G4ElectronIonPair::SetVerbose(G4in    164 inline void G4ElectronIonPair::SetVerbose(G4int val)
165 {                                                 165 {
166   verbose = val;                                  166   verbose = val;
167 }                                                 167 }
168                                                   168 
169 #endif                                            169 #endif
170                                                   170 
171                                                   171