G4INCL::PiNToMultiPionsChannel Class Reference

#include <G4INCLPiNToMultiPionsChannel.hh>

Inheritance diagram for G4INCL::PiNToMultiPionsChannel:

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

Constructor & Destructor Documentation

PiNToMultiPionsChannel()

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

Definition at line 51 of file G4INCLPiNToMultiPionsChannel.cc.

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

~PiNToMultiPionsChannel()

G4INCL::PiNToMultiPionsChannel::~PiNToMultiPionsChannel ( )

virtual

Definition at line 60 of file G4INCLPiNToMultiPionsChannel.cc.

                                                 {
 
  }

Member Function Documentation

fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 64 of file G4INCLPiNToMultiPionsChannel.cc.

                                                            {
 
// assert(npion > 1 && npion < 5);
 
    Particle * nucleon;
    Particle * pion;
    if(particle1->isNucleon()) {
      nucleon = particle1;
      pion = particle2;
    } else {
      nucleon = particle2;
      pion = particle1;
    }
 
      // Erase the parent resonance information of the nucleon and pion
      nucleon->setParentResonancePDGCode(0);
      nucleon->setParentResonanceID(0);
      pion->setParentResonancePDGCode(0);
      pion->setParentResonanceID(0);
    
      G4int ipi=ParticleTable::getIsospin(pion->getType());
      ind2=ParticleTable::getIsospin(nucleon->getType());
 
      ParticleList list;
      list.push_back(nucleon);
      list.push_back(pion);
      fs->addModifiedParticle(nucleon);
      fs->addModifiedParticle(pion);
 
      isospinRepartition(ipi);
 
      const ParticleType tn=ParticleTable::getNucleonType(ind2);
      nucleon->setType(tn);
      ParticleType pionType=ParticleTable::getPionType(isosp[0]);
      pion->setType(pionType);
      const ThreeVector &rcolpion = pion->getPosition();
      const ThreeVector zero;
      for(G4int i=1; i<npion; ++i) {
        pionType=ParticleTable::getPionType(isosp[i]);
        Particle *newPion = new Particle(pionType,zero,rcolpion);
        newPion->setType(pionType);
        list.push_back(newPion);
        fs->addCreatedParticle(newPion);
      }
 
      const G4double sqrtS = KinematicsUtils::totalEnergyInCM(nucleon, pion);
      PhaseSpaceGenerator::generateBiased(sqrtS, list, 0, angularSlope);
 
  }

---

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

• G4INCLPiNToMultiPionsChannel.hh
• G4INCLPiNToMultiPionsChannel.cc