Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/utils/include/G4EmSaturation.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/G4EmSaturation.hh (Version 11.3.0) and /processes/electromagnetic/utils/include/G4EmSaturation.hh (Version 9.4)


  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: G4EmSaturation.hh,v 1.8 2010/08/17 17:36:58 vnivanch Exp $
                                                   >>  27 // GEANT4 tag $Name: geant4-09-04 $
                                                   >>  28 //
 26 //                                                 29 //
 27 #ifndef G4EmSaturation_h                           30 #ifndef G4EmSaturation_h
 28 #define G4EmSaturation_h 1                         31 #define G4EmSaturation_h 1
 29                                                    32 
 30 // -------------------------------------------     33 // -------------------------------------------------------------
 31 //                                                 34 //
 32 // GEANT4 Class header file                        35 // GEANT4 Class header file
 33 //                                                 36 //
 34 //                                                 37 //
 35 // File name:     G4EmSaturation                   38 // File name:     G4EmSaturation
 36 //                                                 39 //
 37 // Author:        Vladimir Ivanchenko              40 // Author:        Vladimir Ivanchenko
 38 //                                                 41 //
 39 // Creation date: 18.02.2008                       42 // Creation date: 18.02.2008
 40 //                                                 43 //
 41 // Modifications:                                  44 // Modifications:
 42 //                                                 45 //
 43 //                                                 46 //
 44 // Class Description:                              47 // Class Description:
 45 //   Compution on saturation effect, which red     48 //   Compution on saturation effect, which reduce visible energy 
 46 //   deposition at the step. Default implement     49 //   deposition at the step. Default implementation takes into 
 47 //   account Birks effect. Birks coefficients      50 //   account Birks effect. Birks coefficients for some materials
 48 //   from G4 database on materials are provide     51 //   from G4 database on materials are provided
 49 //                                             << 
 50 //   This class assumed to be G4ThreadLocal, b << 
 51 //   cache value for material                  << 
 52 //                                                 52 // 
 53 // -------------------------------------------     53 // -------------------------------------------------------------
 54                                                    54 
 55 //....oooOO0OOooo........oooOO0OOooo........oo     55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 56 //....oooOO0OOooo........oooOO0OOooo........oo     56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 57                                                    57 
 58 #include "globals.hh"                              58 #include "globals.hh"
 59 #include "G4Step.hh"                               59 #include "G4Step.hh"
 60 #include "G4ParticleDefinition.hh"                 60 #include "G4ParticleDefinition.hh"
 61 #include <vector>                                  61 #include <vector>
 62                                                    62 
 63 //....oooOO0OOooo........oooOO0OOooo........oo     63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 64                                                    64 
                                                   >>  65 class G4LossTableManager;
 65 class G4NistManager;                               66 class G4NistManager;
 66 class G4MaterialCutsCouple;                        67 class G4MaterialCutsCouple;
 67 class G4Material;                                  68 class G4Material;
 68                                                    69 
 69 class G4EmSaturation                               70 class G4EmSaturation
 70 {                                                  71 {
 71 public:                                            72 public: 
 72                                                    73 
 73   explicit G4EmSaturation(G4int verb);         <<  74   G4EmSaturation();
 74   virtual ~G4EmSaturation();                       75   virtual ~G4EmSaturation();
 75                                                    76 
 76   // this method may be overwritten in the der <<  77   G4double VisibleEnergyDeposition(const G4ParticleDefinition*, 
 77   // which implements alternative algorithm of <<  78            const G4MaterialCutsCouple*,
 78   virtual G4double VisibleEnergyDeposition(con <<  79            G4double length, 
 79              const G4MaterialCutsCouple*,      <<  80            G4double edepTotal,
 80              G4double length,                  <<  81            G4double edepNIEL = 0.0);
 81              G4double edepTotal,               << 
 82              G4double edepNIEL = 0.0) const;   << 
 83                                                    82 
 84   // activate default model                    <<  83   inline G4double VisibleEnergyDeposition(const G4Step*); 
 85   void InitialiseG4Saturation();               << 
 86                                                    84 
 87   // find and Birks coefficient                    85   // find and Birks coefficient 
 88   G4double FindG4BirksCoefficient(const G4Mate     86   G4double FindG4BirksCoefficient(const G4Material*);
 89                                                    87 
                                                   >>  88   inline void SetVerbose(G4int);
                                                   >>  89 
 90   // dump coeffitients used in run time            90   // dump coeffitients used in run time
 91   void DumpBirksCoefficients();                    91   void DumpBirksCoefficients();
 92                                                    92 
 93   // dump G4 list                                  93   // dump G4 list
 94   void DumpG4BirksCoefficients();                  94   void DumpG4BirksCoefficients();
 95                                                    95 
 96   // this method should not be overwritten     <<  96 private:
 97   inline G4double VisibleEnergyDepositionAtASt << 
 98                                                << 
 99   inline void SetVerbose(G4int);               << 
100                                                    97 
101   // hide assignment operator                      98   // hide assignment operator
102   G4EmSaturation & operator=(const G4EmSaturat <<  99   G4EmSaturation & operator=(const G4EmSaturation &right);
103   G4EmSaturation(const G4EmSaturation&) = dele << 100   G4EmSaturation(const G4EmSaturation&);
104                                                << 
105 private:                                       << 
106                                                   101 
107   void InitialiseBirksCoefficient(const G4Mate << 102   G4double FindBirksCoefficient(const G4Material*);
108                                                   103 
109   void InitialiseG4materials();                << 104   void Initialise();
110                                                   105 
111   const G4ParticleDefinition* electron;           106   const G4ParticleDefinition* electron;
112   const G4ParticleDefinition* proton;             107   const G4ParticleDefinition* proton;
                                                   >> 108   G4LossTableManager*         manager;
113   G4NistManager*              nist;               109   G4NistManager*              nist;
114                                                   110 
                                                   >> 111   // cash
                                                   >> 112   const G4Material*           curMaterial;
                                                   >> 113   G4double                    curBirks;
                                                   >> 114   G4double                    curRatio;
                                                   >> 115   G4double                    curChargeSq;
                                                   >> 116 
115   G4int    verbose;                               117   G4int    verbose;             
                                                   >> 118   G4int    nMaterials;
116   G4int    nG4Birks;                              119   G4int    nG4Birks;
117   G4int    nWarnings;                          << 
118                                                << 
119   static size_t nMaterials;                    << 
120                                                   120 
121   // list of materials used in run time           121   // list of materials used in run time
122   static std::vector<G4double> massFactors;    << 122   std::vector<const G4Material*>    matPointers;
123   static std::vector<G4double> effCharges;     << 123   std::vector<G4String>             matNames;
                                                   >> 124   std::vector<G4double>             massFactors;
                                                   >> 125   std::vector<G4double>             effCharges;
124                                                   126 
125   // list of G4 materials                         127   // list of G4 materials 
126   static std::vector<G4double>  g4MatData;     << 128   std::vector<G4double>             g4MatData;
127   static std::vector<G4String>  g4MatNames;    << 129   std::vector<G4String>             g4MatNames;
128 };                                             << 
129                                                   130 
130 //....oooOO0OOooo........oooOO0OOooo........oo << 131 };
131                                                   132 
132 inline void G4EmSaturation::SetVerbose(G4int v    133 inline void G4EmSaturation::SetVerbose(G4int val)
133 {                                                 134 {
134   verbose = val;                                  135   verbose = val;
135 }                                                 136 }
136                                                   137 
137 //....oooOO0OOooo........oooOO0OOooo........oo << 138 inline G4double G4EmSaturation::VisibleEnergyDeposition(
138                                                << 139                 const G4Step* step)
139 inline G4double G4EmSaturation::VisibleEnergyD << 
140                 const G4Step* step) const      << 
141 {                                                 140 {
142   return VisibleEnergyDeposition(step->GetTrac << 141   G4Track* track = step->GetTrack();
143                                  step->GetTrac << 142   return VisibleEnergyDeposition(track->GetParticleDefinition(),
                                                   >> 143                                  track->GetMaterialCutsCouple(),
144          step->GetStepLength(),                   144          step->GetStepLength(),
145                                  step->GetTota    145                                  step->GetTotalEnergyDeposit(),
146                                  step->GetNonI    146                                  step->GetNonIonizingEnergyDeposit());
147 }                                                 147 }
148                                                << 
149 //....oooOO0OOooo........oooOO0OOooo........oo << 
150                                                   148 
151 #endif                                            149 #endif
152                                                   150 
153                                                   151