G4INCL::NNToMultiPionsChannel Class Reference

#include <G4INCLNNToMultiPionsChannel.hh>

Inheritance diagram for G4INCL::NNToMultiPionsChannel:

Public Member Functions
	NNToMultiPionsChannel (const G4int, Particle *, Particle *)
virtual	~NNToMultiPionsChannel ()
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 G4INCLNNToMultiPionsChannel.hh.

Constructor & Destructor Documentation

NNToMultiPionsChannel()

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

Definition at line 51 of file G4INCLNNToMultiPionsChannel.cc.

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

~NNToMultiPionsChannel()

G4INCL::NNToMultiPionsChannel::~NNToMultiPionsChannel ( )

virtual

Definition at line 61 of file G4INCLNNToMultiPionsChannel.cc.

                                               {
 
  }

Member Function Documentation

fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 65 of file G4INCLNNToMultiPionsChannel.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);
      }
 
      const G4double sqrtS = KinematicsUtils::totalEnergyInCM(particle1, particle2);
      G4int biasIndex = ((Random::shoot()<0.5) ? 0 : 1);
      PhaseSpaceGenerator::generateBiased(sqrtS, list, biasIndex, angularSlope);
 
  }

---

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

• G4INCLNNToMultiPionsChannel.hh
• G4INCLNNToMultiPionsChannel.cc