#include <G4INCLNNOmegaToMultiPionsChannel.hh>

Inheritance diagram for G4INCL::NNOmegaToMultiPionsChannel:

Public Member Functions
	NNOmegaToMultiPionsChannel (const G4int, Particle , Particle )

virtual	~NNOmegaToMultiPionsChannel ()

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

Constructor & Destructor Documentation

◆ NNOmegaToMultiPionsChannel()

G4INCL::NNOmegaToMultiPionsChannel::NNOmegaToMultiPionsChannel	(	const G4int	npi,
		Particle *	p1,
		Particle *	p2
	)

Definition at line 51 of file G4INCLNNOmegaToMultiPionsChannel.cc.

    : npion(npi),
    iso1(0),
    iso2(0),
    particle1(p1),
    particle2(p2)
  {
    std::fill(isosp, isosp+4, 0);
  }

◆ ~NNOmegaToMultiPionsChannel()

G4INCL::NNOmegaToMultiPionsChannel::~NNOmegaToMultiPionsChannel ( )

virtual

Definition at line 61 of file G4INCLNNOmegaToMultiPionsChannel.cc.

                                                         {
 
  }

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 65 of file G4INCLNNOmegaToMultiPionsChannel.cc.

                                                                {
// assert(npion > 0 && npion < 5);
 
      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);
 
      isospinRepartition();
 
      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;
      for(G4int i=0; i<npion; ++i) {
        const ParticleType pionType=ParticleTable::getPionType(isosp[i]);
        Particle *pion = new Particle(pionType,zero,rcol);
        list.push_back(pion);
        fs->addCreatedParticle(pion);
      }
      Particle *omega = new Particle(Omega,zero,rcol);
      list.push_back(omega);
      fs->addCreatedParticle(omega);
   
      const G4double sqrtS = KinematicsUtils::totalEnergyInCM(particle1, particle2);
      G4int biasIndex = ((Random::shoot()<0.5) ? 0 : 1);
      PhaseSpaceGenerator::generateBiased(sqrtS, list, biasIndex, angularSlope);
 
  }