Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/xrays/include/G4VXTRenergyLoss.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 //
 28 // base class for 'fast' parametrisation model describing X-ray transition
 29 // created in some G4Envelope. Angular distribuiton is very rough !!! (see DoIt
 30 // method
 31 //
 32 // History:
 33 // 06.10.05  V. Grichine first step to discrete process
 34 // 15.01.02  V. Grichine first version
 35 // 28.07.05, P.Gumplinger add G4ProcessType to constructor
 36 // 28.09.07, V.Ivanchenko general cleanup without change of algorithms
 37 // 19.09.21, V. Grichine, set/get functions for angle anf energy ranges and number of bins
 38 
 39 #ifndef G4VXTRenergyLoss_h
 40 #define G4VXTRenergyLoss_h 1
 41 
 42 #include "globals.hh"
 43 #include "G4Gamma.hh"
 44 #include "G4LogicalVolume.hh"
 45 #include "G4Material.hh"
 46 #include "G4ParticleChange.hh"
 47 #include "G4PhysicsTable.hh"
 48 #include "G4Step.hh"
 49 #include "G4Track.hh"
 50 #include "G4VDiscreteProcess.hh"
 51 
 52 class G4SandiaTable;
 53 class G4VParticleChange;
 54 class G4PhysicsFreeVector;
 55 class G4PhysicsLinearVector;
 56 class G4PhysicsLogVector;
 57 
 58 class G4VXTRenergyLoss : public G4VDiscreteProcess
 59 {
 60  public:
 61   explicit G4VXTRenergyLoss(G4LogicalVolume* anEnvelope, G4Material*,
 62                             G4Material*, G4double, G4double, G4int,
 63                             const G4String& processName = "XTRenergyLoss",
 64                             G4ProcessType type          = fElectromagnetic);
 65   virtual ~G4VXTRenergyLoss();
 66 
 67   virtual void ProcessDescription(std::ostream&) const override;
 68   virtual void DumpInfo() const override { ProcessDescription(G4cout); };
 69 
 70   G4VXTRenergyLoss(G4VXTRenergyLoss&) = delete;
 71   G4VXTRenergyLoss& operator=(const G4VXTRenergyLoss& right) = delete;
 72 
 73   // Virtual methods to be implemented in inherited particular TR radiators
 74   virtual G4double GetStackFactor(G4double energy, G4double gamma,
 75                                   G4double varAngle);
 76 
 77   virtual G4bool IsApplicable(const G4ParticleDefinition&) override;
 78 
 79   virtual G4VParticleChange* PostStepDoIt(const G4Track& aTrack,
 80                                           const G4Step& aStep) override;
 81 
 82   virtual G4double GetMeanFreePath(const G4Track& aTrack,
 83                                    G4double previousStepSize,
 84                                    G4ForceCondition* condition) override;
 85 
 86   virtual void BuildPhysicsTable(const G4ParticleDefinition&) override;
 87   void BuildEnergyTable();
 88   void BuildAngleForEnergyBank();
 89 
 90   void BuildTable(){};
 91   void BuildAngleTable();
 92   void BuildGlobalAngleTable();
 93 
 94   G4complex OneInterfaceXTRdEdx(G4double energy, G4double gamma,
 95                                 G4double varAngle);
 96 
 97   G4double SpectralAngleXTRdEdx(G4double varAngle);
 98 
 99   virtual G4double SpectralXTRdEdx(G4double energy);
100 
101   G4double AngleSpectralXTRdEdx(G4double energy);
102 
103   G4double AngleXTRdEdx(G4double varAngle);
104 
105   G4double OneBoundaryXTRNdensity(G4double energy, G4double gamma,
106                                   G4double varAngle) const;
107 
108   // for photon energy distribution tables
109   G4double XTRNSpectralAngleDensity(G4double varAngle);
110   G4double XTRNSpectralDensity(G4double energy);
111 
112   // for photon angle distribution tables
113   G4double XTRNAngleSpectralDensity(G4double energy);
114   G4double XTRNAngleDensity(G4double varAngle);
115 
116   void GetNumberOfPhotons();
117 
118   // Auxiliary functions for plate/gas material parameters
119   G4double GetPlateFormationZone(G4double, G4double, G4double);
120   G4complex GetPlateComplexFZ(G4double, G4double, G4double);
121   void ComputePlatePhotoAbsCof();
122   G4double GetPlateLinearPhotoAbs(G4double);
123   void GetPlateZmuProduct();
124   G4double GetPlateZmuProduct(G4double, G4double, G4double);
125 
126   G4double GetGasFormationZone(G4double, G4double, G4double);
127   G4complex GetGasComplexFZ(G4double, G4double, G4double);
128   void ComputeGasPhotoAbsCof();
129   G4double GetGasLinearPhotoAbs(G4double);
130   void GetGasZmuProduct();
131   G4double GetGasZmuProduct(G4double, G4double, G4double);
132 
133   G4double GetPlateCompton(G4double);
134   G4double GetGasCompton(G4double);
135   G4double GetComptonPerAtom(G4double, G4double);
136 
137   G4double GetXTRrandomEnergy(G4double scaledTkin, G4int iTkin);
138   G4double GetXTRenergy(G4int iPlace, G4double position, G4int iTransfer);
139 
140   G4double GetRandomAngle(G4double energyXTR, G4int iTkin);
141   G4double GetAngleXTR(G4int iTR, G4double position, G4int iAngle);
142 
143   // set/get methods for class fields
144 
145   void     SetGamma(G4double gamma) { fGamma = gamma; };
146   G4double GetGamma() { return fGamma; };
147   void     SetEnergy(G4double energy) { fEnergy = energy; };
148   G4double GetEnergy() { return fEnergy; };
149   void     SetVarAngle(G4double varAngle) { fVarAngle = varAngle; };
150   G4double GetVarAngle() { return fVarAngle; };
151   void   SetCompton(G4bool pC) { fCompton = pC; };
152   G4bool GetCompton() { return fCompton; };
153 
154   G4int GetKrange(){ return fKrange;};
155   void SetKrange( G4int kk ){ fKrange = kk;};
156 
157 
158   void     SetAlphaGas(G4double ag){ fAlphaGas = ag;};
159   G4double GetAlphaGas() { return fAlphaGas; };  
160   void     SetAlphaPlate(G4double ap){ fAlphaPlate = ap;};
161   G4double GetAlphaPlate() { return fAlphaPlate; };
162 
163   void     SetTheMinEnergyTR(G4double minetr){ fTheMinEnergyTR = minetr;};
164   G4double GetTheMinEnergyTR() { return fTheMinEnergyTR; };  
165   void     SetTheMaxEnergyTR(G4double maxetr){ fTheMaxEnergyTR = maxetr;};
166   G4double GetTheMaxEnergyTR() { return fTheMaxEnergyTR; };  
167 
168   void     SetMinEnergyTR(G4double minetr){ fMinEnergyTR = minetr;};
169   G4double GetMinEnergyTR() { return fMinEnergyTR; };  
170   void     SetMaxEnergyTR(G4double maxetr){ fMaxEnergyTR = maxetr;};
171   G4double GetMaxEnergyTR() { return fMaxEnergyTR; };  
172   
173   void     SetTheMinAngle(G4double minang){ fTheMinAngle = minang;};
174   G4double GetTheMinAngle() { return fTheMinAngle; };  
175   void     SetTheMaxAngle(G4double maxang){ fTheMaxAngle = maxang;};
176   G4double GetTheMaxAngle() { return fTheMaxAngle; };
177 
178   void     SetMinThetaTR(G4double minatr){ fMinThetaTR = minatr;};
179   G4double GetMinThetaTR() { return fMinThetaTR; };  
180   void     SetMaxThetaTR(G4double maxatr){ fMaxThetaTR = maxatr;};
181   G4double GetMaxThetaTR() { return fMaxThetaTR; };
182 
183   // modes of XTR angle distribution
184   
185   void   SetFastAngle(G4bool fatr){ fFastAngle = fatr;};
186   G4bool GetFastAngle() { return fFastAngle; };  
187   void   SetAngleRadDistr(G4bool fatr){ fAngleRadDistr = fatr;};
188   G4bool GetAngleRadDistr() { return fAngleRadDistr; };  
189   
190 
191   
192 
193   G4PhysicsLogVector* GetProtonVector() { return fProtonEnergyVector; };
194   G4int GetTotBin() { return fTotBin; };
195   G4PhysicsFreeVector* GetAngleVector(G4double energy, G4int n);
196 
197  protected:
198   //   min TR energy
199   G4double fTheMinEnergyTR; 
200   //   max TR energy
201   G4double fTheMaxEnergyTR; 
202   G4double fTheMinAngle;  //  min theta of TR quanta
203   G4double fTheMaxAngle;  //  1.e-4;  //  max theta of TR quanta
204 
205   // static const members
206   
207   // min Tkin of proton in tables
208   static constexpr G4double fMinProtonTkin = 100. * CLHEP::GeV;
209   // max Tkin of proton in tables
210   static constexpr G4double fMaxProtonTkin = 100. * CLHEP::TeV;
211   // physical constants for plasma energy
212   static constexpr G4double fPlasmaCof =
213     4. * CLHEP::pi * CLHEP::fine_structure_const * CLHEP::hbarc * CLHEP::hbarc *
214     CLHEP::hbarc / CLHEP::electron_mass_c2;
215   static constexpr G4double fCofTR = CLHEP::fine_structure_const / CLHEP::pi;
216 
217   G4int fTotBin;  //  number of bins in log-gamma scale
218   G4int fBinTR;   //  number of bins in TR energy-angle vectors
219   G4int fKrange; 
220   G4ParticleDefinition* fPtrGamma;  // pointer to TR photon
221 
222   G4double* fGammaCutInKineticEnergy;  // TR photon cut in energy array
223   G4LogicalVolume* fEnvelope;
224   G4PhysicsTable* fAngleDistrTable;
225   G4PhysicsTable* fEnergyDistrTable;
226   G4PhysicsTable* fAngleForEnergyTable;
227   G4PhysicsLogVector* fProtonEnergyVector;
228   G4PhysicsLogVector* fXTREnergyVector;
229   G4SandiaTable* fPlatePhotoAbsCof;
230   G4SandiaTable* fGasPhotoAbsCof;
231 
232   G4ParticleChange fParticleChange;
233   std::vector<G4PhysicsTable*> fAngleBank;
234 
235   G4double fGammaTkinCut;  // Tkin cut of TR photon in current mat.
236   G4double fMinEnergyTR;   //  min TR energy in material
237   G4double fMaxEnergyTR;   //  max TR energy in material
238   G4double fMinThetaTR, fMaxThetaTR;    //  min-max theta of TR quanta
239   G4double fTotalDist;
240   G4double fPlateThick;
241   G4double fGasThick;
242   G4double fAlphaPlate;
243   G4double fAlphaGas;
244   G4double fGamma;     // current Lorentz factor
245   G4double fEnergy;    // energy and
246   G4double fVarAngle;  // angle squared!
247   G4double fLambda;
248   G4double fSigma1;
249   G4double fSigma2;  // plasma energy Sq of matter1/2
250 
251   G4int fMatIndex1;
252   G4int fMatIndex2;
253   G4int fPlateNumber;
254 
255   G4bool fExitFlux;
256   G4bool fFastAngle, fAngleRadDistr;
257   G4bool fCompton;
258 
259   G4int secID = -1;  // creator modelID
260 };
261 
262 #endif
263