Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 // ------------------------------------------------------------------- 28 // GEANT4 Class file 29 // 30 // For information related to this code contact: 31 // 32 // File name: G4VElasticCollision 33 // 34 // Author: Maria Grazia Pia 35 // 36 // Creation date: 15 April 1999 37 // 38 // Modifications: 39 // 40 // ------------------------------------------------------------------- 41 42 #ifndef G4VSCATTERINGCOLLISION_HH 43 #define G4VSCATTERINGCOLLISION_HH 44 45 #include "globals.hh" 46 #include "G4Log.hh" 47 #include "G4VCollision.hh" 48 #include "G4VCrossSectionSource.hh" 49 #include "G4VAngularDistribution.hh" 50 #include "G4KineticTrackVector.hh" 51 52 class G4KineticTrack; 53 54 55 class G4VScatteringCollision : public G4VCollision 56 { 57 58 public: 59 60 G4VScatteringCollision(); 61 void establish_G4MT_TLS_G4VScatteringCollision(); 62 virtual ~G4VScatteringCollision(); 63 64 G4bool operator==(const G4VScatteringCollision &right) const; 65 G4bool operator!=(const G4VScatteringCollision &right) const; 66 67 virtual G4KineticTrackVector* FinalState(const G4KineticTrack& trk1, 68 const G4KineticTrack& trk2) const; 69 virtual const G4VAngularDistribution* GetAngularDistribution() const 70 { 71 return theAngularDistribution; 72 } 73 74 private: 75 G4VScatteringCollision(const G4VScatteringCollision &); 76 G4VScatteringCollision & operator= (const G4VScatteringCollision &); 77 78 protected: 79 80 virtual const std::vector<const G4ParticleDefinition*> & GetOutgoingParticles() const = 0; 81 82 private: 83 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 ); } 86 87 G4double BrWigInt1(const G4double x, const G4double gamma, const G4double m0) const 88 { return 0.5*gamma*gamma*G4Log( (x-m0)*(x-m0)+gamma*gamma/4.0 ) + m0*BrWigInt0(x,gamma,m0); } 89 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; } 92 93 double SampleResonanceMass(const double poleMass, 94 const double width, 95 const double minMass, 96 const double maxMass) const; 97 private: 98 99 G4VAngularDistribution * theAngularDistribution; 100 101 }; 102 103 104 105 106 #endif 107