#include <G4INCLNNToNNEtaChannel.hh>

Inheritance diagram for G4INCL::NNToNNEtaChannel:

Public Member Functions
	NNToNNEtaChannel (Particle , Particle )

virtual	~NNToNNEtaChannel ()

void	fillFinalState (FinalState *fs)

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

virtual	~IChannel ()

FinalState *	getFinalState ()

virtual void	fillFinalState (FinalState *fs)=0

Detailed Description

Definition at line 47 of file G4INCLNNToNNEtaChannel.hh.

Constructor & Destructor Documentation

◆ NNToNNEtaChannel()

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

Definition at line 51 of file G4INCLNNToNNEtaChannel.cc.

    : iso1(0),
    iso2(0),
    particle1(p1),
    particle2(p2)
    {}

◆ ~NNToNNEtaChannel()

G4INCL::NNToNNEtaChannel::~NNToNNEtaChannel ( )

virtual

Definition at line 58 of file G4INCLNNToNNEtaChannel.cc.

                                       {
        
    }

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 62 of file G4INCLNNToNNEtaChannel.cc.

                                                        {
        
        iso1=ParticleTable::getIsospin(particle1->getType());
        iso2=ParticleTable::getIsospin(particle2->getType());
        
        ParticleList list;
        list.push_back(particle1);
        list.push_back(particle2);
        fs->addModifiedParticle(particle1);
        fs->addModifiedParticle(particle2);
                
        const G4double sqrtS = KinematicsUtils::totalEnergyInCM(particle1, particle2);
 
        const ParticleType tn1=ParticleTable::getNucleonType(iso1);
        particle1->setType(tn1);
        const ParticleType tn2=ParticleTable::getNucleonType(iso2);
        particle2->setType(tn2);
        const ThreeVector &rcolnucleon1 = particle1->getPosition();
        const ThreeVector &rcolnucleon2 = particle2->getPosition();
        const ThreeVector rcol = (rcolnucleon1+rcolnucleon2)*0.5;
        const ThreeVector zero;
        ParticleType etaType=Eta;
        Particle *etaCreated = new Particle(etaType,zero,rcol);
        list.push_back(etaCreated);
        fs->addCreatedParticle(etaCreated);
 
        G4int biasIndex = ((Random::shoot()<0.5) ? 0 : 1);
        PhaseSpaceGenerator::generateBiased(sqrtS, list, biasIndex, angularSlope);
        
    }