#include <G4INCLPiNElasticChannel.hh>

Inheritance diagram for G4INCL::PiNElasticChannel:

Public Member Functions
	PiNElasticChannel (Particle , Particle )

virtual	~PiNElasticChannel ()

void	fillFinalState (FinalState *fs)

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

virtual	~IChannel ()

FinalState *	getFinalState ()

virtual void	fillFinalState (FinalState *fs)=0

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);
 
    ParticleType startingNucleonType = nucleon->getType();
    ParticleType startingPionType = pion->getType();
        
        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);
        }
 
    // Erase the parent resonance information if the nucleon or pion changes type
    if ( startingNucleonType != nucleon->getType() || startingPionType != pion->getType() ) {
      nucleon->setParentResonancePDGCode(0);
      nucleon->setParentResonanceID(0);
      pion->setParentResonancePDGCode(0);
      pion->setParentResonanceID(0);
    }
    
        fs->addModifiedParticle(nucleon);
        fs->addModifiedParticle(pion);
    }