#include <G4INCLNSToNSChannel.hh>

Inheritance diagram for G4INCL::NSToNSChannel:

Public Member Functions
	NSToNSChannel (Particle , Particle )

virtual	~NSToNSChannel ()

void	fillFinalState (FinalState *fs)

Public Member Functions inherited from G4INCL::IChannel
	IChannel ()

virtual	~IChannel ()

FinalState *	getFinalState ()

Detailed Description

Definition at line 47 of file G4INCLNSToNSChannel.hh.

Constructor & Destructor Documentation

◆ NSToNSChannel()

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

Definition at line 49 of file G4INCLNSToNSChannel.cc.

50 : particle1(p1), particle2(p2)

51 {}

◆ ~NSToNSChannel()

G4INCL::NSToNSChannel::~NSToNSChannel ( )

virtual

Definition at line 53 of file G4INCLNSToNSChannel.cc.

53{}

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 55 of file G4INCLNSToNSChannel.cc.

                                                     {
        
        Particle *nucleon;
        Particle *sigma;
        
        if(particle1->isNucleon()){
            nucleon = particle1;
            sigma = particle2;
        }
        else{
            nucleon = particle2;
            sigma = particle1;
        }
        
        const G4double sqrtS = KinematicsUtils::totalEnergyInCM(nucleon, sigma);
        const G4int iso = ParticleTable::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(sigma->getType());
// assert(iso == -1 || iso == 1);
        
        if(sigma->getType() == SigmaZero){
            nucleon->setType(ParticleTable::getNucleonType(-iso));
            sigma->setType(ParticleTable::getSigmaType(2*iso));
        }
        else{
            nucleon->setType(ParticleTable::getNucleonType(iso));
            sigma->setType(SigmaZero);
        }
        
        G4double mn=nucleon->getMass();
        G4double my=sigma->getMass();
        
        G4double ey=(sqrtS*sqrtS+my*my-mn*mn)/(2*sqrtS);
        G4double en=std::sqrt(ey*ey-my*my+mn*mn);
        nucleon->setEnergy(en);
        sigma->setEnergy(ey);
        G4double py=std::sqrt(ey*ey-my*my);
        
        ThreeVector mom_hyperon = Random::normVector(py);
 
        sigma->setMomentum(mom_hyperon);
        nucleon->setMomentum(-mom_hyperon);
        
        sigma->adjustEnergyFromMomentum();
        nucleon->adjustEnergyFromMomentum();
        
        fs->addModifiedParticle(nucleon);
        fs->addModifiedParticle(sigma);
                
    }