#include <G4INCLNKToNKChannel.hh>

Inheritance diagram for G4INCL::NKToNKChannel:

Public Member Functions
	NKToNKChannel (Particle , Particle )

virtual	~NKToNKChannel ()

void	fillFinalState (FinalState *fs)

Public Member Functions inherited from G4INCL::IChannel
	IChannel ()

virtual	~IChannel ()

FinalState *	getFinalState ()

Detailed Description

Definition at line 47 of file G4INCLNKToNKChannel.hh.

Constructor & Destructor Documentation

◆ NKToNKChannel()

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

Definition at line 49 of file G4INCLNKToNKChannel.cc.

50 : particle1(p1), particle2(p2)

51 {}

◆ ~NKToNKChannel()

G4INCL::NKToNKChannel::~NKToNKChannel ( )

virtual

Definition at line 53 of file G4INCLNKToNKChannel.cc.

53{}

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 55 of file G4INCLNKToNKChannel.cc.

                                                     {
        
        Particle *nucleon;
        Particle *kaon;
        
        if(particle1->isNucleon()){
            nucleon = particle1;
            kaon = particle2;
        }
        else{
            nucleon = particle2;
            kaon = particle1;
        }
        
// assert(ParticleTable::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(kaon->getType()) == 0);
        
        if(kaon->getType() == KZero){
            nucleon->setType(Neutron);
            kaon->setType(KPlus);
        }
        else{
            nucleon->setType(Proton);
            kaon->setType(KZero);
        }
        
        ThreeVector mom_kaon;
        
//      const G4double pLab = KinematicsUtils::momentumInLab(kaon,nucleon);
 
/*      if(pLab==0) mom_kaon = Random::normVector();
        else{
            const G4double x = kaon->getMomentum().getX();
            const G4double y = kaon->getMomentum().getY();
            const G4double z = kaon->getMomentum().getZ();
            
            const G4double r = std::sqrt(x*x+y*y+z*z);
            const G4double rho = std::sqrt(x*x+y*y);
            
            const G4double b = 12. * pLab/2375.; // correspond to the forward slope description at 2375 MeV/c in K- p elastic
            const G4double cos_theta = std::log(Random::shoot()*(std::exp(b)-std::exp(-b))+std::exp(-b))/b;
            const G4double sin_theta = std::sqrt(1-cos_theta*cos_theta);
            
            const G4double cos_phi = std::cos(Random::shoot()*Math::twoPi);
            const G4double sin_phi = std::sqrt(1-cos_phi*cos_phi);
            
            if(rho == 0) mom_kaon = ThreeVector(sin_theta*cos_phi,sin_theta*sin_phi,cos_theta);
            // Rotation in the direction of the incident kaon
            const G4double px = x/r*cos_theta - y/rho*sin_theta*cos_phi + z/r*x/rho*sin_theta*sin_phi;
            const G4double py = y/r*cos_theta + x/rho*sin_theta*cos_phi + z/r*y/rho*sin_theta*sin_phi;
            const G4double pz = z/r*cos_theta - rho/r*sin_theta*sin_phi;
            
            mom_kaon = ThreeVector(px,py,pz);
        }*/
        
        mom_kaon = Random::normVector();
        
        G4double norm = KinematicsUtils::momentumInCM(kaon,nucleon);
 
        kaon->setMomentum(mom_kaon*norm);
        nucleon->setMomentum(-mom_kaon*norm);
        
        kaon->adjustEnergyFromMomentum();
        nucleon->adjustEnergyFromMomentum();
        
        fs->addModifiedParticle(nucleon);
        fs->addModifiedParticle(kaon);
                
    }