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 // $Id: G4LundStringFragmentation.hh 100828 2016-11-02 15:25:59Z gcosmo $ 27 // 28 // 28 // ------------------------------------------- 29 // ----------------------------------------------------------------------------- 29 // GEANT 4 class implementation file 30 // GEANT 4 class implementation file 30 // 31 // 31 // History: first implementation, Maxim K 32 // History: first implementation, Maxim Komogorov, 10-Jul-1998 32 // ------------------------------------------- 33 // ----------------------------------------------------------------------------- 33 34 34 #ifndef G4LundStringFragmentation_h 35 #ifndef G4LundStringFragmentation_h 35 #define G4LundStringFragmentation_h 1 36 #define G4LundStringFragmentation_h 1 36 37 37 #include "G4VLongitudinalStringDecay.hh" 38 #include "G4VLongitudinalStringDecay.hh" 38 39 39 //******************************************** << 40 //************************************************************************************************************** 40 41 41 class G4LundStringFragmentation: public G4VLon 42 class G4LundStringFragmentation: public G4VLongitudinalStringDecay 42 { 43 { 43 public: 44 public: >> 45 44 G4LundStringFragmentation(); 46 G4LundStringFragmentation(); 45 virtual ~G4LundStringFragmentation(); 47 virtual ~G4LundStringFragmentation(); 46 48 47 virtual G4KineticTrackVector* FragmentStri 49 virtual G4KineticTrackVector* FragmentString(const G4ExcitedString& theString); 48 50 49 private: 51 private: 50 // not implemented to protect/forbid use 52 // not implemented to protect/forbid use 51 G4LundStringFragmentation(const G4LundStri 53 G4LundStringFragmentation(const G4LundStringFragmentation &right); 52 const G4LundStringFragmentation & operator 54 const G4LundStringFragmentation & operator=(const G4LundStringFragmentation &right); 53 G4bool operator==(const G4LundStringFragme << 55 int operator==(const G4LundStringFragmentation &right) const; 54 G4bool operator!=(const G4LundStringFragme << 56 int operator!=(const G4LundStringFragmentation &right) const; 55 57 56 private: 58 private: >> 59 void SetMinimalStringMass(const G4FragmentingString * const string); >> 60 void SetMinimalStringMass2(const G4double aValue); >> 61 57 virtual G4bool StopFragmenting(const G4Fra 62 virtual G4bool StopFragmenting(const G4FragmentingString * const string); 58 virtual G4bool IsItFragmentable(const G4Fr << 63 virtual G4bool IsFragmentable(const G4FragmentingString * const string); 59 64 60 virtual G4bool SplitLast(G4FragmentingStri 65 virtual G4bool SplitLast(G4FragmentingString * string, 61 G4KineticTrackVector * LeftVe 66 G4KineticTrackVector * LeftVector, 62 G4KineticTrackVector * RightV << 67 G4KineticTrackVector * RightVector); 63 68 64 virtual void Sample4Momentum(G4LorentzVect 69 virtual void Sample4Momentum(G4LorentzVector* Mom, G4double Mass, 65 G4LorentzVect 70 G4LorentzVector* AntiMom, G4double AntiMass, 66 G4double Init 71 G4double InitialMass); 67 72 68 virtual G4KineticTrack * Splitup(G4Fragmen << 69 G4Fragmen << 70 << 71 virtual G4LorentzVector * SplitEandP(G4Par 73 virtual G4LorentzVector * SplitEandP(G4ParticleDefinition * pHadron, 72 G4Fra 74 G4FragmentingString * string, 73 G4Fra 75 G4FragmentingString * newString); 74 76 75 virtual G4double GetLightConeZ(G4double zm 77 virtual G4double GetLightConeZ(G4double zmin, G4double zmax, 76 G4int Parto 78 G4int PartonEncoding, 77 G4ParticleD 79 G4ParticleDefinition* pHadron, 78 G4double Px 80 G4double Px, G4double Py); 79 81 80 G4double lambda(G4double s, G4double m1_Sq 82 G4double lambda(G4double s, G4double m1_Sqr, G4double m2_Sqr); 81 83 82 // virtual G4ParticleDefinition * QuarkSpl << 84 virtual G4ParticleDefinition * DiQuarkSplitup(G4ParticleDefinition* decay, 83 // << 84 << 85 virtual G4ParticleDefinition * DiQuarkSpli << 86 85 G4ParticleDefinition *&created); 87 86 88 private: 87 private: 89 G4bool Loop_toFragmentString(const G4Excit << 88 // Internal methods introduced to improve the code structure (AR Nov 2011) >> 89 >> 90 G4bool Loop_toFragmentString(G4ExcitedString * & theStringInCMS, 90 G4KineticTrac 91 G4KineticTrackVector * & LeftVector, 91 G4KineticTrac 92 G4KineticTrackVector * & RightVector); 92 93 93 G4bool Diquark_AntiDiquark_belowThreshold_ 94 G4bool Diquark_AntiDiquark_belowThreshold_lastSplitting(G4FragmentingString * & string, 94 95 G4ParticleDefinition * & LeftHadron, 95 96 G4ParticleDefinition * & RightHadron); 96 97 97 G4bool Diquark_AntiDiquark_aboveThreshold_ 98 G4bool Diquark_AntiDiquark_aboveThreshold_lastSplitting(G4FragmentingString * & string, 98 99 G4ParticleDefinition * & LeftHadron, 99 100 G4ParticleDefinition * & RightHadron); 100 101 101 G4bool Quark_AntiQuark_lastSplitting(G4Fra 102 G4bool Quark_AntiQuark_lastSplitting(G4FragmentingString * & string, 102 G4Par 103 G4ParticleDefinition * & LeftHadron, 103 G4Par 104 G4ParticleDefinition * & RightHadron); 104 105 105 G4bool Quark_Diquark_lastSplitting(G4Fragm 106 G4bool Quark_Diquark_lastSplitting(G4FragmentingString * & string, 106 G4Parti 107 G4ParticleDefinition * & LeftHadron, 107 G4Parti 108 G4ParticleDefinition * & RightHadron ); 108 109 109 G4int SampleState(void); << 110 G4int SampleState(void); 110 << 111 G4double Tmt; // "Temperature" for sampli << 112 111 >> 112 private: >> 113 // ------ For estimation of a minimal string mass --------------- >> 114 G4double Mass_of_light_quark; >> 115 G4double Mass_of_heavy_quark; >> 116 G4double Mass_of_string_junction; >> 117 // ------ An estimated minimal string mass ---------------------- >> 118 G4double MinimalStringMass; >> 119 G4double MinimalStringMass2; >> 120 // ------ Minimal invariant mass used at a string fragmentation - >> 121 G4double WminLUND; >> 122 >> 123 G4int Meson[3][3][6]; >> 124 G4double MesonWeight[3][3][6]; >> 125 >> 126 G4int Baryon[3][3][3][4]; >> 127 G4double BaryonWeight[3][3][3][4]; >> 128 >> 129 G4double Prob_QQbar[3]; >> 130 >> 131 // ------ To improve the code structure >> 132 G4ParticleDefinition * FS_LeftHadron[35], * FS_RightHadron[35]; >> 133 G4double FS_Weight[35]; >> 134 G4int NumberOf_FS; 113 }; 135 }; 114 136 115 //******************************************** << 137 //************************************************************************************************************** 116 // Class G4LundStringFragmentation 138 // Class G4LundStringFragmentation 117 139 118 #endif 140 #endif 119 141 120 142