G4SPBaryon Class Reference

#include <G4SPBaryon.hh>

Collaboration diagram for G4SPBaryon:

Public Member Functions
	G4SPBaryon (G4Proton *aProton)
	G4SPBaryon (G4Neutron *aNeutron)
	G4SPBaryon (G4Lambda *aLambda)
	G4SPBaryon (G4SigmaPlus *aSigmaPlus)
	G4SPBaryon (G4SigmaZero *aSigmaZero)
	G4SPBaryon (G4SigmaMinus *aSigmaMinus)
	G4SPBaryon (G4XiMinus *aXiMinus)
	G4SPBaryon (G4XiZero *aXiZero)
	G4SPBaryon (G4OmegaMinus *anOmegaMinus)
	G4SPBaryon (G4AntiProton *aAntiProton)
	G4SPBaryon (G4AntiNeutron *aAntiNeutron)
	G4SPBaryon (G4AntiLambda *aAntiLambda)
	G4SPBaryon (G4AntiSigmaPlus *aAntiSigmaPlus)
	G4SPBaryon (G4AntiSigmaZero *aAntiSigmaZero)
	G4SPBaryon (G4AntiSigmaMinus *aAntiSigmaMinus)
	G4SPBaryon (G4AntiXiMinus *aAntiXiMinus)
	G4SPBaryon (G4AntiXiZero *aAntiXiZero)
	G4SPBaryon (G4AntiOmegaMinus *anAntiOmegaMinus)
	G4SPBaryon (G4ParticleDefinition *aDefinition)
	~G4SPBaryon ()
G4ParticleDefinition *	GetDefinition ()
void	SampleQuarkAndDiquark (G4int &quark, G4int &diQuark) const
void	FindDiquark (G4int quark, G4int &diQuark) const
G4int	FindQuark (G4int diQuark) const
G4double	GetProbability (G4int diQuark) const
G4int	MatchDiQuarkAndGetQuark (const G4SPBaryon &aBaryon, G4int &aDiQuark) const

Private Member Functions
G4bool	operator== (const G4SPBaryon &aBaryon) const

Private Attributes
G4ParticleDefinition *	theDefinition
std::vector< G4SPPartonInfo * >	thePartonInfo

Detailed Description

Definition at line 56 of file G4SPBaryon.hh.

Constructor & Destructor Documentation

G4SPBaryon() [1/19]

G4SPBaryon::G4SPBaryon

(

G4Proton *

aProton

)

Definition at line 154 of file G4SPBaryon.cc.

{
    theDefinition = aProton;
    thePartonInfo.push_back(new G4SPPartonInfo(2203, 1, 1./3.)); // uu_1, d
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 2, 1./6.)); // ud_1, u
    thePartonInfo.push_back(new G4SPPartonInfo(2101, 2, 1./2.)); // ud_0, u
}

Here is the caller graph for this function:

G4SPBaryon() [2/19]

G4SPBaryon::G4SPBaryon

(

G4Neutron *

aNeutron

)

Definition at line 172 of file G4SPBaryon.cc.

{
    theDefinition = aNeutron;
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 1, 1./6.)); // ud_1, d
    thePartonInfo.push_back(new G4SPPartonInfo(2101, 1, 1./2.)); // ud_0, d
    thePartonInfo.push_back(new G4SPPartonInfo(1103, 2, 1./3.)); // dd_1, u
}

Here is the call graph for this function:

G4SPBaryon() [3/19]

G4SPBaryon::G4SPBaryon

(

G4Lambda *

aLambda

)

Definition at line 190 of file G4SPBaryon.cc.

{
    theDefinition = aLambda;
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 3, 1./3.));  // ud_1, s
    thePartonInfo.push_back(new G4SPPartonInfo(3203, 1, 1./4.));  // su_1, d
    thePartonInfo.push_back(new G4SPPartonInfo(3201, 1, 1./12.)); // su_0, d
    thePartonInfo.push_back(new G4SPPartonInfo(3103, 2, 1./4.));  // sd_1, u
    thePartonInfo.push_back(new G4SPPartonInfo(3101, 2, 1./12.)); // sd_0, u
}

Here is the call graph for this function:

G4SPBaryon() [4/19]

G4SPBaryon::G4SPBaryon

(

G4SigmaPlus *

aSigmaPlus

)

Definition at line 212 of file G4SPBaryon.cc.

{
    theDefinition = aSigmaPlus;
    thePartonInfo.push_back(new G4SPPartonInfo(2203, 3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3203, 2, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(3201, 2, 1./2.));
}

Here is the call graph for this function:

G4SPBaryon() [5/19]

G4SPBaryon::G4SPBaryon

(

G4SigmaZero *

aSigmaZero

)

Definition at line 230 of file G4SPBaryon.cc.

{
    theDefinition = aSigmaZero;
    thePartonInfo.push_back(new G4SPPartonInfo(2103, 3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3203, 1, 1./12.));
    thePartonInfo.push_back(new G4SPPartonInfo(3201, 1, 1./4.));
    thePartonInfo.push_back(new G4SPPartonInfo(3103, 2, 1./12.));
    thePartonInfo.push_back(new G4SPPartonInfo(3101, 2, 1./4.));
}

Here is the call graph for this function:

G4SPBaryon() [6/19]

G4SPBaryon::G4SPBaryon

(

G4SigmaMinus *

aSigmaMinus

)

Definition at line 252 of file G4SPBaryon.cc.

{
    theDefinition = aSigmaMinus;
    thePartonInfo.push_back(new G4SPPartonInfo(1103, 3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(3103, 1, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(3101, 1, 1./2.));
}

Here is the call graph for this function:

G4SPBaryon() [7/19]

G4SPBaryon::G4SPBaryon

(

G4XiMinus *

aXiMinus

)

Definition at line 270 of file G4SPBaryon.cc.

{
    theDefinition = aXiMinus;
    thePartonInfo.push_back(new G4SPPartonInfo(3103, 3, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(3101, 3, 1./2.));
    thePartonInfo.push_back(new G4SPPartonInfo(3303, 1, 1./3.));
}

Here is the call graph for this function:

G4SPBaryon() [8/19]

G4SPBaryon::G4SPBaryon

(

G4XiZero *

aXiZero

)

Definition at line 288 of file G4SPBaryon.cc.

{
    theDefinition = aXiZero;
    thePartonInfo.push_back(new G4SPPartonInfo(3203, 3, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(3201, 3, 1./2.));
    thePartonInfo.push_back(new G4SPPartonInfo(3303, 2, 1./3.));
}

Here is the call graph for this function:

G4SPBaryon() [9/19]

G4SPBaryon::G4SPBaryon

(

G4OmegaMinus *

anOmegaMinus

)

Definition at line 306 of file G4SPBaryon.cc.

{
    theDefinition = anOmegaMinus;
    thePartonInfo.push_back(new G4SPPartonInfo(3303, 3, 1.));
}

Here is the call graph for this function:

G4SPBaryon() [10/19]

G4SPBaryon::G4SPBaryon

(

G4AntiProton *

aAntiProton

)

Definition at line 163 of file G4SPBaryon.cc.

{
    theDefinition = aAntiProton;
    thePartonInfo.push_back(new G4SPPartonInfo(-2203, -1, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -2, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(-2101, -2, 1./2.));
}

Here is the call graph for this function:

G4SPBaryon() [11/19]

G4SPBaryon::G4SPBaryon

(

G4AntiNeutron *

aAntiNeutron

)

Definition at line 181 of file G4SPBaryon.cc.

{
    theDefinition = aAntiNeutron;
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -1, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(-2101, -1, 1./2.));
    thePartonInfo.push_back(new G4SPPartonInfo(-1103, -2, 1./3.));
}

Here is the call graph for this function:

G4SPBaryon() [12/19]

G4SPBaryon::G4SPBaryon

(

G4AntiLambda *

aAntiLambda

)

Definition at line 201 of file G4SPBaryon.cc.

{
    theDefinition = aAntiLambda;
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3203, -1, 1./4.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3201, -1, 1./12.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3103, -2, 1./4.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3101, -2, 1./12.));
}

Here is the call graph for this function:

G4SPBaryon() [13/19]

G4SPBaryon::G4SPBaryon

(

G4AntiSigmaPlus *

aAntiSigmaPlus

)

Definition at line 221 of file G4SPBaryon.cc.

{
    theDefinition = aAntiSigmaPlus;
    thePartonInfo.push_back(new G4SPPartonInfo(-2203, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3203, -2, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3201, -2, 1./2.));
}

Here is the call graph for this function:

G4SPBaryon() [14/19]

G4SPBaryon::G4SPBaryon

(

G4AntiSigmaZero *

aAntiSigmaZero

)

Definition at line 241 of file G4SPBaryon.cc.

{
    theDefinition = aAntiSigmaZero;
    thePartonInfo.push_back(new G4SPPartonInfo(-2103, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3203, -1, 1./12.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3201, -1, 1./4.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3103, -2, 1./12.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3101, -2, 1./4.));
}

Here is the call graph for this function:

G4SPBaryon() [15/19]

G4SPBaryon::G4SPBaryon

(

G4AntiSigmaMinus *

aAntiSigmaMinus

)

Definition at line 261 of file G4SPBaryon.cc.

{
    theDefinition = aAntiSigmaMinus;
    thePartonInfo.push_back(new G4SPPartonInfo(-1103, -3, 1./3.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3103, -1, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3101, -1, 1./2.));
}

Here is the call graph for this function:

G4SPBaryon() [16/19]

G4SPBaryon::G4SPBaryon

(

G4AntiXiMinus *

aAntiXiMinus

)

Definition at line 279 of file G4SPBaryon.cc.

{
    theDefinition = aAntiXiMinus;
    thePartonInfo.push_back(new G4SPPartonInfo(-3103, -3, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3101, -3, 1./2.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3303, -1, 1./3.));
}

Here is the call graph for this function:

G4SPBaryon() [17/19]

G4SPBaryon::G4SPBaryon

(

G4AntiXiZero *

aAntiXiZero

)

Definition at line 297 of file G4SPBaryon.cc.

{
    theDefinition = aAntiXiZero;
    thePartonInfo.push_back(new G4SPPartonInfo(-3203, -3, 1./6.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3201, -3, 1./2.));
    thePartonInfo.push_back(new G4SPPartonInfo(-3303, -2, 1./3.));
}

Here is the call graph for this function:

G4SPBaryon() [18/19]

G4SPBaryon::G4SPBaryon

(

G4AntiOmegaMinus *

anAntiOmegaMinus

)

Definition at line 314 of file G4SPBaryon.cc.

{
    theDefinition = anAntiOmegaMinus;
    thePartonInfo.push_back(new G4SPPartonInfo(-3303, -3, 1.));
}

Here is the call graph for this function:

G4SPBaryon() [19/19]

G4SPBaryon::G4SPBaryon

(

G4ParticleDefinition *

aDefinition

)

Definition at line 322 of file G4SPBaryon.cc.

{
    theDefinition = aDefinition;
    if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(2224))// D++
    {
        thePartonInfo.push_back(new G4SPPartonInfo(2203, 2, 1.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-2224))// anti D++
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-2203, -2, 1.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(2214))// D+
    {
        thePartonInfo.push_back(new G4SPPartonInfo(2203, 1, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(2103, 2, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-2214))// anti D+
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-2203, -1, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-2103, -2, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(2114))// D0
    {
        thePartonInfo.push_back(new G4SPPartonInfo(2103, 1, 2./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(2103, 2, 1./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-2114))// anti D0
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-2103, -1, 2./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-2103, -2, 1./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(1114))// D-
    {
        thePartonInfo.push_back(new G4SPPartonInfo(1103, 1, 1.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-1114))// anti D-
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-1103, -1, 1.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(3224))// S*+
    {
        thePartonInfo.push_back(new G4SPPartonInfo(2203, 3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(3203, 2, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-3224))// anti S*+
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-2203, -3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-3203, -2, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(3214))// S*0
    {
        thePartonInfo.push_back(new G4SPPartonInfo(2103, 3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(3203, 1, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(3103, 2, 1./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-3214))// anti S*0
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-2103, -3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-3203, -1, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-3103, -2, 1./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(3114))// S*-
    {
        thePartonInfo.push_back(new G4SPPartonInfo(1103, 3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(3103, 1, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-3224))// anti S*-
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-1103, -3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-3103, -1, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(3324))// Xi*0
    {
        thePartonInfo.push_back(new G4SPPartonInfo(3203, 3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(3303, 2, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-3324))// anti Xi*0
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-3203, -3, 1./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-3303, -2, 2./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(3314))// Xi*-
    {
        thePartonInfo.push_back(new G4SPPartonInfo(3103, 3, 2./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(3303, 1, 1./3.));
    }
    else if(theDefinition ==
            G4ParticleTable::GetParticleTable()->FindParticle(-3314))// anti Xi*-
    {
        thePartonInfo.push_back(new G4SPPartonInfo(-3103, -3, 2./3.));
        thePartonInfo.push_back(new G4SPPartonInfo(-3303, -1, 1./3.));
    }
}

Here is the call graph for this function:

~G4SPBaryon()

G4SPBaryon::~G4SPBaryon

(

)

Definition at line 433 of file G4SPBaryon.cc.

{
    for(unsigned int i=0;i<thePartonInfo.size(); i++) delete thePartonInfo[i];
}

Member Function Documentation

FindDiquark()

void G4SPBaryon::FindDiquark	(	G4int	quark,
		G4int &	diQuark
	)		const

Definition at line 124 of file G4SPBaryon.cc.

{
    typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
    G4double sum = 0;
    iter i;
    for(i=thePartonInfo.begin() ; i!=thePartonInfo.end(); i++)
    {
        if (std::abs((*i)->GetQuark()) == std::abs(quark))
        {
            sum += (*i)->GetProbability();
        }
    }
    G4double random = G4UniformRand();
    G4double running = 0;
    for(i=thePartonInfo.begin() ; i!=thePartonInfo.end(); i++)
    {
        if (std::abs((*i)->GetQuark()) == std::abs(quark))
        {
            running += (*i)->GetProbability();
            if (running/sum >= random)
            {
                diQuark = (*i)->GetDiQuark();
                break;
            }
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

FindQuark()

G4int G4SPBaryon::FindQuark

(

G4int

diQuark

)

const

Definition at line 32 of file G4SPBaryon.cc.

{
    G4double sum = GetProbability(diQuark);
    G4double random = G4UniformRand();
    G4double running = 0;
    G4int Quark(0);
    typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
    iter i;
    for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
    {
        if (std::abs((*i)->GetDiQuark()) == std::abs(diQuark))
        {
            running += (*i)->GetProbability();
            if (running/sum >= random)
            {
                Quark = (*i)->GetQuark();
                break;
            }
        }
    }
    return Quark;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetDefinition()

G4ParticleDefinition* G4SPBaryon::GetDefinition ( )

inline

Definition at line 86 of file G4SPBaryon.hh.

{return theDefinition;}

Here is the call graph for this function:

GetProbability()

G4double G4SPBaryon::GetProbability

(

G4int

diQuark

)

const

Definition at line 55 of file G4SPBaryon.cc.

{
    G4double sum = 0;
    typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
    iter i;
    for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
    {
        if (std::abs((*i)->GetDiQuark()) == std::abs(diQuark))
        {
            sum += (*i)->GetProbability();
        }
    }
    return sum;
}

Here is the call graph for this function:

Here is the caller graph for this function:

MatchDiQuarkAndGetQuark()

G4int G4SPBaryon::MatchDiQuarkAndGetQuark	(	const G4SPBaryon &	aBaryon,
		G4int &	aDiQuark
	)		const

Definition at line 71 of file G4SPBaryon.cc.

{
    G4int    result=0;
    typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
    iter i;
    G4double running = 0;
    G4double total = 0;
    for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
    {
        total += aBaryon.GetProbability((*i)->GetDiQuark());
    }
    G4double random = G4UniformRand()*total;
    for(i = thePartonInfo.begin(); i!=thePartonInfo.end(); i++)
    {
        running += aBaryon.GetProbability((*i)->GetDiQuark());
        if(random/total<running)
        {
            result = (*i)->GetQuark(); // (diquark annihilated)
            aDiQuark = (*i)->GetDiQuark();
            break;
        }
    }
    return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator==()

G4bool G4SPBaryon::operator== ( const G4SPBaryon &  aBaryon ) const

private

SampleQuarkAndDiquark()

void G4SPBaryon::SampleQuarkAndDiquark	(	G4int &	quark,
		G4int &	diQuark
	)		const

Definition at line 97 of file G4SPBaryon.cc.

{
    typedef std::vector<G4SPPartonInfo *>::const_iterator iter;
    G4double random = G4UniformRand();
    G4double sum = 0;
    iter i;
    for(i=thePartonInfo.begin() ; i!=thePartonInfo.end(); i++)
    {
        sum += (*i)->GetProbability();
        if (sum > random)
        {
            if (theDefinition->GetPDGEncoding() < 0)
            {
                quark = (*i)->GetDiQuark();
                diQuark = (*i)->GetQuark();
            }
            else
            {
                quark = (*i)->GetQuark();
                diQuark = (*i)->GetDiQuark();
            }
            break;
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

theDefinition

G4ParticleDefinition* G4SPBaryon::theDefinition

private

Definition at line 95 of file G4SPBaryon.hh.

thePartonInfo

std::vector<G4SPPartonInfo *> G4SPBaryon::thePartonInfo

private

Definition at line 96 of file G4SPBaryon.hh.

---

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

• G4SPBaryon.hh
• G4SPBaryon.cc