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 6.2.p1)


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