source: trunk/examples/advanced/hadrontherapy/src/EMPositronStandard.cc@ 1305

Last change on this file since 1305 was 807, checked in by garnier, 17 years ago

update

  • Property svn:executable set to *
File size: 4.1 KB
RevLine 
[807]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// HadrontherapyPositronStandard.cc May 2005
27// ----------------------------------------------------------------------------
28// GEANT 4 - Hadrontherapy example
29// ----------------------------------------------------------------------------
30// Code developed by:
31//
32// G.A.P. Cirrone(a)*, F. Di Rosa(a), S. Guatelli(b), G. Russo(a)
33//
34// (a) Laboratori Nazionali del Sud
35// of the National Institute for Nuclear Physics, Catania, Italy
36// (b) National Institute for Nuclear Physics Section of Genova, genova, Italy
37//
38// * cirrone@lns.infn.it
39// ----------------------------------------------------------------------------
40
41#include "EMPositronStandard.hh"
42#include "G4ProcessManager.hh"
43#include "G4ParticleDefinition.hh"
44#include "G4MultipleScattering.hh"
45#include "G4eIonisation.hh"
46#include "G4eBremsstrahlung.hh"
47#include "G4eplusAnnihilation.hh"
48#include "G4StepLimiter.hh"
49#include "G4EmProcessOptions.hh"
50
51
52EMPositronStandard::EMPositronStandard(const G4String& name):
53 G4VPhysicsConstructor(name)
54{
55 G4cout<< "ELECTROMAGNETIC PROCESS(ES): G4MultipleScattering (positron)"
56 << G4endl
57 << " G4eIonisation (positron)"
58 << G4endl
59 << " G4eBremsstrahlung (positron)"
60 << G4endl
61 << " G4eplusAnnihilation (positron)"
62 << G4endl
63 << "APPLIED MODEL(S): -"
64 << G4endl;
65}
66
67EMPositronStandard::~EMPositronStandard()
68{ }
69
70void EMPositronStandard::ConstructProcess()
71{
72
73 // ****************
74 // *** Positron ***
75 // ****************
76
77 G4EmProcessOptions* positronEmProcessOptions = new G4EmProcessOptions();
78 positronEmProcessOptions -> SetDEDXBinning(480);
79
80 G4MultipleScattering* positronMultipScatProcess = new G4MultipleScattering();
81 G4eIonisation* positronIonisationProcess = new G4eIonisation();
82 G4eBremsstrahlung* positronBremsstrProcess = new G4eBremsstrahlung();
83 G4eplusAnnihilation* positronAnnihilationProcess = new G4eplusAnnihilation();
84
85 G4StepLimiter* positronStepLimiter = new G4StepLimiter();
86
87 G4ParticleDefinition* particle = G4Positron::Positron();
88 G4ProcessManager* processManager = particle -> GetProcessManager();
89 processManager -> AddProcess(positronMultipScatProcess, -1, 1, 1);
90 processManager -> AddProcess(positronIonisationProcess, -1, 2, 2);
91 processManager -> AddProcess(positronBremsstrProcess, -1, -1, 3);
92 processManager -> AddProcess(positronAnnihilationProcess, 0, -1, 4);
93 processManager -> AddProcess(positronStepLimiter, -1, -1, 3);
94
95}
Note: See TracBrowser for help on using the repository browser.