include/Geant4/G4InitXscPAI.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // G4InitXscPAI.hh -- header file
0027 //
0028 // History:
0029 //
0030 // 02.04.04, V. Grichine: 1st version based on G4PAIxSection class
0031
0032 #ifndef G4INITXSCPAI_HH
0033 #define G4INITXSCPAI_HH
0034
0035 #include "G4ios.hh"
0036 #include "globals.hh"
0037 #include "Randomize.hh"
0038 #include "G4OrderedTable.hh"
0039 #include "G4PhysicsLogVector.hh"
0040
0041 class G4MaterialCutsCouple;
0042 class G4SandiaTable;
0043
0044 class G4InitXscPAI
0045 {
0046 public:
0047       // Constructors
0048   explicit G4InitXscPAI( const G4MaterialCutsCouple* matCC);
0049
0050   ~G4InitXscPAI() ;
0051
0052   // Methods
0053   // General control functions
0054
0055   void KillCloseIntervals();
0056
0057   void Normalisation();
0058
0059   // Physical methods
0060   G4double RutherfordIntegral( G4int intervalNumber,
0061                                    G4double limitLow,
0062                        G4double limitHigh     ) ;
0063
0064   G4double IntegralTerm(G4double omega);
0065
0066   G4double ImPartDielectricConst( G4int intervalNumber,
0067                                       G4double energy        ) ;
0068
0069   G4double RePartDielectricConst(G4double energy) ;
0070
0071   G4double ModuleSqDielectricConst( G4int intervalNumber,
0072                                       G4double energy        ) ;
0073
0074   G4double DifPAIxSection( G4double omega ) ;
0075   G4double DifPAIdEdx( G4double omega ) ;
0076
0077   G4double PAIdNdxCherenkov( G4double omega ) ;
0078
0079   G4double PAIdNdxPlasmon( G4double omega ) ;
0080
0081   void     IntegralPAIxSection(G4double bg2, G4double Tmax) ;
0082   void     IntegralCherenkov(G4double bg2, G4double Tmax) ;
0083   void     IntegralPlasmon(G4double bg2, G4double Tmax) ;
0084
0085   void      IntegralPAIdEdx(G4double bg2, G4double Tmax) ;
0086
0087   G4double GetPhotonLambda( G4double omega ) ;
0088
0089   G4double GetStepEnergyLoss( G4double step ) ;
0090   G4double GetStepCerenkovLoss( G4double step ) ;
0091   G4double GetStepPlasmonLoss( G4double step ) ;
0092
0093   // Inline access functions
0094
0095
0096   G4int GetIntervalNumber() const { return fIntervalNumber ; }
0097   G4int GetBinPAI() const { return fPAIbin ; }
0098
0099   G4double GetNormalizationCof() const { return fNormalizationCof ; }
0100
0101   G4double GetMatSandiaMatrix(G4int i, G4int j) const
0102           { return (*(*fMatSandiaMatrix)[i])[j]; }
0103
0104   G4PhysicsLogVector* GetPAIxscVector() const { return fPAIxscVector;}
0105   G4PhysicsLogVector* GetPAIdEdxVector() const { return fPAIdEdxVector;}
0106   G4PhysicsLogVector* GetPAIphotonVector() const { return fPAIphotonVector;}
0107   G4PhysicsLogVector* GetPAIelectronVector() const { return fPAIelectronVector;}
0108   G4PhysicsLogVector* GetChCosSqVector() const { return fChCosSqVector;}
0109   G4PhysicsLogVector* GetChWidthVector() const { return fChWidthVector;}
0110
0111   G4InitXscPAI & operator=(const G4InitXscPAI &right) = delete;
0112   G4InitXscPAI(const G4InitXscPAI&) = delete;
0113
0114 private :
0115
0116   // Local class constants
0117
0118   static const G4double fDelta ; // energy shift from interval border = 0.001
0119   static const G4int fPAIbin;
0120   static const G4double fSolidDensity; // ~the border between gases and solids
0121
0122   G4int    fIntervalNumber;    //  The number of energy intervals
0123   G4double fNormalizationCof;   // Normalization cof for PhotoAbsorptionXsection
0124   G4int    fCurrentInterval;
0125   G4int    fIntervalTmax;
0126   G4double fBetaGammaSq ;        // (beta*gamma)^2
0127   G4double fTmax;
0128   G4double fDensity ;            // Current density
0129   G4double fElectronDensity ;    // Current electron (number) density
0130
0131   // Arrays of Sandia coefficients
0132   G4OrderedTable* fMatSandiaMatrix;
0133   G4SandiaTable*  fSandia;
0134
0135   // vectors of integral cross-sections
0136   G4PhysicsLogVector* fPAIxscVector;
0137   G4PhysicsLogVector* fPAIdEdxVector;
0138   G4PhysicsLogVector* fPAIphotonVector;
0139   G4PhysicsLogVector* fPAIelectronVector;
0140   G4PhysicsLogVector* fChCosSqVector;
0141   G4PhysicsLogVector* fChWidthVector;
0142 };
0143
0144 #endif
0145
0146 //
0147 //
0148 /////////////////   end of G4InitXscPAI header file    //////////////////////