Geant4  10.02.p01
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: G4VLongitudinalStringDecay.hh 84219 2014-10-10 14:55:19Z gcosmo $
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
38 #include "G4DynamicParticle.hh"
39 #include "G4KineticTrack.hh"
40 #include "G4KineticTrackVector.hh"
41 #include "G4HadronBuilder.hh"
42 
44 //**************************************************************************************
45 
47  {
48 public:
51 
52 private:
53  // not implemented to protect/forbid use
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
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
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 
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 
106  G4FragmentingString *&newString);
107 
109  G4ParticleDefinition *&created);
110 
111  virtual G4ParticleDefinition * DiQuarkSplitup(G4ParticleDefinition* decay, // Uzhi June 2014
112  G4ParticleDefinition *&created)=0;
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,
139 
140 // Used for some test purposes ------------------------------------------------
141  void ConstructParticle();
142 
144  G4bool theGivenSpin, G4int theSpin);
145 
146 //-----------------------------------------------------------------------------
147 public:
148 
150 
152  void SetStrangenessSuppression(G4double aValue);
153  void SetDiquarkSuppression(G4double aValue);
155 
156  void SetVectorMesonProbability(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:
171 
172  G4double GetStringTensionParameter() {return Kappa;}; // Uzhi 20 June 08
173 
174 //private:
175 protected:
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
186 
189  std::vector<G4double> vectorMesonMix;
190  std::vector<G4double> scalarMesonMix;
191 
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
int operator!=(const G4VLongitudinalStringDecay &right) const
G4ParticleDefinition *(G4HadronBuilder::* Pcreate)(G4ParticleDefinition *, G4ParticleDefinition *)
G4ParticleDefinition * QuarkSplitup(G4ParticleDefinition *decay, G4ParticleDefinition *&created)
virtual G4double GetLightConeZ(G4double zmin, G4double zmax, G4int PartonEncoding, G4ParticleDefinition *pHadron, G4double Px, G4double Py)=0
CLHEP::Hep3Vector G4ThreeVector
G4ExcitedString * CPExcited(const G4ExcitedString &string)
void SetStrangenessSuppression(G4double aValue)
virtual G4KineticTrackVector * FragmentString(const G4ExcitedString &theString)=0
G4ThreeVector SampleQuarkPt(G4double ptMax=-1.)
std::vector< G4double > vectorMesonMix
G4double FragmentationMass(const G4FragmentingString *const string, Pcreate build=0, pDefPair *pdefs=0)
virtual void SetMassCut(G4double aValue)
int G4int
Definition: G4Types.hh:78
virtual G4bool StopFragmenting(const G4FragmentingString *const string)=0
G4KineticTrackVector * DecayResonans(G4KineticTrackVector *aHadrons)
pDefPair CreatePartonPair(G4int NeedParticle, G4bool AllowDiquarks=true)
void SetStringTensionParameter(G4double aValue)
const G4VLongitudinalStringDecay & operator=(const G4VLongitudinalStringDecay &right)
void SetDiquarkBreakProbability(G4double aValue)
G4ParticleDefinition * CreateHadron(G4int id1, G4int id2, G4bool theGivenSpin, G4int theSpin)
ParticleList decay(Cluster *const c)
Carries out a cluster decay.
G4KineticTrackVector * LightFragmentationTest(const G4ExcitedString *const theString)
std::pair< G4ParticleDefinition *, G4ParticleDefinition * > pDefPair
bool G4bool
Definition: G4Types.hh:79
G4ParticleDefinition * FindParticle(G4int Encoding)
virtual void Sample4Momentum(G4LorentzVector *Mom, G4double Mass, G4LorentzVector *AntiMom, G4double AntiMass, G4double InitialMass)=0
std::vector< G4double > scalarMesonMix
virtual G4bool IsFragmentable(const G4FragmentingString *const string)=0
int operator==(const G4VLongitudinalStringDecay &right) const
virtual G4bool SplitLast(G4FragmentingString *string, G4KineticTrackVector *LeftVector, G4KineticTrackVector *RightVector)=0
void SetVectorMesonProbability(G4double aValue)
void SetVectorMesonMixings(std::vector< G4double > aVector)
void SetSpinThreeHalfBarionProbability(G4double aValue)
virtual G4ParticleDefinition * DiQuarkSplitup(G4ParticleDefinition *decay, G4ParticleDefinition *&created)=0
double G4double
Definition: G4Types.hh:76
void SetScalarMesonMixings(std::vector< G4double > aVector)
void CalculateHadronTimePosition(G4double theInitialStringMass, G4KineticTrackVector *)
virtual G4LorentzVector * SplitEandP(G4ParticleDefinition *pHadron, G4FragmentingString *string, G4FragmentingString *newString)=0
G4KineticTrack * Splitup(G4FragmentingString *string, G4FragmentingString *&newString)
void SetDiquarkSuppression(G4double aValue)
CLHEP::HepLorentzVector G4LorentzVector