Geant4  10.03.p03
 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: G4VLongitudinalStringDecay.hh 100828 2016-11-02 15:25:59Z 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
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
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  // Used by LightFragmentationTest for estimation of lowest possible mass of
78  // given quark system
79  G4double FragmentationMass(const G4FragmentingString * const string,
80  Pcreate build=0, pDefPair * pdefs=0);
81 
83 
84  virtual void Sample4Momentum(G4LorentzVector* Mom, G4double Mass,
85  G4LorentzVector* AntiMom, G4double AntiMass,
86  G4double InitialMass)=0;
87 
88  // For decision on continue or stop string fragmentation
89  virtual G4bool StopFragmenting(const G4FragmentingString * const string)=0;
90  virtual G4bool IsFragmentable(const G4FragmentingString * const string)=0;
91 
92  // If a string can not fragment, make last break into 2 hadrons
93  virtual G4bool SplitLast(G4FragmentingString * string,
94  G4KineticTrackVector * LeftVector,
95  G4KineticTrackVector * RightVector)=0;
96 
97  // If a string fragments, do the following
98 
99  // For transver of a string to its CMS frame
100  G4ExcitedString *CPExcited(const G4ExcitedString& string);
101 
103  G4FragmentingString *&newString);
104 
106  G4ParticleDefinition *&created);
107 
109  G4ParticleDefinition *&created)=0;
110 
111  pDefPair CreatePartonPair(G4int NeedParticle, G4bool AllowDiquarks=true);
112 
113  public:
114  // used by G4VKinkyStringDecy.
115  G4int SampleQuarkFlavor(void);
116  G4ThreeVector SampleQuarkPt(G4double ptMax=-1.); // -1. no limit on maxpt.
117 
118  protected:
119 
120  // For determination of kinematical properties of created hadron
121  // virtual G4LorentzVector * SplitEandP(G4ParticleDefinition * pHadron,
122  // G4FragmentingString * string )=0;
123 
124  virtual G4LorentzVector * SplitEandP(G4ParticleDefinition * pHadron,
125  G4FragmentingString * string,
126  G4FragmentingString * newString )=0;
127 
128  virtual G4double GetLightConeZ(G4double zmin, G4double zmax,
129  G4int PartonEncoding,
130  G4ParticleDefinition* pHadron,
131  G4double Px, G4double Py ) = 0;
132 
133  void CalculateHadronTimePosition(G4double theInitialStringMass,
135 
136  // Used for some test purposes
137  void ConstructParticle();
138 
140  G4bool theGivenSpin, G4int theSpin);
141 
142  public:
143 
145 
147  void SetStrangenessSuppression(G4double aValue);
148  void SetDiquarkSuppression(G4double aValue);
150 
151  void SetVectorMesonProbability(G4double aValue);
153 
154  void SetScalarMesonMixings( std::vector<G4double> aVector);
155  void SetVectorMesonMixings( std::vector<G4double> aVector);
156 
157  void SetStringTensionParameter(G4double aValue);
158 
159  protected:
165 
167 
168  protected:
171  G4double SigmaQT; // sigma_q_t is quark transverse momentum distribution parameter
172  G4double DiquarkSuppress; // is Diquark suppression parameter
173  G4double DiquarkBreakProb; // is Diquark breaking probability
174  G4double SmoothParam; // model parameter
179 
182  std::vector<G4double> vectorMesonMix;
183  std::vector<G4double> scalarMesonMix;
184 
186 
187  G4double Kappa; // String tension parameter
188 
189  //G4double MinFragmentationMass(G4ExcitedString * theString,
190  // G4ParticleDefinition*& Hadron1,
191  // G4ParticleDefinition*& Hadron2);
192 
193 };
194 
195 // Class G4VLongitudinalStringDecay
196 
197 #endif
198 
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
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)
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
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)