Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // >> 23 // $Id: G4EmMultiModel.hh,v 1.2 2004/10/25 13:34:48 vnivanch Exp $ >> 24 // GEANT4 tag $Name: geant4-07-00-patch-01 $ 26 // 25 // 27 // ------------------------------------------- 26 // ------------------------------------------------------------------- 28 // 27 // 29 // GEANT4 Class header file 28 // GEANT4 Class header file 30 // 29 // 31 // 30 // 32 // File name: G4EmMultiModel 31 // File name: G4EmMultiModel 33 // 32 // 34 // Author: Vladimir Ivanchenko 33 // Author: Vladimir Ivanchenko 35 // 34 // 36 // Creation date: 03.05.2004 35 // Creation date: 03.05.2004 37 // 36 // 38 // Modifications: 37 // Modifications: 39 // 15-04-05 optimize internal interface (V.Iva << 38 // 40 // 04-07-10 updated interfaces according to g4 << 39 // 41 // 40 // 42 // Class Description: 41 // Class Description: 43 // 42 // 44 // EM model using several G4VEmModels for the << 43 // Energy loss model using several G4VEmModels 45 44 46 // ------------------------------------------- 45 // ------------------------------------------------------------------- 47 // 46 // 48 47 49 #ifndef G4EmMultiModel_h 48 #ifndef G4EmMultiModel_h 50 #define G4EmMultiModel_h 1 49 #define G4EmMultiModel_h 1 51 50 52 #include "globals.hh" 51 #include "globals.hh" 53 #include "G4VEmModel.hh" 52 #include "G4VEmModel.hh" 54 #include "G4ParticleDefinition.hh" << 55 #include "G4MaterialCutsCouple.hh" << 56 #include "G4Material.hh" << 57 #include <vector> 53 #include <vector> 58 54 >> 55 class G4Region; >> 56 class G4PhysicsTable; 59 class G4DynamicParticle; 57 class G4DynamicParticle; 60 58 61 class G4EmMultiModel : public G4VEmModel 59 class G4EmMultiModel : public G4VEmModel 62 { 60 { >> 61 63 public: 62 public: 64 63 65 explicit G4EmMultiModel(const G4String& nam << 64 G4EmMultiModel(const G4String& nam = "MultiModel"); 66 65 67 ~G4EmMultiModel() override; << 66 virtual ~G4EmMultiModel(); 68 67 69 void AddModel(G4VEmModel*); << 68 void Initialise(const G4ParticleDefinition*, const G4DataVector&); 70 69 71 void Initialise(const G4ParticleDefinition*, << 70 G4double HighEnergyLimit(const G4ParticleDefinition*) {return highKinEnergy;}; 72 71 73 G4double ComputeDEDXPerVolume(const G4Materi << 72 G4double LowEnergyLimit(const G4ParticleDefinition*) {return lowKinEnergy;}; 74 const G4ParticleDefinition*, << 75 G4double kineticEnergy, << 76 G4double cutEnergy) final; << 77 << 78 // main method to compute cross section per << 79 G4double ComputeCrossSectionPerAtom(const G4 << 80 G4double kinEnergy, << 81 G4double Z, << 82 G4double A = 0., /* amu */ << 83 G4double cutEnergy = 0.0, << 84 G4double maxEnergy = DBL_MAX) fi << 85 << 86 void SampleSecondaries(std::vector<G4Dynamic << 87 const G4MaterialCutsCouple*, << 88 const G4DynamicParticle*, << 89 G4double tmin, << 90 G4double tmax) final; << 91 73 92 // hide assignment operator << 74 void SetHighEnergyLimit(G4double e) {highKinEnergy = e;}; 93 G4EmMultiModel & operator=(const G4EmMultiM << 75 94 G4EmMultiModel(const G4EmMultiModel&) = del << 76 void SetLowEnergyLimit(G4double e) {lowKinEnergy = e;}; >> 77 >> 78 G4double MinEnergyCut(const G4ParticleDefinition*, >> 79 const G4MaterialCutsCouple*); >> 80 >> 81 G4bool IsInCharge(const G4ParticleDefinition*); >> 82 >> 83 G4double ComputeDEDX(const G4MaterialCutsCouple*, >> 84 const G4ParticleDefinition*, >> 85 G4double kineticEnergy, >> 86 G4double cutEnergy); >> 87 >> 88 G4double CrossSection(const G4MaterialCutsCouple*, >> 89 const G4ParticleDefinition*, >> 90 G4double kineticEnergy, >> 91 G4double cutEnergy, >> 92 G4double maxEnergy); >> 93 >> 94 G4DynamicParticle* SampleSecondary( >> 95 const G4MaterialCutsCouple*, >> 96 const G4DynamicParticle*, >> 97 G4double tmin, >> 98 G4double tmax); >> 99 >> 100 std::vector<G4DynamicParticle*>* SampleSecondaries( >> 101 const G4MaterialCutsCouple*, >> 102 const G4DynamicParticle*, >> 103 G4double tmin, >> 104 G4double tmax); >> 105 >> 106 G4double MaxSecondaryEnergy(const G4DynamicParticle* dynParticle); >> 107 >> 108 void DefineForRegion(const G4Region*); >> 109 >> 110 void SetDynamicParticle(const G4DynamicParticle*); >> 111 >> 112 void AddModel(G4VEmModel*, G4double tmin, G4double tmax); >> 113 >> 114 protected: >> 115 >> 116 G4double MaxSecondaryEnergy(const G4ParticleDefinition*, >> 117 G4double kineticEnergy); 95 118 96 private: 119 private: 97 120 98 G4int nModels = 0; << 121 // hide assignment operator >> 122 G4EmMultiModel & operator=(const G4EmMultiModel &right); >> 123 G4EmMultiModel(const G4EmMultiModel&); >> 124 >> 125 G4int nModels; 99 std::vector<G4VEmModel*> model; 126 std::vector<G4VEmModel*> model; 100 std::vector<G4double> cross_section; << 127 G4DataVector tsecmin; >> 128 G4DataVector cross_section; >> 129 >> 130 G4double highKinEnergy; >> 131 G4double lowKinEnergy; 101 132 102 }; 133 }; 103 134 104 #endif 135 #endif 105 136 106 137