Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 // Author: Luciano Pandola 28 // 29 // History: 30 // ----------- 31 // 23 Aug 2010 L. Pandola 1st implementati 32 // 24 May 2011 L. Pandola Renamed (make v2 33 // 02 Oct 2013 L. Pandola Migrated to mult 34 // 35 // ------------------------------------------- 36 // 37 // Class description: 38 // Low Energy Electromagnetic Physics, e+ and 39 // with Penelope Model, version 2008 40 // ------------------------------------------- 41 42 #ifndef G4PENELOPEBREMSSTRAHLUNGMODEL_HH 43 #define G4PENELOPEBREMSSTRAHLUNGMODEL_HH 1 44 45 #include "globals.hh" 46 #include "G4VEmModel.hh" 47 #include "G4DataVector.hh" 48 #include "G4ParticleChangeForLoss.hh" 49 50 51 class G4PhysicsFreeVector; 52 class G4PhysicsLogVector; 53 class G4ParticleDefinition; 54 class G4DynamicParticle; 55 class G4MaterialCutsCouple; 56 class G4Material; 57 class G4PenelopeOscillatorManager; 58 class G4PenelopeCrossSection; 59 class G4PenelopeBremsstrahlungFS; 60 class G4PenelopeBremsstrahlungAngular; 61 62 class G4PenelopeBremsstrahlungModel : public G 63 { 64 public: 65 explicit G4PenelopeBremsstrahlungModel(const 66 const G4String& processName ="PenBrem") 67 virtual ~G4PenelopeBremsstrahlungModel(); 68 69 void Initialise(const G4ParticleDefinition*, 70 void InitialiseLocal(const G4ParticleDefinit 71 G4VEmModel*) override; 72 73 //DUMMY METHOD 74 G4double ComputeCrossSectionPerAtom(const G4 75 G4double kinEnergy, 76 G4double Z, 77 G4double A=0, 78 G4double cut=0, 79 G4double emax=DBL_MAX) override; 80 81 82 G4double CrossSectionPerVolume(const G4Mater 83 const G4ParticleDefinition* thePartic 84 G4double kineticEnergy, 85 G4double cutEnergy, 86 G4double maxEnergy = DBL_MAX) overrid 87 88 void SampleSecondaries(std::vector<G4Dynamic 89 const G4MaterialCutsCouple*, 90 const G4DynamicParticle*, 91 G4double tmin, 92 G4double maxEnergy) override; 93 94 G4double ComputeDEDXPerVolume(const G4Materi 95 const G4ParticleDefinition*, 96 G4double kineticEnergy, 97 G4double cutEnergy) override; 98 99 // Min cut in kinetic energy allowed by the 100 G4double MinEnergyCut(const G4ParticleDefini 101 const G4MaterialCutsCouple*) override; 102 103 void SetVerbosityLevel(G4int lev){fVerboseLe 104 G4int GetVerbosityLevel(){return fVerboseLev 105 106 G4PenelopeBremsstrahlungModel & operator= 107 (const G4PenelopeBremsstrahlungModel &right) 108 G4PenelopeBremsstrahlungModel(const G4Penelo 109 110 protected: 111 G4ParticleChangeForLoss* fParticleChange; 112 const G4ParticleDefinition* fParticle; 113 114 private: 115 void ClearTables(); 116 G4PenelopeCrossSection* GetCrossSectionTable 117 const G4Material*,G4double cut 118 void SetParticle(const G4ParticleDefinition* 119 // 120 //Members to handle and store cross section 121 void BuildXSTable(const G4Material* material 122 G4double GetPositronXSCorrection(const G4Mat 123 124 //Helpers 125 G4PenelopeOscillatorManager* fOscManager; 126 G4PenelopeBremsstrahlungFS* fPenelopeFSHelpe 127 G4PenelopeBremsstrahlungAngular* fPenelopeAn 128 129 //This is the main energy grid 130 G4PhysicsLogVector* fEnergyGrid; 131 size_t nBins; 132 //G4PenelopeCrossSection takes care of the l 133 std::map< std::pair<const G4Material*,G4doub 134 std::map< std::pair<const G4Material*,G4doub 135 136 //Intrinsic energy limits of the model: cann 137 G4double fIntrinsicLowEnergyLimit; 138 G4double fIntrinsicHighEnergyLimit; 139 140 G4int fVerboseLevel; 141 G4bool fIsInitialised; 142 //Used only for G4EmCalculator and Unit Test 143 G4bool fLocalTable; 144 }; 145 146 #endif 147 148