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 // Hadronic Process: Phase space decay for the Fermi BreakUp model 28 // by V. Lara 29 // 30 // Modifications: 31 // 01.04.2011 General cleanup by V.Ivanchenko: 32 // - IsotropicVector is inlined 33 // - Momentum computation return zero or positive value 34 // - DumpProblem method is added providing more information 35 // - Reduced usage of exotic std functions 36 // 37 38 #ifndef G4FermiPhaseSpaceDecay_hh 39 #define G4FermiPhaseSpaceDecay_hh 1 40 41 #include "G4LorentzVector.hh" 42 #include "G4ThreeVector.hh" 43 #include "Randomize.hh" 44 45 #include <vector> 46 47 class G4Pow; 48 49 class G4FermiPhaseSpaceDecay 50 { 51 public: 52 53 G4FermiPhaseSpaceDecay(); 54 ~G4FermiPhaseSpaceDecay(); 55 56 std::vector<G4LorentzVector*>* Decay(G4double parent_mass, 57 const std::vector<G4double>& fragment_masses) const; 58 59 private: 60 61 inline G4double PtwoBody(G4double E, G4double P1, G4double P2) const; 62 63 G4double BetaKopylov(G4int, CLHEP::HepRandomEngine*) const; 64 65 std::vector<G4LorentzVector*> * 66 KopylovNBodyDecay(G4double, const std::vector<G4double>&) const; 67 68 G4FermiPhaseSpaceDecay(const G4FermiPhaseSpaceDecay&) = delete; 69 const G4FermiPhaseSpaceDecay & operator= 70 (const G4FermiPhaseSpaceDecay &) = delete; 71 G4bool operator==(const G4FermiPhaseSpaceDecay&) = delete; 72 G4bool operator!=(const G4FermiPhaseSpaceDecay&) = delete; 73 74 G4Pow* g4calc; 75 }; 76 77 inline G4double 78 G4FermiPhaseSpaceDecay::PtwoBody(G4double E, G4double P1, G4double P2) const 79 { 80 G4double P = (E+P1+P2)*(E+P1-P2)*(E-P1+P2)*(E-P1-P2)/(4.0*E*E); 81 return (P>0.0) ? std::sqrt(P) : 0.0; 82 } 83 84 #endif 85