Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4VHadDecayAlgorithm.hh
Go to the documentation of this file.
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// Abstract base class for multibody "phase space" generators. Subclasses
28// implement a specific algorithm, such as Kopylov, GENBOD, or Makoto's
29// NBody. Subclasses are used by G4HadPhaseSpaceGenerator.
30//
31// Author: Michael Kelsey (SLAC) <[email protected]>
32
33#ifndef G4VHadDecayAlgorithm_HH
34#define G4VHadDecayAlgorithm_HH 1
35
36#include "globals.hh"
37#include "G4LorentzVector.hh"
38#include "G4ThreeVector.hh"
39#include <vector>
40#include <iosfwd>
41
43public:
44 G4VHadDecayAlgorithm(const G4String& algName, G4int verbose=0)
45 : name(algName), verboseLevel(verbose) {;}
47
48 // Initial state (rest mass) and list of final masses
49 void Generate(G4double initialMass,
50 const std::vector<G4double>& masses,
51 std::vector<G4LorentzVector>& finalState);
52
53 // Enable (or disable if 0) diagnostic messages (subclass may overload)
54 virtual void SetVerboseLevel(G4int verbose) { verboseLevel = verbose; }
55 G4int GetVerboseLevel() const { return verboseLevel; }
56 const G4String& GetName() const { return name; }
57
58protected:
59 // Subclasses MUST implement these functions
60 virtual void GenerateTwoBody(G4double initialMass,
61 const std::vector<G4double>& masses,
62 std::vector<G4LorentzVector>& finalState) = 0;
63
64 virtual void GenerateMultiBody(G4double initialMass,
65 const std::vector<G4double>& masses,
66 std::vector<G4LorentzVector>& finalState) = 0;
67
68 // Validate kinematics (e.g., limit number of final state particles)
69 // Subclasses may override or call back to this function
70 virtual G4bool IsDecayAllowed(G4double initialMass,
71 const std::vector<G4double>& masses) const;
72
73 // Two-body momentum function (c.f. PDK from CERNLIB W505)
75
76 // Convenience functions for uniform angular distributions
77 G4double UniformTheta() const;
78 G4double UniformPhi() const;
79
80 // Utility to dump vector contents to line of output
81 void PrintVector(const std::vector<G4double>& v, const G4String& name,
82 std::ostream& os) const;
83
84private:
85 G4String name;
86 G4int verboseLevel;
87};
88
89#endif /* G4VHadDecayAlgorithm_HH */
double G4double
Definition: G4Types.hh:83
bool G4bool
Definition: G4Types.hh:86
int G4int
Definition: G4Types.hh:85
virtual void GenerateTwoBody(G4double initialMass, const std::vector< G4double > &masses, std::vector< G4LorentzVector > &finalState)=0
G4double UniformTheta() const
void Generate(G4double initialMass, const std::vector< G4double > &masses, std::vector< G4LorentzVector > &finalState)
const G4String & GetName() const
virtual G4bool IsDecayAllowed(G4double initialMass, const std::vector< G4double > &masses) const
void PrintVector(const std::vector< G4double > &v, const G4String &name, std::ostream &os) const
G4VHadDecayAlgorithm(const G4String &algName, G4int verbose=0)
G4double TwoBodyMomentum(G4double M0, G4double M1, G4double M2) const
virtual void SetVerboseLevel(G4int verbose)
virtual void GenerateMultiBody(G4double initialMass, const std::vector< G4double > &masses, std::vector< G4LorentzVector > &finalState)=0