Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/xrays/include/G4ForwardXrayTR.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 ]

Diff markup

Differences between /processes/electromagnetic/xrays/include/G4ForwardXrayTR.hh (Version 11.3.0) and /processes/electromagnetic/xrays/include/G4ForwardXrayTR.hh (Version 1.0)


                                                   >>   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 //                                                  7 //
 26 // G4ForwardXrayTR                             <<   8 // $Id: G4ForwardXrayTR.hh,v 1.2.8.1 1999/12/07 20:51:24 gunter Exp $
                                                   >>   9 // GEANT4 tag $Name: geant4-01-00 $
 27 //                                                 10 //
 28 // Class for description                       <<  11 // G4ForwardXrayTR -- header file
 29 //                                                 12 //
 30 // Class for forward X-ray transition radiatio <<  13 // Class for description of forward X-ray transition radiation generated
 31 // by relativistic charged particle crossed in     14 // by relativistic charged particle crossed interface between material 1
 32 // and material 2 (1 -> 2)                         15 // and material 2 (1 -> 2)
 33                                                <<  16 //
                                                   >>  17 // GEANT 4 class header file --- Copyright CERN 1995
                                                   >>  18 // CERB Geneva Switzerland
                                                   >>  19 //
                                                   >>  20 // for information related to this code, please, contact
                                                   >>  21 // CERN, CN Division, ASD Group
 34 // History:                                        22 // History:
 35 // 22.09.97, V. Grichine (Vladimir.Grichine@ce     23 // 22.09.97, V. Grichine (Vladimir.Grichine@cern.ch)
 36 // 26.01.00, V.Grichine, new constructor and p <<  24 
 37 // 10.03.03, V.Ivanchenko migrade to "cut per  << 
 38 // 03.06.03, V.Ivanchenko fix compilation warn << 
 39                                                    25 
 40 #ifndef G4FORWARDXRAYTR_H                          26 #ifndef G4FORWARDXRAYTR_H
 41 #define G4FORWARDXRAYTR_H                          27 #define G4FORWARDXRAYTR_H
 42                                                    28 
                                                   >>  29 
 43 #include "globals.hh"                              30 #include "globals.hh"
 44 #include "G4Track.hh"                          <<  31 #include "templates.hh"
 45 #include "G4TransitionRadiation.hh"            <<  32 #include "geomdefs.hh"
 46 #include "G4VParticleChange.hh"                <<  33 #include "Randomize.hh"
                                                   >>  34 #include "G4Step.hh"
                                                   >>  35 #include "G4VDiscreteProcess.hh"
                                                   >>  36 #include "G4DynamicParticle.hh"
                                                   >>  37 #include "G4Material.hh"
                                                   >>  38 #include "G4LogicalBorderSurface.hh"
                                                   >>  39 #include "G4LogicalSkinSurface.hh"
                                                   >>  40 #include "G4OpticalSurface.hh"
                                                   >>  41 #include "G4OpticalPhoton.hh"
                                                   >>  42 #include "G4TransportationManager.hh"
 47                                                    43 
 48 class G4ParticleDefinition;                    <<  44 #include "G4TransitionRadiation.hh"
 49 class G4PhysicsTable;                          <<  45 #include "G4PhysicsTable.hh"
 50 class G4PhysicsLogVector;                      <<  46 #include "G4Gamma.hh"
 51                                                    47 
 52 class G4ForwardXrayTR : public G4TransitionRad     48 class G4ForwardXrayTR : public G4TransitionRadiation
 53 {                                                  49 {
 54  public:                                       <<  50 public:
 55   explicit G4ForwardXrayTR(const G4String& mat <<  51 
 56                            const G4String& pro <<  52 // Constructors
                                                   >>  53 
                                                   >>  54   //     G4ForwardXrayTR() ;
                                                   >>  55 
                                                   >>  56 
                                                   >>  57   G4ForwardXrayTR(  const G4String& matName1,    //  G4Material* pMat1,
                                                   >>  58         const G4String& matName2,   //  G4Material* pMat2,     
                                                   >>  59                     const G4String& processName="XrayTR"     ) ;
                                                   >>  60 
                                                   >>  61 //  G4ForwardXrayTR(const G4ForwardXrayTR& right) ;
                                                   >>  62 
                                                   >>  63 // Destructor
                                                   >>  64    ~G4ForwardXrayTR() ;
                                                   >>  65 
                                                   >>  66 // Operators
                                                   >>  67 // G4ForwardXrayTR& operator=(const G4ForwardXrayTR& right) ;
                                                   >>  68 // G4int operator==(const G4ForwardXrayTR& right)const ;
                                                   >>  69 // G4int operator!=(const G4ForwardXrayTR& right)const ;
                                                   >>  70 
                                                   >>  71 ///////////////////////    Methods    /////////////////////////////////
                                                   >>  72 
                                                   >>  73         void BuildXrayTRtables();
                                                   >>  74 
                                                   >>  75   G4double GetMeanFreePath(const G4Track& aTrack,
                                                   >>  76          G4double previousStepSize,
                                                   >>  77          G4ForceCondition* condition)
                                                   >>  78         {
                                                   >>  79           *condition = Forced;
                                                   >>  80     return DBL_MAX;      // so TR doesn't limit mean free path
                                                   >>  81         }
                                                   >>  82 
                                                   >>  83   G4VParticleChange* PostStepDoIt( const G4Track& aTrack,
                                                   >>  84                  const G4Step&  aStep    )   ;
                                                   >>  85 
                                                   >>  86         G4double GetEnergyTR(G4int iMat, G4int jMat, G4int iTkin) const ;
                                                   >>  87        
                                                   >>  88         G4double GetThetaTR(G4int iMat, G4int jMat, G4int iTkin) const ;     
                                                   >>  89 
                                                   >>  90 
                                                   >>  91 ///////////////////// Angle distribution  /////////////////////////////
                                                   >>  92 //
                                                   >>  93 
                                                   >>  94 G4double SpectralAngleTRdensity( G4double energy,
                                                   >>  95                                  G4double varAngle ) const;
                                                   >>  96 
                                                   >>  97 
                                                   >>  98 G4double AngleDensity( G4double energy,
                                                   >>  99                        G4double varAngle ) const  ;
 57                                                   100 
 58   explicit G4ForwardXrayTR(const G4String& pro << 101 G4double EnergyInterval( G4double energy1,
                                                   >> 102                          G4double energy2,
                                                   >> 103                          G4double varAngle ) const ;
 59                                                   104 
 60   ~G4ForwardXrayTR();                          << 105 G4double AngleSum( G4double varAngle1,
                                                   >> 106                    G4double varAngle2     )   const ;
 61                                                   107 
 62   G4ForwardXrayTR(const G4ForwardXrayTR& right << 108 /////////////////////////  Energy distribution ///////////////////////////////
 63   G4ForwardXrayTR& operator=(const G4ForwardXr << 
 64                                                   109 
 65   ///////////////////////    Methods    ////// << 110 G4double SpectralDensity( G4double energy,
                                                   >> 111                           G4double x       ) const ;
 66                                                   112 
 67   void ProcessDescription(std::ostream&) const << 113 G4double AngleInterval( G4double energy,
 68   void DumpInfo() const override { ProcessDesc << 114                         G4double varAngle1,
                                                   >> 115                         G4double varAngle2   ) const ;
 69                                                   116 
 70   void BuildXrayTRtables();                    << 117 G4double EnergySum( G4double energy1,
                                                   >> 118                     G4double energy2     )   const ;
 71                                                   119 
 72   G4double GetMeanFreePath(const G4Track&, G4d << 
 73                            G4ForceCondition* c << 
 74                                                   120 
 75   G4VParticleChange* PostStepDoIt(const G4Trac << 121 ///////////////////////////   Access functions  ////////////////////////////
 76                                   const G4Step << 
 77                                                   122 
 78   G4double GetEnergyTR(G4int iMat, G4int jMat, << 123   G4PhysicsTable* GetAngleDistrTable()  { return fAngleDistrTable  ; } ;
                                                   >> 124   G4PhysicsTable* GetEnergyDistrTable() { return fEnergyDistrTable ; } ;
 79                                                   125 
 80   G4double GetThetaTR(G4int iMat, G4int jMat,  << 126   static G4int    GetSympsonNumber()       { return fSympsonNumber ;  } ;
                                                   >> 127   static G4int            GetBinTR()       { return fBinTR ;       } ;
 81                                                   128 
 82   ///////////////////// Angle distribution  // << 129   static G4double GetMinProtonTkin()       { return fMinProtonTkin ; } ;
                                                   >> 130   static G4double GetMaxProtonTkin()       { return fMaxProtonTkin ; } ;
                                                   >> 131   static G4int           GetTotBin()       { return fTotBin ;       } ;
 83                                                   132 
 84   G4double SpectralAngleTRdensity(G4double ene << 
 85                                   G4double var << 
 86                                                   133 
 87   G4double AngleDensity(G4double energy, G4dou << 134 protected :
 88                                                   135 
 89   G4double EnergyInterval(G4double energy1, G4 << 136 private :  ///////////////  Data members   ///////////////////////////
 90                           G4double varAngle) c << 
 91                                                   137 
 92   G4double AngleSum(G4double varAngle1, G4doub << 138 G4Gamma* fPtrGamma ;  // pointer to TR photon
 93                                                   139 
 94   /////////////////////////  Energy distributi << 140 G4double* fGammaCutInKineticEnergy ; // TR photon cut in energy array
                                                   >> 141 G4double  fGammaTkinCut ;            // Tkin cut of TR photon in current mat.
 95                                                   142 
 96   G4double SpectralDensity(G4double energy, G4 << 143 G4PhysicsTable* fAngleDistrTable ;
                                                   >> 144 G4PhysicsTable* fEnergyDistrTable ;
 97                                                   145 
 98   G4double AngleInterval(G4double energy, G4do << 146 static G4PhysicsLogVector* fProtonEnergyVector ;
 99                          G4double varAngle2) c << 
100                                                   147 
101   G4double EnergySum(G4double energy1, G4doubl << 148 static G4int fSympsonNumber ;                // Accuracy of Sympson integration 
102                                                   149 
103   ///////////////////////////   Access functio << 150 static G4double fTheMinEnergyTR ;            //  static min TR energy
                                                   >> 151 static G4double fTheMaxEnergyTR ;            //  static max TR energy
                                                   >> 152        G4double fMinEnergyTR ;               //  min TR energy in material
                                                   >> 153        G4double fMaxEnergyTR ;               //  max TR energy in material
                                                   >> 154 static G4double  fTheMaxAngle ;               //  max theta of TR quanta
                                                   >> 155 static G4double  fTheMinAngle ;               //  max theta of TR quanta
                                                   >> 156        G4double  fMaxThetaTR ;               //  max theta of TR quanta
                                                   >> 157 static G4int          fBinTR ;               //  number of bins in TR vectors
104                                                   158 
105   G4PhysicsTable* GetAngleDistrTable();        << 159 static G4double fMinProtonTkin ;             // min Tkin of proton in tables
106   G4PhysicsTable* GetEnergyDistrTable();       << 160 static G4double fMaxProtonTkin ;             // max Tkin of proton in tables
                                                   >> 161 static G4int    fTotBin        ;             // number of bins in log scale
                                                   >> 162        G4double fGamma         ;             // current Lorentz factor
107                                                   163 
108   static G4int GetSympsonNumber();             << 164 static G4double fPlasmaCof ;               // physical consts for plasma energy
109   static G4int GetBinTR();                     << 165 static G4double fCofTR ;
110                                                   166 
111   static G4double GetMinProtonTkin();          << 167 G4double fSigma1 ;                       // plasma energy Sq of matter1
112   static G4double GetMaxProtonTkin();          << 168 G4double fSigma2 ;                       // plasma energy Sq of matter2
113   static G4int GetTotBin();                    << 
114                                                   169 
115  protected:  // for access from X-ray TR fast  << 
116   static constexpr G4double fTheMinEnergyTR =  << 
117     1. * CLHEP::keV;  //  static min TR energy << 
118   static constexpr G4double fTheMaxEnergyTR =  << 
119     100. * CLHEP::keV;                         << 
120   static constexpr G4double fTheMaxAngle = 1.0 << 
121   static constexpr G4double fTheMinAngle = 5.0 << 
122   static constexpr G4double fMinProtonTkin =   << 
123     100. * CLHEP::GeV;  // min Tkin of proton  << 
124   static constexpr G4double fMaxProtonTkin =   << 
125     100. * CLHEP::TeV;  // max Tkin of proton  << 
126   static constexpr G4double fPlasmaCof =       << 
127     4.0 * CLHEP::pi * CLHEP::fine_structure_co << 
128     CLHEP::hbarc * CLHEP::hbarc /              << 
129     CLHEP::electron_mass_c2;  // physical cons << 
130   static constexpr G4double fCofTR = CLHEP::fi << 
131                                                   170 
132   static constexpr G4int fSympsonNumber =      << 171 } ;    // end of G4ForwardXrayTR class ---------------------------
133     100;                                // Acc << 
134   static constexpr G4int fBinTR  = 50;  //  nu << 
135   static constexpr G4int fTotBin = 50;  // num << 
136                                                   172 
137   const std::vector<G4double>* fGammaCutInKine << 
138   // TR photon cut in energy array             << 
139                                                   173 
140   G4ParticleDefinition* fPtrGamma;  // pointer << 
141                                                   174 
142   G4PhysicsTable* fAngleDistrTable;            << 
143   G4PhysicsTable* fEnergyDistrTable;           << 
144                                                   175 
145   G4PhysicsLogVector* fProtonEnergyVector;     << 176 #endif   // G4FORWARDXRAYTR_H
146                                                   177 
147   G4double fMinEnergyTR;   //  min TR energy i << 
148   G4double fMaxEnergyTR;   //  max TR energy i << 
149   G4double fMaxThetaTR;    //  max theta of TR << 
150   G4double fGamma;         // current Lorentz  << 
151   G4double fGammaTkinCut;  // Tkin cut of TR p << 
152   G4double fSigma1;        // plasma energy Sq << 
153   G4double fSigma2;        // plasma energy Sq << 
154                                                   178 
155   G4int secID = -1;  // creator modelID        << 
156 };                                             << 
157                                                   179 
158 #endif  // G4FORWARDXRAYTR_H                   << 
159                                                   180