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

#include <G4INCLNKbToL2piChannel.hh>

Inheritance diagram for G4INCL::NKbToL2piChannel:
G4INCL::IChannel

Public Member Functions

void fillFinalState (FinalState *fs)
 
FinalStategetFinalState ()
 
 NKbToL2piChannel (Particle *, Particle *)
 
virtual ~NKbToL2piChannel ()
 

Private Member Functions

 INCL_DECLARE_ALLOCATION_POOL (NKbToL2piChannel)
 

Private Attributes

Particleparticle1
 
Particleparticle2
 

Static Private Attributes

static const G4double angularSlope = 4.
 

Detailed Description

Definition at line 47 of file G4INCLNKbToL2piChannel.hh.

Constructor & Destructor Documentation

◆ NKbToL2piChannel()

G4INCL::NKbToL2piChannel::NKbToL2piChannel ( Particle p1,
Particle p2 
)

Definition at line 51 of file G4INCLNKbToL2piChannel.cc.

◆ ~NKbToL2piChannel()

G4INCL::NKbToL2piChannel::~NKbToL2piChannel ( )
virtual

Definition at line 55 of file G4INCLNKbToL2piChannel.cc.

55{}

Member Function Documentation

◆ fillFinalState()

void G4INCL::NKbToL2piChannel::fillFinalState ( FinalState fs)
virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLNKbToL2piChannel.cc.

57 {
58
59 // p K0b -> L pi+ pi0 (1)
60 //
61 // p K- -> L pi+ pi- (1)
62 // p K- -> L pi0 pi0 (1/4)
63
64 Particle *nucleon;
65 Particle *kaon;
66
67 if(particle1->isNucleon()){
69 kaon = particle2;
70 }
71 else{
73 kaon = particle1;
74 }
75
77
78 const G4int iso = ParticleTable::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(kaon->getType());
79 const G4int iso_n = ParticleTable::getIsospin(nucleon->getType());
80 G4double rdm = Random::shoot();
81
82 ParticleType PionType;
83
84 if(iso == 2 || iso == -2){
85 PionType = ParticleTable::getPionType(iso);
86 kaon->setType(PiZero);
87 }
88 else{
89 if(rdm*5. < 4.){
90 PionType = ParticleTable::getPionType(2*iso_n);
91 kaon->setType(ParticleTable::getPionType(-2*iso_n));
92 }
93 else{
94 PionType = PiZero;
95 kaon->setType(PiZero);
96 }
97 }
98
99 nucleon->setType(Lambda);
100
101 ParticleList list;
102 list.push_back(nucleon);
103 list.push_back(kaon);
104 const ThreeVector &rcol = nucleon->getPosition();
105 const ThreeVector zero;
106 Particle *pion = new Particle(PionType,zero,rcol);
107 list.push_back(pion);
108
110
111 fs->addModifiedParticle(nucleon);
112 fs->addModifiedParticle(kaon);
113 fs->addCreatedParticle(pion);
114
115 }
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
static const G4double angularSlope
G4bool isNucleon() const
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
ParticleType getPionType(const G4int isosp)
Get the type of pion.
void generateBiased(const G4double sqrtS, ParticleList &particles, const size_t index, const G4double slope)
Generate a biased event in the CM system.
G4double shoot()
Definition: G4INCLRandom.cc:93
G4bool pion(G4int ityp)
G4bool nucleon(G4int ityp)
static const G4LorentzVector zero(0., 0., 0., 0.)

References G4INCL::FinalState::addCreatedParticle(), G4INCL::FinalState::addModifiedParticle(), angularSlope, G4INCL::PhaseSpaceGenerator::generateBiased(), G4INCL::ParticleTable::getIsospin(), G4INCL::ParticleTable::getPionType(), G4INCL::Particle::getType(), G4INCL::Particle::isNucleon(), G4INCL::Lambda, G4InuclParticleNames::nucleon(), particle1, particle2, G4InuclParticleNames::pion(), 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;
53 return fs;
54 }
virtual void fillFinalState(FinalState *fs)=0

References G4INCL::IChannel::fillFinalState().

◆ INCL_DECLARE_ALLOCATION_POOL()

G4INCL::NKbToL2piChannel::INCL_DECLARE_ALLOCATION_POOL ( NKbToL2piChannel  )
private

Field Documentation

◆ angularSlope

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

Definition at line 57 of file G4INCLNKbToL2piChannel.hh.

Referenced by fillFinalState().

◆ particle1

Particle* G4INCL::NKbToL2piChannel::particle1
private

Definition at line 55 of file G4INCLNKbToL2piChannel.hh.

Referenced by fillFinalState().

◆ particle2

Particle * G4INCL::NKbToL2piChannel::particle2
private

Definition at line 55 of file G4INCLNKbToL2piChannel.hh.

Referenced by fillFinalState().


The documentation for this class was generated from the following files: