#include <G4INCLNNbarCEXChannel.hh>

Inheritance diagram for G4INCL::NNbarCEXChannel:

Public Member Functions
	NNbarCEXChannel (Particle , Particle )

virtual	~NNbarCEXChannel ()

void	fillFinalState (FinalState *fs)

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

virtual	~IChannel ()

FinalState *	getFinalState ()

Detailed Description

Definition at line 47 of file G4INCLNNbarCEXChannel.hh.

Constructor & Destructor Documentation

◆ NNbarCEXChannel()

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

Definition at line 49 of file G4INCLNNbarCEXChannel.cc.

50 : particle1(p1), particle2(p2)

51 {}

◆ ~NNbarCEXChannel()

G4INCL::NNbarCEXChannel::~NNbarCEXChannel ( )

virtual

Definition at line 53 of file G4INCLNNbarCEXChannel.cc.

53{}

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 55 of file G4INCLNNbarCEXChannel.cc.

                                                       {
 
        //brief ppbar
        // p pbar -> n nbar (BFMM 204)
        //
        //brief nnbar
        // n nbar -> p pbar (same as BFMM 204, but no threshold)
        //
 
        Particle *nucleon;
        Particle *antinucleon;
        
        if(particle1->isNucleon()){
            nucleon = particle1;
            antinucleon = particle2;
        }
        else{
            nucleon = particle2;
            antinucleon = particle1;
        }
        
        const G4double sqrtS = KinematicsUtils::totalEnergyInCM(nucleon, antinucleon);
 
        //setting types of new particles
        if(nucleon->getType()==Proton){
            if(antinucleon->getType()==antiProton){ //ppbar case
                nucleon->setType(Neutron);
                antinucleon->setType(antiNeutron);
            }
            else{ //pnbar case
                //no CEX for pnbar
                INCL_ERROR("We should not be in this channel " << '\n');
            }
        }
        else{ // neutron
            if(antinucleon->getType()==antiNeutron){ //nnbar case
                nucleon->setType(Proton);
                antinucleon->setType(antiProton);
            }
            else{ //npbar case
                //no CEX for npbar
                INCL_ERROR("We should not be in this channel " << '\n');
            }
        }
        
        G4double mn=nucleon->getMass();
        G4double my=antinucleon->getMass();
        
        G4double ey=(sqrtS*sqrtS+my*my-mn*mn)/(2*sqrtS);
        G4double en=std::sqrt(ey*ey-my*my+mn*mn);
        nucleon->setEnergy(en);
        antinucleon->setEnergy(ey);
        G4double py=std::sqrt(ey*ey-my*my);
        
        
        ThreeVector mom_antinucleon = Random::normVector(py);
 
        antinucleon->setMomentum(mom_antinucleon);
        nucleon->setMomentum(-mom_antinucleon);
        
        fs->addModifiedParticle(nucleon);
        fs->addModifiedParticle(antinucleon);
                
    }