G4INCL::SigmaZeroDecayChannel Class Reference

#include <G4INCLSigmaZeroDecayChannel.hh>

Inheritance diagram for G4INCL::SigmaZeroDecayChannel:

Public Member Functions
void	fillFinalState (FinalState *fs)
FinalState *	getFinalState ()
	SigmaZeroDecayChannel (Particle *, ThreeVector const &)
virtual	~SigmaZeroDecayChannel ()

Static Public Member Functions
static G4double	computeDecayTime (Particle *p)

Private Member Functions
	INCL_DECLARE_ALLOCATION_POOL (SigmaZeroDecayChannel)
void	sampleAngles (G4double *, G4double *, G4double *)

Private Attributes
ThreeVector const	incidentDirection
Particle *	theParticle

Detailed Description

Definition at line 49 of file G4INCLSigmaZeroDecayChannel.hh.

Constructor & Destructor Documentation

SigmaZeroDecayChannel()

G4INCL::SigmaZeroDecayChannel::SigmaZeroDecayChannel	(	Particle *	p,
		ThreeVector const &	dir
	)

Definition at line 46 of file G4INCLSigmaZeroDecayChannel.cc.

    :theParticle(p), incidentDirection(dir)
    {}

~SigmaZeroDecayChannel()

G4INCL::SigmaZeroDecayChannel::~SigmaZeroDecayChannel ( )

virtual

Definition at line 50 of file G4INCLSigmaZeroDecayChannel.cc.

{}

Member Function Documentation

computeDecayTime()

G4double G4INCL::SigmaZeroDecayChannel::computeDecayTime ( Particle *  p )

static

Definition at line 53 of file G4INCLSigmaZeroDecayChannel.cc.

                                                                {
// assert(p->getType() == SigmaZero);
        const G4double gamma = std::sqrt(1+std::pow(p->getMomentum().mag()/p->getMass(),2));
        const G4double tau = ParticleTable::getWidth(SigmaZero)*3E8*1E15*gamma; // fm
        const G4double t = -tau * std::log(Random::shoot());
        return t;
    }

References G4INCL::Particle::getMass(), G4INCL::Particle::getMomentum(), G4INCL::ParticleTable::getWidth(), G4INCL::ThreeVector::mag(), G4INCL::Random::shoot(), and G4INCL::SigmaZero.

Referenced by G4INCL::StandardPropagationModel::generateDecays().

fillFinalState()

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

virtual

Implements G4INCL::IChannel.

Definition at line 69 of file G4INCLSigmaZeroDecayChannel.cc.

                                                             {
        
// assert( theParticle->getType() == SigmaZero);
        ParticleType createdType = Photon;
        
        const G4double sqrtS = theParticle->getMass();
        
        theParticle->setType(Lambda);
        G4double phi, c_tet, s_tet;
        sampleAngles(&c_tet, &s_tet, &phi);
        
        G4double c_phi = std::cos(phi);
        G4double s_phi = std::sin(phi);
        G4double beta = incidentDirection.mag();
        
        G4double q1, q2, q3;
        G4double sal=0.0;
        if (beta >= 1.0e-10)
            sal = incidentDirection.perp()/beta;
        if (sal >= 1.0e-6) {
            G4double b1 = incidentDirection.getX();
            G4double b2 = incidentDirection.getY();
            G4double b3 = incidentDirection.getZ();
            G4double cal = b3/beta;
            G4double t1 = c_tet+cal*s_tet*s_phi/sal;
            G4double t2 = s_tet/sal;
            q1=(b1*t1+b2*t2*c_phi)/beta;
            q2=(b2*t1-b1*t2*c_phi)/beta;
            q3=(b3*t1/beta-t2*s_phi);
        } else {
            q1 = s_tet*c_phi;
            q2 = s_tet*s_phi;
            q3 = c_tet;
        }
                    
        G4double xq = KinematicsUtils::momentumInCM(sqrtS,
                                                  theParticle->getMass(),
                                                  ParticleTable::getINCLMass(createdType));
        q1 *= xq;
        q2 *= xq;
        q3 *= xq;
        
        ThreeVector createdMomentum(q1, q2, q3);
        ThreeVector createdPosition(theParticle->getPosition());
        Particle *createdParticle = new Particle(createdType, createdMomentum, createdPosition);
        theParticle->setMomentum(-createdMomentum);
        theParticle->adjustEnergyFromMomentum();
        
        fs->addModifiedParticle(theParticle);
        fs->addCreatedParticle(createdParticle);
        
    }

References G4INCL::FinalState::addCreatedParticle(), G4INCL::FinalState::addModifiedParticle(), G4INCL::Particle::adjustEnergyFromMomentum(), anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, G4INCL::ParticleTable::getINCLMass(), G4INCL::Particle::getMass(), G4INCL::Particle::getPosition(), G4INCL::ThreeVector::getX(), G4INCL::ThreeVector::getY(), G4INCL::ThreeVector::getZ(), incidentDirection, G4INCL::Lambda, G4INCL::ThreeVector::mag(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::ThreeVector::perp(), G4INCL::Photon, sampleAngles(), G4INCL::Particle::setMomentum(), G4INCL::Particle::setType(), and theParticle.

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::SigmaZeroDecayChannel::INCL_DECLARE_ALLOCATION_POOL ( SigmaZeroDecayChannel  )

private

sampleAngles()

void G4INCL::SigmaZeroDecayChannel::sampleAngles ( G4double *  ctet_par, G4double *  stet_par, G4double *  phi_par  )

private

Definition at line 61 of file G4INCLSigmaZeroDecayChannel.cc.

                                                                                                      {
 
        (*ctet_par) = -1.0 + 2.0*Random::shoot();
        if(std::abs(*ctet_par) > 1.0) (*ctet_par) = Math::sign(*ctet_par); // needed?
        (*stet_par) = std::sqrt(1.-(*ctet_par)*(*ctet_par));
        (*phi_par) = Math::twoPi * Random::shoot();
    }

References G4INCL::Random::shoot(), G4INCL::Math::sign(), and G4INCL::Math::twoPi.

Referenced by fillFinalState().

Field Documentation

incidentDirection

ThreeVector const G4INCL::SigmaZeroDecayChannel::incidentDirection

private

Definition at line 61 of file G4INCLSigmaZeroDecayChannel.hh.

Referenced by fillFinalState().

theParticle

Particle* G4INCL::SigmaZeroDecayChannel::theParticle

private

Definition at line 60 of file G4INCLSigmaZeroDecayChannel.hh.

Referenced by fillFinalState().

---

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

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