69 G4cout <<
"G4PhononDownconversion::GetMeanFreePath = " << mfp <<
G4endl;
92 else MakeTTSecondaries(aTrack);
111 inline double G4PhononDownconversion::GetLTDecayProb(
double d,
double x)
const {
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));
120 inline double G4PhononDownconversion::GetTTDecayProb(
double d,
double x)
const {
122 G4double A = 0.5*(1-d*
d)*(fBeta+fLambda+(1+d*d)*(fGamma+fMu));
123 G4double B = fBeta+fLambda+2*d*d*(fGamma+fMu);
124 G4double C = fBeta + fLambda + 2*(fGamma+fMu);
125 G4double D = (1-d*
d)*(2*fBeta+4*fGamma+fLambda+3*fMu);
127 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)));
133 inline double G4PhononDownconversion::MakeLDeviation(
double d,
double x)
const {
136 return std::acos((1+(x*x)-((d*d)*(1-x)*(1-x)))/(2*x));
142 inline double G4PhononDownconversion::MakeTDeviation(
double d,
double x)
const {
145 return std::acos((1-x*x+d*d*(1-x)*(1-x))/(2*d*(1-x)));
151 inline double G4PhononDownconversion::MakeTTDeviation(
double d,
double x)
const {
154 return std::acos((1-d*d*(1-x)*(1-x)+d*d*x*x)/(2*d*x));
162 void G4PhononDownconversion::MakeTTSecondaries(
const G4Track& aTrack) {
175 while(p >= GetTTDecayProb(d, x*d)) {
181 G4double theta1=MakeTTDeviation(d, x);
182 G4double theta2=MakeTTDeviation(d, 1-x);
194 G4double Esec1 = x*E, Esec2 = E-Esec1;
218 void G4PhononDownconversion::MakeLTSecondaries(
const G4Track& aTrack) {
231 while(p >= GetLTDecayProb(d, x)) {
237 G4double thetaL=MakeLDeviation(d, x);
238 G4double thetaT=MakeTDeviation(d, x);
247 G4double Esec1 = x*E, Esec2 = E-Esec1;
G4PhononDownconversion(const G4String &processName="phononDownconversion")
G4double condition(const G4ErrorSymMatrix &m)
G4double GetGamma() const
const G4LatticePhysical * theLattice
G4double GetAnhDecConstant() const
G4double GetVelocity() const
virtual G4double GetMeanFreePath(const G4Track &, G4double, G4ForceCondition *)
G4double GetLambda() const
G4double GetSTDOS() const
virtual ~G4PhononDownconversion()
G4ThreeVector G4RandomDirection()
Definition of the G4PhononPolarization enum.
Definition of the G4PhononDownconversion class.
G4double GetFTDOS() const
G4PhononTrackMap * trackKmap
G4double GetKineticEnergy() const
Definition of the G4PhononTrackMap base class.
G4GLOB_DLL std::ostream G4cout
const G4ThreeVector & GetK(const G4Track *track) const
static G4PhononLong * PhononDefinition()
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
virtual G4bool IsApplicable(const G4ParticleDefinition &)
virtual void Initialize(const G4Track &)
virtual G4int ChoosePolarization(G4double Ldos, G4double STdos, G4double FTdos) const
void SetNumberOfSecondaries(G4int totSecondaries)
Hep3Vector orthogonal() const
void ProposeEnergy(G4double finalEnergy)
G4ParticleChange aParticleChange
Definition of the G4LatticePhysical class.
void AddSecondary(G4Track *aSecondary)
virtual G4Track * CreateSecondary(G4int polarization, const G4ThreeVector &K, G4double energy) const
void ProposeTrackStatus(G4TrackStatus status)
Hep3Vector & rotate(double, const Hep3Vector &)