Geant4-11
Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
G4INCL::NNToNNKKbChannel Class Reference

#include <G4INCLNNToNNKKbChannel.hh>

Inheritance diagram for G4INCL::NNToNNKKbChannel:
G4INCL::IChannel

Public Member Functions

void fillFinalState (FinalState *fs)
 
FinalState * getFinalState ()
 
 NNToNNKKbChannel (Particle *, Particle *)
 
virtual ~NNToNNKKbChannel ()
 

Private Member Functions

 INCL_DECLARE_ALLOCATION_POOL (NNToNNKKbChannel)
 

Private Attributes

Particle * particle1
 
Particle * particle2
 

Static Private Attributes

static const G4double angularSlope = 2.
 

Detailed Description

Definition at line 47 of file G4INCLNNToNNKKbChannel.hh.

Constructor & Destructor Documentation

◆ NNToNNKKbChannel()

G4INCL::NNToNNKKbChannel::NNToNNKKbChannel ( Particle *  p1,
Particle *  p2 
)

Definition at line 51 of file G4INCLNNToNNKKbChannel.cc.

52 : particle1(p1), particle2(p2)
53 {}

◆ ~NNToNNKKbChannel()

G4INCL::NNToNNKKbChannel::~NNToNNKKbChannel ( )
virtual

Definition at line 55 of file G4INCLNNToNNKKbChannel.cc.

55{}

Member Function Documentation

◆ fillFinalState()

void G4INCL::NNToNNKKbChannel::fillFinalState ( FinalState *  fs)
virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLNNToNNKKbChannel.cc.

57 {
58
59 // pp -> pp K+ K- (1)
60 // pp -> pp K0 K0 (1)
61 // pp -> pn K+ K0 (4)
62 // pn -> pp K0 K- (4)
63 // pn -> pn K+ K- (9)
64
65 const G4double sqrtS = KinematicsUtils::totalEnergyInCM(particle1, particle2);
66
67 const G4int iso = ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
68
69 ParticleType Kaon1Type;
70 ParticleType Kaon2Type;
71 G4double rdm = Random::shoot();
72
73 if(iso == 2){
74 if(rdm*6. < 1.){
75 Kaon1Type = KPlus;
76 Kaon2Type = KMinus;
77 }
78 else if(rdm*6. < 2){
79 Kaon1Type = KZero;
80 Kaon2Type = KZeroBar;
81 }
82 else{
83 Kaon1Type = KPlus;
84 Kaon2Type = KZeroBar;
85 particle1->setType(Neutron);
86 }
87 }
88 else if(iso == -2){
89 if(rdm*6. < 1.){
90 Kaon1Type = KPlus;
91 Kaon2Type = KMinus;
92 }
93 else if(rdm*6. < 2){
94 Kaon1Type = KZero;
95 Kaon2Type = KZeroBar;
96 }
97 else{
98 Kaon1Type = KZero;
99 Kaon2Type = KMinus;
100 particle1->setType(Proton);
101 }
102 }
103 else if(rdm*26. < 9.){
104 Kaon1Type = KPlus;
105 Kaon2Type = KMinus;
106 }
107 else if(rdm*26. < 18.){
108 Kaon1Type = KZero;
109 Kaon2Type = KZeroBar;
110 }
111 else if(rdm*26. < 22.){
112 Kaon1Type = KZero;
113 Kaon2Type = KMinus;
114 particle1->setType(Proton);
115 particle2->setType(Proton);
116 }
117 else{
118 Kaon1Type = KPlus;
119 Kaon2Type = KZeroBar;
120 particle1->setType(Neutron);
121 particle2->setType(Neutron);
122 }
123
124 ParticleList list;
125 list.push_back(particle1);
126 list.push_back(particle2);
127 const ThreeVector &rcol1 = particle1->getPosition();
128 const ThreeVector &rcol2 = particle2->getPosition();
129 const ThreeVector zero;
130 Particle *kaon = new Particle(Kaon1Type,zero,rcol1);
131 Particle *antikaon = new Particle(Kaon2Type,zero,rcol2);
132 list.push_back(kaon);
133 list.push_back(antikaon);
134
135 if(Random::shoot()<0.5) PhaseSpaceGenerator::generateBiased(sqrtS, list, 0, angularSlope);
136 else PhaseSpaceGenerator::generateBiased(sqrtS, list, 1, angularSlope);
137
138 fs->addModifiedParticle(particle1);
139 fs->addModifiedParticle(particle2);
140 fs->addCreatedParticle(kaon);
141 fs->addCreatedParticle(antikaon);
142
143 }
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4INCL::NNToNNKKbChannel::angularSlope
static const G4double angularSlope
Definition: G4INCLNNToNNKKbChannel.hh:57
G4INCL::Particle::getPosition
const G4INCL::ThreeVector & getPosition() const
Definition: G4INCLParticle.hh:816
G4INCL::Particle::getType
G4INCL::ParticleType getType() const
Definition: G4INCLParticle.hh:178
G4INCL::Particle::setType
void setType(ParticleType t)
Definition: G4INCLParticle.hh:187
G4INCL::KinematicsUtils::totalEnergyInCM
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
Definition: G4INCLKinematicsUtils.cc:94
G4INCL::ParticleTable::getIsospin
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
Definition: G4INCLParticleTable.cc:478
G4INCL::PhaseSpaceGenerator::generateBiased
void generateBiased(const G4double sqrtS, ParticleList &particles, const size_t index, const G4double slope)
Generate a biased event in the CM system.
Definition: G4INCLPhaseSpaceGenerator.cc:98
G4INCL::Random::shoot
G4double shoot()
Definition: G4INCLRandom.cc:93
anonymous_namespace{G4CascadeDeexciteBase.cc}::zero
static const G4LorentzVector zero(0., 0., 0., 0.)

References G4INCL::FinalState::addCreatedParticle(), G4INCL::FinalState::addModifiedParticle(), angularSlope, G4INCL::PhaseSpaceGenerator::generateBiased(), G4INCL::ParticleTable::getIsospin(), G4INCL::Particle::getPosition(), G4INCL::Particle::getType(), G4INCL::KMinus, G4INCL::KPlus, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::Neutron, particle1, particle2, G4INCL::Proton, G4INCL::Particle::setType(), G4INCL::Random::shoot(), G4INCL::KinematicsUtils::totalEnergyInCM(), and anonymous_namespace{G4CascadeDeexciteBase.cc}::zero.

◆ getFinalState()

FinalState * G4INCL::IChannel::getFinalState ( )
inherited

Definition at line 50 of file G4INCLIChannel.cc.

50 {
51 FinalState *fs = new FinalState;
52 fillFinalState(fs);
53 return fs;
54 }
G4INCL::IChannel::fillFinalState
virtual void fillFinalState(FinalState *fs)=0

References G4INCL::IChannel::fillFinalState().

◆ INCL_DECLARE_ALLOCATION_POOL()

G4INCL::NNToNNKKbChannel::INCL_DECLARE_ALLOCATION_POOL ( NNToNNKKbChannel  )
private

Field Documentation

◆ angularSlope

const G4double G4INCL::NNToNNKKbChannel::angularSlope = 2.
staticprivate

Definition at line 57 of file G4INCLNNToNNKKbChannel.hh.

Referenced by fillFinalState().

◆ particle1

Particle* G4INCL::NNToNNKKbChannel::particle1
private

Definition at line 55 of file G4INCLNNToNNKKbChannel.hh.

Referenced by fillFinalState().

◆ particle2

Particle * G4INCL::NNToNNKKbChannel::particle2
private

Definition at line 55 of file G4INCLNNToNNKKbChannel.hh.

Referenced by fillFinalState().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3