Geant4
9.6.p02
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
geant4_9_6_p02
source
processes
hadronic
models
inclxx
incl_physics
src
G4INCLPionNucleonChannel.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
// INCL++ intra-nuclear cascade model
27
// Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
28
// Davide Mancusi, CEA
29
// Alain Boudard, CEA
30
// Sylvie Leray, CEA
31
// Joseph Cugnon, University of Liege
32
//
33
#define INCLXX_IN_GEANT4_MODE 1
34
35
#include "
globals.hh
"
36
37
#include "
G4INCLPionNucleonChannel.hh
"
38
#include "
G4INCLKinematicsUtils.hh
"
39
#include "
G4INCLBinaryCollisionAvatar.hh
"
40
#include "
G4INCLRandom.hh
"
41
#include "
G4INCLGlobals.hh
"
42
#include "
G4INCLLogger.hh
"
43
44
namespace
G4INCL {
45
46
PionNucleonChannel::PionNucleonChannel
(
Particle
*p1,
Particle
*p2,
Nucleus
*nucleus,
const
G4bool
localE)
47
: theNucleus(nucleus), particle1(p1), particle2(p2), locE(localE)
48
{
49
50
}
51
52
PionNucleonChannel::~PionNucleonChannel
(){
53
54
}
55
56
FinalState
*
PionNucleonChannel::getFinalState
() {
57
FinalState
*fs =
new
FinalState
;
58
59
Particle
* nucleon;
60
Particle
* pion;
61
if
(particle1->
isNucleon
()) {
62
nucleon = particle1;
63
pion = particle2;
64
}
else
{
65
nucleon = particle2;
66
pion = particle1;
67
}
68
69
ParticleType
deltaType =
DeltaZero
;
70
if
(
ParticleConfig::isPair
(particle1, particle2,
Proton
,
PiPlus
)) {
71
deltaType =
DeltaPlusPlus
;
72
}
else
if
(
ParticleConfig::isPair
(particle1, particle2,
Neutron
,
PiPlus
)) {
73
deltaType =
DeltaPlus
;
74
}
else
if
(
ParticleConfig::isPair
(particle1, particle2,
Proton
,
PiZero
)) {
75
deltaType =
DeltaPlus
;
76
}
else
if
(
ParticleConfig::isPair
(particle1, particle2,
Neutron
,
PiZero
)) {
77
deltaType =
DeltaZero
;
78
}
else
if
(
ParticleConfig::isPair
(particle1, particle2,
Proton
,
PiMinus
)) {
79
deltaType =
DeltaZero
;
80
}
else
if
(
ParticleConfig::isPair
(particle1, particle2,
Neutron
,
PiMinus
)) {
81
deltaType =
DeltaMinus
;
82
}
else
{
83
ERROR
(
"Unknown particle pair in Pi-N collision."
<< std::endl);
84
}
85
86
G4double
deltaEnergy = nucleon->
getEnergy
() - nucleon->
getPotentialEnergy
()
87
+ pion->
getEnergy
() - pion->
getPotentialEnergy
();
88
89
nucleon->
setType
(deltaType);
// nucleon becomes the delta
90
deltaEnergy += theNucleus->
getPotential
()->
computePotentialEnergy
(nucleon);
91
nucleon->
setEnergy
(deltaEnergy);
// set the energy of the delta
92
93
ThreeVector
deltaMomentum = nucleon->
getMomentum
() + pion->
getMomentum
();
94
nucleon->
setMomentum
(deltaMomentum);
95
96
const
G4double
deltaMass = std::sqrt(deltaEnergy*deltaEnergy - deltaMomentum.
mag2
());
97
nucleon->
setMass
(deltaMass);
98
theNucleus->
updatePotentialEnergy
(nucleon);
99
100
fs->
addModifiedParticle
(nucleon);
// nucleon became a delta
101
fs->
addDestroyedParticle
(pion);
// pion was removed
102
return
fs;
103
}
104
105
}
Generated on Sat May 25 2013 14:33:59 for Geant4 by
1.8.4