G4INCL::NKElasticChannel Class Reference

#include <G4INCLNKElasticChannel.hh>

Inheritance diagram for G4INCL::NKElasticChannel:

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

Private Member Functions
	INCL_DECLARE_ALLOCATION_POOL (NKElasticChannel)

Private Attributes
Particle *	particle1
Particle *	particle2

Detailed Description

Definition at line 47 of file G4INCLNKElasticChannel.hh.

Constructor & Destructor Documentation

NKElasticChannel()

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

Definition at line 49 of file G4INCLNKElasticChannel.cc.

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

~NKElasticChannel()

G4INCL::NKElasticChannel::~NKElasticChannel ( )

virtual

Definition at line 53 of file G4INCLNKElasticChannel.cc.

{}

Member Function Documentation

fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 55 of file G4INCLNKElasticChannel.cc.

                                                        {
        
        Particle *nucleon;
        Particle *kaon;
        
        if(particle1->isNucleon()){
            nucleon = particle1;
            kaon = particle2;
        }
        else{
            nucleon = particle2;
            kaon = particle1;
        }
        
        const G4double norm = KinematicsUtils::momentumInCM(nucleon, kaon);
        
        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();
        
        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::isNucleon(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::Random::normVector(), G4InuclParticleNames::nucleon(), particle1, particle2, and G4INCL::Particle::setMomentum().

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::NKElasticChannel::INCL_DECLARE_ALLOCATION_POOL ( NKElasticChannel  )

private

Field Documentation

particle1

Particle* G4INCL::NKElasticChannel::particle1

private

Definition at line 55 of file G4INCLNKElasticChannel.hh.

Referenced by fillFinalState().

particle2

Particle * G4INCL::NKElasticChannel::particle2

private

Definition at line 55 of file G4INCLNKElasticChannel.hh.

Referenced by fillFinalState().

---

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

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