#include <G4INCLEtaOrOmegaNToSKChannel.hh>

Inheritance diagram for G4INCL::EtaOrOmegaNToSKChannel:

Public Member Functions
	EtaOrOmegaNToSKChannel (Particle , Particle )
virtual	~EtaOrOmegaNToSKChannel ()
void	fillFinalState (FinalState *fs)
Public Member Functions inherited from G4INCL::IChannel
	IChannel ()
virtual	~IChannel ()
FinalState *	getFinalState ()

Detailed Description

Definition at line 47 of file G4INCLEtaOrOmegaNToSKChannel.hh.

Constructor & Destructor Documentation

◆ EtaOrOmegaNToSKChannel()

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

Definition at line 47 of file G4INCLEtaOrOmegaNToSKChannel.cc.

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

◆ ~EtaOrOmegaNToSKChannel()

G4INCL::EtaOrOmegaNToSKChannel::~EtaOrOmegaNToSKChannel ( )

virtual

Definition at line 53 of file G4INCLEtaOrOmegaNToSKChannel.cc.

                                                   {
 
    }

Member Function Documentation

◆ fillFinalState()

void G4INCL::EtaOrOmegaNToSKChannel::fillFinalState ( FinalState * fs )

virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLEtaOrOmegaNToSKChannel.cc.

                                                              {
        Particle * nucleon;
        Particle * meson;
 
        if(particle1->isNucleon()) {
            nucleon = particle1;
            meson = particle2;
        } else {
            nucleon = particle2;
            meson = particle1;
        }
 
        const G4double r2 = Random::shoot();
 
        if (nucleon->getType() == Neutron) {
            if (r2*3. < 2.) {
                nucleon->setType(SigmaZero);
                meson->setType(KZero);
            }
            else {
                nucleon->setType(SigmaMinus);
                meson->setType(KPlus);
            }
        }
        else {
            if (r2*3. < 2.) {
                nucleon->setType(SigmaPlus);
                meson->setType(KZero);
            }
            else {
                nucleon->setType(SigmaZero);
                meson->setType(KPlus);
            }
        }
 
        G4double sh=nucleon->getEnergy()+meson->getEnergy();
        G4double mn=nucleon->getMass();
        G4double me=meson->getMass();
        G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
        nucleon->setEnergy(en);
        G4double ee=std::sqrt(en*en-mn*mn+me*me);
        meson->setEnergy(ee);
        G4double pn=std::sqrt(en*en-mn*mn);
 
        ThreeVector mom_nucleon = Random::normVector(pn);
        
        nucleon->setMomentum(mom_nucleon);
        meson->setMomentum(-mom_nucleon);
        
        fs->addModifiedParticle(nucleon);
        fs->addModifiedParticle(meson);
    }

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