| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 // Abstract base class for multibody "phase sp 27 // Abstract base class for multibody "phase space" generators. Subclasses
28 // implement a specific algorithm, such as Kop 28 // implement a specific algorithm, such as Kopylov, GENBOD, or Makoto's
29 // NBody. Subclasses are used by G4HadPhaseSp 29 // NBody. Subclasses are used by G4HadPhaseSpaceGenerator.
30 // 30 //
31 // Author: Michael Kelsey (SLAC) <kelsey@slac 31 // Author: Michael Kelsey (SLAC) <kelsey@slac.stanford.edu>
32 32
33 #ifndef G4VHadDecayAlgorithm_HH 33 #ifndef G4VHadDecayAlgorithm_HH
34 #define G4VHadDecayAlgorithm_HH 1 34 #define G4VHadDecayAlgorithm_HH 1
35 35
36 #include "globals.hh" 36 #include "globals.hh"
37 #include "G4LorentzVector.hh" 37 #include "G4LorentzVector.hh"
38 #include "G4ThreeVector.hh" 38 #include "G4ThreeVector.hh"
39 #include <vector> 39 #include <vector>
40 #include <iosfwd> 40 #include <iosfwd>
41 41
42 class G4VHadDecayAlgorithm { 42 class G4VHadDecayAlgorithm {
43 public: 43 public:
44 G4VHadDecayAlgorithm(const G4String& algName 44 G4VHadDecayAlgorithm(const G4String& algName, G4int verbose=0)
45 : name(algName), verboseLevel(verbose) {;} 45 : name(algName), verboseLevel(verbose) {;}
46 virtual ~G4VHadDecayAlgorithm() {;} 46 virtual ~G4VHadDecayAlgorithm() {;}
47 47
48 // Initial state (rest mass) and list of fin 48 // Initial state (rest mass) and list of final masses
49 void Generate(G4double initialMass, 49 void Generate(G4double initialMass,
50 const std::vector<G4double>& masses, 50 const std::vector<G4double>& masses,
51 std::vector<G4LorentzVector>& finalState); 51 std::vector<G4LorentzVector>& finalState);
52 52
53 // Enable (or disable if 0) diagnostic messa 53 // Enable (or disable if 0) diagnostic messages (subclass may overload)
54 virtual void SetVerboseLevel(G4int verbose) 54 virtual void SetVerboseLevel(G4int verbose) { verboseLevel = verbose; }
55 G4int GetVerboseLevel() const { return verbo 55 G4int GetVerboseLevel() const { return verboseLevel; }
56 const G4String& GetName() const { return nam 56 const G4String& GetName() const { return name; }
57 57
58 protected: 58 protected:
59 // Subclasses MUST implement these functions 59 // Subclasses MUST implement these functions
60 virtual void GenerateTwoBody(G4double initia 60 virtual void GenerateTwoBody(G4double initialMass,
61 const std::vector<G4double>& mass 61 const std::vector<G4double>& masses,
62 std::vector<G4LorentzVector>& fin 62 std::vector<G4LorentzVector>& finalState) = 0;
63 63
64 virtual void GenerateMultiBody(G4double init 64 virtual void GenerateMultiBody(G4double initialMass,
65 const std::vector<G4double>& masses, 65 const std::vector<G4double>& masses,
66 std::vector<G4LorentzVector>& finalSt 66 std::vector<G4LorentzVector>& finalState) = 0;
67 67
68 // Validate kinematics (e.g., limit number o 68 // Validate kinematics (e.g., limit number of final state particles)
69 // Subclasses may override or call back to t 69 // Subclasses may override or call back to this function
70 virtual G4bool IsDecayAllowed(G4double initi 70 virtual G4bool IsDecayAllowed(G4double initialMass,
71 const std::vector<G4double>& masses) c 71 const std::vector<G4double>& masses) const;
72 72
73 // Two-body momentum function (c.f. PDK from 73 // Two-body momentum function (c.f. PDK from CERNLIB W505)
74 G4double TwoBodyMomentum(G4double M0, G4doub 74 G4double TwoBodyMomentum(G4double M0, G4double M1, G4double M2) const;
75 75
76 // Convenience functions for uniform angular 76 // Convenience functions for uniform angular distributions
77 G4double UniformTheta() const; 77 G4double UniformTheta() const;
78 G4double UniformPhi() const; 78 G4double UniformPhi() const;
79 79
80 // Utility to dump vector contents to line o 80 // Utility to dump vector contents to line of output
81 void PrintVector(const std::vector<G4double> 81 void PrintVector(const std::vector<G4double>& v, const G4String& name,
82 std::ostream& os) const; 82 std::ostream& os) const;
83 83
84 private: 84 private:
85 G4String name; 85 G4String name;
86 G4int verboseLevel; 86 G4int verboseLevel;
87 }; 87 };
88 88
89 #endif /* G4VHadDecayAlgorithm_HH */ 89 #endif /* G4VHadDecayAlgorithm_HH */
90 90