47 G4String text =
"G4Parton::GetDefinition(): Encoding not in particle table";
72 G4String text =
"G4Parton::GetDefinition(): Particle is not a parton";
87 if (thisPDGiIsospin == 0) {
98 if (thisPDGiSpin == 0) {
146 aLightConeMomentum*=
theX;
147 G4double TransverseMass2 =
sqr(a4Momentum.px()) +
sqr(a4Momentum.py()) +
sqr(Mass);
148 a4Momentum.setPz(0.5*(aLightConeMomentum - TransverseMass2/aLightConeMomentum)*(aDirection? 1: -1));
149 a4Momentum.setE( 0.5*(aLightConeMomentum + TransverseMass2/aLightConeMomentum));
157 aLightConeMomentum*=
theX;
159 G4double TransverseMass2 =
sqr(a4Momentum.px()) +
sqr(a4Momentum.py()) +
sqr(Mass);
160 a4Momentum.setPz(0.5*(aLightConeMomentum - aLightConeE - TransverseMass2/aLightConeMomentum)*(aDirection? 1: -1));
161 a4Momentum.setE( 0.5*(aLightConeMomentum + aLightConeE + TransverseMass2/aLightConeMomentum));
G4int GetPDGiSpin() const
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
G4ThreeVector thePosition
G4double GetPDGIsospin3() const
const G4LorentzVector & Get4Momentum() const
G4Parton & operator=(const G4Parton &right)
void Set4Momentum(const G4LorentzVector &aMomentum)
void DefineMomentumInZ(G4double aLightConeMomentum, G4bool aDirection)
G4GLOB_DLL std::ostream G4cout
const G4String & GetParticleType() const
G4ParticleDefinition * theDefinition
const G4ThreeVector & GetPosition() const
static G4ParticleTable * GetParticleTable()
G4int GetPDGiIsospin() const
G4LorentzVector theMomentum
CLHEP::HepLorentzVector G4LorentzVector