#include <G4INCLPiNElasticChannel.hh>

Inheritance diagram for G4INCL::PiNElasticChannel:

Public Member Functions
void	fillFinalState (FinalState *fs)

FinalState *	getFinalState ()

	PiNElasticChannel (Particle , Particle )

virtual	~PiNElasticChannel ()

Private Attributes
Particle *	particle1

Particle *	particle2

Detailed Description

Definition at line 47 of file G4INCLPiNElasticChannel.hh.

Constructor & Destructor Documentation

◆ PiNElasticChannel()

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

Definition at line 47 of file G4INCLPiNElasticChannel.cc.

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

◆ ~PiNElasticChannel()

G4INCL::PiNElasticChannel::~PiNElasticChannel ( )

virtual

Definition at line 53 of file G4INCLPiNElasticChannel.cc.

                                         {
 
    }

Member Function Documentation

◆ fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLPiNElasticChannel.cc.

                                                         {
        Particle * nucleon;
        Particle * pion;
        if(particle1->isNucleon()) {
            nucleon = particle1;
            pion = particle2;
        } else {
            nucleon = particle2;
            pion = particle1;
        }
        
        G4double bpn=8e-6;  // B-parameter of exp(Bt) - (MeV/c)^-2
        G4double px_nucleon=nucleon->getMomentum().getX();
        G4double py_nucleon=nucleon->getMomentum().getY();
        G4double pz_nucleon=nucleon->getMomentum().getZ();
        G4double cnorm1=px_nucleon*px_nucleon+py_nucleon*py_nucleon;
        G4double cnorm2=std::sqrt(cnorm1);
        G4double tnorm=cnorm1+pz_nucleon*pz_nucleon;
        G4double cnorm=std::sqrt(tnorm);
        G4double btm=std::exp(-4.*tnorm*bpn);
        G4double rndm=Random::shoot();
        G4double yrn=1.-rndm*(1.-btm);
        G4double tt=std::log(yrn)/bpn;
        G4double ctheta=1.+0.5*tt/tnorm;
        G4double stheta=std::sqrt(1.-ctheta*ctheta);
        G4double t7=1.-2.*Random::shoot();
        G4double t8=std::sqrt(1-t7*t7);
        G4double t1=-py_nucleon/cnorm2;
        G4double t2=px_nucleon/cnorm2;
        G4double t4=t2*pz_nucleon/cnorm;
        G4double t5=-t1*pz_nucleon/cnorm;
        G4double t6=-cnorm2/cnorm;
        ThreeVector mom_nucleon(
                                px_nucleon*ctheta+cnorm*stheta*(t1*t7+t4*t8),
                                py_nucleon*ctheta+cnorm*stheta*(t2*t7+t5*t8),
                                pz_nucleon*ctheta+cnorm*stheta*t6*t8
                                );
        nucleon->setMomentum(mom_nucleon);
        pion->setMomentum(-mom_nucleon);
 
        G4int iso=ParticleTable::getIsospin(nucleon->getType())+ParticleTable::getIsospin(pion->getType());
        if (iso == 1 || iso == -1) {
            rndm=3*Random::shoot();
            if (rndm < 1.) {
                ParticleType nucleonType=ParticleTable::getNucleonType(-iso);
                nucleon->setType(nucleonType);
                ParticleType pionType=ParticleTable::getPionType(2*iso);
                pion->setType(pionType);
            }
            else {
                ParticleType nucleonType=ParticleTable::getNucleonType(iso);
                nucleon->setType(nucleonType);
                pion->setType(PiZero);
            }
        }
        else { // Useless
            ParticleType nucleonType=ParticleTable::getNucleonType(iso/3);
            nucleon->setType(nucleonType);
            ParticleType pionType=ParticleTable::getPionType(2*iso/3);
            pion->setType(pionType);
        }
 
        fs->addModifiedParticle(nucleon);
        fs->addModifiedParticle(pion);
    }

References G4INCL::FinalState::addModifiedParticle(), G4INCL::ParticleTable::getIsospin(), G4INCL::ParticleTable::getNucleonType(), G4INCL::ParticleTable::getPionType(), G4INCL::Particle::isNucleon(), G4InuclParticleNames::nucleon(), particle1, particle2, G4InuclParticleNames::pion(), G4INCL::PiZero, and G4INCL::Random::shoot().

◆ 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

◆ particle1

Particle* G4INCL::PiNElasticChannel::particle1

private

Definition at line 55 of file G4INCLPiNElasticChannel.hh.

Referenced by fillFinalState().

◆ particle2

Particle * G4INCL::PiNElasticChannel::particle2