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: G4hhIonisation.hh,v 1.1 2005/10/30 15:40:05 vnivanch Exp $ >> 24 // GEANT4 tag $Name: geant4-08-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: G4hhIonisation 31 // File name: G4hhIonisation 33 // 32 // 34 // Author: Vladimir Ivanchenko 33 // Author: Vladimir Ivanchenko 35 // 34 // 36 // Creation date: 30.09.2005 35 // Creation date: 30.09.2005 37 // 36 // 38 // Modifications: 37 // Modifications: 39 // 38 // 40 // 39 // 41 // Class Description: 40 // Class Description: 42 // 41 // 43 // This class manages the ionisation process f 42 // This class manages the ionisation process for exsotic hadrons 44 // without any delta electrons production 43 // without any delta electrons production 45 // 44 // 46 45 47 // ------------------------------------------- 46 // ------------------------------------------------------------------- 48 // 47 // 49 48 50 #ifndef G4hhIonisation_h 49 #ifndef G4hhIonisation_h 51 #define G4hhIonisation_h 1 50 #define G4hhIonisation_h 1 52 51 53 #include "G4VEnergyLossProcess.hh" 52 #include "G4VEnergyLossProcess.hh" 54 #include "globals.hh" 53 #include "globals.hh" 55 #include "G4VEmModel.hh" 54 #include "G4VEmModel.hh" 56 55 57 class G4Material; 56 class G4Material; 58 class G4VEmFluctuationModel; 57 class G4VEmFluctuationModel; 59 58 60 class G4hhIonisation : public G4VEnergyLossPro 59 class G4hhIonisation : public G4VEnergyLossProcess 61 { 60 { 62 61 63 public: 62 public: 64 63 65 explicit G4hhIonisation(const G4String& name << 64 G4hhIonisation(const G4String& name = "hhIoni"); 66 65 67 ~G4hhIonisation() override; << 66 virtual ~G4hhIonisation(); 68 67 69 G4bool IsApplicable(const G4ParticleDefiniti << 68 G4bool IsApplicable(const G4ParticleDefinition& p); 70 69 71 G4double MinPrimaryEnergy(const G4ParticleDe 70 G4double MinPrimaryEnergy(const G4ParticleDefinition* p, 72 const G4Material*, G4double cut) ove << 71 const G4Material*, G4double cut); 73 72 74 // print description in html << 73 // Print out of the class parameters 75 void ProcessDescription(std::ostream&) const << 74 virtual void PrintInfo(); 76 << 77 // hide assignment operator << 78 G4hhIonisation & operator=(const G4hhIonisat << 79 G4hhIonisation(const G4hhIonisation&) = dele << 80 75 81 protected: 76 protected: 82 77 83 void InitialiseEnergyLossProcess(const G4Par << 78 std::vector<G4DynamicParticle*>* SecondariesPostStep( 84 const G4Par << 79 G4VEmModel*, >> 80 const G4MaterialCutsCouple*, >> 81 const G4DynamicParticle*, >> 82 G4double&); >> 83 >> 84 virtual void InitialiseEnergyLossProcess(const G4ParticleDefinition*, >> 85 const G4ParticleDefinition*); 85 86 86 private: 87 private: 87 88 88 const G4ParticleDefinition* theParticle = nu << 89 // hide assignment operator 89 G4VEmFluctuationModel* flucModel = nullptr; << 90 G4hhIonisation & operator=(const G4hhIonisation &right); >> 91 G4hhIonisation(const G4hhIonisation&); >> 92 >> 93 G4double mass; >> 94 G4double ratio; >> 95 G4double minKinEnergy; >> 96 >> 97 const G4ParticleDefinition* theParticle; >> 98 const G4ParticleDefinition* theBaseParticle; >> 99 G4VEmFluctuationModel* flucModel; 90 100 91 G4double mass = 0.0; << 101 G4bool isInitialised; 92 G4double ratio = 0.0; << 102 >> 103 G4double eth; 93 104 94 G4bool isInitialised = false; << 95 }; 105 }; >> 106 >> 107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 108 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 109 >> 110 inline G4bool G4hhIonisation::IsApplicable(const G4ParticleDefinition& p) >> 111 { >> 112 return (p.GetPDGCharge() != 0.0 && p.GetPDGMass() > 100.0*MeV && >> 113 !p.IsShortLived()); >> 114 } >> 115 >> 116 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 117 >> 118 inline G4double G4hhIonisation::MinPrimaryEnergy(const G4ParticleDefinition*, >> 119 const G4Material*, >> 120 G4double cut) >> 121 { >> 122 G4double x = 0.5*cut/electron_mass_c2; >> 123 G4double y = electron_mass_c2/mass; >> 124 G4double g = x*y + std::sqrt((1. + x)*(1. + x*y*y)); >> 125 return mass*(g - 1.0); >> 126 } >> 127 >> 128 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 129 >> 130 inline std::vector<G4DynamicParticle*>* G4hhIonisation::SecondariesPostStep( >> 131 G4VEmModel*, >> 132 const G4MaterialCutsCouple*, >> 133 const G4DynamicParticle*, >> 134 G4double&) >> 135 { >> 136 return 0; >> 137 } 96 138 97 //....oooOO0OOooo........oooOO0OOooo........oo 139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 98 140 99 #endif 141 #endif 100 142