include/Geant4/G4ecpssrBaseKxsModel.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 #ifndef G4ecpssrBaseKxsModel_hh
0028 #define G4ecpssrBaseKxsModel_hh 1
0029
0030 #include "G4VecpssrKModel.hh"
0031 #include "globals.hh"
0032 #include <map>
0033 #include <vector>
0034
0035 #include "G4CrossSectionDataSet.hh"
0036
0037 class G4ecpssrBaseKxsModel : public G4VecpssrKModel
0038 {
0039 public:
0040   explicit G4ecpssrBaseKxsModel();
0041   ~G4ecpssrBaseKxsModel();
0042
0043   G4double CalculateCrossSection(G4int, G4double, G4double) override;//according to W.Brandt and G.Lapicki, Phys.Rev.A23(1981)
0044   G4double  ExpIntFunction(G4int n,G4double x);//Exponential Integral Function
0045
0046   G4ecpssrBaseKxsModel(const G4ecpssrBaseKxsModel&) = delete;
0047   G4ecpssrBaseKxsModel & operator = (const G4ecpssrBaseKxsModel &right) = delete;
0048
0049 private:
0050   G4double FunctionFK(G4double k, G4double theta);
0051   G4double LogLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
0052   G4double LinLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
0053
0054   G4double QuadInterpolator(G4double e11,
0055                     G4double e12,
0056                 G4double e21,
0057                 G4double e22,
0058                 G4double x11,
0059                 G4double x12,
0060                 G4double x21,
0061                 G4double x22,
0062                 G4double t1,
0063                 G4double t2,
0064                 G4double t,
0065                 G4double e);
0066
0067   typedef std::map<G4double, std::map<G4double, G4double> > TriDimensionMap;
0068   TriDimensionMap FKData;
0069   std::vector<G4double> dummyVec;
0070
0071   typedef std::map<G4double, std::vector<G4double> > VecMap;
0072   VecMap aVecMap;
0073
0074   G4CrossSectionDataSet* tableC1;
0075   G4CrossSectionDataSet* tableC2;
0076   G4CrossSectionDataSet* tableC3;
0077
0078   G4int verboseLevel;
0079 };
0080
0081 #endif