42#ifndef G4VSCATTERINGCOLLISION_HH
43#define G4VSCATTERINGCOLLISION_HH
84 double BrWigInt0(
const double x,
const double gamma,
const double m0)
const
85 {
return 2.0*gamma*std::atan( 2.0 * (x-m0)/ gamma ); }
88 {
return 0.5*gamma*gamma*
G4Log( (x-m0)*(x-m0)+gamma*gamma/4.0 ) + m0*
BrWigInt0(x,gamma,m0); }
90 double BrWigInv(
const double x,
const double gamma,
const double m0)
const
91 {
return 0.5*gamma*std::tan( 0.5*x/gamma )+m0; }
96 const double maxMass)
const;
G4double G4Log(G4double x)
G4double BrWigInt1(const G4double x, const G4double gamma, const G4double m0) const
G4VAngularDistribution * theAngularDistribution
G4bool operator!=(const G4VScatteringCollision &right) const
double BrWigInt0(const double x, const double gamma, const double m0) const
void establish_G4MT_TLS_G4VScatteringCollision()
virtual const G4VAngularDistribution * GetAngularDistribution() const
G4bool operator==(const G4VScatteringCollision &right) const
G4VScatteringCollision & operator=(const G4VScatteringCollision &)
double SampleResonanceMass(const double poleMass, const double width, const double minMass, const double maxMass) const
virtual ~G4VScatteringCollision()
G4VScatteringCollision(const G4VScatteringCollision &)
double BrWigInv(const double x, const double gamma, const double m0) const
virtual G4KineticTrackVector * FinalState(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual const std::vector< const G4ParticleDefinition * > & GetOutgoingParticles() const =0