1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
|
#include <globals.hh>
#include <G4ParticleDefinition.hh>
#include <G4ParticleTable.hh>
#include <G4Gamma.hh>
#include <G4Electron.hh>
#include <G4Positron.hh>
#include <G4NeutrinoE.hh>
#include <G4AntiNeutrinoE.hh>
#include <G4ProcessManager.hh>
#include "G4BeamTestEMPhysics.h"
G4BeamTestEMPhysics::G4BeamTestEMPhysics()
: G4VPhysicsConstructor("standard")
{}
G4BeamTestEMPhysics::~G4BeamTestEMPhysics()
{}
void G4BeamTestEMPhysics::ConstructParticle()
{
// gamma
G4Gamma::GammaDefinition();
// electron
G4Electron::ElectronDefinition();
G4Positron::PositronDefinition();
G4NeutrinoE::NeutrinoEDefinition();
G4AntiNeutrinoE::AntiNeutrinoEDefinition();
}
void G4BeamTestEMPhysics::ConstructProcess()
{
G4ProcessManager *pManager = 0;
// Gamma Physics
pManager = G4Gamma::Gamma()->GetProcessManager();
pManager->AddDiscreteProcess(&photoEffect);
pManager->AddDiscreteProcess(&comptonEffect);
pManager->AddDiscreteProcess(&pairProduction);
// Electron Physics
pManager = G4Electron::Electron()->GetProcessManager();
pManager->AddProcess(&electronMultipleScattering, -1, 1, 1);
pManager->AddProcess(&electronIonisation, -1, 2, 2);
pManager->AddProcess(&electronBremsStrahlung, -1, 3, 3);
//Positron Physics
pManager = G4Positron::Positron()->GetProcessManager();
pManager->AddProcess(&positronMultipleScattering, -1, 1, 1);
pManager->AddProcess(&positronIonisation, -1, 2, 2);
pManager->AddProcess(&positronBremsStrahlung, -1, 3, 3);
pManager->AddProcess(&annihilation, 0,-1, 4);
}
|
