82 if(etrans <= 0.0) {
return result; }
97 if(ne < 0) { ne = 0; }
100 CLHEP::electron_mass_c2 - bond_energy);
104 G4double sinTheta = std::sqrt(1. - cosTheta * cosTheta);
116 G4double energy = 0.5*((ecm - mass)*(ecm + mass) + emass*emass)/ecm;
119 G4double mom = std::sqrt((energy - emass)*(energy + emass));
122 mom * sinTheta * std::sin(phi),
123 mom * cosTheta, energy);
138 if(P2 > 0.0) { p = std::sqrt(P2); }
140 lv.set(v.x()*p, v.y()*p, v.z()*p, E);
CLHEP::Hep3Vector G4ThreeVector
G4int GetNumberOfElectrons() const
void SetNumberOfElectrons(G4int value)
G4PolarizationTransition * fPolarization
G4PolarizedGammaTransition()
static G4double GetBindingEnergy(G4int Z, G4int SubshellNb)
const G4LorentzVector & GetMomentum() const
void SetMomentum(const G4LorentzVector &value)
static const double twopi
G4double GetGroundStateMass() const
virtual ~G4PolarizedGammaTransition()
G4double GetPDGMass() const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double energy(const ThreeVector &p, const G4double m)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static const G4double * P2[nN]
static G4Electron * Electron()
virtual G4Fragment * SampleTransition(G4Fragment *nucleus, G4double newExcEnergy, G4int deltaS, size_t shell, G4bool isGamma, G4bool isLongLived)
G4double GetExcitationEnergy() const
CLHEP::HepLorentzVector G4LorentzVector
static G4int GetNumberOfShells(G4int Z)