Geant4
10.03.p01
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
G4INCLOmegaNElasticChannel.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
// Alain Boudard, CEA-Saclay, France
28
// Joseph Cugnon, University of Liege, Belgium
29
// Jean-Christophe David, CEA-Saclay, France
30
// Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
31
// Sylvie Leray, CEA-Saclay, France
32
// Davide Mancusi, CEA-Saclay, France
33
//
34
#define INCLXX_IN_GEANT4_MODE 1
35
36
#include "
globals.hh
"
37
38
#include "
G4INCLOmegaNElasticChannel.hh
"
39
#include "
G4INCLKinematicsUtils.hh
"
40
#include "
G4INCLBinaryCollisionAvatar.hh
"
41
#include "
G4INCLRandom.hh
"
42
#include "
G4INCLGlobals.hh
"
43
#include "
G4INCLLogger.hh
"
44
45
namespace
G4INCL {
46
47
OmegaNElasticChannel::OmegaNElasticChannel
(
Particle
*p1,
Particle
*p2)
48
: particle1(p1), particle2(p2)
49
{
50
51
}
52
53
OmegaNElasticChannel::~OmegaNElasticChannel
(){
54
55
}
56
57
void
OmegaNElasticChannel::fillFinalState
(
FinalState
*fs) {
58
Particle
*
nucleon
;
59
Particle
* omega;
60
if
(particle1->
isNucleon
()) {
61
nucleon = particle1;
62
omega = particle2;
63
}
else
{
64
nucleon = particle2;
65
omega = particle1;
66
}
67
68
G4double
sh=nucleon->
getEnergy
()+omega->
getEnergy
();
69
G4double
mn=nucleon->
getMass
();
70
G4double
me=omega->
getMass
();
71
G4double
en=(sh*sh+mn*mn-me*me)/(2*sh);
72
nucleon->
setEnergy
(en);
73
G4double
ee=std::sqrt(en*en-mn*mn+me*me);
74
omega->
setEnergy
(ee);
75
G4double
pn
=std::sqrt(en*en-mn*mn);
76
77
ThreeVector
mom_nucleon;
78
79
// Isotropy
80
mom_nucleon =
Random::normVector
(pn);
81
82
83
nucleon->
setMomentum
(mom_nucleon);
84
omega->
setMomentum
(-mom_nucleon);
85
86
fs->
addModifiedParticle
(nucleon);
87
fs->
addModifiedParticle
(omega);
88
89
}
90
}
G4INCLOmegaNElasticChannel.hh
G4INCL::Particle::getMass
G4double getMass() const
Get the cached particle mass.
Definition:
G4INCLParticle.hh:348
G4INCLRandom.hh
G4INCL::Particle::getEnergy
G4double getEnergy() const
Definition:
G4INCLParticle.hh:554
G4INCLGlobals.hh
G4InuclParticleNames::pn
Definition:
G4InuclParticleNames.hh:67
G4INCL::Random::normVector
ThreeVector normVector(G4double norm=1.)
Definition:
G4INCLRandom.cc:134
G4INCLBinaryCollisionAvatar.hh
G4INCL::Particle::setEnergy
void setEnergy(G4double energy)
Definition:
G4INCLParticle.hh:570
G4INCL::OmegaNElasticChannel::fillFinalState
void fillFinalState(FinalState *fs)
Definition:
G4INCLOmegaNElasticChannel.cc:57
G4INCL::OmegaNElasticChannel::OmegaNElasticChannel
OmegaNElasticChannel(Particle *, Particle *)
Definition:
G4INCLOmegaNElasticChannel.cc:47
G4INCL::FinalState
Definition:
G4INCLFinalState.hh:64
G4InuclParticleNames::nucleon
G4bool nucleon(G4int ityp)
Definition:
G4InuclParticleNames.hh:95
globals.hh
G4INCL::Particle
Definition:
G4INCLParticle.hh:73
G4INCL::ThreeVector
Definition:
G4INCLThreeVector.hh:54
G4INCL::Particle::isNucleon
G4bool isNucleon() const
Definition:
G4INCLParticle.hh:232
G4double
double G4double
Definition:
G4Types.hh:76
G4INCLKinematicsUtils.hh
G4INCL::OmegaNElasticChannel::~OmegaNElasticChannel
virtual ~OmegaNElasticChannel()
Definition:
G4INCLOmegaNElasticChannel.cc:53
G4INCL::FinalState::addModifiedParticle
void addModifiedParticle(Particle *p)
Definition:
G4INCLFinalState.cc:60
G4INCLLogger.hh
G4INCL::Particle::setMomentum
virtual void setMomentum(const G4INCL::ThreeVector &momentum)
Definition:
G4INCLParticle.hh:592
geant4.10.03.p01
source
processes
hadronic
models
inclxx
incl_physics
src
G4INCLOmegaNElasticChannel.cc
Generated on Thu Mar 16 2017 22:38:10 for Geant4 by
1.8.5