G4QGSMSplitableHadron Class Reference

#include <G4QGSMSplitableHadron.hh>

Inheritance diagram for G4QGSMSplitableHadron:

Collaboration diagram for G4QGSMSplitableHadron:

Public Member Functions
	G4QGSMSplitableHadron ()
	G4QGSMSplitableHadron (const G4ReactionProduct &aPrimary)
	G4QGSMSplitableHadron (const G4ReactionProduct &aPrimary, G4bool Direction)
	G4QGSMSplitableHadron (const G4Nucleon &aNucleon)
	G4QGSMSplitableHadron (const G4Nucleon &aNucleon, G4bool Direction)
virtual	~G4QGSMSplitableHadron ()
virtual void	SplitUp ()
virtual void	SetFirstParton (G4int PDGcode)
virtual void	SetSecondParton (G4int PDGcode)
virtual G4Parton *	GetNextParton ()
virtual G4Parton *	GetNextAntiParton ()
Public Member Functions inherited from G4VSplitableHadron
	G4VSplitableHadron ()
	G4VSplitableHadron (const G4ReactionProduct &aPrimary)
	G4VSplitableHadron (const G4Nucleon &aNucleon)
	G4VSplitableHadron (const G4VKineticNucleon *aNucleon)
virtual	~G4VSplitableHadron ()
int	operator== (const G4VSplitableHadron &right) const
int	operator!= (const G4VSplitableHadron &right) const
void	Set4Momentum (const G4LorentzVector &a4Momentum)
const G4LorentzVector &	Get4Momentum () const
void	SetDefinition (const G4ParticleDefinition *aDefinition)
const G4ParticleDefinition *	GetDefinition () const
void	IncrementCollisionCount (G4int aCount)
void	SetCollisionCount (G4int aCount)
void	SetTimeOfCreation (G4double aTime)
G4double	GetTimeOfCreation ()
void	SetPosition (const G4ThreeVector &aPosition)
const G4ThreeVector &	GetPosition () const
void	SetStatus (const G4int aStatus)
G4int	GetStatus ()
G4bool	IsSplit ()
G4int	GetSoftCollisionCount ()

Private Member Functions
const G4QGSMSplitableHadron &	operator= (const G4QGSMSplitableHadron &right)
void	InitParameters ()
void	DiffractiveSplitUp ()
void	SoftSplitUp ()
G4ThreeVector	GaussianPt (G4double widthSquare, G4double maxPtSquare)
void	GetValenceQuarkFlavors (const G4ParticleDefinition *aPart, G4Parton *&Parton1, G4Parton *&Parton2)
G4Parton *	BuildSeaQuark (G4bool isAntiQuark, G4int aPDGCode, G4int nSeaPair)
G4double	SampleX (G4double anXmin, G4int nSea, G4int theTotalSea, G4double aBeta)

Private Attributes
G4bool	Direction
std::deque< G4Parton * >	Color
std::deque< G4Parton * >	AntiColor
G4MesonSplitter	theMesonSplitter
G4BaryonSplitter	theBaryonSplitter
G4double	alpha
G4double	beta
G4double	theMinPz
G4double	StrangeSuppress
G4double	sigmaPt
G4double	widthOfPtSquare
G4double	minTransverseMass

Additional Inherited Members
Protected Member Functions inherited from G4VSplitableHadron
void	Splitting ()

Detailed Description

Definition at line 41 of file G4QGSMSplitableHadron.hh.

Constructor & Destructor Documentation

G4QGSMSplitableHadron() [1/5]

G4QGSMSplitableHadron::G4QGSMSplitableHadron

(

)

Definition at line 76 of file G4QGSMSplitableHadron.cc.

{
    InitParameters();
}

Here is the call graph for this function:

G4QGSMSplitableHadron() [2/5]

G4QGSMSplitableHadron::G4QGSMSplitableHadron

(

const G4ReactionProduct &

aPrimary

)

Definition at line 89 of file G4QGSMSplitableHadron.cc.

:  G4VSplitableHadron(aPrimary)
{
    InitParameters();
}

Here is the call graph for this function:

G4QGSMSplitableHadron() [3/5]

G4QGSMSplitableHadron::G4QGSMSplitableHadron	(	const G4ReactionProduct &	aPrimary,
		G4bool	Direction
	)

Definition at line 81 of file G4QGSMSplitableHadron.cc.

:G4VSplitableHadron(aPrimary)
{
    InitParameters();
    Direction = aDirection;
}

Here is the call graph for this function:

G4QGSMSplitableHadron() [4/5]

G4QGSMSplitableHadron::G4QGSMSplitableHadron

(

const G4Nucleon &

aNucleon

)

Definition at line 95 of file G4QGSMSplitableHadron.cc.

:  G4VSplitableHadron(aNucleon)
{
    InitParameters();
}

Here is the call graph for this function:

G4QGSMSplitableHadron() [5/5]

G4QGSMSplitableHadron::G4QGSMSplitableHadron	(	const G4Nucleon &	aNucleon,
		G4bool	Direction
	)

Definition at line 101 of file G4QGSMSplitableHadron.cc.

:  G4VSplitableHadron(aNucleon)
{
    InitParameters();
    Direction = aDirection;
}

Here is the call graph for this function:

~G4QGSMSplitableHadron()

G4QGSMSplitableHadron::~G4QGSMSplitableHadron ( )

virtual

Definition at line 108 of file G4QGSMSplitableHadron.cc.

{}

Member Function Documentation

BuildSeaQuark()

G4Parton * G4QGSMSplitableHadron::BuildSeaQuark ( G4bool  isAntiQuark, G4int  aPDGCode, G4int  nSeaPair  )

private

Definition at line 398 of file G4QGSMSplitableHadron.cc.

{
    if (isAntiQuark) aPDGCode*=-1;
    G4Parton* result = new G4Parton(aPDGCode);
    result->SetPosition(GetPosition());
    G4ThreeVector aPtVector = GaussianPt(sigmaPt, DBL_MAX);
    G4LorentzVector a4Momentum(aPtVector, 0);
    result->Set4Momentum(a4Momentum);
    return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DiffractiveSplitUp()

void G4QGSMSplitableHadron::DiffractiveSplitUp ( )

private

Definition at line 129 of file G4QGSMSplitableHadron.cc.

{   
    // take the particle definitions and get the partons HPW
    G4Parton * Left = NULL;
    G4Parton * Right = NULL;
    GetValenceQuarkFlavors(GetDefinition(), Left, Right);
    Left->SetPosition(GetPosition());
    Right->SetPosition(GetPosition());

    G4LorentzVector HadronMom = Get4Momentum();
    //std::cout << "DSU 1 - "<<HadronMom<<std::endl;

    // momenta of string ends
    G4double pt2 = HadronMom.perp2();
    G4double transverseMass2 = HadronMom.plus()*HadronMom.minus();
    G4double maxAvailMomentum2 = sqr(std::sqrt(transverseMass2) - std::sqrt(pt2));
    G4ThreeVector pt(minTransverseMass, minTransverseMass, 0);
    if(maxAvailMomentum2/widthOfPtSquare>0.01) pt = GaussianPt(widthOfPtSquare, maxAvailMomentum2);
    //std::cout << "DSU 1.1 - "<< maxAvailMomentum2<< pt <<std::endl;

    G4LorentzVector LeftMom(pt, 0.);
    G4LorentzVector RightMom;
    RightMom.setPx(HadronMom.px() - pt.x());
    RightMom.setPy(HadronMom.py() - pt.y());
    //std::cout << "DSU 2 - "<<RightMom<<" "<< LeftMom <<std::endl;

    G4double Local1 = HadronMom.minus() + (RightMom.perp2() - LeftMom.perp2())/HadronMom.plus();
    G4double Local2 = std::sqrt(std::max(0., sqr(Local1) - 4.*RightMom.perp2()*HadronMom.minus()/HadronMom.plus()));
    //std::cout << "DSU 3 - "<< Local1 <<" "<< Local2 <<std::endl;
    if (Direction) Local2 = -Local2;
    G4double RightMinus   = 0.5*(Local1 + Local2);
    G4double LeftMinus = HadronMom.minus() - RightMinus;
    //std::cout << "DSU 4 - "<< RightMinus <<" "<< LeftMinus << " "<<HadronMom.minus() <<std::endl;

    G4double LeftPlus  = LeftMom.perp2()/LeftMinus;
    G4double RightPlus = HadronMom.plus() - LeftPlus;
    //std::cout << "DSU 5 - "<< RightPlus <<" "<< LeftPlus <<std::endl;
    LeftMom.setPz(0.5*(LeftPlus - LeftMinus));
    LeftMom.setE (0.5*(LeftPlus + LeftMinus));
    RightMom.setPz(0.5*(RightPlus - RightMinus));
    RightMom.setE (0.5*(RightPlus + RightMinus));
    //std::cout << "DSU 6 - "<< LeftMom <<" "<< RightMom <<std::endl;
    Left->Set4Momentum(LeftMom);
    Right->Set4Momentum(RightMom);
    Color.push_back(Left);
    AntiColor.push_back(Right);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GaussianPt()

G4ThreeVector G4QGSMSplitableHadron::GaussianPt ( G4double  widthSquare, G4double  maxPtSquare  )

private

Definition at line 384 of file G4QGSMSplitableHadron.cc.

{
    G4double R;
        const G4int maxNumberOfLoops = 1000;
        G4int loopCounter = -1;
    while ( ((R = -widthSquare*G4Log(G4UniformRand())) > maxPtSquare) &&  /* Loop checking, 26.10.2015, A.Ribon */
                ++loopCounter < maxNumberOfLoops ) {;}
        if ( loopCounter >= maxNumberOfLoops ) R = 0.0;
    R = std::sqrt(R);
    G4double phi = twopi*G4UniformRand();
    return G4ThreeVector (R*std::cos(phi), R*std::sin(phi), 0.);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetNextAntiParton()

G4Parton * G4QGSMSplitableHadron::GetNextAntiParton ( )

inlinevirtual

Implements G4VSplitableHadron.

Definition at line 103 of file G4QGSMSplitableHadron.hh.

{
    if(AntiColor.size() == 0) return 0;
    G4Parton * result = AntiColor.front();
    AntiColor.pop_front();
    return result;
}

GetNextParton()

G4Parton * G4QGSMSplitableHadron::GetNextParton ( )

inlinevirtual

Implements G4VSplitableHadron.

Definition at line 95 of file G4QGSMSplitableHadron.hh.

{
    if(Color.size()==0) return 0;
    G4Parton * result = Color.back();
    Color.pop_back();
    return result;
}

GetValenceQuarkFlavors()

void G4QGSMSplitableHadron::GetValenceQuarkFlavors ( const G4ParticleDefinition *  aPart, G4Parton *&  Parton1, G4Parton *&  Parton2  )

private

Definition at line 321 of file G4QGSMSplitableHadron.cc.

{
    // Note! convention aEnd = q or (qq)bar and bEnd = qbar or qq.
    G4int aEnd;
    G4int bEnd;
    G4int HadronEncoding = aPart->GetPDGEncoding();
    if (aPart->GetBaryonNumber() == 0)
    {
        theMesonSplitter.SplitMeson(HadronEncoding, &aEnd, &bEnd);
    }
    else
    {
        theBaryonSplitter.SplitBarion(HadronEncoding, &aEnd, &bEnd);
    }

    Parton1 = new G4Parton(aEnd);
    Parton1->SetPosition(GetPosition());

    //  G4cerr << "G4QGSMSplitableHadron::GetValenceQuarkFlavors()" << G4endl;
    //  G4cerr << "Parton 1: "
    //         << " PDGcode: "  << aEnd
    //         << " - Name: "   << Parton1->GetDefinition()->GetParticleName()
    //         << " - Type: "   << Parton1->GetDefinition()->GetParticleType()
    //         << " - Spin-3: " << Parton1->GetSpinZ()
    //         << " - Colour: " << Parton1->GetColour() << G4endl;

    Parton2 = new G4Parton(bEnd);
    Parton2->SetPosition(GetPosition());

    //  G4cerr << "Parton 2: "
    //         << " PDGcode: "  << bEnd
    //         << " - Name: "   << Parton2->GetDefinition()->GetParticleName()
    //         << " - Type: "   << Parton2->GetDefinition()->GetParticleType()
    //         << " - Spin-3: " << Parton2->GetSpinZ()
    //         << " - Colour: " << Parton2->GetColour() << G4endl;
    //  G4cerr << "... now checking for color and spin conservation - yielding: " << G4endl;

    // colour of parton 1 choosen at random by G4Parton(aEnd)
    // colour of parton 2 is the opposite:

    Parton2->SetColour(-(Parton1->GetColour()));

    // isospin-3 of both partons is handled by G4Parton(PDGCode)

    // spin-3 of parton 1 and 2 choosen at random by G4Parton(aEnd)
    // spin-3 of parton 2 may be constrained by spin of original particle:

    if ( std::abs(Parton1->GetSpinZ() + Parton2->GetSpinZ()) > aPart->GetPDGSpin())
    {
        Parton2->SetSpinZ(-(Parton2->GetSpinZ()));    
    }

    //  G4cerr << "Parton 2: "
    //         << " PDGcode: "  << bEnd
    //         << " - Name: "   << Parton2->GetDefinition()->GetParticleName()
    //         << " - Type: "   << Parton2->GetDefinition()->GetParticleType()
    //         << " - Spin-3: " << Parton2->GetSpinZ()
    //         << " - Colour: " << Parton2->GetColour() << G4endl;
    //  G4cerr << "------------" << G4endl;

}

Here is the call graph for this function:

Here is the caller graph for this function:

InitParameters()

void G4QGSMSplitableHadron::InitParameters ( )

private

Definition at line 55 of file G4QGSMSplitableHadron.cc.

{
    // changing rapidity distribution for all
    alpha = -0.5; // Note that this number is still assumed in the algorithm
    // needs to be generalized.
    // changing rapidity distribution for projectile like
    beta = 2.5;// Note that this number is still assumed in the algorithm
    // needs to be generalized.
    theMinPz = 0.5*G4PionMinus::PionMinus()->GetPDGMass();
    //  theMinPz = 0.1*G4PionMinus::PionMinus()->GetPDGMass();
    //  theMinPz = G4PionMinus::PionMinus()->GetPDGMass();
    // as low as possible, otherwise, we have unphysical boundary conditions in the sampling.
    StrangeSuppress = 0.48;
    sigmaPt = 0.*GeV; // widens eta slightly, if increased to 1.7,
    // but Maxim's algorithm breaks energy conservation
    // to be revised.
    widthOfPtSquare = 0.01*GeV*GeV;
    Direction = FALSE;
    minTransverseMass = 1*keV;
} 

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

const G4QGSMSplitableHadron& G4QGSMSplitableHadron::operator= ( const G4QGSMSplitableHadron &  right )

private

SampleX()

G4double G4QGSMSplitableHadron::SampleX ( G4double  anXmin, G4int  nSea, G4int  theTotalSea, G4double  aBeta  )

private

Definition at line 410 of file G4QGSMSplitableHadron.cc.

{
    G4double result;
    G4double x1, x2;
    G4double ymax = 0;
    for(G4int ii=1; ii<100; ii++)
    {
        G4double y = G4Pow::GetInstance()->powA(1./G4double(ii), alpha);
        y *= G4Pow::GetInstance()->powN( G4Pow::GetInstance()->powA(1-anXmin-totalSea*anXmin, alpha+1) - G4Pow::GetInstance()->powA(anXmin, alpha+1), nSea);
        y *= G4Pow::GetInstance()->powA(1-anXmin-totalSea*anXmin, aBeta+1) - G4Pow::GetInstance()->powA(anXmin, aBeta+1);
        if(y>ymax) ymax = y;
    }
    G4double y;
    G4double xMax=1-(totalSea+1)*anXmin;
    if(anXmin > xMax)
    {
        G4cout << "anXmin = "<<anXmin<<" nSea = "<<nSea<<" totalSea = "<< totalSea<<G4endl;
        throw G4HadronicException(__FILE__, __LINE__, "G4QGSMSplitableHadron - Fatal: Cannot sample parton densities under these constraints.");
    }
        const G4int maxNumberOfLoops = 10000;
        G4int loopCounter = -1;
    do
    {
        x1 = G4RandFlat::shoot(anXmin, xMax);
        y = G4Pow::GetInstance()->powA(x1, alpha);
        y *= G4Pow::GetInstance()->powN( G4Pow::GetInstance()->powA(1-x1-totalSea*anXmin, alpha+1) - G4Pow::GetInstance()->powA(anXmin, alpha+1), nSea);
        y *= G4Pow::GetInstance()->powA(1-x1-totalSea*anXmin, aBeta+1) - G4Pow::GetInstance()->powA(anXmin, aBeta+1);
        x2 = ymax*G4UniformRand();
    }
    while ( (x2>y) && ++loopCounter < maxNumberOfLoops );  /* Loop checking, 26.10.2015, A.Ribon */
        if ( loopCounter >= maxNumberOfLoops ) {
          G4ExceptionDescription ed;
          ed << " Failed sampling after maxNumberOfLoops attempts : forced exit! " << G4endl;
          G4Exception( "G4QGSMSplitableHadron::SampleX ", "HAD_QGS_002", JustWarning, ed );
        }
    result = x1;
    return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetFirstParton()

void G4QGSMSplitableHadron::SetFirstParton ( G4int  PDGcode )

inlinevirtual

Implements G4VSplitableHadron.

Definition at line 111 of file G4QGSMSplitableHadron.hh.

{PDGcode++;}

SetSecondParton()

void G4QGSMSplitableHadron::SetSecondParton ( G4int  PDGcode )

inlinevirtual

Implements G4VSplitableHadron.

Definition at line 113 of file G4QGSMSplitableHadron.hh.

{PDGcode++;}

SoftSplitUp()

void G4QGSMSplitableHadron::SoftSplitUp ( )

private

Definition at line 178 of file G4QGSMSplitableHadron.cc.

{
    //... sample transversal momenta for sea and valence quarks
    G4double phi, pts;
    G4double SumPy = 0.;
    G4double SumPx = 0.;
    G4ThreeVector Pos    = GetPosition();
    G4int nSeaPair = GetSoftCollisionCount()-1;

    // here the condition,to ensure viability of splitting, also in cases
    // where difractive excitation occured together with soft scattering.
    //   G4double LightConeMomentum = (Direction)? Get4Momentum().plus() : Get4Momentum().minus();
    //   G4double Xmin = theMinPz/LightConeMomentum;
    G4double Xmin = theMinPz/( Get4Momentum().e() - GetDefinition()->GetPDGMass() );
    while(Xmin>=1-(2*nSeaPair+1)*Xmin) Xmin*=0.95;  /* Loop checking, 26.10.2015, A.Ribon */

    G4int aSeaPair;
    for (aSeaPair = 0; aSeaPair < nSeaPair; aSeaPair++)
    {
        //  choose quark flavour, d:u:s = 1:1:(1/StrangeSuppress-2)

        G4int aPDGCode = 1 + (G4int)(G4UniformRand()/StrangeSuppress);

        //  BuildSeaQuark() determines quark spin, isospin and colour
        //  via parton-constructor G4Parton(aPDGCode)

        G4Parton * aParton = BuildSeaQuark(false, aPDGCode, nSeaPair);

        //      G4cerr << "G4QGSMSplitableHadron::SoftSplitUp()" << G4endl;

        //      G4cerr << "Parton 1: "
        //             << " PDGcode: "  << aPDGCode
        //             << " - Name: "   << aParton->GetDefinition()->GetParticleName()
        //             << " - Type: "   << aParton->GetDefinition()->GetParticleType()
        //             << " - Spin-3: " << aParton->GetSpinZ()
        //             << " - Colour: " << aParton->GetColour() << G4endl;

        // save colour a spin-3 for anti-quark

        G4int firstPartonColour = aParton->GetColour();
        G4double firstPartonSpinZ = aParton->GetSpinZ();

        SumPx += aParton->Get4Momentum().px();
        SumPy += aParton->Get4Momentum().py();
        Color.push_back(aParton);

        // create anti-quark

        aParton = BuildSeaQuark(true, aPDGCode, nSeaPair);
        aParton->SetSpinZ(-firstPartonSpinZ);
        aParton->SetColour(-firstPartonColour);

        //      G4cerr << "Parton 2: "
        //             << " PDGcode: "  << -aPDGCode
        //             << " - Name: "   << aParton->GetDefinition()->GetParticleName()
        //             << " - Type: "   << aParton->GetDefinition()->GetParticleType()
        //             << " - Spin-3: " << aParton->GetSpinZ()
        //             << " - Colour: " << aParton->GetColour() << G4endl;
        //      G4cerr << "------------" << G4endl;

        SumPx += aParton->Get4Momentum().px();
        SumPy += aParton->Get4Momentum().py();
        AntiColor.push_back(aParton);
    }
    // Valence quark
    G4Parton* pColorParton = NULL;
    G4Parton* pAntiColorParton = NULL;
    GetValenceQuarkFlavors(GetDefinition(), pColorParton, pAntiColorParton);
    G4int ColorEncoding = pColorParton->GetPDGcode();

    pts   =  sigmaPt*std::sqrt(-G4Log(G4UniformRand()));
    phi   = 2.*pi*G4UniformRand();
    G4double Px = pts*std::cos(phi);
    G4double Py = pts*std::sin(phi);
    SumPx += Px;
    SumPy += Py;

    if (ColorEncoding < 0) // use particle definition
    {
        G4LorentzVector ColorMom(-SumPx, -SumPy, 0, 0);
        pColorParton->Set4Momentum(ColorMom);
        G4LorentzVector AntiColorMom(Px, Py, 0, 0);
        pAntiColorParton->Set4Momentum(AntiColorMom);
    }
    else
    {
        G4LorentzVector ColorMom(Px, Py, 0, 0);
        pColorParton->Set4Momentum(ColorMom);
        G4LorentzVector AntiColorMom(-SumPx, -SumPy, 0, 0);
        pAntiColorParton->Set4Momentum(AntiColorMom);
    }
    Color.push_back(pColorParton);
    AntiColor.push_back(pAntiColorParton);

    // Sample X
    G4int nAttempt = 0;
    G4double SumX = 0;
    G4double aBeta = beta;
    G4double ColorX, AntiColorX;
    if (GetDefinition() == G4PionMinus::PionMinusDefinition()) aBeta = 1.;
    if (GetDefinition() == G4Gamma::GammaDefinition()) aBeta = 1.;
    if (GetDefinition() == G4PionPlus::PionPlusDefinition()) aBeta = 1.;
    if (GetDefinition() == G4PionZero::PionZeroDefinition()) aBeta = 1.;
    if (GetDefinition() == G4KaonPlus::KaonPlusDefinition()) aBeta = 0.;
    if (GetDefinition() == G4KaonMinus::KaonMinusDefinition()) aBeta = 0.;
        const G4int maxNumberOfAttempts = 1000;
    do
    {
        SumX = 0;
        nAttempt++;
        G4int    NumberOfUnsampledSeaQuarks = 2*nSeaPair;
        ColorX = SampleX(Xmin, NumberOfUnsampledSeaQuarks, 2*nSeaPair, aBeta);
        Color.back()->SetX(SumX = ColorX);// this is the valenz quark.
        for(G4int aPair = 0; aPair < nSeaPair; aPair++)
        {
            NumberOfUnsampledSeaQuarks--;
            ColorX = SampleX(Xmin, NumberOfUnsampledSeaQuarks, 2*nSeaPair, aBeta);
            Color[aPair]->SetX(ColorX);
            SumX += ColorX;
            NumberOfUnsampledSeaQuarks--;
            AntiColorX = SampleX(Xmin, NumberOfUnsampledSeaQuarks, 2*nSeaPair, aBeta);
            AntiColor[aPair]->SetX(AntiColorX); // the 'sea' partons
            SumX += AntiColorX;
            if (1. - SumX <= Xmin)  break;
        }
    }
    while ( (1. - SumX <= Xmin) && nAttempt < maxNumberOfAttempts );  /* Loop checking, 26.10.2015, A.Ribon */   
        if ( nAttempt >= maxNumberOfAttempts ) return;

    (*(AntiColor.end()-1))->SetX(1. - SumX); // the di-quark takes the rest, then go to momentum
    G4double lightCone  = ((!Direction) ? Get4Momentum().minus() : Get4Momentum().plus());
    G4double lightCone2 = ((!Direction) ? Get4Momentum().plus() : Get4Momentum().minus());
    for(aSeaPair = 0; aSeaPair < nSeaPair+1; aSeaPair++)
    {
        G4Parton* aParton = Color[aSeaPair];
        aParton->DefineMomentumInZ(lightCone, lightCone2, Direction);

        aParton = AntiColor[aSeaPair];
        aParton->DefineMomentumInZ(lightCone, lightCone2, Direction);
    }
    return;
} 

Here is the call graph for this function:

Here is the caller graph for this function:

SplitUp()

void G4QGSMSplitableHadron::SplitUp ( )

virtual

Implements G4VSplitableHadron.

Definition at line 114 of file G4QGSMSplitableHadron.cc.

{  
    if (IsSplit()) return;
    Splitting();
    if (Color.size()!=0) return;
    if (GetSoftCollisionCount() == 0)
    {
        DiffractiveSplitUp();
    }
    else
    {
        SoftSplitUp();
    }
}

Here is the call graph for this function:

Member Data Documentation

alpha

G4double G4QGSMSplitableHadron::alpha

private

Definition at line 86 of file G4QGSMSplitableHadron.hh.

AntiColor

std::deque<G4Parton *> G4QGSMSplitableHadron::AntiColor

private

Definition at line 78 of file G4QGSMSplitableHadron.hh.

beta

G4double G4QGSMSplitableHadron::beta

private

Definition at line 87 of file G4QGSMSplitableHadron.hh.

Color

std::deque<G4Parton *> G4QGSMSplitableHadron::Color

private

Definition at line 77 of file G4QGSMSplitableHadron.hh.

Direction

G4bool G4QGSMSplitableHadron::Direction

private

Definition at line 75 of file G4QGSMSplitableHadron.hh.

minTransverseMass

G4double G4QGSMSplitableHadron::minTransverseMass

private

Definition at line 92 of file G4QGSMSplitableHadron.hh.

sigmaPt

G4double G4QGSMSplitableHadron::sigmaPt

private

Definition at line 90 of file G4QGSMSplitableHadron.hh.

StrangeSuppress

G4double G4QGSMSplitableHadron::StrangeSuppress

private

Definition at line 89 of file G4QGSMSplitableHadron.hh.

theBaryonSplitter

G4BaryonSplitter G4QGSMSplitableHadron::theBaryonSplitter

private

Definition at line 82 of file G4QGSMSplitableHadron.hh.

theMesonSplitter

G4MesonSplitter G4QGSMSplitableHadron::theMesonSplitter

private

Definition at line 81 of file G4QGSMSplitableHadron.hh.

theMinPz

G4double G4QGSMSplitableHadron::theMinPz

private

Definition at line 88 of file G4QGSMSplitableHadron.hh.

widthOfPtSquare

G4double G4QGSMSplitableHadron::widthOfPtSquare

private

Definition at line 91 of file G4QGSMSplitableHadron.hh.

---

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

• G4QGSMSplitableHadron.hh
• G4QGSMSplitableHadron.cc