69 G4cout <<
"G4PhononDownconversion::GetMeanFreePath = " << mfp <<
G4endl;
113 return (1/(x*x))*(1-x*
x)*(1-x*x)*((1+
x)*(1+x)-d*d*((1-
x)*(1-x)))*(1+x*x-d*d*(1-
x)*(1-x))*(1+x*x-d*d*(1-x)*(1-
x));
126 return (A+B*d*x-B*x*x)*(A+B*d*x-B*x*
x)+(C*x*(d-x)-D/(d-
x)*(x-d-(1-d*d)/(4*
x)))*(C*x*(d-x)-D/(d-
x)*(x-d-(1-d*d)/(4*
x)));
135 return std::acos((1+(x*x)-((d*d)*(1-x)*(1-x)))/(2*x));
144 return std::acos((1-x*x+d*d*(1-x)*(1-x))/(2*d*(1-x)));
153 return std::acos((1-d*d*(1-x)*(1-x)+d*d*x*x)/(2*d*x));
G4double GetTTDecayProb(G4double, G4double) const
G4double GetKineticEnergy() const
virtual G4double GetMeanFreePath(const G4Track &, G4double, G4ForceCondition *)
void MakeLTSecondaries(const G4Track &)
G4double MakeTTDeviation(G4double, G4double) const
Definition of the G4LatticePhysical class.
void MakeTTSecondaries(const G4Track &)
G4ThreeVector G4RandomDirection()
G4double condition(const G4ErrorSymMatrix &m)
static constexpr double h_Planck
void AddSecondary(G4Track *aSecondary)
G4ParticleChange aParticleChange
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
G4double GetAnhDecConstant() const
virtual void Initialize(const G4Track &)
virtual G4int ChoosePolarization(G4double Ldos, G4double STdos, G4double FTdos) const
G4PhononDownconversion(const G4String &processName="phononDownconversion")
G4double GetLambda() const
static constexpr double twopi
double A(double temperature)
const G4LatticePhysical * theLattice
G4double GetSTDOS() const
virtual G4bool IsApplicable(const G4ParticleDefinition &)
void ProposeEnergy(G4double finalEnergy)
G4double MakeLDeviation(G4double, G4double) const
G4double GetVelocity() const
G4double GetLTDecayProb(G4double, G4double) const
const G4ThreeVector & GetK(const G4Track *track) const
Hep3Vector & rotate(double, const Hep3Vector &)
static G4PhononLong * PhononDefinition()
G4double GetFTDOS() const
G4PhononTrackMap * trackKmap
virtual ~G4PhononDownconversion()
G4GLOB_DLL std::ostream G4cout
virtual G4Track * CreateSecondary(G4int polarization, const G4ThreeVector &K, G4double energy) const
Hep3Vector orthogonal() const
G4double MakeTDeviation(G4double, G4double) const
double B(double temperature)
void ProposeTrackStatus(G4TrackStatus status)
G4double GetGamma() const
void SetNumberOfSecondaries(G4int totSecondaries)