Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | List of all members
G4VPartonStringModel Class Referenceabstract

#include <G4VPartonStringModel.hh>

+ Inheritance diagram for G4VPartonStringModel:

Public Member Functions

 G4VPartonStringModel (const G4String &modelName="Parton String Model")
 
 ~G4VPartonStringModel () override
 
 G4VPartonStringModel (const G4VPartonStringModel &right)=delete
 
const G4VPartonStringModel & operator= (const G4VPartonStringModel &right)=delete
 
G4bool operator== (const G4VPartonStringModel &right) const =delete
 
G4bool operator!= (const G4VPartonStringModel &right) const =delete
 
void SetFragmentationModel (G4VStringFragmentation *aModel)
 
G4KineticTrackVector * Scatter (const G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary) override
 
void ModelDescription (std::ostream &outFile) const override
 
G4V3DNucleus * GetProjectileNucleus () const override
 
- Public Member Functions inherited from G4VHighEnergyGenerator
 G4VHighEnergyGenerator (const G4String &modelName="HighEnergyGenerator")
 
 ~G4VHighEnergyGenerator () override
 
 G4VHighEnergyGenerator (const G4VHighEnergyGenerator &right)=delete
 
const G4VHighEnergyGenerator & operator= (const G4VHighEnergyGenerator &right)=delete
 
G4bool operator== (const G4VHighEnergyGenerator &right) const =delete
 
G4bool operator!= (const G4VHighEnergyGenerator &right) const =delete
 
virtual G4V3DNucleus * GetWoundedNucleus () const =0
 
virtual G4V3DNucleus * GetProjectileNucleus () const
 
virtual G4KineticTrackVector * Scatter (const G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)=0
 
void ModelDescription (std::ostream &) const override
 
- 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
 

Protected Member Functions

virtual void Init (const G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)=0
 
virtual G4ExcitedStringVector * GetStrings ()=0
 
G4bool EnergyAndMomentumCorrector (G4KineticTrackVector *Output, G4LorentzVector &TotalCollisionMomentum)
 
- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 

Additional Inherited Members

- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 50 of file G4VPartonStringModel.hh.

Constructor & Destructor Documentation

◆ G4VPartonStringModel() [1/2]

G4VPartonStringModel::G4VPartonStringModel ( const G4String &  modelName = "Parton String Model")

Definition at line 47 of file G4VPartonStringModel.cc.

48 : G4VHighEnergyGenerator(modelName),
49 stringFragmentationModel(nullptr)
50{
51 // Make shure Shotrylived particles are constructed.
52 // VI: should not instantiate particles by any model
53 /*
54 G4ShortLivedConstructor ShortLived;
55 ShortLived.ConstructParticle();
56 */
57}

◆ ~G4VPartonStringModel()

G4VPartonStringModel::~G4VPartonStringModel ( )
override

Definition at line 59 of file G4VPartonStringModel.cc.

60{
61}

◆ G4VPartonStringModel() [2/2]

G4VPartonStringModel::G4VPartonStringModel ( const G4VPartonStringModel &  right)
delete

Member Function Documentation

◆ EnergyAndMomentumCorrector()

G4bool G4VPartonStringModel::EnergyAndMomentumCorrector ( G4KineticTrackVector *  Output,
G4LorentzVector &  TotalCollisionMomentum 
)
protected

◆ GetProjectileNucleus()

G4V3DNucleus * G4VPartonStringModel::GetProjectileNucleus ( ) const
overridevirtual

Reimplemented from G4VHighEnergyGenerator.

Reimplemented in G4QGSModel< ParticipantType >, and G4QGSModel< G4GammaParticipants >.

Definition at line 456 of file G4VPartonStringModel.cc.

457{ return nullptr; }

Referenced by Scatter().

◆ GetStrings()

virtual G4ExcitedStringVector * G4VPartonStringModel::GetStrings ( )
protectedpure virtual

Implemented in G4FTFModel, G4QGSModel< ParticipantType >, and G4QGSModel< G4GammaParticipants >.

Referenced by Scatter().

◆ Init()

virtual void G4VPartonStringModel::Init ( const G4Nucleus &  theNucleus,
const G4DynamicParticle &  thePrimary 
)
protectedpure virtual

Implemented in G4FTFModel, G4QGSModel< ParticipantType >, and G4QGSModel< G4GammaParticipants >.

Referenced by Scatter().

◆ ModelDescription()

void G4VPartonStringModel::ModelDescription ( std::ostream &  outFile) const
overridevirtual

Reimplemented from G4VHighEnergyGenerator.

Reimplemented in G4QGSModel< ParticipantType >, and G4QGSModel< G4GammaParticipants >.

Definition at line 451 of file G4VPartonStringModel.cc.

452{
453 outFile << GetModelName() << " has no description yet.\n";
454}
G4HadronicInteraction::GetModelName
const G4String & GetModelName() const
Definition: G4HadronicInteraction.hh:115

◆ operator!=()

G4bool G4VPartonStringModel::operator!= ( const G4VPartonStringModel &  right) const
delete

◆ operator=()

const G4VPartonStringModel & G4VPartonStringModel::operator= ( const G4VPartonStringModel &  right)
delete

◆ operator==()

G4bool G4VPartonStringModel::operator== ( const G4VPartonStringModel &  right) const
delete

◆ Scatter()

G4KineticTrackVector * G4VPartonStringModel::Scatter ( const G4Nucleus &  theNucleus,
const G4DynamicParticle &  thePrimary 
)
overridevirtual

Implements G4VHighEnergyGenerator.

Definition at line 63 of file G4VPartonStringModel.cc.

65{
66 G4ExcitedStringVector * strings = nullptr;
67 G4DynamicParticle thePrimary=aPrimary;
68 G4LorentzVector SumStringMom(0.,0.,0.,0.);
69 G4KineticTrackVector * theResult = 0;
70 G4Nucleon * theNuclNucleon(nullptr);
71
72 #ifdef debug_PartonStringModel
73 G4cout<<G4endl;
74 G4cout<<"-----------------------Parton-String model is runnung ------------"<<G4endl;
75 G4cout<<"Projectile Name Mass "<<thePrimary.GetDefinition()->GetParticleName()<<" "
76 <<thePrimary.GetMass()<<G4endl;
77 G4cout<<" Momentum "<<thePrimary.Get4Momentum()<<G4endl;
78 G4cout<<"Target nucleus A Z "<<theNucleus.GetA_asInt()<<" "
79 <<theNucleus.GetZ_asInt()<<G4endl<<G4endl;
80 G4int Bsum=thePrimary.GetDefinition()->GetBaryonNumber() + theNucleus.GetA_asInt();
81 G4int Qsum=thePrimary.GetDefinition()->GetPDGCharge() + theNucleus.GetZ_asInt();
82 G4cout<<"Initial baryon number "<<Bsum<<G4endl;
83 G4cout<<"Initial charge "<<Qsum<<G4endl;
84 G4cout<<"-------------- Parton-String model: Generation of strings -------"<<G4endl<<G4endl;
85 Bsum -= theNucleus.GetA_asInt(); Qsum -= theNucleus.GetZ_asInt();
86 if(GetProjectileNucleus()) {
87 Bsum -= thePrimary.GetDefinition()->GetBaryonNumber();
88 Qsum -= thePrimary.GetDefinition()->GetPDGCharge();
89 }
90 G4int QsumSec(0), BsumSec(0);
91 G4LorentzVector SumPsecondr(0.,0.,0.,0.);
92 #endif
93
94 G4LorentzRotation toZ;
95 G4LorentzVector Ptmp=thePrimary.Get4Momentum();
96 toZ.rotateZ(-1*Ptmp.phi());
97 toZ.rotateY(-1*Ptmp.theta());
98 thePrimary.Set4Momentum(toZ*Ptmp);
99 G4LorentzRotation toLab(toZ.inverse());
100
101 G4bool Success=true;
102 G4int attempts = 0, maxAttempts=1000;
103 do
104 {
105 if (attempts++ > maxAttempts )
106 {
107 Init(theNucleus,thePrimary); // To put a nucleus into ground state
108 // But marks of hitted nucleons are left. They must be erased.
109 G4V3DNucleus * ResNucleus = GetWoundedNucleus();
110 theNuclNucleon = ResNucleus->StartLoop() ? ResNucleus->GetNextNucleon() : nullptr;
111 while( theNuclNucleon )
112 {
113 if(theNuclNucleon->AreYouHit()) theNuclNucleon->Hit(nullptr);
114 theNuclNucleon = ResNucleus->GetNextNucleon();
115 }
116
117 G4V3DNucleus * ProjResNucleus = GetProjectileNucleus();
118 if(ProjResNucleus != 0)
119 {
120 theNuclNucleon = ProjResNucleus->StartLoop() ? ProjResNucleus->GetNextNucleon() : nullptr;
121 while( theNuclNucleon )
122 {
123 if(theNuclNucleon->AreYouHit()) theNuclNucleon->Hit(nullptr);
124 theNuclNucleon = ProjResNucleus->GetNextNucleon();
125 }
126 }
127
128 G4ExceptionDescription ed;
129 ed << "Projectile Name Mass " <<thePrimary.GetDefinition()->GetParticleName()
130 << " " << thePrimary.GetMass()<< G4endl;
131 ed << "           Momentum  " << thePrimary.Get4Momentum() <<G4endl;
132 ed << "Target nucleus   A Z " << theNucleus.GetA_asInt() << " "
133 << theNucleus.GetZ_asInt() <<G4endl;
134 ed << "Initial states of projectile and target nucleus will be returned!"<<G4endl;
135 G4Exception( "G4VPartonStringModel::Scatter(): fails to generate or fragment strings ",
136 "HAD_PARTON_STRING_001", JustWarning, ed );
137
138 G4ThreeVector Position(0.,0.,2*ResNucleus->GetOuterRadius());
139 G4KineticTrack* Hadron = new G4KineticTrack(aPrimary.GetParticleDefinition(), 0.,
140 Position, aPrimary.Get4Momentum());
141 if(theResult == nullptr) theResult = new G4KineticTrackVector();
142 theResult->push_back(Hadron);
143 return theResult;
144 }
145
146 Success=true;
147
148 Init(theNucleus,thePrimary);
149
150 strings = GetStrings();
151
152 if (strings->empty()) { Success=false; continue; }
153
154 G4double stringEnergy(0);
155 SumStringMom=G4LorentzVector(0.,0.,0.,0.);
156
157 #ifdef debug_PartonStringModel
158 G4cout<<"------------ Parton-String model: Number of produced strings ---- "<<strings->size()<<G4endl;
159 #endif
160
161 #ifdef debug_heavyHadrons
162 // Check charm and bottom numbers of the projectile:
163 G4int count_charm_projectile = thePrimary.GetDefinition()->GetQuarkContent( 4 ) -
164 thePrimary.GetDefinition()->GetAntiQuarkContent( 4 );
165 G4int count_bottom_projectile = thePrimary.GetDefinition()->GetQuarkContent( 5 ) -
166 thePrimary.GetDefinition()->GetAntiQuarkContent( 5 );
167 G4int count_charm_strings = 0, count_bottom_strings = 0;
168 G4int count_charm_hadrons = 0, count_bottom_hadrons = 0;
169 #endif
170
171 for ( unsigned int astring=0; astring < strings->size(); astring++)
172 {
173 // rotate string to lab frame, models have it aligned to z
174 if((*strings)[astring]->IsExcited())
175 {
176 stringEnergy += (*strings)[astring]->GetLeftParton()->Get4Momentum().t();
177 stringEnergy += (*strings)[astring]->GetRightParton()->Get4Momentum().t();
178 (*strings)[astring]->LorentzRotate(toLab);
179 SumStringMom+=(*strings)[astring]->Get4Momentum();
180 #ifdef debug_PartonStringModel
181 G4cout<<"String No "<<astring+1<<" "<<(*strings)[astring]->Get4Momentum()<<" "
182 <<(*strings)[astring]->Get4Momentum().mag()
183 <<" Partons "<<(*strings)[astring]->GetLeftParton()->GetDefinition()->GetPDGEncoding()
184 <<" "<<(*strings)[astring]->GetRightParton()->GetDefinition()->GetPDGEncoding()<<G4endl;
185 #endif
186
187 #ifdef debug_heavyHadrons
188 G4int left_charm = (*strings)[astring]->GetLeftParton()->GetDefinition()->GetQuarkContent( 4 );
189 G4int left_anticharm = (*strings)[astring]->GetLeftParton()->GetDefinition()->GetAntiQuarkContent( 4 );
190 G4int right_charm = (*strings)[astring]->GetRightParton()->GetDefinition()->GetQuarkContent( 4 );
191 G4int right_anticharm = (*strings)[astring]->GetRightParton()->GetDefinition()->GetAntiQuarkContent( 4 );
192 G4int left_bottom = (*strings)[astring]->GetLeftParton()->GetDefinition()->GetQuarkContent( 5 );
193 G4int left_antibottom = (*strings)[astring]->GetLeftParton()->GetDefinition()->GetAntiQuarkContent( 5 );
194 G4int right_bottom = (*strings)[astring]->GetRightParton()->GetDefinition()->GetQuarkContent( 5 );
195 G4int right_antibottom = (*strings)[astring]->GetRightParton()->GetDefinition()->GetAntiQuarkContent( 5 );
196 if ( left_charm != 0 || left_anticharm != 0 || right_charm != 0 || right_anticharm != 0 ||
197 left_bottom != 0 || left_antibottom != 0 || right_bottom != 0 || right_antibottom != 0 ) {
198 count_charm_strings += left_charm - left_anticharm + right_charm - right_anticharm;
199 count_bottom_strings += left_bottom - left_antibottom + right_bottom - right_antibottom;
200 G4cout << "G4VPartonStringModel::Scatter : string #" << astring << " ("
201 << (*strings)[astring]->GetLeftParton()->GetDefinition()->GetParticleName() << " , "
202 << (*strings)[astring]->GetRightParton()->GetDefinition()->GetParticleName() << ")" << G4endl;
203 }
204 #endif
205 }
206 else
207 {
208 stringEnergy += (*strings)[astring]->GetKineticTrack()->Get4Momentum().t();
209 (*strings)[astring]->LorentzRotate(toLab);
210 SumStringMom+=(*strings)[astring]->GetKineticTrack()->Get4Momentum();
211 #ifdef debug_PartonStringModel
212 G4cout<<"A track No "<<astring+1<<" "
213 <<(*strings)[astring]->GetKineticTrack()->Get4Momentum()<<" "
214 <<(*strings)[astring]->GetKineticTrack()->Get4Momentum().mag()<<" "
215 <<(*strings)[astring]->GetKineticTrack()->GetDefinition()->GetParticleName()<<G4endl;
216 #endif
217
218 #ifdef debug_heavyHadrons
219 G4int charm = (*strings)[astring]->GetKineticTrack()->GetDefinition()->GetQuarkContent( 4 );
220 G4int anticharm = (*strings)[astring]->GetKineticTrack()->GetDefinition()->GetAntiQuarkContent( 4 );
221 G4int bottom = (*strings)[astring]->GetKineticTrack()->GetDefinition()->GetQuarkContent( 5 );
222 G4int antibottom = (*strings)[astring]->GetKineticTrack()->GetDefinition()->GetAntiQuarkContent( 5 );
223 if ( charm != 0 || anticharm != 0 || bottom != 0 || antibottom != 0 ) {
224 count_charm_strings += charm - anticharm;
225 count_bottom_strings += bottom - antibottom;
226 G4cout << "G4VPartonStringModel::Scatter : track #" << astring << "\t"
227 << (*strings)[astring]->GetKineticTrack()->GetDefinition()->GetParticleName() << G4endl;
228 }
229 #endif
230 }
231 }
232
233 #ifdef debug_heavyHadrons
234 if ( count_charm_projectile != count_charm_strings ) {
235 G4cout << "G4VPartonStringModel::Scatter : CHARM VIOLATION in String formation ! #projectile="
236 << count_charm_projectile << " ; #strings=" << count_charm_strings << G4endl;
237 }
238 if ( count_bottom_projectile != count_bottom_strings ) {
239 G4cout << "G4VPartonStringModel::Scatter : BOTTOM VIOLATION in String formation ! #projectile="
240 << count_bottom_projectile << " ; #strings=" << count_bottom_strings << G4endl;
241 }
242 #endif
243
244 #ifdef debug_PartonStringModel
245 G4cout<<G4endl<<"SumString4Mom "<<SumStringMom<<G4endl;
246 G4LorentzVector TargetResidual4Momentum(0.,0.,0.,0.);
247 G4LorentzVector ProjectileResidual4Momentum(0.,0.,0.,0.);
248 G4int hitsT(0), charged_hitsT(0);
249 G4int hitsP(0), charged_hitsP(0);
250 G4double ExcitationEt(0.), ExcitationEp(0.);
251 #endif
252
253 G4V3DNucleus * ProjResNucleus = GetProjectileNucleus();
254
255 G4int numberProtonProjectileResidual( 0 ), numberNeutronProjectileResidual( 0 );
256 if(ProjResNucleus != 0)
257 {
258 theNuclNucleon = ProjResNucleus->StartLoop() ? ProjResNucleus->GetNextNucleon() : nullptr;
259 G4int numberProtonProjectileHits( 0 ), numberNeutronProjectileHits( 0 );
260 while( theNuclNucleon )
261 {
262 if(theNuclNucleon->AreYouHit())
263 {
264 G4LorentzVector tmp=toLab*theNuclNucleon->Get4Momentum();
265 #ifdef debug_PartonStringModel
266 ProjectileResidual4Momentum += tmp;
267 hitsP++;
268 if ( theNuclNucleon->GetDefinition() == G4Proton::Proton() ) ++charged_hitsP;
269 ExcitationEp +=theNuclNucleon->GetBindingEnergy();
270 #endif
271 theNuclNucleon->SetMomentum(tmp);
272 if ( theNuclNucleon->GetDefinition() == G4Proton::Proton() ) ++numberProtonProjectileHits;
273 if ( theNuclNucleon->GetDefinition() == G4Neutron::Neutron() ) ++numberNeutronProjectileHits;
274 }
275 theNuclNucleon = ProjResNucleus->GetNextNucleon();
276 }
277 #ifdef debug_PartonStringModel
278 G4cout<<"Projectile residual A, Z and E* "
279 <<thePrimary.GetDefinition()->GetBaryonNumber() - hitsP<<" "
280 <<thePrimary.GetDefinition()->GetPDGCharge() - charged_hitsP<<" "
281 <<ExcitationEp<<G4endl;
282 G4cout<<"Projectile residual 4 momentum "<<ProjectileResidual4Momentum<<G4endl;
283 #endif
284 numberProtonProjectileResidual = thePrimary.GetDefinition()->GetPDGCharge() - numberProtonProjectileHits;
285 numberNeutronProjectileResidual = thePrimary.GetDefinition()->GetBaryonNumber()
286 - thePrimary.GetDefinition()->GetPDGCharge() - numberNeutronProjectileHits;
287 }
288
289 G4V3DNucleus * ResNucleus = GetWoundedNucleus();
290
291 // loop over wounded nucleus
292 theNuclNucleon = ResNucleus->StartLoop() ? ResNucleus->GetNextNucleon() : nullptr;
293 G4int numberProtonTargetHits( 0 ), numberNeutronTargetHits( 0 );
294 while( theNuclNucleon )
295 {
296 if(theNuclNucleon->AreYouHit())
297 {
298 G4LorentzVector tmp=toLab*theNuclNucleon->Get4Momentum();
299 #ifdef debug_PartonStringModel
300 TargetResidual4Momentum += tmp;
301 hitsT++;
302 if ( theNuclNucleon->GetDefinition() == G4Proton::Proton() ) ++charged_hitsT;
303 ExcitationEt +=theNuclNucleon->GetBindingEnergy();
304 #endif
305 theNuclNucleon->SetMomentum(tmp);
306 if ( theNuclNucleon->GetDefinition() == G4Proton::Proton() ) ++numberProtonTargetHits;
307 if ( theNuclNucleon->GetDefinition() == G4Neutron::Neutron() ) ++numberNeutronTargetHits;
308 }
309 theNuclNucleon = ResNucleus->GetNextNucleon();
310 }
311
312 #ifdef debug_PartonStringModel
313 G4cout<<"Target residual A, Z and E* "
314 <<theNucleus.GetA_asInt() - hitsT<<" "
315 <<theNucleus.GetZ_asInt() - charged_hitsT<<" "
316 <<ExcitationEt<<G4endl;
317 G4cout<<"Target residual 4 momentum "<<TargetResidual4Momentum<<G4endl;
318 Bsum+=( hitsT + hitsP);
319 Qsum+=(charged_hitsT + charged_hitsP);
320 G4cout<<"Hitted # of nucleons of projectile and target "<<hitsP<<" "<<hitsT<<G4endl;
321 G4cout<<"Hitted # of protons of projectile and target "
322 <<charged_hitsP<<" "<<charged_hitsT<<G4endl<<G4endl;
323 G4cout<<"Bsum Qsum "<<Bsum<<" "<<Qsum<<G4endl<<G4endl;
324 #endif
325
326 // Re-sample in the case of unphysical nuclear residual:
327 // 1 (H), 2 (2He), and 3 (3Li) protons alone without neutrons can exist, but not more;
328 // no bound states of 2 or more neutrons without protons can exist.
329 G4int numberProtonTargetResidual = theNucleus.GetZ_asInt() - numberProtonTargetHits;
330 G4int numberNeutronTargetResidual = theNucleus.GetA_asInt() - theNucleus.GetZ_asInt() - numberNeutronTargetHits;
331 G4bool unphysicalResidual = false;
332 if ( ( numberProtonTargetResidual > 3 && numberNeutronTargetResidual == 0 ) ||
333 ( numberProtonTargetResidual == 0 && numberNeutronTargetResidual > 1 ) ) {
334 unphysicalResidual = true;
335 //G4cout << "***UNPHYSICAL TARGET RESIDUAL*** Z=" << numberProtonTargetResidual
336 // << " ; N=" << numberNeutronTargetResidual;
337 }
338 if ( ( numberProtonProjectileResidual > 3 && numberNeutronProjectileResidual == 0 ) ||
339 ( numberProtonProjectileResidual == 0 && numberNeutronProjectileResidual > 1 ) ) {
340 unphysicalResidual = true;
341 //G4cout << "***UNPHYSICAL PROJECTILE RESIDUAL*** Z=" << numberProtonProjectileResidual
342 // << " ; N=" << numberNeutronProjectileResidual;
343 }
344 if ( unphysicalResidual ) {
345 //G4cout << " -> REJECTING COLLISION because of unphysical residual !" << G4endl;
346 Success = false;
347 continue;
348 }
349
350 //=========================================================================================
351 // Fragment strings
352 #ifdef debug_PartonStringModel
353 G4cout<<"---------------- Attempt to fragment strings ------------- "<<attempts<<G4endl;
354 #endif
355
356 G4double InvMass=SumStringMom.mag();
357 G4double SumMass(0.);
358
359 #ifdef debug_PartonStringModel
360 QsumSec=0; BsumSec=0;
361 SumPsecondr=G4LorentzVector(0.,0.,0.,0.);
362 #endif
363
364 if(theResult != nullptr)
365 {
366 std::for_each(theResult->begin(), theResult->end(), DeleteKineticTrack());
367 delete theResult;
368 }
369
370 theResult = stringFragmentationModel->FragmentStrings(strings);
371
372 #ifdef debug_PartonStringModel
373 G4cout<<"String fragmentation success (OK - #0, Not OK - 0) : "<<theResult<<G4endl;
374 #endif
375
376 if(theResult == 0) {Success=false; continue;}
377
378 #ifdef debug_PartonStringModel
379 G4cout<<"Parton-String model: Number of produced particles "<<theResult->size()<<G4endl;
380 SumPsecondr = G4LorentzVector(0.,0.,0.,0.);
381 QsumSec = 0; BsumSec = 0;
382 #endif
383
384 SumMass=0.;
385 for ( unsigned int i=0; i < theResult->size(); i++)
386 {
387 SumMass+=(*theResult)[i]->Get4Momentum().mag();
388 #ifdef debug_PartonStringModel
389 G4cout<<i<<" : "<<(*theResult)[i]->GetDefinition()->GetParticleName()<<" "
390 <<(*theResult)[i]->Get4Momentum()<<" "
391 <<(*theResult)[i]->Get4Momentum().mag()<<" "
392 <<(*theResult)[i]->GetDefinition()->GetPDGMass()<<G4endl;
393 SumPsecondr+=(*theResult)[i]->Get4Momentum();
394 BsumSec += (*theResult)[i]->GetDefinition()->GetBaryonNumber();
395 QsumSec += (*theResult)[i]->GetDefinition()->GetPDGCharge();
396 #endif
397
398 #ifdef debug_heavyHadrons
399 G4int charm = (*theResult)[i]->GetDefinition()->GetQuarkContent( 4 );
400 G4int anticharm = (*theResult)[i]->GetDefinition()->GetAntiQuarkContent( 4 );
401 G4int bottom = (*theResult)[i]->GetDefinition()->GetQuarkContent( 5 );
402 G4int antibottom = (*theResult)[i]->GetDefinition()->GetAntiQuarkContent( 5 );
403 if ( charm != 0 || anticharm != 0 || bottom != 0 || antibottom != 0 ) {
404 count_charm_hadrons += charm - anticharm;
405 count_bottom_hadrons += bottom - antibottom;
406 G4cout << "G4VPartonStringModel::Scatter : hadron #" << i << "\t"
407 << (*theResult)[i]->GetDefinition()->GetParticleName() << G4endl;
408 }
409 #endif
410 }
411
412 #ifdef debug_heavyHadrons
413 if ( count_charm_projectile != count_charm_hadrons ) {
414 G4cout << "G4VPartonStringModel::Scatter : CHARM VIOLATION in String hadronization ! #projectile="
415 << count_charm_projectile << " ; #hadrons=" << count_charm_hadrons << G4endl;
416 }
417 if ( count_bottom_projectile != count_bottom_hadrons ) {
418 G4cout << "G4VPartonStringModel::Scatter : BOTTOM VIOLATION in String hadronization ! #projectile="
419 << count_bottom_projectile << " ; #hadrons=" << count_bottom_hadrons << G4endl;
420 }
421 #endif
422
423 #ifdef debug_PartonStringModel
424 G4cout<<G4endl<<"-----------------------Parton-String model: balances -------------"<<G4endl;
425 if(Qsum != QsumSec) {
426 G4cout<<"Charge is not conserved!!! ----"<<G4endl;
427 G4cout<<" Qsum != QsumSec "<<Qsum<<" "<<QsumSec<<G4endl;
428 }
429 if(Bsum != BsumSec) {
430 G4cout<<"Baryon number is not conserved!!!"<<G4endl;
431 G4cout<<" Bsum != BsumSec "<<Bsum<<" "<<BsumSec<<G4endl;
432 }
433 #endif
434
435 if((SumMass > InvMass)||(SumMass == 0.)) {Success=false;}
436 std::for_each(strings->begin(), strings->end(), DeleteString() );
437 delete strings;
438
439 } while(!Success);
440
441 #ifdef debug_PartonStringModel
442 G4cout <<"Baryon number balance "<<Bsum-BsumSec<<G4endl;
443 G4cout <<"Charge balance "<<Qsum-QsumSec<<G4endl;
444 G4cout <<"4 momentum balance "<<SumStringMom-SumPsecondr<<G4endl;
445 G4cout<<"---------------------End of Parton-String model work -------------"<<G4endl<<G4endl;
446 #endif
447
448 return theResult;
449}
JustWarning
@ JustWarning
Definition: G4ExceptionSeverity.hh:65
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
G4ExceptionDescription
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
G4ExcitedStringVector
std::vector< G4ExcitedString * > G4ExcitedStringVector
Definition: G4ExcitedStringVector.hh:41
G4LorentzVector
CLHEP::HepLorentzVector G4LorentzVector
Definition: G4LorentzVector.hh:34
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::HepLorentzRotation::rotateY
HepLorentzRotation & rotateY(double delta)
Definition: LorentzRotation.cc:158
CLHEP::HepLorentzRotation::rotateZ
HepLorentzRotation & rotateZ(double delta)
Definition: LorentzRotation.cc:170
CLHEP::HepLorentzRotation::inverse
HepLorentzRotation inverse() const
CLHEP::HepLorentzVector::theta
double theta() const
G4DynamicParticle::GetMass
G4double GetMass() const
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4DynamicParticle::Get4Momentum
G4LorentzVector Get4Momentum() const
G4DynamicParticle::Set4Momentum
void Set4Momentum(const G4LorentzVector &momentum)
Definition: G4DynamicParticle.cc:410
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:103
G4Nucleon::Hit
void Hit(G4VSplitableHadron *aHit)
Definition: G4Nucleon.hh:89
G4Nucleus::GetA_asInt
G4int GetA_asInt() const
Definition: G4Nucleus.hh:109
G4Nucleus::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Nucleus.hh:115
G4ParticleDefinition::GetQuarkContent
G4int GetQuarkContent(G4int flavor) const
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:112
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:108
G4ParticleDefinition::GetAntiQuarkContent
G4int GetAntiQuarkContent(G4int flavor) const
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:92
G4V3DNucleus::GetOuterRadius
virtual G4double GetOuterRadius()=0
G4V3DNucleus::GetNextNucleon
virtual G4Nucleon * GetNextNucleon()=0
G4V3DNucleus::StartLoop
virtual G4bool StartLoop()=0
G4VHighEnergyGenerator::GetWoundedNucleus
virtual G4V3DNucleus * GetWoundedNucleus() const =0
G4VPartonStringModel::Init
virtual void Init(const G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)=0
G4VPartonStringModel::GetStrings
virtual G4ExcitedStringVector * GetStrings()=0
G4VPartonStringModel::GetProjectileNucleus
G4V3DNucleus * GetProjectileNucleus() const override
Definition: G4VPartonStringModel.cc:456
G4VStringFragmentation::FragmentStrings
virtual G4KineticTrackVector * FragmentStrings(const G4ExcitedStringVector *theStrings)=0
field_utils::Value3D::Position
@ Position

◆ SetFragmentationModel()

void G4VPartonStringModel::SetFragmentationModel ( G4VStringFragmentation *  aModel)
inline

Definition at line 78 of file G4VPartonStringModel.hh.

79{
80 if(aModel != stringFragmentationModel) {
81 stringFragmentationModel = aModel;
82 }
83}

Referenced by G4BertiniElectroNuclearBuilder::Build(), G4FTFBuilder::BuildModel(), G4QGSBuilder::BuildModel(), LBE::ConstructHad(), G4ElectroVDNuclearModel::G4ElectroVDNuclearModel(), G4FTFBinaryKaonBuilder::G4FTFBinaryKaonBuilder(), G4FTFBinaryNeutronBuilder::G4FTFBinaryNeutronBuilder(), G4FTFBinaryPiKBuilder::G4FTFBinaryPiKBuilder(), G4FTFBinaryPionBuilder::G4FTFBinaryPionBuilder(), G4FTFBinaryProtonBuilder::G4FTFBinaryProtonBuilder(), G4FTFPAntiBarionBuilder::G4FTFPAntiBarionBuilder(), G4FTFPKaonBuilder::G4FTFPKaonBuilder(), G4FTFPNeutronBuilder::G4FTFPNeutronBuilder(), G4FTFPPiKBuilder::G4FTFPPiKBuilder(), G4FTFPPionBuilder::G4FTFPPionBuilder(), G4FTFPProtonBuilder::G4FTFPProtonBuilder(), G4HadronicAbsorptionFritiof::G4HadronicAbsorptionFritiof(), G4HadronicAbsorptionFritiofWithBinaryCascade::G4HadronicAbsorptionFritiofWithBinaryCascade(), G4HyperonFTFPBuilder::G4HyperonFTFPBuilder(), G4HyperonQGSPBuilder::G4HyperonQGSPBuilder(), G4MuonVDNuclearModel::G4MuonVDNuclearModel(), G4QGSBinaryKaonBuilder::G4QGSBinaryKaonBuilder(), G4QGSBinaryNeutronBuilder::G4QGSBinaryNeutronBuilder(), G4QGSBinaryPiKBuilder::G4QGSBinaryPiKBuilder(), G4QGSBinaryPionBuilder::G4QGSBinaryPionBuilder(), G4QGSBinaryProtonBuilder::G4QGSBinaryProtonBuilder(), G4QGSPAntiBarionBuilder::G4QGSPAntiBarionBuilder(), G4QGSPKaonBuilder::G4QGSPKaonBuilder(), G4QGSPLundStrFragmProtonBuilder::G4QGSPLundStrFragmProtonBuilder(), G4QGSPNeutronBuilder::G4QGSPNeutronBuilder(), G4QGSPPiKBuilder::G4QGSPPiKBuilder(), G4QGSPPionBuilder::G4QGSPPionBuilder(), and G4QGSPProtonBuilder::G4QGSPProtonBuilder().

The documentation for this class was generated from the following files: