G4INCL::NNOmegaToMultiPionsChannel Class Reference

#include <G4INCLNNOmegaToMultiPionsChannel.hh>

Inheritance diagram for G4INCL::NNOmegaToMultiPionsChannel:

Collaboration 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 ()

Detailed Description

Definition at line 47 of file G4INCLNNOmegaToMultiPionsChannel.hh.

Constructor & Destructor Documentation

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);
  }

G4INCL::NNOmegaToMultiPionsChannel::~NNOmegaToMultiPionsChannel ( )

virtual

Definition at line 61 of file G4INCLNNOmegaToMultiPionsChannel.cc.

                                                         {

  }

Member Function Documentation

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);

  }

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLNNOmegaToMultiPionsChannel.hh
• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLNNOmegaToMultiPionsChannel.cc