include/Geant4/G4ecpssrBaseLixsModel.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 //
0027 // Author: Haifa Ben Abdelouahed
0028 //
0029 //
0030 // History:
0031 // -----------
0032 //  23 Apr 2008   H. Ben Abdelouahed   1st implementation
0033 //  28 Apr 2008   MGP        Major revision according to a design iteration
0034 //  29 Apr 2009   ALF Updated Desing for Integration
0035 //  01 Sep 2009   ALF Updated to G4AnalyticalEcpssrLiCrossSection
0036 //  28 Oct 2011   ALF Changed name G4AnalyticalEcpssrLiCrossSection to G4ecpssrBaseLixsModel
0037 //
0038 // -------------------------------------------------------------------
0039
0040 // Class description:
0041 // Low Energy Electromagnetic Physics, Cross section, p and alpha ionisation, L shell
0042 // -------------------------------------------------------------------
0043
0044 #ifndef G4ecpssrBaseLixsModel_hh
0045 #define G4ecpssrBaseLixsModel_hh 1
0046
0047 #include "G4VecpssrLiModel.hh"
0048 #include "globals.hh"
0049 #include <map>
0050 #include <vector>
0051
0052 class G4ecpssrBaseLixsModel : public G4VecpssrLiModel
0053 {
0054 public:
0055   explicit G4ecpssrBaseLixsModel();
0056
0057   ~G4ecpssrBaseLixsModel();
0058
0059   G4double CalculateL1CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident) override;//according to W.Brandt and G.Lapicki, Phys.Rev.A23(1981)
0060
0061   G4double CalculateL2CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident) override;//according to W.Brandt and G.Lapicki, Phys.Rev.A23(1981)
0062
0063   G4double CalculateL3CrossSection(G4int zTarget,G4double massIncident, G4double energyIncident) override;//according to W.Brandt and G.Lapicki, Phys.Rev.A23(1981)
0064
0065   G4double CalculateVelocity(G4int subShell, G4int zTarget,G4double massIncident, G4double energyIncident);
0066   G4double  ExpIntFunction(G4int n,G4double x);//Exponential Integral Function
0067   G4ecpssrBaseLixsModel(const G4ecpssrBaseLixsModel&) = delete;
0068   G4ecpssrBaseLixsModel & operator = (const G4ecpssrBaseLixsModel &right) = delete;
0069
0070 private:
0071   G4double FunctionFL1(G4double k, G4double theta);
0072   G4double FunctionFL2(G4double k, G4double theta);
0073   G4double LogLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
0074   G4double LinLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
0075   G4double LinLinInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
0076   G4double QuadInterpolator(G4double e11,
0077                     G4double e12,
0078                 G4double e21,
0079                 G4double e22,
0080                 G4double x11,
0081                 G4double x12,
0082                 G4double x21,
0083                 G4double x22,
0084                 G4double t1,
0085                 G4double t2,
0086                 G4double t,
0087                 G4double e);
0088
0089   typedef std::map<G4double, std::map<G4double, G4double> > TriDimensionMap;
0090   TriDimensionMap FL1Data;
0091   TriDimensionMap FL2Data;
0092   std::vector<G4double> dummyVec1;
0093   std::vector<G4double> dummyVec2;
0094
0095   typedef std::map<G4double, std::vector<G4double> > VecMap;
0096   VecMap aVecMap1;
0097   VecMap aVecMap2;
0098
0099   G4int verboseLevel;
0100
0101 };
0102
0103 #endif