G4INCL::PiNToOmegaChannel Class Reference

#include <G4INCLPiNToOmegaChannel.hh>

Inheritance diagram for G4INCL::PiNToOmegaChannel:

Collaboration diagram for G4INCL::PiNToOmegaChannel:

Public Member Functions
	PiNToOmegaChannel (Particle *, Particle *)
virtual	~PiNToOmegaChannel ()
void	fillFinalState (FinalState *fs)
Public Member Functions inherited from G4INCL::IChannel
	IChannel ()
virtual	~IChannel ()
FinalState *	getFinalState ()

Detailed Description

Definition at line 47 of file G4INCLPiNToOmegaChannel.hh.

Constructor & Destructor Documentation

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

Definition at line 47 of file G4INCLPiNToOmegaChannel.cc.

    : particle1(p1), particle2(p2)
    {

    }

G4INCL::PiNToOmegaChannel::~PiNToOmegaChannel ( )

virtual

Definition at line 53 of file G4INCLPiNToOmegaChannel.cc.

                                         {

    }

Member Function Documentation

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

virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLPiNToOmegaChannel.cc.

                                                         {
        Particle * nucleon;
        Particle * pion;
        if(particle1->isNucleon()) {
            nucleon = particle1;
            pion = particle2;
        } else {
            nucleon = particle2;
            pion = particle1;
        }

        
        G4int iso=ParticleTable::getIsospin(nucleon->getType())+ParticleTable::getIsospin(pion->getType());
// assert(iso == 1 || iso == -1);
        if (iso == 1) {
            nucleon->setType(Proton);
        }
        else if (iso == -1) {
            nucleon->setType(Neutron);
        }
        pion->setType(Omega);
//        nucleon->setEnergy(std::sqrt((nucleon->getMass())*(nucleon->getMass())+(mom_nucleon.mag()*mom_nucleon.mag())));
//        pion->setEnergy(std::sqrt((pion->getMass())*(pion->getMass())+(mom_nucleon.mag()*mom_nucleon.mag())));
  G4double sh=nucleon->getEnergy()+pion->getEnergy();
        G4double mn=nucleon->getMass();
        G4double me=pion->getMass();
        G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
        nucleon->setEnergy(en);
        G4double ee=std::sqrt(en*en-mn*mn+me*me);
        pion->setEnergy(ee);
        G4double pn=std::sqrt(en*en-mn*mn);

        ThreeVector mom_nucleon = Random::normVector(pn);
        
        nucleon->setMomentum(mom_nucleon);
        pion->setMomentum(-mom_nucleon);
        
  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/G4INCLPiNToOmegaChannel.hh
• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLPiNToOmegaChannel.cc