68 G4cout <<
"G4PhononDownconversion::GetMeanFreePath = " << mfp <<
G4endl;
91 else MakeTTSecondaries(aTrack);
112 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));
120 G4double A = 0.5*(1-d*d)*(fBeta+fLambda+(1+d*d)*(fGamma+fMu));
121 G4double B = fBeta+fLambda+2*d*d*(fGamma+fMu);
122 G4double C = fBeta + fLambda + 2*(fGamma+fMu);
123 G4double D = (1-d*d)*(2*fBeta+4*fGamma+fLambda+3*fMu);
125 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)));
134 return std::acos((1+(x*x)-((d*d)*(1-x)*(1-x)))/(2*x));
143 return std::acos((1-x*x+d*d*(1-x)*(1-x))/(2*d*(1-x)));
152 return std::acos((1-d*d*(1-x)*(1-x)+d*d*x*x)/(2*d*x));
160void G4PhononDownconversion::MakeTTSecondaries(
const G4Track& aTrack) {
173 while(p >= GetTTDecayProb(d, x*d)) {
179 G4double theta1=MakeTTDeviation(d, x);
180 G4double theta2=MakeTTDeviation(d, 1-x);
192 G4double Esec1 = x*E, Esec2 = E-Esec1;
216void G4PhononDownconversion::MakeLTSecondaries(
const G4Track& aTrack) {
229 while(p >= GetLTDecayProb(d, x)) {
235 G4double thetaL=MakeLDeviation(d, x);
236 G4double thetaT=MakeTDeviation(d, x);
245 G4double Esec1 = x*E, Esec2 = E-Esec1;
double B(double temperature)
double A(double temperature)
G4double condition(const G4ErrorSymMatrix &m)
Definition of the G4LatticePhysical class.
G4ThreeVector G4RandomDirection()
G4GLOB_DLL std::ostream G4cout
Hep3Vector orthogonal() const
Hep3Vector & rotate(double, const Hep3Vector &)
G4double GetFTDOS() const
G4double GetAnhDecConstant() const
G4double GetGamma() const
G4double GetSTDOS() const
G4double GetLambda() const
void AddSecondary(G4Track *aSecondary)
void ProposeEnergy(G4double finalEnergy)
virtual void Initialize(const G4Track &)
virtual G4double GetMeanFreePath(const G4Track &, G4double, G4ForceCondition *)
virtual G4bool IsApplicable(const G4ParticleDefinition &)
G4PhononDownconversion(const G4String &processName="phononDownconversion")
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
virtual ~G4PhononDownconversion()
static G4PhononLong * PhononDefinition()
const G4ThreeVector & GetK(const G4Track *track) const
G4double GetVelocity() const
G4double GetKineticEnergy() const
void ProposeTrackStatus(G4TrackStatus status)
void SetNumberOfSecondaries(G4int totSecondaries)
virtual G4int ChoosePolarization(G4double Ldos, G4double STdos, G4double FTdos) const
virtual G4Track * CreateSecondary(G4int polarization, const G4ThreeVector &K, G4double energy) const
const G4LatticePhysical * theLattice
G4PhononTrackMap * trackKmap
G4ParticleChange aParticleChange