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

#include <G4INCLNpiToSKpiChannel.hh>

Inheritance diagram for G4INCL::NpiToSKpiChannel:
G4INCL::IChannel

Public Member Functions

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

Private Member Functions

 INCL_DECLARE_ALLOCATION_POOL (NpiToSKpiChannel)
 

Private Attributes

Particle * particle1
 
Particle * particle2
 

Static Private Attributes

static const G4double angularSlope = 6.
 

Detailed Description

Definition at line 47 of file G4INCLNpiToSKpiChannel.hh.

Constructor & Destructor Documentation

◆ NpiToSKpiChannel()

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

Definition at line 51 of file G4INCLNpiToSKpiChannel.cc.

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

◆ ~NpiToSKpiChannel()

G4INCL::NpiToSKpiChannel::~NpiToSKpiChannel ( )
virtual

Definition at line 55 of file G4INCLNpiToSKpiChannel.cc.

55{}

Member Function Documentation

◆ fillFinalState()

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

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLNpiToSKpiChannel.cc.

57 {
58
59 // p pi+ -> S+ pi+ K0 (5/4)
60 // p pi+ -> S+ pi0 K+ (3/4)
61 // p pi+ -> S0 pi+ K+ (1/4)
62 //
63 // p pi0 -> S+ pi0 K0 (1/2)
64 // p pi0 -> S+ pi- K+ (1/2)
65 // p pi0 -> S0 pi+ K0 (3/4)
66 // p pi0 -> S0 pi0 K+ (3/8)
67 // p pi0 -> S- pi+ K+ (1/2)
68 //
69 // p pi- -> S+ pi- K0 (3/8)
70 // p pi- -> S0 pi0 K0 (5/8)
71 // p pi- -> S0 pi- K+ (5/8)
72 // p pi- -> S- pi+ K0 (1)
73 // p pi- -> S- pi0 K+ (3/8)
74
75 Particle *nucleon;
76 Particle *pion;
77
78
79 if(particle1->isNucleon()){
80 nucleon = particle1;
81 pion = particle2;
82 }
83 else{
84 nucleon = particle2;
85 pion = particle1;
86 }
87
88 const G4double sqrtS = KinematicsUtils::totalEnergyInCM(nucleon, pion);
89
90 const G4int iso = ParticleTable::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(pion->getType());
91 G4double rdm = Random::shoot();
92
93 ParticleType KaonType;
94
95 if(iso == 3 || iso == -3){
96 if(rdm*9. < 5.){
97 KaonType = ParticleTable::getKaonType(-iso/3);
98 nucleon->setType(ParticleTable::getSigmaType(iso/3*2));
99 }
100 else if(rdm*9. < 8.){
101 KaonType = ParticleTable::getKaonType(iso/3);
102 pion->setType(PiZero);
103 nucleon->setType(ParticleTable::getSigmaType(iso/3*2));
104 }
105 else{
106 KaonType = ParticleTable::getKaonType(iso/3);
107 nucleon->setType(SigmaZero);
108 }
109 }
110 else if(pion->getType() == PiZero){
111 if(rdm*21. < 4.){
112 KaonType = ParticleTable::getKaonType(-iso);
113 nucleon->setType(ParticleTable::getSigmaType(iso*2));
114 }
115 else if(rdm*21. < 8.){
116 KaonType = ParticleTable::getKaonType(iso);
117 pion->setType(ParticleTable::getPionType(-2*iso));
118 nucleon->setType(ParticleTable::getSigmaType(iso*2));
119 }
120 else if(rdm*21. < 14.){
121 KaonType = ParticleTable::getKaonType(-iso);
122 pion->setType(ParticleTable::getPionType(2*iso));
123 nucleon->setType(SigmaZero);
124 }
125 else if(rdm*21. < 17.){
126 KaonType = ParticleTable::getKaonType(iso);
127 nucleon->setType(SigmaZero);
128 }
129 else{
130 KaonType = ParticleTable::getKaonType(iso);
131 pion->setType(ParticleTable::getPionType(2*iso));
132 nucleon->setType(ParticleTable::getSigmaType(-iso*2));
133 }
134 }
135 else{
136 if(rdm*24. < 3.){
137 KaonType = ParticleTable::getKaonType(iso);
138 nucleon->setType(ParticleTable::getSigmaType(-iso*2));
139 }
140 else if(rdm*24. < 8.){
141 KaonType = ParticleTable::getKaonType(iso);
142 pion->setType(PiZero);
143 nucleon->setType(SigmaZero);
144 }
145 else if(rdm*24. < 13.){
146 KaonType = ParticleTable::getKaonType(-iso);
147 nucleon->setType(SigmaZero);
148 }
149 else if(rdm*24. < 21.){
150 KaonType = ParticleTable::getKaonType(iso);
151 pion->setType(ParticleTable::getPionType(-2*iso));
152 nucleon->setType(ParticleTable::getSigmaType(iso*2));
153 }
154 else{
155 KaonType = ParticleTable::getKaonType(-iso);
156 pion->setType(PiZero);
157 nucleon->setType(ParticleTable::getSigmaType(iso*2));
158 }
159 }
160
161 ParticleList list;
162 list.push_back(nucleon);
163 list.push_back(pion);
164 const ThreeVector &rcol = nucleon->getPosition();
165 const ThreeVector zero;
166 Particle *kaon = new Particle(KaonType,zero,rcol);
167 list.push_back(kaon);
168
169 PhaseSpaceGenerator::generateBiased(sqrtS, list, 0, angularSlope);
170
171 INCL_DEBUG("NpiToSKpi " << (kaon->getMomentum().theta()) * 180. / G4INCL::Math::pi << '\n');
172
173 fs->addModifiedParticle(nucleon);
174 fs->addModifiedParticle(pion);
175 fs->addCreatedParticle(kaon);
176
177 }
INCL_DEBUG
#define INCL_DEBUG(x)
Definition: G4INCLLogger.hh:240
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4INCL::NpiToSKpiChannel::angularSlope
static const G4double angularSlope
Definition: G4INCLNpiToSKpiChannel.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::Math::pi
const G4double pi
Definition: G4INCLGlobals.hh:68
G4INCL::ParticleTable::getKaonType
ParticleType getKaonType(const G4int isosp)
Get the type of kaon.
Definition: G4INCLParticleTable.cc:1293
G4INCL::ParticleTable::getSigmaType
ParticleType getSigmaType(const G4int isosp)
Get the type of sigma.
Definition: G4INCLParticleTable.cc:1280
G4INCL::ParticleTable::getIsospin
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
Definition: G4INCLParticleTable.cc:478
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::pion
G4bool pion(G4int ityp)
Definition: G4InuclParticleNames.hh:90
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::Particle::getMomentum(), G4INCL::ParticleTable::getPionType(), G4INCL::ParticleTable::getSigmaType(), INCL_DEBUG, G4INCL::Particle::isNucleon(), G4InuclParticleNames::nucleon(), particle1, particle2, G4INCL::Math::pi, G4InuclParticleNames::pion(), G4INCL::PiZero, G4INCL::Random::shoot(), G4INCL::SigmaZero, G4INCL::ThreeVector::theta(), 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::NpiToSKpiChannel::INCL_DECLARE_ALLOCATION_POOL ( NpiToSKpiChannel  )
private

Field Documentation

◆ angularSlope

const G4double G4INCL::NpiToSKpiChannel::angularSlope = 6.
staticprivate

Definition at line 57 of file G4INCLNpiToSKpiChannel.hh.

Referenced by fillFinalState().

◆ particle1

Particle* G4INCL::NpiToSKpiChannel::particle1
private

Definition at line 55 of file G4INCLNpiToSKpiChannel.hh.

Referenced by fillFinalState().

◆ particle2

Particle * G4INCL::NpiToSKpiChannel::particle2
private

Definition at line 55 of file G4INCLNpiToSKpiChannel.hh.

Referenced by fillFinalState().

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