49 : polarFlag(false), fDirection(0.,0.,0.), fTwoJMAX(10), fVerbose(0)
82 G4cout <<
"G4GammaTransition::GenerateGamma - Etrans(MeV)= "
83 << etrans <<
" Eexnew= " << newExcEnergy
84 <<
" Ebond= " << bond_energy <<
G4endl;
87 etrans += bond_energy;
121 G4double energy = 0.5*((ecm - mass)*(ecm + mass) + emass*emass)/ecm;
146 G4cout <<
"G4GammaTransition::SampleTransition : " << *result <<
G4endl;
158 G4cout <<
"G4GammaTransition::SampleDirection : 2J1= " << twoJ1
159 <<
" 2J2= " << twoJ2 <<
" ratio= " << ratio
160 <<
" mp= " << mp <<
G4endl;
177 G4int L0 = 0, Lp = 0;
189 G4double sinTheta = std::sqrt((1.-cosTheta)*(1.+cosTheta));
190 fDirection.
set(sinTheta*std::cos(phi),sinTheta*std::sin(phi),cosTheta);
G4ThreeVector G4RandomDirection()
G4GLOB_DLL std::ostream G4cout
void set(double x, double y, double z)
Hep3Vector boostVector() const
HepLorentzVector & boost(double, double, double)
void set(double x, double y, double z, double t)
static G4double GetBindingEnergy(G4int Z, G4int SubshellNb)
static G4int GetNumberOfShells(G4int Z)
static G4Electron * Electron()
G4double GetGroundStateMass() const
G4double GetExcitationEnergy() const
const G4LorentzVector & GetMomentum() const
void SetNumberOfElectrons(G4int value)
void SetExcEnergyAndMomentum(G4double eexc, const G4LorentzVector &)
G4int GetNumberOfElectrons() const
virtual void SampleDirection(G4Fragment *nuc, G4double ratio, G4int twoJ1, G4int twoJ2, G4int mp)
virtual ~G4GammaTransition()
G4PolarizationTransition fPolTrans
virtual G4Fragment * SampleTransition(G4Fragment *nucleus, G4double newExcEnergy, G4double mpRatio, G4int JP1, G4int JP2, G4int MP, G4int shell, G4bool isDiscrete, G4bool isGamma)
G4double GetPDGMass() const
void SampleGammaTransition(G4NuclearPolarization *np, G4int twoJ1, G4int twoJ2, G4int L0, G4int Lp, G4double mpRatio, G4double &cosTheta, G4double &phi)
static constexpr double electron_mass_c2
static constexpr double twopi
G4double energy(const ThreeVector &p, const G4double m)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
T min(const T t1, const T t2)
brief Return the smallest of the two arguments