Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // >> 26 // $Id: G4PenelopeIonisationModel.hh,v 1.1 2010-07-28 07:12:13 pandola Exp $ >> 27 // GEANT4 tag $Name: not supported by cvs2svn $ 26 // 28 // 27 // Author: Luciano Pandola 29 // Author: Luciano Pandola 28 // 30 // 29 // History: 31 // History: 30 // ----------- 32 // ----------- 31 // 30 Mar 2010 L. Pandola 1st implementati 33 // 30 Mar 2010 L. Pandola 1st implementation. 32 // 25 May 2011 L. Pandola Renamed (make v2 << 34 // 25 May 2011 L. Pandola Renamed (make v2008 as default Penelope) 33 // 09 Mar 2012 L. Pandola Moved the manage << 34 // cross sections t << 35 // 07 Oct 2013 L. Pandola Migration to MT << 36 // 23 Jun 2015 L. Pandola Added private me << 37 // 35 // 38 // ------------------------------------------- 36 // ------------------------------------------------------------------- 39 // 37 // 40 // Class description: 38 // Class description: 41 // Low Energy Electromagnetic Physics, e+ and 39 // Low Energy Electromagnetic Physics, e+ and e- ionisation 42 // with Penelope Model, version 2008 40 // with Penelope Model, version 2008 43 // ------------------------------------------- 41 // ------------------------------------------------------------------- 44 42 45 #ifndef G4PENELOPEIONISATIONMODEL_HH 43 #ifndef G4PENELOPEIONISATIONMODEL_HH 46 #define G4PENELOPEIONISATIONMODEL_HH 1 44 #define G4PENELOPEIONISATIONMODEL_HH 1 47 45 48 #include "globals.hh" 46 #include "globals.hh" 49 #include "G4VEmModel.hh" 47 #include "G4VEmModel.hh" 50 #include "G4DataVector.hh" 48 #include "G4DataVector.hh" 51 #include "G4ParticleChangeForLoss.hh" 49 #include "G4ParticleChangeForLoss.hh" 52 #include "G4VAtomDeexcitation.hh" << 50 #include "G4AtomicDeexcitation.hh" 53 51 54 class G4PhysicsFreeVector; 52 class G4PhysicsFreeVector; 55 class G4PhysicsLogVector; 53 class G4PhysicsLogVector; 56 class G4ParticleDefinition; 54 class G4ParticleDefinition; 57 class G4DynamicParticle; 55 class G4DynamicParticle; 58 class G4MaterialCutsCouple; 56 class G4MaterialCutsCouple; 59 class G4Material; 57 class G4Material; 60 class G4PenelopeOscillatorManager; 58 class G4PenelopeOscillatorManager; 61 class G4PenelopeOscillator; 59 class G4PenelopeOscillator; 62 class G4PenelopeCrossSection; 60 class G4PenelopeCrossSection; 63 class G4PenelopeIonisationXSHandler; << 64 61 65 class G4PenelopeIonisationModel : public G4VEm 62 class G4PenelopeIonisationModel : public G4VEmModel 66 { 63 { >> 64 67 public: 65 public: 68 explicit G4PenelopeIonisationModel(const G4P << 66 >> 67 G4PenelopeIonisationModel(const G4ParticleDefinition* p=0, 69 const G4String& processName ="PenIon 68 const G4String& processName ="PenIoni"); >> 69 70 virtual ~G4PenelopeIonisationModel(); 70 virtual ~G4PenelopeIonisationModel(); 71 71 72 void Initialise(const G4ParticleDefinition*, << 72 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&); 73 void InitialiseLocal(const G4ParticleDefinit << 74 G4VEmModel*) override; << 75 73 76 //*This is a dummy method. Never inkoved by 74 //*This is a dummy method. Never inkoved by the tracking, it just issues 77 //*a warning if one tries to get Cross Secti 75 //*a warning if one tries to get Cross Sections per Atom via the 78 //*G4EmCalculator. 76 //*G4EmCalculator. 79 G4double ComputeCrossSectionPerAtom(const G4 << 77 virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*, 80 G4double, << 78 G4double, 81 G4double, << 79 G4double, 82 G4double, << 80 G4double, 83 G4double, << 81 G4double, 84 G4double) override; << 82 G4double); 85 << 83 86 G4double CrossSectionPerVolume(const G4Mater << 84 virtual G4double CrossSectionPerVolume(const G4Material* material, 87 const G4ParticleDefinition* << 85 const G4ParticleDefinition* theParticle, 88 theParticle, << 86 G4double kineticEnergy, 89 G4double kineticEnergy, << 87 G4double cutEnergy, 90 G4double cutEnergy, << 88 G4double maxEnergy = DBL_MAX); 91 G4double maxEnergy = DBL_MAX) overrid << 92 89 93 void SampleSecondaries(std::vector<G4Dynamic << 90 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*, 94 const G4MaterialCutsCouple*, << 91 const G4MaterialCutsCouple*, 95 const G4DynamicParticle*, << 92 const G4DynamicParticle*, 96 G4double tmin, << 93 G4double tmin, 97 G4double maxEnergy) override; << 94 G4double maxEnergy); 98 95 99 G4double ComputeDEDXPerVolume(const G4Materi << 96 virtual G4double ComputeDEDXPerVolume(const G4Material*, 100 const G4ParticleDefinition*, << 97 const G4ParticleDefinition*, 101 G4double kineticEnergy, << 98 G4double kineticEnergy, 102 G4double cutEnergy) override; << 99 G4double cutEnergy); 103 100 104 // Min cut in kinetic energy allowed by the 101 // Min cut in kinetic energy allowed by the model 105 G4double MinEnergyCut(const G4ParticleDefini << 102 virtual G4double MinEnergyCut(const G4ParticleDefinition*, 106 const G4MaterialCutsCouple*) override; << 103 const G4MaterialCutsCouple*); 107 104 108 void SetVerbosityLevel(G4int lev){fVerboseLe << 105 void SetVerbosityLevel(G4int lev){verboseLevel = lev;}; 109 G4int GetVerbosityLevel(){return fVerboseLev << 106 G4int GetVerbosityLevel(){return verboseLevel;}; 110 << 107 111 G4PenelopeIonisationModel & operator=(const << 108 void ActivateAuger(G4bool); 112 G4PenelopeIonisationModel(const G4PenelopeIo << 109 >> 110 G4double GetDensityCorrection(const G4Material*,G4double energy); 113 111 114 protected: 112 protected: 115 G4ParticleChangeForLoss* fParticleChange; 113 G4ParticleChangeForLoss* fParticleChange; 116 const G4ParticleDefinition* fParticle; << 117 114 118 private: 115 private: 119 void SetParticle(const G4ParticleDefinition* << 116 void ClearTables(); 120 void SampleFinalStateElectron(const G4Materi << 117 121 G4double cutEnergy, << 118 G4PenelopeIonisationModel & operator=(const G4PenelopeIonisationModel &right); 122 G4double kineticEnergy); << 119 G4PenelopeIonisationModel(const G4PenelopeIonisationModel&); 123 void SampleFinalStatePositron(const G4Materi << 120 124 G4double cutEnergy, << 121 G4PenelopeCrossSection* GetCrossSectionTableForCouple(const G4ParticleDefinition*, 125 G4double kineticEnergy); << 122 const G4Material*,G4double cut); 126 << 123 127 G4PenelopeOscillatorManager* fOscManager; << 124 void BuildXSTable(const G4Material*,G4double cut, 128 G4PenelopeIonisationXSHandler* fCrossSection << 125 const G4ParticleDefinition*); 129 G4VAtomDeexcitation* fAtomDeexcitation; << 126 130 << 127 131 G4double fKineticEnergy1; << 128 void BuildDeltaTable(const G4Material*); 132 G4double fCosThetaPrimary; << 129 133 G4double fEnergySecondary; << 130 G4DataVector* ComputeShellCrossSectionsElectron(G4PenelopeOscillator* , 134 G4double fCosThetaSecondary; << 131 G4double energy,G4double cut, >> 132 G4double delta); >> 133 >> 134 G4DataVector* ComputeShellCrossSectionsPositron(G4PenelopeOscillator* , >> 135 G4double energy,G4double cut, >> 136 G4double delta); >> 137 >> 138 void SampleFinalStateElectron(const G4Material*,G4double cutEnergy,G4double kineticEnergy); >> 139 void SampleFinalStatePositron(const G4Material*,G4double cutEnergy,G4double kineticEnergy); 135 140 136 //Intrinsic energy limits of the model: cann 141 //Intrinsic energy limits of the model: cannot be extended by the parent process 137 G4double fIntrinsicLowEnergyLimit; 142 G4double fIntrinsicLowEnergyLimit; 138 G4double fIntrinsicHighEnergyLimit; 143 G4double fIntrinsicHighEnergyLimit; 139 144 140 G4int fVerboseLevel; << 145 G4int verboseLevel; 141 G4int fTargetOscillator; << 146 142 size_t fNBins; << 147 G4bool isInitialised; 143 G4bool fIsInitialised; << 148 144 G4bool fPIXEflag; << 149 G4AtomicDeexcitation deexcitationManager; 145 //Used only for G4EmCalculator and Unit Test << 150 146 G4bool fLocalTable; << 151 G4double kineticEnergy1; >> 152 G4double cosThetaPrimary; >> 153 G4double energySecondary; >> 154 G4double cosThetaSecondary; >> 155 G4int targetOscillator; >> 156 >> 157 G4PenelopeOscillatorManager* oscManager; >> 158 >> 159 //G4PenelopeCrossSection takes care of the logs >> 160 std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*> *XSTableElectron; >> 161 std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*> *XSTablePositron; >> 162 >> 163 //delta vs. log(energy) >> 164 std::map<const G4Material*,G4PhysicsFreeVector*> *theDeltaTable; >> 165 G4PhysicsLogVector* energyGrid; >> 166 >> 167 size_t nBins; >> 168 >> 169 147 }; 170 }; 148 171 149 #endif 172 #endif 150 173 151 174