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 // 27 // 28 // ------------------------------------------- 28 // ----------------------------------------------------------------------------- 29 // GEANT 4 class implementation file 29 // GEANT 4 class implementation file 30 // 30 // 31 // History: first implementation, Maxim K 31 // History: first implementation, Maxim Komogorov, 10-Jul-1998 32 // ------------------------------------------- 32 // ----------------------------------------------------------------------------- 33 #ifndef G4QGSMFragmentation_h 33 #ifndef G4QGSMFragmentation_h 34 #define G4QGSMFragmentation_h 1 34 #define G4QGSMFragmentation_h 1 35 35 36 #include "G4VLongitudinalStringDecay.hh" 36 #include "G4VLongitudinalStringDecay.hh" 37 37 38 //******************************************** 38 //****************************************************************************** 39 class G4QGSMFragmentation:public G4VLongitudin 39 class G4QGSMFragmentation:public G4VLongitudinalStringDecay 40 { 40 { 41 public: 41 public: 42 G4QGSMFragmentation(); 42 G4QGSMFragmentation(); 43 ~G4QGSMFragmentation(); 43 ~G4QGSMFragmentation(); 44 virtual G4KineticTrackVector* FragmentStri 44 virtual G4KineticTrackVector* FragmentString(const G4ExcitedString& theString); 45 45 46 private: 46 private: 47 // not implemented to protect/forbid use 47 // not implemented to protect/forbid use 48 G4QGSMFragmentation(const G4QGSMFragmentat 48 G4QGSMFragmentation(const G4QGSMFragmentation &right); 49 const G4QGSMFragmentation & operator=(cons 49 const G4QGSMFragmentation & operator=(const G4QGSMFragmentation &right); 50 G4bool operator==(const G4QGSMFragmentatio << 50 int operator==(const G4QGSMFragmentation &right) const; 51 G4bool operator!=(const G4QGSMFragmentatio << 51 int operator!=(const G4QGSMFragmentation &right) const; 52 52 53 private: 53 private: 54 virtual G4bool StopFragmenting(const G4Fra << 54 void SetMinMasses(); 55 virtual G4bool IsItFragmentable(const G4Fr << 55 void SetMinimalStringMass(const G4FragmentingString * const string); >> 56 void SetMinimalStringMass2(const G4double aValue); >> 57 >> 58 virtual G4bool StopFragmenting(const G4FragmentingString * const string); >> 59 virtual G4bool IsFragmentable(const G4FragmentingString * const string); 56 60 57 virtual G4bool SplitLast(G4FragmentingStri 61 virtual G4bool SplitLast(G4FragmentingString * string, 58 G4KineticTrackVector * LeftVe 62 G4KineticTrackVector * LeftVector, 59 G4KineticTrackVec 63 G4KineticTrackVector * RightVector); 60 64 61 virtual void Sample4Momentum(G4LorentzVect 65 virtual void Sample4Momentum(G4LorentzVector* Mom, G4double Mass, G4LorentzVector* AntiMom, 62 G4double Anti 66 G4double AntiMass, G4double InitialMass); 63 67 64 virtual G4KineticTrack * Splitup(G4Fragmen 68 virtual G4KineticTrack * Splitup(G4FragmentingString *string, 65 G4Fragmen 69 G4FragmentingString *&newString); 66 70 67 // The hadron can be producet at QuarkSpli << 68 // virtual G4ParticleDefinition * QuarkSpl << 69 // G4Particl << 70 << 71 virtual G4ParticleDefinition * DiQuarkSpli << 72 << 73 << 74 virtual G4LorentzVector * SplitEandP(G4Par 71 virtual G4LorentzVector * SplitEandP(G4ParticleDefinition * pHadron, 75 G4Fra 72 G4FragmentingString * string, 76 G4Fra 73 G4FragmentingString * newString); 77 74 78 virtual G4double GetLightConeZ(G4double zm 75 virtual G4double GetLightConeZ(G4double zmin, G4double zmax, G4int PartonEncoding, 79 G4ParticleD 76 G4ParticleDefinition* pHadron, G4double Px, G4double Py); 80 77 81 private: << 78 virtual G4ParticleDefinition * DiQuarkSplitup(G4ParticleDefinition* decay, 82 // model parameters << 79 G4ParticleDefinition *&created); 83 G4double arho; << 84 G4double aphi; << 85 G4double aJPs; // alpha_J/Psi << 86 G4double aUps; // alpha_Y << 87 G4double an; << 88 G4double ala; << 89 << 90 G4double alaC; << 91 G4double alaB; << 92 G4double aXi; << 93 G4double aXiC; << 94 G4double aXiB; << 95 G4double aXiCC; << 96 G4double aXiCB; << 97 G4double aXiBB; << 98 << 99 G4double aksi; << 100 G4double alft; << 101 << 102 G4double FFq2q[5][5][2]; << 103 G4double FFq2qq[5][15][2]; << 104 G4double FFqq2q[15][5][2]; << 105 G4double FFqq2qq[15][5][2]; << 106 80 107 private: 81 private: 108 void SetFFq2q(); << 82 // ------ For estimation of a minimal string mass --------------- 109 void SetFFq2qq(); << 83 G4double Mass_of_light_quark; 110 void SetFFqq2qq(); << 84 G4double Mass_of_heavy_quark; 111 void SetFFqq2q(); << 85 G4double Mass_of_string_junction; >> 86 >> 87 G4double minMassQQbarStr[3][3]; >> 88 G4double minMassQDiQStr[3][3][3]; >> 89 >> 90 // ------ An estimated minimal string mass ---------------------- >> 91 G4double MinimalStringMass; >> 92 G4double MinimalStringMass2; >> 93 >> 94 // model parameters >> 95 const G4double arho; >> 96 const G4double aphi; >> 97 const G4double an; >> 98 const G4double ala; >> 99 const G4double aksi; >> 100 const G4double alft; 112 101 113 G4int IndexDiQ[5][5]; << 114 }; 102 }; 115 103 116 // Class G4QGSMFragmentation 104 // Class G4QGSMFragmentation 117 << 118 #endif 105 #endif 119 106 120 107