Geant4-11
G4DNARuddAngle.cc
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// -------------------------------------------------------------------
28//
29// GEANT4 Class file
30//
31//
32// File name: G4DNARuddAngle
33//
34// Author: Vladimir Ivantcheko
35//
36// Creation date: 23 August 2013
37//
38// Modifications:
39//
40// Class Description:
41//
42// Delta-electron Angular Distribution Generation
43//
44// Class Description: End
45//
46// -------------------------------------------------------------------
47//
48
49#include "G4DNARuddAngle.hh"
51#include "Randomize.hh"
53#include "G4Electron.hh"
54#include "G4SystemOfUnits.hh"
55
56using namespace std;
57
59 : G4VEmAngularDistribution("deltaRudd")
60{
62}
63
65{}
66
69 G4double secKinetic, G4int,
70 G4int,
71 const G4Material*)
72{
73 G4double k = dp->GetKineticEnergy();
74 G4double cosTheta = 1.0;
75
76 const G4ParticleDefinition* particle = dp->GetDefinition();
77 G4double mass = particle->GetPDGMass();
78
79 G4double maximumEnergyTransfer = k;
80 if(particle == fElectron) { maximumEnergyTransfer *= 0.5; }
81 else if(mass > MeV) {
82 G4double ratio = electron_mass_c2/mass;
83 G4double tau = k/mass;
84 maximumEnergyTransfer = 2.0*electron_mass_c2*tau*(tau + 2.) /
85 (1. + 2.0*(tau + 1.)*ratio + ratio*ratio);
86
87 }
88
89 if (secKinetic>100*eV && secKinetic <= maximumEnergyTransfer) {
90 cosTheta = std::sqrt(secKinetic / maximumEnergyTransfer);
91 } else {
92 cosTheta = (2.*G4UniformRand())-1.;
93 }
94
95 G4double sint = sqrt((1.0 - cosTheta)*(1.0 + cosTheta));
97
98 fLocalDirection.set(sint*cos(phi), sint*sin(phi), cosTheta);
100
101 return fLocalDirection;
102}
103
106 G4double secEkin, G4int Z,
107 const G4Material* mat)
108{
109 return SampleDirectionForShell(dp, secEkin, Z, 0, mat);
110}
111
113{}
static constexpr double twopi
Definition: G4SIunits.hh:56
static constexpr double eV
Definition: G4SIunits.hh:201
static constexpr double MeV
Definition: G4SIunits.hh:200
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
const G4int Z[17]
#define G4UniformRand()
Definition: Randomize.hh:52
void set(double x, double y, double z)
Hep3Vector & rotateUz(const Hep3Vector &)
Definition: ThreeVector.cc:33
~G4DNARuddAngle() override
G4ThreeVector & SampleDirectionForShell(const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, G4int shellIdx, const G4Material *mat=nullptr) override
void PrintGeneratorInformation() const override
G4DNARuddAngle(const G4String &name="")
const G4ParticleDefinition * fElectron
G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, const G4Material *mat=nullptr) override
const G4ThreeVector & GetMomentumDirection() const
G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
static G4Electron * Electron()
Definition: G4Electron.cc:93
float electron_mass_c2
Definition: hepunit.py:273