Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4SynchrotronRadiationInMat.hh
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // $Id: G4SynchrotronRadiationInMat.hh 97385 2016-06-02 09:59:53Z gcosmo $
28 //
29 // ------------------------------------------------------------
30 // GEANT 4 class header file
31 // CERN Geneva Switzerland
32 //
33 //
34 // History:
35 // 21-5-98 1 version , V. Grichine
36 // 28-05-01, V.Ivanchenko minor changes to provide ANSI -wall compilation
37 // 19-05-06, V.Ivanchenko rename from G4SynchrotronRadiation
38 //
39 //
40 // ------------------------------------------------------------
41 
42 #ifndef G4SynchrotronRadiationInMat_h
43 #define G4SynchrotronRadiationInMat_h 1
44 
45 #include "G4ios.hh"
46 #include "globals.hh"
47 #include "Randomize.hh"
48 #include "G4VDiscreteProcess.hh"
50 #include "G4FieldManager.hh"
51 #include "G4Field.hh"
52 #include "G4ThreeVector.hh"
53 #include "G4PropagatorInField.hh"
54 
55 #include "G4Track.hh"
56 #include "G4Step.hh"
57 
58 
59 #include "G4Gamma.hh"
60 #include "G4Electron.hh"
61 #include "G4Positron.hh"
62 
63 
64 #include "G4PhysicsTable.hh"
65 #include "G4PhysicsLogVector.hh"
66 
67 
69 {
70 public:
71 
72  explicit G4SynchrotronRadiationInMat(const G4String& processName =
73  "SynchrotronRadiation",
75 
77 
78 private:
79 
81  operator=(const G4SynchrotronRadiationInMat &right) = delete;
83 
84 public:
85 
86  G4double GetMeanFreePath( const G4Track& track,
87  G4double previousStepSize,
88  G4ForceCondition* condition ) override;
89 
90  G4VParticleChange *PostStepDoIt( const G4Track& track,
91  const G4Step& Step ) override;
92 
93  G4double GetPhotonEnergy( const G4Track& trackData,
94  const G4Step& stepData );
95 
97 
100 
103 
107 
108  G4bool IsApplicable(const G4ParticleDefinition&) override;
109 
110  static G4double GetLambdaConst(){ return fLambdaConst; };
111  static G4double GetEnergyConst(){ return fEnergyConst; };
112 
113  void SetRootNumber(G4int rn){ fRootNumber = rn; };
114  void SetVerboseLevel(G4int v){ fVerboseLevel = v; };
115  void SetKsi(G4double ksi){ fKsi = ksi; };
116  void SetEta(G4double eta){ fEta = eta; };
117  void SetPsiGamma(G4double psg){ fPsiGamma = psg; };
118  void SetOrderAngleK(G4double ord){ fOrderAngleK = ord; }; // should be 1/3 or 2/3
119 
120 private:
121 
122  static const G4double fLambdaConst;
123 
124  static const G4double fEnergyConst;
125 
126  static const G4double fIntegralProbabilityOfSR[200];
127 
128  const G4double
129  LowestKineticEnergy; // low energy limit of the cross-section formula
130 
131  //const G4double
132  //HighestKineticEnergy; // high energy limit of the cross-section formula
133 
134  //G4int TotBin; // number of bins in the tables
135 
136  G4double CutInRange;
137 
138  const G4ParticleDefinition* theGamma;
139  const G4ParticleDefinition* theElectron;
140  const G4ParticleDefinition* thePositron;
141 
142  //const G4double* GammaCutInKineticEnergy;
143  //const G4double* ElectronCutInKineticEnergy;
144  //const G4double* PositronCutInKineticEnergy;
145  //const G4double* ParticleCutInKineticEnergy;
146 
147 
148  G4double GammaCutInKineticEnergyNow;
149  G4double ElectronCutInKineticEnergyNow;
150  G4double PositronCutInKineticEnergyNow;
151  G4double ParticleCutInKineticEnergyNow;
152 
153  G4double fAlpha;
154  G4int fRootNumber;
155  G4double fKsi; // omega/omega_c
156  G4double fPsiGamma; // Psi-angle*gamma
157  G4double fEta; //
158  G4double fOrderAngleK; // 1/3 or 2/3
159 
160 
161  G4int fVerboseLevel;
162  G4PropagatorInField* fFieldPropagator;
163 
164 };
165 
166 #endif // end of G4SynchrotronRadiationInMat.hh
167 
G4double condition(const G4ErrorSymMatrix &m)
G4VParticleChange * PostStepDoIt(const G4Track &track, const G4Step &Step) override
G4double GetMeanFreePath(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
int G4int
Definition: G4Types.hh:78
G4bool IsApplicable(const G4ParticleDefinition &) override
G4SynchrotronRadiationInMat(const G4String &processName="SynchrotronRadiation", G4ProcessType type=fElectromagnetic)
G4double GetPhotonEnergy(const G4Track &trackData, const G4Step &stepData)
bool G4bool
Definition: G4Types.hh:79
Definition: G4Step.hh:76
tuple v
Definition: test.py:18
Definition: Step.hh:41
double G4double
Definition: G4Types.hh:76
G4ForceCondition
G4double GetRandomEnergySR(G4double, G4double)
G4ProcessType