#include <G4INCLRecombinationChannel.hh>

Inheritance diagram for G4INCL::RecombinationChannel:

Public Member Functions
void	fillFinalState (FinalState *fs)

FinalState *	getFinalState ()

	RecombinationChannel (Particle p1, Particle p2)

virtual	~RecombinationChannel ()

Private Attributes
Particle *	theDelta

Particle *	theNucleon

Detailed Description

Definition at line 55 of file G4INCLRecombinationChannel.hh.

Constructor & Destructor Documentation

◆ RecombinationChannel()

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

Definition at line 56 of file G4INCLRecombinationChannel.cc.

  {
    if(p1->isDelta()) {
// assert(p2->isNucleon());
      theDelta = p1;
      theNucleon = p2;
    } else {
// assert(p1->isNucleon());
      theDelta = p2;
      theNucleon = p1;
    }
  }

References G4INCL::Particle::isDelta(), theDelta, and theNucleon.

◆ ~RecombinationChannel()

G4INCL::RecombinationChannel::~RecombinationChannel ( )

virtual

Definition at line 69 of file G4INCLRecombinationChannel.cc.

70 {

71 }

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 73 of file G4INCLRecombinationChannel.cc.

  {
    // Compute the total available energy in the CM
    const G4double sqrts = KinematicsUtils::totalEnergyInCM(theDelta, theNucleon);
 
    // Assign the types of the final-state particles
    switch(theDelta->getType()) {
      case DeltaPlusPlus:
// assert(theNucleon->getType()!=Proton);
        theDelta->setType(Proton);
        theNucleon->setType(Proton);
        break;
      case DeltaPlus:
        theDelta->setType(Proton);
        break;
      case DeltaZero:
        theDelta->setType(Neutron);
        break;
      case DeltaMinus:
// assert(theNucleon->getType()!=Neutron);
        theDelta->setType(Neutron);
        theNucleon->setType(Neutron);
        break;
      default:
        INCL_ERROR("Unknown particle type in RecombinationChannel" << '\n');
        break;
    }
 
    // Calculate the momenta of the nucleons in the final state
    const G4double pCM = KinematicsUtils::momentumInCM(sqrts, theDelta->getMass(), theNucleon->getMass());
 
    // The angular distribution of final-state nucleons is isotropic
    ThreeVector momentum = Random::normVector(pCM);
 
    // Assign the momenta
    theDelta->setMomentum(momentum);
    theNucleon->setMomentum(-momentum);
 
    // Update the kinetic energies
    theDelta->adjustEnergyFromMomentum();
    theNucleon->adjustEnergyFromMomentum();
 
    // Create the final state
    fs->addModifiedParticle(theDelta);
    fs->addModifiedParticle(theNucleon);
 
  }

References G4INCL::FinalState::addModifiedParticle(), G4INCL::Particle::adjustEnergyFromMomentum(), G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::Particle::getMass(), G4INCL::Particle::getType(), INCL_ERROR, G4INCL::KinematicsUtils::momentumInCM(), G4INCL::Neutron, G4INCL::Random::normVector(), G4INCL::Proton, G4INCL::Particle::setMomentum(), G4INCL::Particle::setType(), theDelta, theNucleon, and G4INCL::KinematicsUtils::totalEnergyInCM().

◆ getFinalState()

FinalState * G4INCL::IChannel::getFinalState ( )

inherited

Definition at line 50 of file G4INCLIChannel.cc.

                                      {
    FinalState *fs = new FinalState;
    fillFinalState(fs);
    return fs;
  }

References G4INCL::IChannel::fillFinalState().

Field Documentation

◆ theDelta

Particle * G4INCL::RecombinationChannel::theDelta

private

Definition at line 64 of file G4INCLRecombinationChannel.hh.

Referenced by fillFinalState(), and RecombinationChannel().

◆ theNucleon

Particle* G4INCL::RecombinationChannel::theNucleon