Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4ExcitedStringDecay Class Reference

#include <G4ExcitedStringDecay.hh>

+ Inheritance diagram for G4ExcitedStringDecay:

Public Member Functions

 G4ExcitedStringDecay (G4VLongitudinalStringDecay *aStringDecay=nullptr)
 
virtual ~G4ExcitedStringDecay ()
 
virtual G4KineticTrackVector * FragmentStrings (const G4ExcitedStringVector *theStrings)
 
- Public Member Functions inherited from G4VStringFragmentation
 G4VStringFragmentation (const G4String &name="StringFragmentation")
 
 ~G4VStringFragmentation () override
 
 G4VStringFragmentation (const G4VStringFragmentation &right)=delete
 
const G4VStringFragmentation & operator= (const G4VStringFragmentation &right)=delete
 
G4bool operator== (const G4VStringFragmentation &right) const =delete
 
G4bool operator!= (const G4VStringFragmentation &right) const =delete
 
virtual G4KineticTrackVector * FragmentStrings (const G4ExcitedStringVector *theStrings)=0
 
- Public Member Functions inherited from G4HadronicInteraction
 G4HadronicInteraction (const G4String &modelName="HadronicModel")
 
virtual ~G4HadronicInteraction ()
 
virtual G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
 
virtual G4double SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
 
virtual G4bool IsApplicable (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
 
G4double GetMinEnergy () const
 
G4double GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMinEnergy (G4double anEnergy)
 
void SetMinEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
 
G4double GetMaxEnergy () const
 
G4double GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMaxEnergy (const G4double anEnergy)
 
void SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
 
G4int GetVerboseLevel () const
 
void SetVerboseLevel (G4int value)
 
const G4String & GetModelName () const
 
void DeActivateFor (const G4Material *aMaterial)
 
void ActivateFor (const G4Material *aMaterial)
 
void DeActivateFor (const G4Element *anElement)
 
void ActivateFor (const G4Element *anElement)
 
G4bool IsBlocked (const G4Material *aMaterial) const
 
G4bool IsBlocked (const G4Element *anElement) const
 
void SetRecoilEnergyThreshold (G4double val)
 
G4double GetRecoilEnergyThreshold () const
 
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels () const
 
virtual std::pair< G4double, G4double > GetEnergyMomentumCheckLevels () const
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
virtual void ModelDescription (std::ostream &outFile) const
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void InitialiseModel ()
 
 G4HadronicInteraction (const G4HadronicInteraction &right)=delete
 
const G4HadronicInteraction & operator= (const G4HadronicInteraction &right)=delete
 
G4bool operator== (const G4HadronicInteraction &right) const =delete
 
G4bool operator!= (const G4HadronicInteraction &right) const =delete
 

Additional Inherited Members

- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 38 of file G4ExcitedStringDecay.hh.

Constructor & Destructor Documentation

◆ G4ExcitedStringDecay()

G4ExcitedStringDecay::G4ExcitedStringDecay ( G4VLongitudinalStringDecay *  aStringDecay = nullptr)

Definition at line 39 of file G4ExcitedStringDecay.cc.

40 : G4VStringFragmentation(), theStringDecay(ptr)
41{
42 if(!ptr) {
43 G4HadronicInteraction* p =
44 G4HadronicInteractionRegistry::Instance()->FindModel("LundStringFragmentation");
45 theStringDecay = static_cast<G4VLongitudinalStringDecay*>(p);
46 if(!theStringDecay) { theStringDecay = new G4LundStringFragmentation(); }
47 }
48 SetModelName(theStringDecay->GetModelName());
49}
G4HadronicInteractionRegistry::FindModel
G4HadronicInteraction * FindModel(const G4String &name)
Definition: G4HadronicInteractionRegistry.cc:110
G4HadronicInteractionRegistry::Instance
static G4HadronicInteractionRegistry * Instance()
Definition: G4HadronicInteractionRegistry.cc:36
G4HadronicInteraction::SetModelName
void SetModelName(const G4String &nam)
Definition: G4HadronicInteraction.hh:166
G4HadronicInteraction::GetModelName
const G4String & GetModelName() const
Definition: G4HadronicInteraction.hh:115

◆ ~G4ExcitedStringDecay()

G4ExcitedStringDecay::~G4ExcitedStringDecay ( )
virtual

Definition at line 51 of file G4ExcitedStringDecay.cc.

52{}

Member Function Documentation

◆ FragmentStrings()

G4KineticTrackVector * G4ExcitedStringDecay::FragmentStrings ( const G4ExcitedStringVector *  theStrings)
virtual

Implements G4VStringFragmentation.

Definition at line 59 of file G4ExcitedStringDecay.cc.

60{
61 G4LorentzVector KTsum(0.,0.,0.,0.);
62
63 #ifdef debug_G4ExcitedStringDecay
64 G4cout<<G4endl;
65 G4cout<<"--------------------------- G4ExcitedStringDecay ----------------------"<<G4endl;
66 G4cout<<"Hadronization of Excited Strings: theStrings->size() "<<theStrings->size()<<G4endl;
67 #endif
68
69 for ( unsigned int astring=0; astring < theStrings->size(); astring++)
70 // for ( unsigned int astring=0; astring < 1; astring++)
71 {
72 // G4cout<<"theStrings->operator[](astring)->IsExcited() "<<" "<<astring<<" "<<theStrings->operator[](astring)->IsExcited()<<G4endl;
73 if ( theStrings->operator[](astring)->IsExcited() )
74 {KTsum+= theStrings->operator[](astring)->Get4Momentum();}
75 else {KTsum+=theStrings->operator[](astring)->GetKineticTrack()->Get4Momentum();}
76 }
77
78 G4LorentzRotation toCms( -1 * KTsum.boostVector() );
79 G4LorentzRotation toLab(toCms.inverse());
80 G4LorentzVector Ptmp;
81 KTsum=G4LorentzVector(0.,0.,0.,0.);
82
83 for ( unsigned int astring=0; astring < theStrings->size(); astring++)
84 // for ( unsigned int astring=0; astring < 1; astring++)
85 {
86 if ( theStrings->operator[](astring)->IsExcited() )
87 {
88 Ptmp=toCms * theStrings->operator[](astring)->GetLeftParton()->Get4Momentum();
89 theStrings->operator[](astring)->GetLeftParton()->Set4Momentum(Ptmp);
90
91 Ptmp=toCms * theStrings->operator[](astring)->GetRightParton()->Get4Momentum();
92 theStrings->operator[](astring)->GetRightParton()->Set4Momentum(Ptmp);
93
94 KTsum+= theStrings->operator[](astring)->Get4Momentum();
95 }
96 else
97 {
98 Ptmp=toCms * theStrings->operator[](astring)->GetKineticTrack()->Get4Momentum();
99 theStrings->operator[](astring)->GetKineticTrack()->Set4Momentum(Ptmp);
100 KTsum+= theStrings->operator[](astring)->GetKineticTrack()->Get4Momentum();
101 }
102 }
103
104 G4SampleResonance BrW;
105 const G4ParticleDefinition* TrackDefinition=0;
106
107 G4KineticTrackVector * theResult = new G4KineticTrackVector;
108 G4int attempts(0);
109 G4bool success=false;
110 G4bool NeedEnergyCorrector=false;
111 do {
112 #ifdef debug_G4ExcitedStringDecay
113 G4cout<<"New try No "<<attempts<<" to hadronize strings"<<G4endl;
114 #endif
115
116 std::for_each(theResult->begin() , theResult->end() , DeleteKineticTrack());
117 theResult->clear();
118
119 attempts++;
120
121 G4LorentzVector KTsecondaries(0.,0.,0.,0.);
122 NeedEnergyCorrector=false;
123
124 for ( unsigned int astring=0; astring < theStrings->size(); astring++)
125 // for ( unsigned int astring=0; astring < 1; astring++) // Proj. Last Str. Decay for FTF
126 // for ( unsigned int astring=1; astring < 2; astring++) // Proj. Last Str. Decay for QGS
127 // for ( unsigned int astring=0; astring < 1; astring++) // For testing purposes
128 {
129 #ifdef debug_G4ExcitedStringDecay
130 G4cout<<"String No "<<astring+1<<" Excited? "<<theStrings->operator[](astring)->IsExcited()<<G4endl;
131
132 G4cout<<"String No "<<astring+1<<" 4Momentum "<<theStrings->operator[](astring)->Get4Momentum()
133 <<" "<<theStrings->operator[](astring)->Get4Momentum().mag()<<G4endl;
134 #endif
135
136 G4KineticTrackVector * generatedKineticTracks = NULL;
137 if ( theStrings->operator[](astring)->IsExcited() )
138 {
139 #ifdef debug_G4ExcitedStringDecay
140 G4cout<<"Fragment String with partons: "
141 <<theStrings->operator[](astring)->GetLeftParton()->GetPDGcode() <<" "
142 <<theStrings->operator[](astring)->GetRightParton()->GetPDGcode()<<" "
143 <<"Direction "<<theStrings->operator[](astring)->GetDirection()<<G4endl;
144 #endif
145 generatedKineticTracks=FragmentString(*theStrings->operator[](astring));
146 #ifdef debug_G4ExcitedStringDecay
147 G4cout<<"(G4ExcitedStringDecay) Number of produced hadrons = "
148 <<generatedKineticTracks->size()<<G4endl;
149 #endif
150 } else {
151 #ifdef debug_G4ExcitedStringDecay
152 G4cout<<" GetTrack from the String"<<G4endl;
153 #endif
154 G4LorentzVector Mom=theStrings->operator[](astring)->GetKineticTrack()->Get4Momentum();
155 G4KineticTrack * aTrack= new G4KineticTrack(
156 theStrings->operator[](astring)->GetKineticTrack()->GetDefinition(),
157 theStrings->operator[](astring)->GetKineticTrack()->GetFormationTime(),
158 G4ThreeVector(0), Mom);
159
160 aTrack->SetPosition(theStrings->operator[](astring)->GetKineticTrack()->GetPosition());
161
162 #ifdef debug_G4ExcitedStringDecay
163 G4cout<<" A particle stored in the track is "<<aTrack->GetDefinition()->GetParticleName()<<G4endl;
164 #endif
165
166 generatedKineticTracks = new G4KineticTrackVector;
167 generatedKineticTracks->push_back(aTrack);
168 }
169
170 if (generatedKineticTracks == nullptr || generatedKineticTracks->size() == 0)
171 {
172 // G4cerr << "G4VPartonStringModel:No KineticTracks produced" << G4endl;
173 continue;
174 }
175
176 G4LorentzVector KTsum1(0.,0.,0.,0.);
177 for ( unsigned int aTrack=0; aTrack<generatedKineticTracks->size();aTrack++)
178 {
179 #ifdef debug_G4ExcitedStringDecay
180 G4cout<<"Prod part No. "<<aTrack+1<<" "
181 <<(*generatedKineticTracks)[aTrack]->GetDefinition()->GetParticleName()<<" "
182 <<(*generatedKineticTracks)[aTrack]->Get4Momentum()
183 <<(*generatedKineticTracks)[aTrack]->Get4Momentum().mag()<<G4endl;
184 #endif
185 // --------------- Sampling mass of unstable hadronic resonances ----------------
186 TrackDefinition = (*generatedKineticTracks)[aTrack]->GetDefinition();
187
188 if (TrackDefinition->IsShortLived())
189 {
190 G4double NewTrackMass =
191 BrW.SampleMass( TrackDefinition->GetPDGMass(), TrackDefinition->GetPDGWidth(),
192 BrW.GetMinimumMass( TrackDefinition ) + 10.0*MeV,
193 TrackDefinition->GetPDGMass() + 5.0*TrackDefinition->GetPDGWidth() );
194 G4LorentzVector Tmp=G4LorentzVector((*generatedKineticTracks)[aTrack]->Get4Momentum());
195 Tmp.setE(std::sqrt(sqr(NewTrackMass) + Tmp.vect().mag2()));
196
197 (*generatedKineticTracks)[aTrack]->Set4Momentum(Tmp);
198
199 #ifdef debug_G4ExcitedStringDecay
200 G4cout<<"Resonance *** "<<aTrack+1<<" "
201 <<(*generatedKineticTracks)[aTrack]->GetDefinition()->GetParticleName()<<" "
202 <<(*generatedKineticTracks)[aTrack]->Get4Momentum()
203 <<(*generatedKineticTracks)[aTrack]->Get4Momentum().mag()<<G4endl;
204 #endif
205 }
206 //-------------------------------------------------------------------------------
207
208 theResult->push_back(generatedKineticTracks->operator[](aTrack));
209 KTsum1+= (*generatedKineticTracks)[aTrack]->Get4Momentum();
210 }
211 KTsecondaries+=KTsum1;
212
213 #ifdef debug_G4ExcitedStringDecay
214 G4cout << "String secondaries(" <<generatedKineticTracks->size()<< ")"<<G4endl
215 <<"Init string momentum: "<< theStrings->operator[](astring)->Get4Momentum()<<G4endl
216 <<"Final hadrons momentum: "<< KTsum1 << G4endl;
217 #endif
218
219 if ( KTsum1.e() > 0 &&
220 std::abs((KTsum1.e()-theStrings->operator[](astring)->Get4Momentum().e()) / KTsum1.e()) > perMillion )
221 {
222 NeedEnergyCorrector=true;
223 }
224
225 #ifdef debug_G4ExcitedStringDecay
226 G4cout<<"NeedEnergyCorrection yes/no "<<NeedEnergyCorrector<<G4endl;
227 #endif
228
229 // clean up
230 delete generatedKineticTracks;
231 success=true;
232 }
233
234 if ( NeedEnergyCorrector ) success=EnergyAndMomentumCorrector(theResult, KTsum);
235 } while (!success && (attempts < 100)); /* Loop checking, 07.08.2015, A.Ribon */
236
237 for ( unsigned int aTrack=0; aTrack<theResult->size();aTrack++)
238 {
239 Ptmp=(*theResult)[aTrack]->Get4Momentum();
240 Ptmp.transform( toLab);
241 (*theResult)[aTrack]->Set4Momentum(Ptmp);
242 }
243
244 #ifdef debug_G4ExcitedStringDecay
245 G4cout<<"End of the strings fragmentation (G4ExcitedStringDecay)"<<G4endl;
246
247 G4LorentzVector KTsum1(0.,0.,0.,0.);
248
249 for ( unsigned int aTrack=0; aTrack<theResult->size();aTrack++)
250 {
251 G4cout << " corrected tracks .. " << (*theResult)[aTrack]->GetDefinition()->GetParticleName()
252 <<" " << (*theResult)[aTrack]->Get4Momentum()
253 <<" " << (*theResult)[aTrack]->Get4Momentum().mag()<< G4endl;
254 KTsum1+= (*theResult)[aTrack]->Get4Momentum();
255 }
256
257 G4cout << "Needcorrector/success " << NeedEnergyCorrector << "/" << success
258 << ", Corrected total 4 momentum " << KTsum1 << G4endl;
259 if ( ! success ) G4cout << "failed to correct E/p" << G4endl;
260
261 G4cout<<"End of the Hadronization (G4ExcitedStringDecay)"<<G4endl;
262 #endif
263
264 if (!success)
265 {
266 if (theResult->size() != 0)
267 {
268 std::for_each(theResult->begin() , theResult->end() , DeleteKineticTrack());
269 theResult->clear();
270 delete theResult; theResult=0;
271 }
272 for ( unsigned int astring=0; astring < theStrings->size(); astring++)
273 // for ( unsigned int astring=0; astring < 1; astring++) // Need more correct. For testing purposes.
274 {
275 if ( theStrings->operator[](astring)->IsExcited() )
276 {
277 Ptmp=theStrings->operator[](astring)->GetLeftParton()->Get4Momentum();
278 Ptmp.transform( toLab);
279 theStrings->operator[](astring)->GetLeftParton()->Set4Momentum(Ptmp);
280
281 Ptmp=theStrings->operator[](astring)->GetRightParton()->Get4Momentum();
282 Ptmp.transform( toLab);
283 theStrings->operator[](astring)->GetRightParton()->Set4Momentum(Ptmp);
284 }
285 else
286 {
287 Ptmp=theStrings->operator[](astring)->GetKineticTrack()->Get4Momentum();
288 Ptmp.transform( toLab);
289 theStrings->operator[](astring)->GetKineticTrack()->Set4Momentum(Ptmp);
290 }
291 }
292 }
293 return theResult;
294}
G4LorentzVector
CLHEP::HepLorentzVector G4LorentzVector
Definition: G4LorentzVector.hh:33
G4double
double G4double
Definition: G4Types.hh:83
G4bool
bool G4bool
Definition: G4Types.hh:86
G4int
int G4int
Definition: G4Types.hh:85
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
CLHEP::Hep3Vector::mag2
double mag2() const
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
CLHEP::HepLorentzVector::transform
HepLorentzVector & transform(const HepRotation &)
G4KineticTrack::SetPosition
void SetPosition(const G4ThreeVector aPosition)
Definition: G4KineticTrack.hh:246
G4KineticTrack::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4KineticTrack.hh:221
G4ParticleDefinition::GetPDGWidth
G4double GetPDGWidth() const
Definition: G4ParticleDefinition.hh:112
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:109
G4SampleResonance::SampleMass
G4double SampleMass(const G4double poleMass, const G4double gamma, const G4double minMass, const G4double maxMass) const
Definition: G4SampleResonance.cc:138
G4SampleResonance::GetMinimumMass
G4double GetMinimumMass(const G4ParticleDefinition *p) const
Definition: G4SampleResonance.cc:46
sqr
T sqr(const T &x)
Definition: templates.hh:128
The documentation for this class was generated from the following files: