43 #ifndef G4VANNIHILATIONCOLLISION_HH
44 #define G4VANNIHILATIONCOLLISION_HH
72 return theAngularDistribution;
82 double BrWigInt0(
const double x,
const double gamma,
const double m0)
const
83 {
return 2.0*gamma*std::atan( 2.0 * (x-m0)/ gamma ); }
86 {
return 0.5*gamma*gamma*
G4Log( (x-m0)*(x-m0)+gamma*gamma/4.0 ) + m0*BrWigInt0(x,gamma,m0); }
88 double BrWigInv(
const double x,
const double gamma,
const double m0)
const
89 {
return 0.5*gamma*std::tan( 0.5*x/gamma )+m0; }
91 double SampleResonanceMass(
const double poleMass,
94 const double maxMass)
const;
virtual G4KineticTrackVector * FinalState(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual const G4ParticleDefinition * GetOutgoingParticle(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0
virtual const G4VAngularDistribution * GetAngularDistribution() const
G4bool operator!=(const G4VAnnihilationCollision &right) const
G4VAnnihilationCollision()
G4double G4Log(G4double x)
virtual ~G4VAnnihilationCollision()
G4bool operator==(const G4VAnnihilationCollision &right) const