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
59
60
61
62
63
64
65
66
67
68
69
70
71
|
/*
* copyright (C) 2010
* The Icecube Collaboration
*
* $Id: G4BeamTestEMPhysics.cxx 86420 2012-03-20 16:00:37Z kislat $
*
* @version $Revision: 86420 $
* @date $LastChangedDate: 2012-03-20 16:00:37 +0000 (Tue, 20 Mar 2012) $
* @author Fabian Kislat <fabian.kislat@desy.de>, Peter Nießen (tanktop)
* Last changed by: $LastChangedBy: kislat $
*/
#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);
}
|
