include/Geant4/G4DNABornExcitationModel2.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
0028 #ifndef G4DNABornExcitationModel2_h
0029 #define G4DNABornExcitationModel2_h 1
0030
0031 #include "G4VEmModel.hh"
0032 #include "G4ParticleChangeForGamma.hh"
0033 #include "G4ProductionCutsTable.hh"
0034
0035 //#include "G4DNACrossSectionDataSet.hh"
0036 #include "G4LogLogInterpolation.hh"
0037 #include "G4Electron.hh"
0038 #include "G4Proton.hh"
0039 #include "G4DNAWaterExcitationStructure.hh"
0040 #include "G4NistManager.hh"
0041
0042 class G4PhysicsVector;
0043
0044 class G4DNABornExcitationModel2: public G4VEmModel
0045 {
0046 public:
0047   G4DNABornExcitationModel2(const G4ParticleDefinition* p = nullptr,
0048                            const G4String& nam = "DNABornExcitationModel");
0049
0050   ~G4DNABornExcitationModel2() override;
0051
0052   G4DNABornExcitationModel2 & operator=(const  G4DNABornExcitationModel2 &right) = delete;
0053   G4DNABornExcitationModel2(const  G4DNABornExcitationModel2&) = delete;
0054
0055   void Initialise(const G4ParticleDefinition*,
0056                           const G4DataVector& = *(new G4DataVector())) override;
0057
0058   G4double CrossSectionPerVolume(const G4Material* material,
0059                                          const G4ParticleDefinition* p,
0060                                          G4double ekin,
0061                                          G4double emin,
0062                                          G4double emax) override;
0063
0064   G4double GetPartialCrossSection(const G4Material*,
0065                                           G4int level,
0066                                           const G4ParticleDefinition*,
0067                                           G4double kineticEnergy) override;
0068
0069   void SampleSecondaries(std::vector<G4DynamicParticle*>*,
0070                                  const G4MaterialCutsCouple*,
0071                                  const G4DynamicParticle*,
0072                                  G4double tmin,
0073                                  G4double maxEnergy) override;
0074
0075   inline void SelectStationary(G4bool input);
0076
0077 protected:
0078
0079   G4ParticleChangeForGamma* fParticleChangeForGamma;
0080
0081 private:
0082
0083   G4bool statCode;
0084
0085   // Water density table
0086   const std::vector<G4double>* fpMolWaterDensity;
0087
0088   G4bool isInitialised{false};
0089   G4int verboseLevel;
0090   const G4ParticleDefinition* fParticleDefinition;
0091
0092   G4double fLowEnergy;
0093   G4double fHighEnergy;
0094
0095   G4PhysicsTable* fTableData{nullptr};
0096   G4PhysicsVector* fTotalXS;
0097   size_t fLastBinCallForFinalXS;
0098
0099   // Partial cross section
0100   G4int RandomSelect(G4double energy);
0101
0102   G4DNAWaterExcitationStructure waterStructure;
0103
0104   //
0105 };
0106
0107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
0108
0109 inline void G4DNABornExcitationModel2::SelectStationary (G4bool input)
0110 {
0111     statCode = input;
0112 }
0113
0114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
0115
0116 #endif