| 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 $
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