G4INCL::PiNElasticChannel Class Reference

#include <G4INCLPiNElasticChannel.hh>

Inheritance diagram for G4INCL::PiNElasticChannel:

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

Detailed Description

Definition at line 47 of file G4INCLPiNElasticChannel.hh.

Constructor & Destructor Documentation

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

Definition at line 47 of file G4INCLPiNElasticChannel.cc.

    : particle1(p1), particle2(p2)
    {

    }

G4INCL::PiNElasticChannel::~PiNElasticChannel ( )

virtual

Definition at line 53 of file G4INCLPiNElasticChannel.cc.

                                         {

    }

Member Function Documentation

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

        /*    ParticleType deltaType = DeltaZero;
         if(ParticleConfig::isPair(particle1, particle2, Proton, PiPlus)) {
         deltaType = DeltaPlusPlus;
         } else if(ParticleConfig::isPair(particle1, particle2, Neutron, PiPlus)) {
         deltaType = DeltaPlus;
         } else if(ParticleConfig::isPair(particle1, particle2, Proton, PiZero)) {
         deltaType = DeltaPlus;
         } else if(ParticleConfig::isPair(particle1, particle2, Neutron, PiZero)) {
         deltaType = DeltaZero;
         } else if(ParticleConfig::isPair(particle1, particle2, Proton, PiMinus)) {
         deltaType = DeltaZero;
         } else if(ParticleConfig::isPair(particle1, particle2, Neutron, PiMinus)) {
         deltaType = DeltaMinus;
         } else {
         INCL_ERROR("Unknown particle pair in Pi-N collision." << '\n');
         }

         G4double deltaEnergy = nucleon->getEnergy()+ pion->getEnergy();

         nucleon->setType(deltaType); // nucleon becomes the delta
         nucleon->setEnergy(deltaEnergy); // set the energy of the delta

         ThreeVector deltaMomentum = nucleon->getMomentum() + pion->getMomentum();
         nucleon->setMomentum(deltaMomentum);

         const G4double deltaMass = std::sqrt(deltaEnergy*deltaEnergy - deltaMomentum.mag2());
         nucleon->setMass(deltaMass);
         nucleon->setHelicity(1.0);

         fs->addModifiedParticle(nucleon); // nucleon became a delta
         fs->addDestroyedParticle(pion);  // pion was removed
         fs->setDeltaFixed(true);
         */
        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 {
            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);
    }

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLPiNElasticChannel.hh
• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLPiNElasticChannel.cc