G4INCL::NKToNKChannel Class Reference

#include <G4INCLNKToNKChannel.hh>

Inheritance diagram for G4INCL::NKToNKChannel:

Public Member Functions
void	fillFinalState (FinalState *fs)
FinalState *	getFinalState ()
	NKToNKChannel (Particle *, Particle *)
virtual	~NKToNKChannel ()

Private Member Functions
	INCL_DECLARE_ALLOCATION_POOL (NKToNKChannel)

Private Attributes
Particle *	particle1
Particle *	particle2

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.

        : particle1(p1), particle2(p2)
        {}

~NKToNKChannel()

G4INCL::NKToNKChannel::~NKToNKChannel ( )

virtual

Definition at line 53 of file G4INCLNKToNKChannel.cc.

{}

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

References G4INCL::FinalState::addModifiedParticle(), G4INCL::Particle::adjustEnergyFromMomentum(), G4INCL::Particle::getType(), G4INCL::Particle::isNucleon(), G4INCL::KPlus, G4INCL::KZero, G4INCL::KinematicsUtils::momentumInCM(), G4INCL::Neutron, G4INCL::Random::normVector(), G4InuclParticleNames::nucleon(), particle1, particle2, G4INCL::Proton, G4INCL::Particle::setMomentum(), and G4INCL::Particle::setType().

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

INCL_DECLARE_ALLOCATION_POOL()

G4INCL::NKToNKChannel::INCL_DECLARE_ALLOCATION_POOL ( NKToNKChannel  )

private

Field Documentation

particle1

Particle* G4INCL::NKToNKChannel::particle1

private

Definition at line 55 of file G4INCLNKToNKChannel.hh.

Referenced by fillFinalState().

particle2

Particle * G4INCL::NKToNKChannel::particle2

private

Definition at line 55 of file G4INCLNKToNKChannel.hh.

Referenced by fillFinalState().

---

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

• source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLNKToNKChannel.hh
• source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLNKToNKChannel.cc