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 #ifndef G4QGSMSplitableHadron_h 27 #define G4QGSMSplitableHadron_h 1 28 29 #include "G4VSplitableHadron.hh" 30 #include "G4PartonVector.hh" 31 #include "G4MesonSplitter.hh" 32 #include "G4BaryonSplitter.hh" 33 #include "Randomize.hh" 34 #include <deque> 35 36 // based on prototype by Maxim Komogorov 37 // Splitting into methods, and centralizing of model parameters HPW Feb 1999 38 // continued clean-up of interfaces and algorithms HPW 1999. 39 // Redesign of data structures and algorithms HPW Feb 1999 40 41 class G4QGSMSplitableHadron : public G4VSplitableHadron 42 { 43 public: 44 G4QGSMSplitableHadron(); 45 G4QGSMSplitableHadron(const G4ReactionProduct & aPrimary); 46 G4QGSMSplitableHadron(const G4ReactionProduct & aPrimary, G4bool Direction); 47 G4QGSMSplitableHadron(const G4Nucleon & aNucleon); 48 G4QGSMSplitableHadron(const G4Nucleon & aNucleon, G4bool Direction); 49 50 virtual ~G4QGSMSplitableHadron(); 51 52 private: 53 const G4QGSMSplitableHadron & operator=(const G4QGSMSplitableHadron &right); 54 55 public: 56 virtual void SplitUp(); 57 virtual void SetFirstParton(G4int PDGcode); 58 virtual void SetSecondParton(G4int PDGcode); 59 virtual G4Parton * GetNextParton(); 60 virtual G4Parton * GetNextAntiParton(); 61 62 private: 63 void InitParameters(); 64 void DiffractiveSplitUp(); 65 void SoftSplitUp(); 66 67 G4ThreeVector GaussianPt(G4double widthSquare, G4double maxPtSquare); 68 void GetValenceQuarkFlavors(const G4ParticleDefinition * aPart, 69 G4Parton *& Parton1, G4Parton *& Parton2); 70 G4Parton * BuildSeaQuark(G4bool isAntiQuark, G4int aPDGCode, G4int nSeaPair); 71 G4double SampleX(G4double anXmin, G4int nSea, G4int theTotalSea, G4double aBeta); 72 73 private: 74 // aggregated data 75 G4bool Direction; // FALSE is target. - candidate for more detailed design. 76 77 std::deque<G4Parton *> Color; 78 std::deque<G4Parton *> AntiColor; 79 unsigned int iP; 80 unsigned int iAP; 81 82 private: 83 // associated classes 84 G4MesonSplitter theMesonSplitter; 85 G4BaryonSplitter theBaryonSplitter; 86 87 private: 88 // model parameters 89 G4double alpha; 90 G4double beta; 91 G4double theMinPz; 92 G4double StrangeSuppress; 93 G4double sigmaPt; 94 G4double widthOfPtSquare; 95 G4double minTransverseMass; 96 }; 97 98 inline G4Parton* G4QGSMSplitableHadron::GetNextParton() 99 { 100 if(Color.size()==0) return 0; 101 102 G4Parton * result = Color.operator[](iP); 103 iP++; if(iP == Color.size()) iP=0; 104 return result; 105 } 106 107 inline G4Parton* G4QGSMSplitableHadron::GetNextAntiParton() 108 { 109 if(AntiColor.size() == 0) return 0; 110 111 G4Parton * result = AntiColor.operator[](iAP); 112 iAP++; if(iAP == AntiColor.size()) iAP=0; 113 return result; 114 } 115 116 inline void G4QGSMSplitableHadron::SetFirstParton(G4int) 117 {} 118 inline void G4QGSMSplitableHadron::SetSecondParton(G4int) 119 {} 120 #endif 121 122