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

#include <G4INCLNKToNK2piChannel.hh>

Inheritance diagram for G4INCL::NKToNK2piChannel:
G4INCL::IChannel

Public Member Functions

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

Private Member Functions

 INCL_DECLARE_ALLOCATION_POOL (NKToNK2piChannel)
 

Private Attributes

Particle * particle1
 
Particle * particle2
 

Static Private Attributes

static const G4double angularSlope = 4.
 

Detailed Description

Definition at line 47 of file G4INCLNKToNK2piChannel.hh.

Constructor & Destructor Documentation

◆ NKToNK2piChannel()

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

Definition at line 51 of file G4INCLNKToNK2piChannel.cc.

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

◆ ~NKToNK2piChannel()

G4INCL::NKToNK2piChannel::~NKToNK2piChannel ( )
virtual

Definition at line 55 of file G4INCLNKToNK2piChannel.cc.

55{}

Member Function Documentation

◆ fillFinalState()

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

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLNKToNK2piChannel.cc.

57 {
58
59 // p K+ -> p K+ pi+ pi- (1)
60 // p K+ -> p K+ pi0 pi0 (1/8)
61 // p K+ -> p K0 pi+ pi0 (1)
62 // p K+ -> n K+ pi+ pi0 (1/2)
63 // p K+ -> n K0 pi+ pi+ (1/4)
64 //
65 // p K0 -> p K0 pi+ pi- (1)
66 // p K0 -> p K0 pi0 pi0 (1/8)
67 // p K0 -> p K+ pi0 pi- (1)
68 // p K0 -> n K+ pi+ pi- (1/4)
69 // p K0 -> n K+ pi0 pi0 (1/4)
70 // p K0 -> n K0 pi+ pi0 (1/2)
71
72 Particle *nucleon;
73 Particle *kaon;
74
75 if(particle1->isNucleon()){
76 nucleon = particle1;
77 kaon = particle2;
78 }
79 else{
80 nucleon = particle2;
81 kaon = particle1;
82 }
83
84 const G4double sqrtS = KinematicsUtils::totalEnergyInCM(nucleon, kaon);
85
86 const G4int iso = ParticleTable::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(kaon->getType());
87 const G4int iso_n = ParticleTable::getIsospin(nucleon->getType());
88 G4double rdm = Random::shoot();
89
90 ParticleType Pion1Type;
91 ParticleType Pion2Type;
92
93 if(iso == 2 || iso == -2){
94 if(rdm*23. < 8.){
95 Pion1Type = PiPlus;
96 Pion2Type = PiMinus;
97 }
98 else if(rdm*23. < 9.){
99 Pion1Type = PiZero;
100 Pion2Type = PiZero;
101 }
102 else if(rdm*23. < 17.){
103 Pion1Type = ParticleTable::getPionType(iso);
104 Pion2Type = PiZero;
105 kaon->setType(ParticleTable::getKaonType(-iso/2));
106 }
107 else if(rdm*23. < 21.){
108 Pion1Type = ParticleTable::getPionType(iso);
109 Pion2Type = PiZero;
110 nucleon->setType(ParticleTable::getNucleonType(-iso/2));
111 }
112 else{
113 Pion1Type = ParticleTable::getPionType(iso);
114 Pion2Type = ParticleTable::getPionType(iso);
115 kaon->setType(ParticleTable::getKaonType(-iso/2));
116 nucleon->setType(ParticleTable::getNucleonType(-iso/2));
117 }
118 }
119 else{
120 if(rdm*25. < 8.){
121 Pion1Type = PiPlus;
122 Pion2Type = PiMinus;
123 }
124 else if(rdm*25. < 9.){
125 Pion1Type = PiZero;
126 Pion2Type = PiZero;
127 }
128 else if(rdm*25. < 17.){
129 Pion1Type = ParticleTable::getPionType(-2*iso_n);
130 Pion2Type = PiZero;
131 kaon->setType(ParticleTable::getKaonType(iso_n));
132 }
133 else if(rdm*25. < 19.){
134 Pion1Type = PiPlus;
135 Pion2Type = PiMinus;
136 kaon->setType(ParticleTable::getKaonType(iso_n));
137 nucleon->setType(ParticleTable::getNucleonType(-iso_n));
138 }
139 else if(rdm*25. < 21.){
140 Pion1Type = PiZero;
141 Pion2Type = PiZero;
142 kaon->setType(ParticleTable::getKaonType(iso_n));
143 nucleon->setType(ParticleTable::getNucleonType(-iso_n));
144 }
145 else{
146 Pion1Type = ParticleTable::getPionType(2*iso_n);
147 Pion2Type = PiZero;
148 nucleon->setType(ParticleTable::getNucleonType(-iso_n));
149 }
150 }
151
152 ParticleList list;
153 list.push_back(nucleon);
154 list.push_back(kaon);
155 const ThreeVector &rcol1 = nucleon->getPosition();
156 const ThreeVector &rcol2 = kaon->getPosition();
157 const ThreeVector zero;
158 Particle *pion1 = new Particle(Pion1Type,zero,rcol1);
159 Particle *pion2 = new Particle(Pion2Type,zero,rcol2);
160 list.push_back(pion1);
161 list.push_back(pion2);
162
163 PhaseSpaceGenerator::generateBiased(sqrtS, list, 0, angularSlope);
164
165 fs->addModifiedParticle(nucleon);
166 fs->addModifiedParticle(kaon);
167 fs->addCreatedParticle(pion1);
168 fs->addCreatedParticle(pion2);
169
170 }
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4INCL::NKToNK2piChannel::angularSlope
static const G4double angularSlope
Definition: G4INCLNKToNK2piChannel.hh:57
G4INCL::Particle::isNucleon
G4bool isNucleon() const
Definition: G4INCLParticle.hh:303
G4INCL::KinematicsUtils::totalEnergyInCM
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
Definition: G4INCLKinematicsUtils.cc:94
G4INCL::ParticleTable::getKaonType
ParticleType getKaonType(const G4int isosp)
Get the type of kaon.
Definition: G4INCLParticleTable.cc:1293
G4INCL::ParticleTable::getIsospin
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
Definition: G4INCLParticleTable.cc:478
G4INCL::ParticleTable::getNucleonType
ParticleType getNucleonType(const G4int isosp)
Get the type of nucleon.
Definition: G4INCLParticleTable.cc:1253
G4INCL::ParticleTable::getPionType
ParticleType getPionType(const G4int isosp)
Get the type of pion.
Definition: G4INCLParticleTable.cc:1240
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
G4InuclParticleNames::nucleon
G4bool nucleon(G4int ityp)
Definition: G4InuclParticleNames.hh:94
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::ParticleTable::getKaonType(), G4INCL::ParticleTable::getNucleonType(), G4INCL::ParticleTable::getPionType(), G4INCL::Particle::getPosition(), G4INCL::Particle::getType(), G4INCL::Particle::isNucleon(), G4InuclParticleNames::nucleon(), particle1, particle2, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, 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::NKToNK2piChannel::INCL_DECLARE_ALLOCATION_POOL ( NKToNK2piChannel  )
private

Field Documentation

◆ angularSlope

const G4double G4INCL::NKToNK2piChannel::angularSlope = 4.
staticprivate

Definition at line 57 of file G4INCLNKToNK2piChannel.hh.

Referenced by fillFinalState().

◆ particle1

Particle* G4INCL::NKToNK2piChannel::particle1
private

Definition at line 55 of file G4INCLNKToNK2piChannel.hh.

Referenced by fillFinalState().

◆ particle2

Particle * G4INCL::NKToNK2piChannel::particle2
private

Definition at line 55 of file G4INCLNKToNK2piChannel.hh.

Referenced by fillFinalState().

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