78 etrans += bond_energy;
95 if(ne < 0) { ne = 0; }
98 CLHEP::electron_mass_c2 - bond_energy);
102 G4double sinTheta = std::sqrt(1. - cosTheta * cosTheta);
114 G4double energy = 0.5*((ecm - mass)*(ecm + mass) + emass*emass)/ecm;
117 G4double mom = std::sqrt((energy - emass)*(energy + emass));
120 mom * sinTheta * std::sin(phi),
121 mom * cosTheta, energy);
136 if(P2 > 0.0) { p = std::sqrt(P2); }
138 lv.set(v.x()*p, v.y()*p, v.z()*p, E);
CLHEP::Hep3Vector G4ThreeVector
G4int GetNumberOfElectrons() const
void SetNumberOfElectrons(G4int value)
virtual ~G4GammaTransition()
static G4double GetBindingEnergy(G4int Z, G4int SubshellNb)
const G4LorentzVector & GetMomentum() const
void SetMomentum(const G4LorentzVector &value)
static const double twopi
G4double GetGroundStateMass() const
G4double GetPDGMass() const
virtual G4Fragment * SampleTransition(G4Fragment *nucleus, G4double newExcEnergy, G4int deltaS, size_t shell, G4bool isGamma, G4bool isLongLived)
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()
G4double GetExcitationEnergy() const
CLHEP::HepLorentzVector G4LorentzVector
static G4int GetNumberOfShells(G4int Z)