Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4VLongitudinalStringDecay.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 // $Id$
28 // Maxim Komogorov
29 //
30 // -----------------------------------------------------------------------------
31 // GEANT 4 class implementation file
32 //
33 // History: first implementation, Maxim Komogorov, 1-Jul-1998
34 // -----------------------------------------------------------------------------
35 #ifndef G4VLongitudinalStringDecay_h
36 #define G4VLongitudinalStringDecay_h 1
37 #include "G4VStringFragmentation.hh"
38 #include "G4DynamicParticle.hh"
39 #include "G4KineticTrack.hh"
40 #include "G4KineticTrackVector.hh"
41 #include "G4HadronBuilder.hh"
42 
43 class G4FragmentingString;
44 //**************************************************************************************
45 
46 class G4VLongitudinalStringDecay
47  {
48 public:
49  G4VLongitudinalStringDecay();
50  virtual ~G4VLongitudinalStringDecay();
51 
52 private:
53  // not implemented to protect/forbid use
54  G4VLongitudinalStringDecay(const G4VLongitudinalStringDecay &right);
55  const G4VLongitudinalStringDecay & operator=(const G4VLongitudinalStringDecay &right);
56  int operator==(const G4VLongitudinalStringDecay &right) const;
57  int operator!=(const G4VLongitudinalStringDecay &right) const;
58 
59 public:
60  virtual G4KineticTrackVector* FragmentString(const G4ExcitedString& theString)=0;
61 
62 protected:
63 
64 // For changing Mass Cut used for selection of very small mass strings
65  virtual void SetMassCut(G4double aValue);
66 
67 // For handling a string with very low mass
68  G4KineticTrackVector * LightFragmentationTest(const G4ExcitedString * const theString);
69 
70 // To store created quarks or 2 last hadrons
71  typedef std::pair<G4ParticleDefinition*, G4ParticleDefinition*> pDefPair;
72 
73 // For creation of hadrons from given quark pair
74  typedef G4ParticleDefinition * (G4HadronBuilder::*Pcreate)
75  (G4ParticleDefinition*, G4ParticleDefinition*);
76 
77 //-----------------------------------------------------------------------------
78 // Used by LightFragmentationTest for estimation of lowest possible mass of
79 // given quark system
80  G4double FragmentationMass(const G4FragmentingString * const string,
81  Pcreate build=0,
82  pDefPair * pdefs=0);
83 
84  G4ParticleDefinition* FindParticle(G4int Encoding);
85 
86  virtual void Sample4Momentum(G4LorentzVector* Mom, G4double Mass,
87  G4LorentzVector* AntiMom, G4double AntiMass,
88  G4double InitialMass)=0;
89 //-----------------------------------------------------------------------------
90 // For decision on continue or stop string fragmentation
91  virtual G4bool StopFragmenting(const G4FragmentingString * const string)=0;
92  virtual G4bool IsFragmentable(const G4FragmentingString * const string)=0;
93 
94 // If a string can not fragment, make last break into 2 hadrons
95  virtual G4bool SplitLast(G4FragmentingString * string,
96  G4KineticTrackVector * LeftVector,
97  G4KineticTrackVector * RightVector)=0;
98 //-----------------------------------------------------------------------------
99 
100 // If a string fragments, do the following
101 
102 // For transver of a string to its CMS frame
103  G4ExcitedString *CPExcited(const G4ExcitedString& string);
104 
105  G4KineticTrack * Splitup(G4FragmentingString *string,
106  G4FragmentingString *&newString);
107 
108  G4ParticleDefinition * QuarkSplitup(G4ParticleDefinition* decay,
109  G4ParticleDefinition *&created);
110 
111  G4ParticleDefinition * DiQuarkSplitup(G4ParticleDefinition* decay,
112  G4ParticleDefinition *&created);
113 
114  pDefPair CreatePartonPair(G4int NeedParticle, G4bool AllowDiquarks=true);
115 
116 public:
117 // used by G4VKinkyStringDecy..
118  G4int SampleQuarkFlavor(void);
119  G4ThreeVector SampleQuarkPt(G4double ptMax=-1.); // -1. no limit on maxpt.
120 
121 protected:
122 
123 //-----------------------------------------------------------------------------
124 // For determination of kinematical properties of created hadron
125 // virtual G4LorentzVector * SplitEandP(G4ParticleDefinition * pHadron, // Uzhi
126 // G4FragmentingString * string )=0; // Uzhi
127 
128  virtual G4LorentzVector * SplitEandP(G4ParticleDefinition * pHadron, // Uzhi
129  G4FragmentingString * string, // Uzhi
130  G4FragmentingString * newString )=0;// Uzhi
131 
132  virtual G4double GetLightConeZ(G4double zmin, G4double zmax,
133  G4int PartonEncoding,
134  G4ParticleDefinition* pHadron,
135  G4double Px, G4double Py ) = 0;
136 
137  void CalculateHadronTimePosition(G4double theInitialStringMass,
138  G4KineticTrackVector *);
139 
140 // Used for some test purposes ------------------------------------------------
141  void ConstructParticle();
142 
143  G4ParticleDefinition* CreateHadron(G4int id1, G4int id2,
144  G4bool theGivenSpin, G4int theSpin);
145 
146 //-----------------------------------------------------------------------------
147 public:
148 
149  G4KineticTrackVector* DecayResonans (G4KineticTrackVector* aHadrons);
150 
151  void SetSigmaTransverseMomentum(G4double aQT);
152  void SetStrangenessSuppression(G4double aValue);
153  void SetDiquarkSuppression(G4double aValue);
154  void SetDiquarkBreakProbability(G4double aValue);
155 
156  void SetVectorMesonProbability(G4double aValue);
157  void SetSpinThreeHalfBarionProbability(G4double aValue);
158 
159  void SetScalarMesonMixings( std::vector<G4double> aVector);
160  void SetVectorMesonMixings( std::vector<G4double> aVector);
161 
162  void SetStringTensionParameter(G4double aValue); // Uzhi 20 June 08
163 
164 //private:
165 protected:
166  G4double GetDiquarkSuppress() {return DiquarkSuppress;};
167  G4double GetDiquarkBreakProb() {return DiquarkBreakProb;};
168  G4double GetStrangeSuppress() {return StrangeSuppress;};
169  G4double GetClusterMass() {return ClusterMass;};
170  G4int GetClusterLoopInterrupt() {return ClusterLoopInterrupt;};
171 
172  G4double GetStringTensionParameter() {return Kappa;}; // Uzhi 20 June 08
173 
174 //private:
175 protected:
176  G4double MassCut;
177  G4double ClusterMass;
178  G4double SigmaQT; // sigma_q_t is quark transverse momentum distribution parameter
179  G4double DiquarkSuppress; // is Diquark suppression parameter
180  G4double DiquarkBreakProb; // is Diquark breaking probability
181  G4double SmoothParam; // model parameter
182  G4double StrangeSuppress ;
183  G4int StringLoopInterrupt;
184  G4int ClusterLoopInterrupt;
185  G4HadronBuilder *hadronizer;
186 
187  G4double pspin_meson;
188  G4double pspin_barion;
189  std::vector<G4double> vectorMesonMix;
190  std::vector<G4double> scalarMesonMix;
191 
192  G4bool PastInitPhase;
193 
194  G4double Kappa; // String tension parameter // Uzhi 20 June 08
195 
196 // G4double MinFragmentationMass(G4ExcitedString * theString,
197 // G4ParticleDefinition*& Hadron1,
198 // G4ParticleDefinition*& Hadron2);
199 
200 };
201 
202 //*************************************************************************************
203 // Class G4VLongitudinalStringDecay
204 #endif
Generated on Sat May 25 2013 14:34:05 for Geant4 by   doxygen 1.8.4