| Geant4 Cross Reference |
>> 1 // This code implementation is the intellectual property of
>> 2 // the GEANT4 collaboration.
1 // 3 //
2 // ******************************************* << 4 // By copying, distributing or modifying the Program (or any work
3 // * License and Disclaimer << 5 // based on the Program) you indicate your acceptance of this statement,
4 // * << 6 // and all its terms.
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 // ------------------------------------------- <<
28 // <<
29 // GEANT4 Class header file <<
30 // <<
31 // <<
32 // File name: G4eplusAnnihilation <<
33 // <<
34 // Author: Vladimir Ivanchenko on base <<
35 // <<
36 // Creation date: 02.08.2004 <<
37 // <<
38 // Modifications: <<
39 // 08-11-04 Migration to new interface of Stor <<
40 // 15-03-05 Update interface according to chan <<
41 // in G4VEmProcess (V.Ivanchenko) <<
42 // 08-04-05 Major optimisation of internal int <<
43 // 04-05-05 Make class to be default (V.Ivanch <<
44 // 04-12-05 SetProposedKineticEnergy(0.) for a <<
45 // 09-08-06 add SetModel(G4VEmModel*) (mma) <<
46 // 12-09-06, move SetModel(G4VEmModel*) in G4V <<
47 // <<
48 // <<
49 // Class Description: <<
50 // <<
51 // This class manages the process of e+ annihi <<
52 // It is possible to enable ApplyCuts and Enta <<
53 // G4EmParameters class using UI commands or C <<
54 // EM splitting or Russian roulette are allowe <<
55 // are enabled. <<
56 // 7 //
57 // ------------------------------------------- << 8 // $Id: G4eplusAnnihilation.hh,v 1.3 1999/12/17 18:26:12 maire Exp $
>> 9 // GEANT4 tag $Name: geant4-03-00 $
58 // 10 //
>> 11 // ------------ G4eplusAnnihilation process ------
>> 12 // by Michel Maire, 7 july 1996
>> 13 //
>> 14 // 10-01-97, crossection table + meanfreepath table, M.Maire
>> 15 // 17-03-97, merge 'in fly' and 'at rest', M.Maire
>> 16 // 31-08-98, new methods SetBining() and PrintInfo()
>> 17 // ------------------------------------------------------------
59 18
60 #ifndef G4eplusAnnihilation_h << 19 // class description
61 #define G4eplusAnnihilation_h 1 << 20 //
62 << 21 // e+ e- ---> gamma gamma
63 #include "G4VEmProcess.hh" << 22 // inherit from G4VRestDiscreteProcess
64 << 23 //
65 class G4ParticleDefinition; <<
66 class G4VPositronAtRestModel; <<
67 <<
68 class G4eplusAnnihilation : public G4VEmProces <<
69 { <<
70 <<
71 public: <<
72 <<
73 explicit G4eplusAnnihilation(const G4String& <<
74 <<
75 ~G4eplusAnnihilation() override; <<
76 <<
77 G4bool IsApplicable(const G4ParticleDefiniti <<
78 <<
79 G4VParticleChange* AtRestDoIt( <<
80 const G4Track& tr <<
81 const G4Step& ste <<
82 <<
83 G4double AtRestGetPhysicalInteractionLength( <<
84 const G4Track& tr <<
85 G4ForceCondition* <<
86 ) override; <<
87 <<
88 // print documentation in html format <<
89 void ProcessDescription(std::ostream&) const <<
90 <<
91 G4eplusAnnihilation & operator=(const G4eplu <<
92 G4eplusAnnihilation(const G4eplusAnnihilatio <<
93 24
94 protected: << 25 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 26 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
95 27
96 void InitialiseProcess(const G4ParticleDefin << 28 #ifndef G4eplusAnnihilation_h
>> 29 #define G4eplusAnnihilation_h 1
97 30
98 // Print out of the class parameters << 31 #include "G4ios.hh"
99 void StreamProcessInfo(std::ostream& outFile << 32 #include "globals.hh"
>> 33 #include "Randomize.hh"
>> 34 #include "G4VRestDiscreteProcess.hh"
>> 35 #include "G4PhysicsTable.hh"
>> 36 #include "G4PhysicsLogVector.hh"
>> 37 #include "G4ElementTable.hh"
>> 38 #include "G4Gamma.hh"
>> 39 #include "G4Positron.hh"
>> 40 #include "G4Step.hh"
100 41
101 private: << 42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
102 43
103 G4VPositronAtRestModel* f2GammaAtRestModel{n << 44 class G4eplusAnnihilation : public G4VRestDiscreteProcess
104 G4VPositronAtRestModel* f3GammaAtRestModel{n << 45
105 G4int fEntanglementModelID; << 46 {
106 G4bool isInitialised{false}; << 47 public: // with description
107 G4bool fEntangled{false}; << 48
108 G4bool fApplyCuts{false}; << 49 G4eplusAnnihilation(const G4String& processName ="annihil");
>> 50
>> 51 ~G4eplusAnnihilation();
>> 52
>> 53 G4bool IsApplicable(const G4ParticleDefinition&);
>> 54 // true for positron only.
>> 55
>> 56 void SetPhysicsTableBining(G4double lowE, G4double highE, G4int nBins);
>> 57 // Allows to define the binning of the PhysicsTables,
>> 58 // before to build them.
>> 59
>> 60 void BuildPhysicsTable(const G4ParticleDefinition& PositronType);
>> 61 // It builds the total CrossSectionPerAtom table, for e+,
>> 62 // and for every element contained in the elementTable.
>> 63 // It builds the MeanFreePath table, for e+,
>> 64 // and for every material contained in the materialTable.
>> 65 // This function overloads a virtual function of the base class.
>> 66 // It is invoked by the G4ParticleWithCuts::SetCut() method.
>> 67
>> 68 void PrintInfoDefinition();
>> 69 // Print few lines of informations about the process: validity range,
>> 70 // origine ..etc..
>> 71 // Invoked by BuildPhysicsTable().
>> 72
>> 73 G4double GetMeanFreePath(const G4Track& aTrack,
>> 74 G4double previousStepSize,
>> 75 G4ForceCondition* condition);
>> 76 // It returns the MeanFreePath of the process for the current track :
>> 77 // (energy, material)
>> 78 // The previousStepSize and G4ForceCondition* are not used.
>> 79 // This function overloads a virtual function of the base class.
>> 80 // It is invoked by the ProcessManager of the Particle.
>> 81
>> 82 G4double GetCrossSectionPerAtom(G4DynamicParticle* aDynamicPositron,
>> 83 G4Element* anElement);
>> 84 // It returns the total CrossSectionPerAtom of the process,
>> 85 // for the current DynamicPositron (energy), in anElement.
>> 86
>> 87 G4VParticleChange* PostStepDoIt(const G4Track& aTrack,
>> 88 const G4Step& aStep);
>> 89 // It computes the final state of the process (at end of step),
>> 90 // returned as a ParticleChange object.
>> 91 // This function overloads a virtual function of the base class.
>> 92 // It is invoked by the ProcessManager of the Particle.
>> 93
>> 94 G4double GetMeanLifeTime(const G4Track& aTrack,
>> 95 G4ForceCondition* condition);
>> 96 // It is invoked by the ProcessManager of the Positron if this
>> 97 // e+ has a kinetic energy null. Then it return 0 to force the
>> 98 // call of AtRestDoIt.
>> 99 // This function overloads a virtual function of the base class.
>> 100
>> 101 G4VParticleChange* AtRestDoIt(const G4Track& aTrack,
>> 102 const G4Step& aStep);
>> 103 // It computes the final state of the process:
>> 104 // e+ (at rest) e- (at rest) ---> gamma gamma,
>> 105 // returned as a ParticleChange object.
>> 106 // This function overloads a virtual function of the base class.
>> 107 // It is invoked by the ProcessManager of the Particle.
>> 108
>> 109 protected:
>> 110
>> 111 virtual G4double ComputeCrossSectionPerAtom(G4double PositKinEnergy,
>> 112 G4double AtomicNumber);
>> 113
>> 114 virtual G4double ComputeMeanFreePath(G4double PositKinEnergy,
>> 115 G4Material* aMaterial);
>> 116
>> 117 private:
>> 118
>> 119 // hide assignment operator as private
>> 120 G4eplusAnnihilation& operator=(const G4eplusAnnihilation &right);
>> 121 G4eplusAnnihilation(const G4eplusAnnihilation& );
>> 122
>> 123 private:
>> 124
>> 125 G4PhysicsTable* theCrossSectionTable; // table for crossection
>> 126 G4PhysicsTable* theMeanFreePathTable;
>> 127
>> 128 G4double LowestEnergyLimit ; // low energy limit of the crossection formula
>> 129 G4double HighestEnergyLimit ; // high energy limit of the crossection formula
>> 130 G4int NumbBinTable ; // number of bins in the crossection table
109 }; 131 };
110 132
111 //....oooOO0OOooo........oooOO0OOooo........oo << 133 #include "G4eplusAnnihilation.icc"
112 << 134
113 #endif 135 #endif
>> 136
114 137