#include <G4INCLPiNToOmegaChannel.hh>

Inheritance diagram for G4INCL::PiNToOmegaChannel:

Public Member Functions
	PiNToOmegaChannel (Particle , Particle )

virtual	~PiNToOmegaChannel ()

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 G4INCLPiNToOmegaChannel.hh.

Constructor & Destructor Documentation

◆ PiNToOmegaChannel()

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

Definition at line 47 of file G4INCLPiNToOmegaChannel.cc.

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

◆ ~PiNToOmegaChannel()

G4INCL::PiNToOmegaChannel::~PiNToOmegaChannel ( )

virtual

Definition at line 53 of file G4INCLPiNToOmegaChannel.cc.

                                         {
 
    }

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLPiNToOmegaChannel.cc.

                                                         {
        Particle * nucleon;
        Particle * pion;
        if(particle1->isNucleon()) {
            nucleon = particle1;
            pion = particle2;
        } else {
            nucleon = particle2;
            pion = particle1;
        }
 
        G4int iso=ParticleTable::getIsospin(nucleon->getType())+ParticleTable::getIsospin(pion->getType());
// assert(iso == 1 || iso == -1);
        if (iso == 1) {
            nucleon->setType(Proton);
        }
        else if (iso == -1) {
            nucleon->setType(Neutron);
        }
        pion->setType(Omega);
 
                // Erase the parent resonance information of the nucleon and pion
            nucleon->setParentResonancePDGCode(0);
            nucleon->setParentResonanceID(0);
                pion->setParentResonancePDGCode(0);
                pion->setParentResonanceID(0);
        
//        nucleon->setEnergy(std::sqrt((nucleon->getMass())*(nucleon->getMass())+(mom_nucleon.mag()*mom_nucleon.mag())));
//        pion->setEnergy(std::sqrt((pion->getMass())*(pion->getMass())+(mom_nucleon.mag()*mom_nucleon.mag())));
  G4double sh=nucleon->getEnergy()+pion->getEnergy();
        G4double mn=nucleon->getMass();
        G4double me=pion->getMass();
        G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
        nucleon->setEnergy(en);
        G4double ee=std::sqrt(en*en-mn*mn+me*me);
        pion->setEnergy(ee);
        G4double pn=std::sqrt(en*en-mn*mn);
 
        ThreeVector mom_nucleon = Random::normVector(pn);
        
        nucleon->setMomentum(mom_nucleon);
        pion->setMomentum(-mom_nucleon);
        
  fs->addModifiedParticle(nucleon);
  fs->addModifiedParticle(pion);
  }