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 // Multibody "phase space" generator using GENBOD (CERNLIB W515) method. 28 // 29 // Author: Michael Kelsey (SLAC) <kelsey@slac.stanford.edu> 30 31 #ifndef G4HadPhaseSpaceGenbod_HH 32 #define G4HadPhaseSpaceGenbod_HH 1 33 34 #include "G4VHadPhaseSpaceAlgorithm.hh" 35 36 37 class G4HadPhaseSpaceGenbod : public G4VHadPhaseSpaceAlgorithm { 38 public: 39 G4HadPhaseSpaceGenbod(G4int verbose=0); 40 virtual ~G4HadPhaseSpaceGenbod() {;} 41 42 protected: 43 virtual void GenerateMultiBody(G4double initialMass, 44 const std::vector<G4double>& masses, 45 std::vector<G4LorentzVector>& finalState); 46 47 protected: 48 void Initialize(G4double initialMass, 49 const std::vector<G4double>& masses); 50 51 void FillRandomBuffer(); 52 53 void ComputeWeightScale(const std::vector<G4double>& masses); 54 55 void FillEnergySteps(G4double initialMass, 56 const std::vector<G4double>& masses); 57 58 void GenerateMomenta(const std::vector<G4double>& masses, 59 std::vector<G4LorentzVector>& finalState); 60 61 void AccumulateFinalState(size_t i, 62 const std::vector<G4double>& masses, 63 std::vector<G4LorentzVector>& finalState); 64 65 G4bool AcceptEvent() const; // Use accept-reject to generate distribution 66 G4double ComputeWeight() const; 67 68 private: 69 size_t nFinal; // Number of final state particles 70 G4double totalMass; // Sum of final state masses 71 G4double massExcess; // Available kinetic energy 72 G4double weightMax; // Maximum possible weight 73 G4int nTrials; // Accept/reject cycles taken 74 75 std::vector<G4double> msum; // Cumulative sum of masses 76 std::vector<G4double> msq; // Final state squared masses 77 std::vector<G4double> rndm; // Random sequence for effective masses 78 std::vector<G4double> meff; // Random final-state effective masses 79 std::vector<G4double> pd; // Random momentum magnitudes 80 }; 81 82 #endif /* G4HadPhaseSpaceGenbod_HH */ 83