G4Parton Class Reference

#include <G4Parton.hh>

Public Member Functions
	G4Parton ()
	G4Parton (G4int PDGencoding)
	G4Parton (const G4Parton &right)
	~G4Parton ()
G4Parton &	operator= (const G4Parton &right)
int	operator== (const G4Parton &right) const
int	operator!= (const G4Parton &right) const
G4int	GetPDGcode () const
G4ParticleDefinition *	GetDefinition ()
void	DefineMomentumInZ (G4double aLightConeMomentum, G4bool aDirection)
void	DefineMomentumInZ (G4double aLightConeMomentum, G4double aLightConeE, G4bool aDirection)
const G4ThreeVector &	GetPosition () const
void	SetPosition (const G4ThreeVector &aPosition)
const G4LorentzVector &	Get4Momentum () const
void	Set4Momentum (const G4LorentzVector &aMomentum)
void	SetX (G4double anX)
G4double	GetX ()
void	SetColour (G4int aColour)
G4int	GetColour ()
void	SetIsoSpinZ (G4double anIsoSpinZ)
G4double	GetIsoSpinZ ()
void	SetSpinZ (G4double aSpinZ)
G4double	GetSpinZ ()

Detailed Description

Definition at line 47 of file G4Parton.hh.

Constructor & Destructor Documentation

G4Parton::G4Parton ( )

inline

Definition at line 50 of file G4Parton.hh.

References G4UniformRand.

      {
        // CAUTION: 
        // this is a preliminary definition yielding u and d quarks only!
        //
        PDGencoding=(G4int)(2.*G4UniformRand()); 
        theColour = (G4int)(3.*G4UniformRand())+1;
        theIsoSpinZ = ((G4int)(G4UniformRand()))-0.5;
        theSpinZ = ((G4int)(G4UniformRand()))-0.5;
      }

G4Parton::G4Parton

(

G4int

PDGencoding

)

Definition at line 39 of file G4Parton.cc.

References G4ParticleTable::FindParticle(), G4cout, G4endl, G4UniformRand, G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetParticleType(), G4ParticleDefinition::GetPDGiIsospin(), G4ParticleDefinition::GetPDGIsospin3(), and G4ParticleDefinition::GetPDGiSpin().

{
    PDGencoding=PDGcode;
    theX = 0;
    theDefinition=G4ParticleTable::GetParticleTable()->FindParticle(PDGencoding);
    if (theDefinition == NULL)
    {
      G4cout << "Encoding = "<<PDGencoding<<G4endl;
      G4String text = "G4Parton::GetDefinition(): Encoding not in particle table";
      throw G4HadronicException(__FILE__, __LINE__, text);
    }
    //
    // colour by random in (1,2,3)=(R,G,B) for quarks and 
    //                  in (-1,-2,-3)=(Rbar,Gbar,Bbar) for anti-quarks:
  //
    if (theDefinition->GetParticleType() == "quarks") {
        theColour = ((G4int)(3.*G4UniformRand())+1)*(std::abs(PDGencoding)/PDGencoding) ;
    }
    //
    // colour by random in (-1,-2,-3)=(Rbar,Gbar,Bbar)=(GB,RB,RG) for di-quarks and
    //                  in (1,2,3)=(R,G,B)=(GB,RB,RG) for anti-di-quarks:
  //
    else if (theDefinition->GetParticleType() == "diquarks") {
        theColour = -((G4int)(3.*G4UniformRand())+1)*(std::abs(PDGencoding)/PDGencoding);
    }
    //
    // colour by random in (-11,-12,...,-33)=(RRbar,RGbar,RBbar,...,BBbar) for gluons:
  //
    else if (theDefinition->GetParticleType() == "gluons") {
        theColour = -(((G4int)(3.*G4UniformRand())+1)*10 + ((G4int)(3.*G4UniformRand())+1));
    }
    else {
      G4cout << "Encoding = "<<PDGencoding<<G4endl;
      G4String text = "G4Parton::GetDefinition(): Particle is not a parton";
      throw G4HadronicException(__FILE__, __LINE__, text);
    }
    //  
    // isospin-z from PDG-encoded isospin-z for 
    // quarks, anti-quarks, di-quarks, and anti-di-quarks:
  //
    if ((theDefinition->GetParticleType() == "quarks") || (theDefinition->GetParticleType() == "diquarks")){
        theIsoSpinZ = theDefinition->GetPDGIsospin3();
    }
    //
  // isospin-z choosen at random from PDG-encoded isospin for gluons (should be zero):
    //
    else {
        G4int thisPDGiIsospin=theDefinition->GetPDGiIsospin();
        if (thisPDGiIsospin == 0) {
            theIsoSpinZ = 0;
        }
        else {
            theIsoSpinZ = ((G4int)((thisPDGiIsospin+1)*G4UniformRand()))-thisPDGiIsospin*0.5;
        }
    }
    //
    // spin-z choosen at random from PDG-encoded spin:
    //
    G4int thisPDGiSpin=theDefinition->GetPDGiSpin();
    if (thisPDGiSpin == 0) {
        theSpinZ = 0;
    }
    else {
        G4int rand=((G4int)((thisPDGiSpin+1)*G4UniformRand()));
        theSpinZ = rand-thisPDGiSpin*0.5;;
    }
}

G4Parton::G4Parton

(

const G4Parton &

right

)

Definition at line 107 of file G4Parton.cc.

{
    PDGencoding = right.PDGencoding;
    theMomentum = right.theMomentum;
    thePosition = right.thePosition;
    theX = right.theX;
    theDefinition = right.theDefinition;
    theColour = right.theColour;
    theIsoSpinZ = right.theIsoSpinZ;
    theSpinZ = right.theSpinZ;
}

G4Parton::~G4Parton

(

)

Definition at line 136 of file G4Parton.cc.

{
//  cout << "G4Parton::~G4Parton(): this = "<<this <<endl;
//  cout << "break here"<<this <<endl;
}

Member Function Documentation

void G4Parton::DefineMomentumInZ	(	G4double	aLightConeMomentum,
		G4bool	aDirection
	)

Definition at line 142 of file G4Parton.cc.

References Get4Momentum(), CLHEP::HepLorentzVector::px(), CLHEP::HepLorentzVector::py(), Set4Momentum(), CLHEP::HepLorentzVector::setE(), CLHEP::HepLorentzVector::setPz(), and sqr().

{
    G4double Mass = GetMass();
    G4LorentzVector a4Momentum = Get4Momentum();
    aLightConeMomentum*=theX;
    G4double TransverseMass2 = sqr(a4Momentum.px()) + sqr(a4Momentum.py()) + sqr(Mass);
    a4Momentum.setPz(0.5*(aLightConeMomentum - TransverseMass2/aLightConeMomentum)*(aDirection? 1: -1)); 
    a4Momentum.setE( 0.5*(aLightConeMomentum + TransverseMass2/aLightConeMomentum));
    Set4Momentum(a4Momentum);
}  

void G4Parton::DefineMomentumInZ	(	G4double	aLightConeMomentum,
		G4double	aLightConeE,
		G4bool	aDirection
	)

Definition at line 153 of file G4Parton.cc.

References Get4Momentum(), CLHEP::HepLorentzVector::px(), CLHEP::HepLorentzVector::py(), Set4Momentum(), CLHEP::HepLorentzVector::setE(), CLHEP::HepLorentzVector::setPz(), and sqr().

{
    G4double Mass = GetMass();
    G4LorentzVector a4Momentum = Get4Momentum();
    aLightConeMomentum*=theX;
    aLightConeE*=theX;
    G4double TransverseMass2 = sqr(a4Momentum.px()) + sqr(a4Momentum.py()) + sqr(Mass);
    a4Momentum.setPz(0.5*(aLightConeMomentum - aLightConeE - TransverseMass2/aLightConeMomentum)*(aDirection? 1: -1)); 
    a4Momentum.setE( 0.5*(aLightConeMomentum + aLightConeE + TransverseMass2/aLightConeMomentum));
    Set4Momentum(a4Momentum);
}  

const G4LorentzVector & G4Parton::Get4Momentum ( ) const

inline

Definition at line 140 of file G4Parton.hh.

Referenced by G4SoftStringBuilder::BuildString(), G4DiffractiveExcitation::CreateStrings(), DefineMomentumInZ(), G4VKinkyStringDecay::FragmentString(), G4FragmentingString::G4FragmentingString(), G4PartonPair::G4PartonPair(), operator=(), G4QGSParticipants::PerformDiffractiveCollisions(), G4QGSParticipants::PerformSoftCollisions(), and G4QGSDiffractiveExcitation::String().

{
    return theMomentum;
}

G4int G4Parton::GetColour ( )

inline

Definition at line 89 of file G4Parton.hh.

Referenced by G4PartonPair::G4PartonPair().

{return theColour;}

G4ParticleDefinition * G4Parton::GetDefinition ( )

inline

Definition at line 158 of file G4Parton.hh.

Referenced by G4DiffractiveExcitation::CreateStrings(), G4FragmentingString::G4FragmentingString(), and G4PartonPair::G4PartonPair().

{
    return theDefinition;
}

G4double G4Parton::GetIsoSpinZ ( )

inline

Definition at line 92 of file G4Parton.hh.

{return theIsoSpinZ;}

G4int G4Parton::GetPDGcode ( ) const

inline

Definition at line 124 of file G4Parton.hh.

Referenced by G4SoftStringBuilder::BuildString(), G4DiffractiveExcitation::CreateStrings(), G4ExcitedString::GetAntiColorParton(), G4ExcitedString::GetColorParton(), operator=(), G4QGSParticipants::PerformDiffractiveCollisions(), G4QGSParticipants::PerformSoftCollisions(), and G4QGSDiffractiveExcitation::String().

{
    return PDGencoding;
}

const G4ThreeVector & G4Parton::GetPosition ( void  ) const

inline

Definition at line 129 of file G4Parton.hh.

Referenced by G4VKinkyStringDecay::FragmentString(), G4ExcitedString::G4ExcitedString(), and operator=().

{
    return thePosition;
}

G4double G4Parton::GetSpinZ ( )

inline

Definition at line 95 of file G4Parton.hh.

Referenced by G4PartonPair::G4PartonPair().

{return theSpinZ;}

G4double G4Parton::GetX ( )

inline

Definition at line 86 of file G4Parton.hh.

Referenced by G4SoftStringBuilder::BuildString(), G4QGSParticipants::PerformDiffractiveCollisions(), and G4QGSParticipants::PerformSoftCollisions().

{return theX;}

int G4Parton::operator!= ( const G4Parton &  right ) const

inline

Definition at line 119 of file G4Parton.hh.

References right.

{
    return this!=&right;
}

G4Parton & G4Parton::operator=

(

const G4Parton &

right

)

Definition at line 119 of file G4Parton.cc.

References Get4Momentum(), GetPDGcode(), and GetPosition().

{
   if (this != &right)
   {
      PDGencoding=right.GetPDGcode();
      theMomentum=right.Get4Momentum();
      thePosition=right.GetPosition();
      theX = right.theX;
      theDefinition = right.theDefinition;
      theColour = right.theColour;
      theIsoSpinZ = right.theIsoSpinZ;
      theSpinZ = right.theSpinZ;
   }

    return *this;
}

int G4Parton::operator== ( const G4Parton &  right ) const

inline

Definition at line 114 of file G4Parton.hh.

References right.

{
    return this==&right;
}   

void G4Parton::Set4Momentum ( const G4LorentzVector &  aMomentum )

inline

Definition at line 145 of file G4Parton.hh.

Referenced by G4DiffractiveExcitation::CreateStrings(), DefineMomentumInZ(), G4VKinkyStringDecay::FragmentString(), and G4QGSDiffractiveExcitation::String().

{
    theMomentum=aMomentum;
}

void G4Parton::SetColour ( G4int  aColour )

inline

Definition at line 88 of file G4Parton.hh.

{theColour = aColour;}

void G4Parton::SetIsoSpinZ ( G4double  anIsoSpinZ )

inline

Definition at line 91 of file G4Parton.hh.

{theIsoSpinZ = anIsoSpinZ;}

void G4Parton::SetPosition ( const G4ThreeVector &  aPosition )

inline

Definition at line 134 of file G4Parton.hh.

Referenced by G4VKinkyStringDecay::FragmentString().

{
    thePosition=aPosition;
}

void G4Parton::SetSpinZ ( G4double  aSpinZ )

inline

Definition at line 94 of file G4Parton.hh.

{theSpinZ = aSpinZ;}

void G4Parton::SetX ( G4double  anX )

inline

Definition at line 85 of file G4Parton.hh.

{ theX = anX; }

---

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

• G4Parton.hh
• G4Parton.cc