Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4ScreeningMottCrossSection.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 //  G4ScreeningMottCrossSection.hh
 27 //-------------------------------------------------------------------
 28 //
 29 // GEANT4 Class header file
 30 //
 31 // File name:    G4ScreeningMottCrossSection
 32 //
 33 // Author:      Cristina Consolandi
 34 //
 35 // Creation date: 20.10.2011
 36 //
 37 // Modifications:
 38 //
 39 // Class Description:
 40 //      Computation of electron Coulomb Scattering Cross Section.
 41 //      Suitable for high energy electrons and light target materials.
 42 //
 43 //      Reference:
 44 //      M.J. Boschini et al.
 45 //     "Non Ionizing Energy Loss induced by Electrons in the Space Environment"
 46 //      Proc. of the 13th Int. Conf. on Particle Physics and Advanced Technology
 47 //      (13th ICPPAT, Como 3-7/10/2011), World Scientific (Singapore).
 48 //      Available at: http://arxiv.org/abs/1111.4042v4
 49 //
 50 //      1) Mott Differential Cross Section Approximation:
 51 //         For Target material up to Z=92 (U):
 52 //         As described in http://arxiv.org/abs/1111.4042v4
 53 //         par. 2.1 , eq. (16)-(17)
 54 //         Else (Z>92):
 55 //         W. A. McKinley and H. Fashbach, Phys. Rev. 74, (1948) 1759.
 56 //      2) Screening coefficient:
 57 //      vomn G. Moliere, Z. Naturforsh A2 (1947), 133-145; A3 (1948), 78.
 58 //      3) Nuclear Form Factor:
 59 //      A.V. Butkevich et al. Nucl. Instr. Meth. A488 (2002), 282-294.
 60 //
 61 // -----------------------------------------------------------------------------
 62 
 63 //
 64 #ifndef G4ScreeningMottCrossSection_h
 65 #define G4ScreeningMottCrossSection_h 1
 66 
 67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 68 
 69 #include "globals.hh"
 70 #include "G4ParticleDefinition.hh"
 71 #include <vector>
 72 
 73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 74 
 75 static const G4int DIMMOTT = 750;
 76 
 77 class G4NistManager;
 78 class G4Pow;
 79 
 80 class G4ScreeningMottCrossSection
 81 {
 82 
 83 public:
 84 
 85   explicit G4ScreeningMottCrossSection();
 86 
 87   ~G4ScreeningMottCrossSection();
 88 
 89   void Initialise(const G4ParticleDefinition*, G4double cosThetaLim);
 90 
 91   void SetupKinematic(G4double kinEnergy, G4int Z);
 92 
 93   G4double NuclearCrossSection(G4int form, G4int fast);
 94   G4double GetScatteringAngle(G4int form, G4int fast);
 95 
 96   G4double RatioMottRutherford(G4double tet);
 97   G4double RatioMottRutherfordCosT(G4double sin2t2);
 98 
 99   G4double McFcorrection(G4double sin2t2);
100   inline void SetupParticle(const G4ParticleDefinition*);
101 
102   G4ScreeningMottCrossSection & operator=
103   (const G4ScreeningMottCrossSection &right) = delete;
104   G4ScreeningMottCrossSection(const G4ScreeningMottCrossSection&) = delete;
105 
106 private:
107 
108   G4double ComputeAngle(G4int idx, G4double& rand);
109 
110   G4double FormFactor2ExpHof(G4double sin2t2);
111   G4double FormFactor2Gauss(G4double sin2t2);
112   G4double FormFactor2UniformHelm(G4double sin2t2);
113   G4double DifferentialXSection(G4int idx, G4int form);
114 
115   G4double  GetTransitionRandom();
116 
117   G4NistManager*  fNistManager;
118   G4Pow*          fG4pow;
119 
120   const G4ParticleDefinition* particle;
121 
122   G4double          fTotalCross;
123   //cost - min - max
124   G4double              cosThetaMin;// def 1.0
125   G4double              cosThetaMax;// def -1.0
126 
127   G4double        cosTetMinNuc;
128   G4double          cosTetMaxNuc;
129 
130   //energy cut
131   G4double              ecut;
132   G4double              etag;
133 
134   G4double              spin;
135   G4double              mass;
136 
137   //lab of incedent particle
138   G4double              tkinLab;
139   G4double              momLab2;
140   G4double              invbetaLab2;
141 
142   //relative system with nucleus
143   G4double    mu_rel;
144   G4double              tkin;
145   G4double              mom2;
146   G4double              invbeta2;
147   G4double    beta;
148   G4double    gamma;
149 
150   //constants
151   G4double              alpha;
152   G4double              htc2;
153   G4double              e2;
154 
155   // target nucleus
156   G4double              targetMass;
157   G4double    As;
158   G4int                 targetZ;
159   G4int           targetA;
160 
161   // working array
162   std::vector<G4double> cross;
163 };
164 
165 
166 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
167 
168 inline 
169 void G4ScreeningMottCrossSection::SetupParticle(const G4ParticleDefinition* p)
170 {
171   particle = p;
172   mass = particle->GetPDGMass();
173   spin = particle->GetPDGSpin();
174   if(0.0 != spin) { spin = 0.5; }
175   tkin = 0.0;
176 }
177 
178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
179 
180 #endif
181