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: G4MuBremsstrahlung.hh,v 1.26 2005/04/11 10:40:47 vnivanch Exp $ >> 24 // GEANT4 tag $Name: geant4-07-01-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: G4MuBremsstrahlung 31 // File name: G4MuBremsstrahlung 33 // 32 // 34 // Author: Laszlo Urban 33 // Author: Laszlo Urban 35 // 34 // 36 // Creation date: 30.09.1997 35 // Creation date: 30.09.1997 37 // 36 // 38 // Modifications: 37 // Modifications: 39 // 38 // 40 // 10/02/00 modifications , new e.m. structure 39 // 10/02/00 modifications , new e.m. structure, L.Urban 41 // 10-08-01 new methods Store/Retrieve Physics 40 // 10-08-01 new methods Store/Retrieve PhysicsTable (mma) 42 // 29-10-01 all static functions no more inlin 41 // 29-10-01 all static functions no more inlined (mma) 43 // 10-05-02 V.Ivanchenko update to new design 42 // 10-05-02 V.Ivanchenko update to new design 44 // 26-12-02 secondary production moved to deri 43 // 26-12-02 secondary production moved to derived classes (VI) 45 // 24-01-03 Make models region aware (V.Ivanch 44 // 24-01-03 Make models region aware (V.Ivanchenko) 46 // 05-02-03 Fix compilation warnings (V.Ivanch 45 // 05-02-03 Fix compilation warnings (V.Ivanchenko) 47 // 08-08-03 STD substitute standard (V.Ivanch 46 // 08-08-03 STD substitute standard (V.Ivanchenko) 48 // 12-11-03 G4EnergyLossSTD -> G4EnergyLossPro 47 // 12-11-03 G4EnergyLossSTD -> G4EnergyLossProcess (V.Ivanchenko) 49 // 21-01-04 Migrade to G4ParticleChangeForLoss 48 // 21-01-04 Migrade to G4ParticleChangeForLoss (V.Ivanchenko) 50 // 10-02-04 Add lowestKinEnergy (V.Ivanchenko) 49 // 10-02-04 Add lowestKinEnergy (V.Ivanchenko) 51 // 17-08-04 Rename the process "Mu" -> "mu" (V 50 // 17-08-04 Rename the process "Mu" -> "mu" (V.Ivanchenko) 52 // 08-11-04 Migration to new interface of Stor 51 // 08-11-04 Migration to new interface of Store/Retrieve tables (V.Ivantchenko) 53 // 08-04-05 Major optimisation of internal int 52 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko) 54 // 53 // 55 // Class Description: 54 // Class Description: 56 // 55 // 57 // This class manages the Bremsstrahlung proce 56 // This class manages the Bremsstrahlung process for muons. 58 // it inherites from G4VContinuousDiscreteProc 57 // it inherites from G4VContinuousDiscreteProcess via G4VEnergyLossProcess. 59 // 58 // 60 59 61 // ------------------------------------------- 60 // ------------------------------------------------------------------- 62 // 61 // 63 62 64 #ifndef G4MuBremsstrahlung_h 63 #ifndef G4MuBremsstrahlung_h 65 #define G4MuBremsstrahlung_h 1 64 #define G4MuBremsstrahlung_h 1 66 65 67 //....oooOO0OOooo........oooOO0OOooo........oo 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 68 67 69 #include "globals.hh" 68 #include "globals.hh" 70 #include "G4VEnergyLossProcess.hh" 69 #include "G4VEnergyLossProcess.hh" 71 #include "G4VEmModel.hh" 70 #include "G4VEmModel.hh" 72 71 73 //....oooOO0OOooo........oooOO0OOooo........oo 72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 74 73 75 class G4ParticleDefinition; 74 class G4ParticleDefinition; 76 75 77 class G4MuBremsstrahlung : public G4VEnergyLos 76 class G4MuBremsstrahlung : public G4VEnergyLossProcess 78 77 79 { 78 { 80 public: 79 public: 81 80 82 explicit G4MuBremsstrahlung(const G4String& << 81 G4MuBremsstrahlung(const G4String& processName = "muBrems"); 83 82 84 ~G4MuBremsstrahlung() override = default; << 83 virtual ~G4MuBremsstrahlung(); 85 84 86 G4bool IsApplicable(const G4ParticleDefiniti << 85 G4bool IsApplicable(const G4ParticleDefinition& p); 87 86 88 G4double MinPrimaryEnergy(const G4ParticleDe 87 G4double MinPrimaryEnergy(const G4ParticleDefinition* p, 89 const G4Material*, G4double cut) ove << 88 const G4Material*, >> 89 G4double cut); 90 90 91 inline void SetLowestKineticEnergy(G4double << 91 std::vector<G4DynamicParticle*>* SecondariesPostStep( >> 92 G4VEmModel*, >> 93 const G4MaterialCutsCouple*, >> 94 const G4DynamicParticle*, >> 95 G4double& tcut); 92 96 93 // print description in html << 97 // Print out of the class parameters 94 void ProcessDescription(std::ostream&) const << 98 void PrintInfo(); 95 << 96 G4MuBremsstrahlung & operator=(const G4MuBre << 97 G4MuBremsstrahlung(const G4MuBremsstrahlung& << 98 99 99 protected: 100 protected: 100 101 101 void InitialiseEnergyLossProcess(const G4Par 102 void InitialiseEnergyLossProcess(const G4ParticleDefinition*, 102 const G4Par << 103 const G4ParticleDefinition*); >> 104 >> 105 private: >> 106 >> 107 G4MuBremsstrahlung & operator=(const G4MuBremsstrahlung &right); >> 108 G4MuBremsstrahlung(const G4MuBremsstrahlung&); >> 109 >> 110 const G4ParticleDefinition* theParticle; >> 111 const G4ParticleDefinition* theBaseParticle; 103 112 104 G4double lowestKinEnergy; 113 G4double lowestKinEnergy; 105 G4bool isInitialised = false; << 114 G4bool isInitialised; 106 115 107 }; 116 }; 108 117 109 //....oooOO0OOooo........oooOO0OOooo........oo 118 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 119 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 120 >> 121 inline G4bool G4MuBremsstrahlung::IsApplicable(const G4ParticleDefinition& p) >> 122 { >> 123 return (p.GetPDGCharge() != 0.0 && p.GetPDGMass() > 10.0*MeV); >> 124 } >> 125 >> 126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 127 >> 128 inline G4double G4MuBremsstrahlung::MinPrimaryEnergy(const G4ParticleDefinition*, >> 129 const G4Material*, >> 130 G4double) >> 131 { >> 132 return lowestKinEnergy; >> 133 } >> 134 >> 135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 110 136 111 inline void G4MuBremsstrahlung::SetLowestKinet << 137 inline std::vector<G4DynamicParticle*>* G4MuBremsstrahlung::SecondariesPostStep( >> 138 G4VEmModel* model, >> 139 const G4MaterialCutsCouple* couple, >> 140 const G4DynamicParticle* dp, >> 141 G4double& tcut) 112 { 142 { 113 lowestKinEnergy = e; << 143 return model->SampleSecondaries(couple, dp, tcut); 114 } 144 } 115 145 116 //....oooOO0OOooo........oooOO0OOooo........oo 146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 117 147 118 #endif 148 #endif 119 149