G4INCL::OmegaNToPiNChannel Class Reference

#include <G4INCLOmegaNToPiNChannel.hh>

Inheritance diagram for G4INCL::OmegaNToPiNChannel:

Public Member Functions
void	fillFinalState (FinalState *fs)
FinalState *	getFinalState ()
	OmegaNToPiNChannel (Particle *, Particle *)
virtual	~OmegaNToPiNChannel ()

Private Member Functions
	INCL_DECLARE_ALLOCATION_POOL (OmegaNToPiNChannel)

Private Attributes
Particle *	particle1
Particle *	particle2

Detailed Description

Definition at line 47 of file G4INCLOmegaNToPiNChannel.hh.

Constructor & Destructor Documentation

OmegaNToPiNChannel()

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

Definition at line 47 of file G4INCLOmegaNToPiNChannel.cc.

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

~OmegaNToPiNChannel()

G4INCL::OmegaNToPiNChannel::~OmegaNToPiNChannel ( )

virtual

Definition at line 53 of file G4INCLOmegaNToPiNChannel.cc.

                                           {
 
    }

Member Function Documentation

fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 57 of file G4INCLOmegaNToPiNChannel.cc.

                                                          {
        Particle * nucleon;
        Particle * omega;
        if(particle1->isNucleon()) {
            nucleon = particle1;
            omega = particle2;
        } else {
            nucleon = particle2;
            omega = particle1;
        }
 
        const G4double r2 = Random::shoot();
        if (nucleon->getType() == Neutron) {
            if (r2*3. < 2.) {
                nucleon->setType(Proton);
                omega->setType(PiMinus);
            }
            else {
                nucleon->setType(Neutron);
                omega->setType(PiZero);
            }
        }
        else {
            if (r2*3. < 2.) {
                nucleon->setType(Neutron);
                omega->setType(PiPlus);
            }
            else {
                nucleon->setType(Proton);
                omega->setType(PiZero);
            }
        }
        
        G4double sh=nucleon->getEnergy()+omega->getEnergy();
        G4double mn=nucleon->getMass();
        G4double me=omega->getMass();
        G4double en=(sh*sh+mn*mn-me*me)/(2*sh);
        nucleon->setEnergy(en);
        G4double ee=std::sqrt(en*en-mn*mn+me*me);
        omega->setEnergy(ee);
        G4double pn=std::sqrt(en*en-mn*mn);
/*// test isotropy
        const G4double pi=std::acos(-1.0);
        G4double x1;
        G4double u1;
        G4double fteta;
        G4double teta;
        G4double fi;
 
        G4int passe1=0;
        while (passe1==0) {
            // Sample x from 0 to pi/2
            x1=(pi/2.)*Random::shoot();
            // Sample u from 0 to 1
            u1=Random::shoot();
            fteta=std::sin(x1);
            // The condition
            if (u1 < fteta) {
                teta=x1;
                passe1=1;
                if (Random::shoot() < 0.5) {
                    teta=pi-teta;
                }
            }
        }
        fi=(2.0*pi)*Random::shoot();
        
// end test     
        ThreeVector mom_nucleon(
                                pn*std::sin(teta)*std::cos(fi),
                                pn*std::sin(teta)*std::sin(fi),
                                pn*std::cos(teta)
                                );
        if (Random::shoot() < 0.5) {
            nucleon->setMomentum(mom_nucleon);
            omega->setMomentum(-mom_nucleon);
        }
        else {
            omega->setMomentum(mom_nucleon);
            nucleon->setMomentum(-mom_nucleon);
        }
*/
        ThreeVector mom_nucleon = Random::normVector(pn);
        
        nucleon->setMomentum(mom_nucleon);
        omega->setMomentum(-mom_nucleon);
        
        fs->addModifiedParticle(nucleon);
        fs->addModifiedParticle(omega);
    }

References G4INCL::FinalState::addModifiedParticle(), G4INCL::Particle::getEnergy(), G4INCL::Particle::getMass(), G4INCL::Particle::isNucleon(), G4INCL::Neutron, G4INCL::Random::normVector(), G4InuclParticleNames::nucleon(), particle1, particle2, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4InuclParticleNames::pn, G4INCL::Proton, G4INCL::Particle::setEnergy(), G4INCL::Particle::setMomentum(), G4INCL::Particle::setType(), 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().

INCL_DECLARE_ALLOCATION_POOL()

G4INCL::OmegaNToPiNChannel::INCL_DECLARE_ALLOCATION_POOL ( OmegaNToPiNChannel  )

private

Field Documentation

particle1

Particle* G4INCL::OmegaNToPiNChannel::particle1

private

Definition at line 55 of file G4INCLOmegaNToPiNChannel.hh.

Referenced by fillFinalState().

particle2

Particle * G4INCL::OmegaNToPiNChannel::particle2

private

Definition at line 55 of file G4INCLOmegaNToPiNChannel.hh.

Referenced by fillFinalState().

---

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

• source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLOmegaNToPiNChannel.hh
• source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLOmegaNToPiNChannel.cc