Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4InitXscPAI.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 ]

  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 // G4InitXscPAI.hh -- header file
 27 //
 28 // History:
 29 //
 30 // 02.04.04, V. Grichine: 1st version based on G4PAIxSection class
 31 
 32 #ifndef G4INITXSCPAI_HH
 33 #define G4INITXSCPAI_HH
 34 
 35 #include "G4ios.hh"
 36 #include "globals.hh"
 37 #include "Randomize.hh"
 38 #include "G4OrderedTable.hh"
 39 #include "G4PhysicsLogVector.hh"
 40 
 41 class G4MaterialCutsCouple;
 42 class G4SandiaTable;
 43 
 44 class G4InitXscPAI
 45 {
 46 public:
 47     // Constructors
 48   explicit G4InitXscPAI( const G4MaterialCutsCouple* matCC);
 49 
 50   ~G4InitXscPAI() ;
 51 
 52   // Methods
 53   // General control functions
 54 
 55   void KillCloseIntervals();
 56 
 57   void Normalisation();
 58 
 59   // Physical methods
 60   G4double RutherfordIntegral( G4int intervalNumber,
 61                                  G4double limitLow,
 62                G4double limitHigh     ) ;
 63 
 64   G4double IntegralTerm(G4double omega);
 65 
 66   G4double ImPartDielectricConst( G4int intervalNumber,
 67                                     G4double energy        ) ;
 68 
 69   G4double RePartDielectricConst(G4double energy) ;
 70 
 71   G4double ModuleSqDielectricConst( G4int intervalNumber,
 72                                     G4double energy        ) ;
 73 
 74   G4double DifPAIxSection( G4double omega ) ;
 75   G4double DifPAIdEdx( G4double omega ) ;
 76 
 77   G4double PAIdNdxCherenkov( G4double omega ) ;
 78 
 79   G4double PAIdNdxPlasmon( G4double omega ) ;
 80 
 81   void     IntegralPAIxSection(G4double bg2, G4double Tmax) ;
 82   void     IntegralCherenkov(G4double bg2, G4double Tmax) ;
 83   void     IntegralPlasmon(G4double bg2, G4double Tmax) ;
 84 
 85   void      IntegralPAIdEdx(G4double bg2, G4double Tmax) ;
 86 
 87   G4double GetPhotonLambda( G4double omega ) ;
 88 
 89   G4double GetStepEnergyLoss( G4double step ) ;
 90   G4double GetStepCerenkovLoss( G4double step ) ;
 91   G4double GetStepPlasmonLoss( G4double step ) ;
 92 
 93   // Inline access functions
 94 
 95 
 96   G4int GetIntervalNumber() const { return fIntervalNumber ; }
 97   G4int GetBinPAI() const { return fPAIbin ; }
 98 
 99   G4double GetNormalizationCof() const { return fNormalizationCof ; }
100 
101   G4double GetMatSandiaMatrix(G4int i, G4int j) const
102           { return (*(*fMatSandiaMatrix)[i])[j]; }
103 
104   G4PhysicsLogVector* GetPAIxscVector() const { return fPAIxscVector;}
105   G4PhysicsLogVector* GetPAIdEdxVector() const { return fPAIdEdxVector;}
106   G4PhysicsLogVector* GetPAIphotonVector() const { return fPAIphotonVector;}
107   G4PhysicsLogVector* GetPAIelectronVector() const { return fPAIelectronVector;}
108   G4PhysicsLogVector* GetChCosSqVector() const { return fChCosSqVector;}
109   G4PhysicsLogVector* GetChWidthVector() const { return fChWidthVector;}
110 
111   G4InitXscPAI & operator=(const G4InitXscPAI &right) = delete;
112   G4InitXscPAI(const G4InitXscPAI&) = delete;
113 
114 private :
115 
116   // Local class constants
117 
118   static const G4double fDelta ; // energy shift from interval border = 0.001
119   static const G4int fPAIbin;
120   static const G4double fSolidDensity; // ~the border between gases and solids
121 
122   G4int    fIntervalNumber;    //  The number of energy intervals
123   G4double fNormalizationCof;   // Normalization cof for PhotoAbsorptionXsection
124   G4int    fCurrentInterval;
125   G4int    fIntervalTmax;
126   G4double fBetaGammaSq ;        // (beta*gamma)^2
127   G4double fTmax;
128   G4double fDensity ;            // Current density
129   G4double fElectronDensity ;    // Current electron (number) density
130 
131   // Arrays of Sandia coefficients
132   G4OrderedTable* fMatSandiaMatrix;
133   G4SandiaTable*  fSandia;
134 
135   // vectors of integral cross-sections
136   G4PhysicsLogVector* fPAIxscVector;
137   G4PhysicsLogVector* fPAIdEdxVector;
138   G4PhysicsLogVector* fPAIphotonVector;
139   G4PhysicsLogVector* fPAIelectronVector;
140   G4PhysicsLogVector* fChCosSqVector;
141   G4PhysicsLogVector* fChWidthVector;
142 };    
143 
144 #endif   
145 
146 //
147 //
148 /////////////////   end of G4InitXscPAI header file    //////////////////////
149