include/Geant4/G4DNARelativisticIonisationModel.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 // $Id: G4DNARelativisticIonisationModel.hh 90057 2015-05-11 22:25:50Z matkara $
0027 //
0028
0029 #ifndef G4DNARelativisticIonisationModel_h
0030 #define G4DNARelativisticIonisationModel_h 1
0031
0032 #include "G4VEmModel.hh"
0033 #include "G4ParticleChangeForGamma.hh"
0034 #include "G4ProductionCutsTable.hh"
0035 #include "G4VAtomDeexcitation.hh"
0036 #include "G4PhysicsFreeVector.hh"
0037
0038 #include "G4LogLogInterpolation.hh"
0039 #include "G4Electron.hh"
0040 #include "G4Proton.hh"
0041 #include "G4NistManager.hh"
0042
0043 #include "G4DNACrossSectionDataSet.hh"
0044 //#include "G4DNAWaterExcitationStructure.hh"
0045
0046 class G4DNARelativisticIonisationModel: public G4VEmModel
0047 {
0048 public:
0049   G4DNARelativisticIonisationModel(const G4ParticleDefinition* p = nullptr,
0050                            const G4String& nam = "DNARelativisticIonisationModel");
0051
0052   ~G4DNARelativisticIonisationModel() override;
0053
0054   G4DNARelativisticIonisationModel & operator
0055           =(const  G4DNARelativisticIonisationModel &right) = delete;
0056   G4DNARelativisticIonisationModel(const  G4DNARelativisticIonisationModel&) = delete;
0057
0058   void Initialise(const G4ParticleDefinition*,
0059                           const G4DataVector& = *(new G4DataVector())) override;
0060
0061   G4double CrossSectionPerVolume(const G4Material* material,
0062                                          const G4ParticleDefinition* p,
0063                                          G4double ekin,
0064                                          G4double emin,
0065                                          G4double emax) override;
0066
0067   virtual G4double GetTotalCrossSection  (const G4Material* material,
0068                                           const G4ParticleDefinition*,
0069                                           G4double kineticEnergy);
0070   G4double GetPartialCrossSection(const G4Material* material,
0071                                           G4int level,
0072                                           const G4ParticleDefinition*,
0073                                           G4double kineticEnergy) override;
0074   virtual G4double GetDifferentialCrossSection(const G4Material* material,
0075                                           const G4ParticleDefinition* particle,
0076                                           G4double kineticEnergy,
0077                                           G4double secondaryEnergy,
0078                                           G4int  level);
0079
0080   void SampleSecondaries(std::vector<G4DynamicParticle*>*,
0081                                  const G4MaterialCutsCouple*,
0082                                  const G4DynamicParticle*,
0083                                  G4double tmin,
0084                                  G4double maxEnergy) override;
0085   virtual void LoadAtomicStates(G4int z, const char *path);
0086   inline  void SelectStationary       (G4bool input){statCode = input;};
0087   inline  void SelectFasterComputation(G4bool input){fasterCode = input;};
0088
0089
0090 protected:
0091
0092   G4ParticleChangeForGamma* fParticleChangeForGamma;
0093
0094 private:
0095
0096   std::vector <G4int   >  iState    [99];
0097   std::vector <G4int   >  iShell    [99];
0098   std::vector <G4int   >  iSubShell [99];
0099   std::vector <G4double>  Nelectrons[99];
0100   std::vector <G4double>  Ebinding  [99];
0101   std::vector <G4double>  Ekinetic  [99];
0102
0103   std::map <G4int, std::vector<G4double> > eVecEZ;
0104
0105   using DeauxDimensionVecMapZ = std::map<G4int, std::map<G4double, std::vector<G4double>>>;
0106   DeauxDimensionVecMapZ eVecEjeEZ;
0107
0108   using TriDimensionVecMapZ = std::map<G4int, std::map<G4int, std::map<G4double, std::vector<G4double>>>>;
0109   TriDimensionVecMapZ eProbaShellMapZ;
0110
0111   using QuadDimensionMapZ = std::map<G4int, std::map<G4int, std::map<G4double, std::map<G4double, G4double>>>>;
0112   QuadDimensionMapZ eDiffCrossSectionDataZ;
0113   QuadDimensionMapZ eEjectedEnergyDataZ;
0114
0115
0116   G4double     fLowEnergyLimit=0.;
0117   G4double     fHighEnergyLimit=0.;
0118
0119   G4bool       isInitialised=false;
0120   G4bool       statCode=false;
0121   G4bool       fasterCode=false;
0122   G4int        verboseLevel=0;
0123
0124   const std::vector<G4double>*  fMaterialDensity=nullptr;
0125   const  G4ParticleDefinition*  fParticleDefinition=nullptr;
0126   G4VAtomDeexcitation*          fAtomDeexcitation=nullptr;
0127
0128   G4int RandomSelect(const G4Material* material,
0129                      const G4ParticleDefinition*,
0130                      G4double kineticEnergy);
0131
0132   G4double       GetEjectedElectronEnergy   (
0133                  const G4Material* material,
0134                  const G4ParticleDefinition* ,
0135                  G4double energy,
0136                  G4int shell          );
0137   G4ThreeVector  GetEjectedElectronDirection(
0138                  const G4ParticleDefinition* ,
0139                  G4double energy,G4double secondaryenergy);
0140
0141   G4double Interpolate     (G4double e1 ,G4double e2 ,G4double e  ,
0142                             G4double xs1, G4double xs2);
0143   G4double QuadInterpolator(G4double e11,G4double e12,G4double e21,G4double e22,
0144                             G4double x11,G4double x12,G4double x21,G4double x22,
0145                             G4double t1 ,G4double t2 ,G4double t ,G4double e);
0146
0147 };
0148
0149 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
0150
0151 #endif