Geant4-11
Public Member Functions | Protected Attributes | Private Attributes
G4FTFParticipants Class Reference

#include <G4FTFParticipants.hh>

Inheritance diagram for G4FTFParticipants:
G4VParticipants

Public Member Functions

void Clean ()
 
 G4FTFParticipants ()
 
 G4FTFParticipants (const G4FTFParticipants &right)=delete
 
G4double GetBmax2 () const
 
G4double GetBmin2 () const
 
G4double GetImpactParameter () const
 
G4InteractionContentGetInteraction ()
 
void GetList (const G4ReactionProduct &thePrimary, G4FTFParameters *theParameters)
 
virtual G4V3DNucleusGetProjectileNucleus () const
 
virtual G4V3DNucleusGetWoundedNucleus () const
 
virtual void Init (G4int theZ, G4int theA)
 
virtual void InitProjectileNucleus (G4int theZ, G4int theA, G4int numberOfLambdasOrAntiLambdas=0)
 
G4bool Next ()
 
G4bool operator!= (const G4FTFParticipants &right) const =delete
 
G4bool operator!= (const G4VParticipants &right) const =delete
 
const G4FTFParticipantsoperator= (const G4FTFParticipants &right)=delete
 
G4bool operator== (const G4FTFParticipants &right) const =delete
 
G4bool operator== (const G4VParticipants &right) const =delete
 
G4bool SampleBinInterval () const
 
void SetBminBmax (const G4double bmin_value, const G4double bmax_value)
 
void SetImpactParameter (const G4double b_value)
 
virtual void SetNucleus (G4V3DNucleus *aNucleus)
 
virtual void SetProjectileNucleus (G4V3DNucleus *aNucleus)
 
void ShiftInteractionTime ()
 
void SortInteractionsIncT ()
 
void StartLoop ()
 
 ~G4FTFParticipants ()
 

Protected Attributes

G4V3DNucleustheNucleus
 
G4V3DNucleustheProjectileNucleus
 

Private Attributes

G4double Bimpact
 
G4bool BinInterval
 
G4double Bmax2
 
G4double Bmin2
 
G4int currentInteraction
 
std::vector< G4InteractionContent * > theInteractions
 

Detailed Description

Definition at line 50 of file G4FTFParticipants.hh.

Constructor & Destructor Documentation

◆ G4FTFParticipants() [1/2]

G4FTFParticipants::G4FTFParticipants ( )

Definition at line 62 of file G4FTFParticipants.cc.

62 : Bimpact( 0.0 ), BinInterval( false ),
63 Bmin2( -1.0 ), Bmax2( -1.0 ),
65{}

◆ ~G4FTFParticipants()

G4FTFParticipants::~G4FTFParticipants ( )

Definition at line 70 of file G4FTFParticipants.cc.

70{}

◆ G4FTFParticipants() [2/2]

G4FTFParticipants::G4FTFParticipants ( const G4FTFParticipants right)
delete

Member Function Documentation

◆ Clean()

void G4FTFParticipants::Clean ( )

Definition at line 384 of file G4FTFParticipants.cc.

384 {
385 for ( size_t i = 0; i < theInteractions.size(); i++ ) {
386 if ( theInteractions[ i ] ) {
387 delete theInteractions[ i ];
388 theInteractions[ i ] = 0;
389 }
390 }
391 theInteractions.clear();
393}
std::vector< G4InteractionContent * > theInteractions

References currentInteraction, and theInteractions.

Referenced by G4FTFModel::GetStrings(), and G4FTFModel::Init().

◆ GetBmax2()

G4double G4FTFParticipants::GetBmax2 ( ) const
inline

Definition at line 119 of file G4FTFParticipants.hh.

119 {
120 return Bmax2;
121}

References Bmax2.

Referenced by GetList().

◆ GetBmin2()

G4double G4FTFParticipants::GetBmin2 ( ) const
inline

Definition at line 115 of file G4FTFParticipants.hh.

115 {
116 return Bmin2;
117}

References Bmin2.

Referenced by GetList().

◆ GetImpactParameter()

G4double G4FTFParticipants::GetImpactParameter ( ) const
inline

Definition at line 99 of file G4FTFParticipants.hh.

99 {
100 return Bimpact;
101}

References Bimpact.

Referenced by G4FTFModel::GetStrings().

◆ GetInteraction()

G4InteractionContent & G4FTFParticipants::GetInteraction ( )
inline

◆ GetList()

void G4FTFParticipants::GetList ( const G4ReactionProduct thePrimary,
G4FTFParameters theParameters 
)

Definition at line 74 of file G4FTFParticipants.cc.

75 {
76
77 #ifdef debugFTFparticipant
78 G4cout << "Participants::GetList" << G4endl
79 << "thePrimary " << thePrimary.GetMomentum() << G4endl << G4endl;
80 #endif
81
82 G4double betta_z = thePrimary.GetMomentum().z() / thePrimary.GetTotalEnergy();
83 if ( betta_z < 1.0e-10 ) betta_z = 1.0e-10;
84
85 StartLoop(); // reset Loop over Interactions
86
87 for ( unsigned int i = 0; i < theInteractions.size(); i++ ) delete theInteractions[i];
88 theInteractions.clear();
89
90 G4double deltaxy = 2.0 * fermi; // Extra nuclear radius
91
92 if ( theProjectileNucleus == nullptr ) { // Hadron-nucleus or anti-baryon-nucleus interactions
93
94 G4double impactX( 0.0 ), impactY( 0.0 );
95 G4double B( 0.0 ), B2( 0.0 );
96
97 G4VSplitableHadron* primarySplitable = new G4DiffractiveSplitableHadron( thePrimary );
98
99 #ifdef debugFTFparticipant
100 G4cout << "Hadron-nucleus or anti-baryon-nucleus interactions" << G4endl;
101 #endif
102
103 G4double xyradius = theNucleus->GetOuterRadius() + deltaxy; // Range of impact parameter sampling
104
105 const G4int maxNumberOfLoops = 1000;
106 G4int loopCounter = 0;
107 do {
108
109 std::pair< G4double, G4double > theImpactParameter;
110 if ( SampleBinInterval() ) {
111 B2 = GetBmin2() + G4UniformRand() * ( GetBmax2() - GetBmin2() );
112 B = B2 > 0.0 ? std::sqrt( B2 ) : 0.0;
113 G4double Phi = twopi * G4UniformRand();
114 impactX = B * std::cos( Phi );
115 impactY = B * std::sin( Phi );
117 } else {
118 theImpactParameter = theNucleus->ChooseImpactXandY( xyradius );
119 impactX = theImpactParameter.first;
120 impactY = theImpactParameter.second;
121 SetImpactParameter( std::sqrt( sqr(impactX) + sqr(impactY) ) );
122 }
123
124 #ifdef debugFTFparticipant
125 G4cout << "New interaction list," << " b[fm]= "
126 << std::sqrt( sqr(impactX ) + sqr( impactY ) )/fermi << G4endl;
127 #endif
128
129 G4ThreeVector thePosition( impactX, impactY, 0.0 );
130 primarySplitable->SetPosition( thePosition );
131
134
135 #ifdef debugFTFparticipant
136 G4int TrN( 0 );
137 #endif
138
139 while ( ( nucleon = theNucleus->GetNextNucleon() ) ) { /* Loop checking, 10.08.2015, A.Ribon */
140
141 G4double impact2 = sqr( impactX - nucleon->GetPosition().x() ) +
142 sqr( impactY - nucleon->GetPosition().y() );
143
144 if ( theParameters->GetProbabilityOfInteraction( impact2/fermi/fermi ) >
145 G4UniformRand() ) {
146 primarySplitable->SetStatus( 1 ); // It takes part in the interaction
147 G4VSplitableHadron* targetSplitable = 0;
148 if ( ! nucleon->AreYouHit() ) {
149 targetSplitable = new G4DiffractiveSplitableHadron( *nucleon );
150 nucleon->Hit( targetSplitable );
151 targetSplitable->SetStatus( 1 ); // It takes part in the interaction
152
153 #ifdef debugFTFparticipant
154 G4cout << "Participated nucleons #, " << TrN << " " << "Splitable Pr* Tr* "
155 << primarySplitable << " " << targetSplitable << G4endl;
156 #endif
157
158 }
159 G4InteractionContent* aInteraction = new G4InteractionContent( primarySplitable );
160 G4Nucleon* PrNucleon = 0;
161 aInteraction->SetProjectileNucleon( PrNucleon );
162 aInteraction->SetTarget( targetSplitable );
163 aInteraction->SetTargetNucleon( nucleon );
164 aInteraction->SetStatus( 1 );
165 aInteraction->SetInteractionTime( ( primarySplitable->GetPosition().z() +
166 nucleon->GetPosition().z() ) / betta_z );
167 theInteractions.push_back( aInteraction );
168 }
169
170 #ifdef debugFTFparticipant
171 TrN++;
172 #endif
173
174 }
175
176 } while ( ( theInteractions.size() == 0 ) &&
177 ++loopCounter < maxNumberOfLoops ); /* Loop checking, 10.08.2015, A.Ribon */
178 if ( loopCounter >= maxNumberOfLoops ) {
179 #ifdef debugFTFparticipant
180 G4cout << "BAD situation: forced exit from the while loop!" << G4endl;
181 #endif
182 return;
183 }
184
185 #ifdef debugFTFparticipant
186 G4cout << "Number of Hit nucleons " << theInteractions.size() << "\t Bx[fm] " << impactX/fermi
187 << "\t By[fm] " << impactY/fermi << "\t B[fm] "
188 << std::sqrt( sqr( impactX ) + sqr( impactY ) )/fermi << G4endl << G4endl;
189 #endif
190
191 //SortInteractionsIncT(); // Not need because nucleons are sorted in increasing z-coordinates.
192 ShiftInteractionTime(); // To put correct times and z-coordinates
193 return;
194
195 } // end of if ( theProjectileNucleus == 0 )
196
197 // Projectile and target are nuclei
198
199 #ifdef debugFTFparticipant
200 G4cout << "Projectile and target are nuclei" << G4endl;
201 #endif
202
203 //G4cout<<theProjectileNucleus->GetOuterRadius()/fermi<<" "<<theNucleus->GetOuterRadius()/fermi<<" "<<deltaxy/fermi<<G4endl;
204
205 // Range of impact parameter sampling
207
208 G4double impactX( 0.0 ), impactY( 0.0 );
209 G4double B( 0.0 ), B2( 0.0 );
210
211 const G4int maxNumberOfLoops = 1000;
212 G4int loopCounter = 0;
213 do {
214
215 std::pair< G4double, G4double > theImpactParameter;
216 if ( SampleBinInterval() ) {
217 B2 = GetBmin2() + G4UniformRand() * ( GetBmax2() - GetBmin2() );
218 B = B2 > 0.0 ? std::sqrt( B2 ) : 0.0; // In G4 internal units (mm)
219 G4double Phi = twopi * G4UniformRand();
220 impactX = B * std::cos( Phi );
221 impactY = B * std::sin( Phi );
223 } else {
224 theImpactParameter = theNucleus->ChooseImpactXandY( xyradius );
225 impactX = theImpactParameter.first;
226 impactY = theImpactParameter.second;
227 SetImpactParameter( std::sqrt( sqr(impactX) + sqr(impactY) ) );
228 }
229
230 #ifdef debugFTFparticipant
231 G4cout << "New interaction list, " << "b[fm] "
232 << std::sqrt( sqr( impactX ) + sqr( impactY ) )/fermi << G4endl;
233 #endif
234
235 G4ThreeVector theBeamPosition( impactX, impactY, 0.0 );
236
238 G4Nucleon* ProjectileNucleon;
239
240 #ifdef debugFTFparticipant
241 G4int PrNuclN( 0 );
242 #endif
243
244 while ( ( ProjectileNucleon = theProjectileNucleus->GetNextNucleon() ) ) { /* Loop checking, 10.08.2015, A.Ribon */
245
246 G4VSplitableHadron* ProjectileSplitable = 0;
248 G4Nucleon* TargetNucleon = 0;
249
250 #ifdef debugFTFparticipant
251 G4int TrNuclN( 0 );
252 #endif
253
254 while ( ( TargetNucleon = theNucleus->GetNextNucleon() ) ) { /* Loop checking, 10.08.2015, A.Ribon */
255
256 G4double impact2 = sqr( impactX + ProjectileNucleon->GetPosition().x() -
257 TargetNucleon->GetPosition().x() ) +
258 sqr( impactY + ProjectileNucleon->GetPosition().y() -
259 TargetNucleon->GetPosition().y() );
260 G4VSplitableHadron* TargetSplitable = 0;
261 if ( theParameters->GetProbabilityOfInteraction( impact2/fermi/fermi ) >
262 G4UniformRand() ) { // An interaction has happend!
263
264 #ifdef debugFTFparticipant
265 G4cout << G4endl << "An Interaction has happend" << G4endl << "Proj N mom " << PrNuclN
266 << " " << ProjectileNucleon->Get4Momentum() << "-------------" << G4endl
267 << "Targ N mom " << TrNuclN << " " << TargetNucleon->Get4Momentum() << G4endl
268 << "PrN TrN Z coords [fm]" << ProjectileNucleon->GetPosition().z()/fermi
269 << " " << TargetNucleon->GetPosition().z()/fermi
270 << " " << ProjectileNucleon->GetPosition().z()/fermi +
271 TargetNucleon->GetPosition().z()/fermi << G4endl;
272 #endif
273
274 if ( ! ProjectileNucleon->AreYouHit() ) {
275 // Projectile nucleon was not involved until now.
276 ProjectileSplitable = new G4DiffractiveSplitableHadron( *ProjectileNucleon );
277 ProjectileNucleon->Hit( ProjectileSplitable );
278 ProjectileSplitable->SetStatus( 1 ); // It takes part in the interaction
279 } else { // Projectile nucleon was involved before.
280 ProjectileSplitable = ProjectileNucleon->GetSplitableHadron();
281 }
282
283 if ( ! TargetNucleon->AreYouHit() ) { // Target nucleon was not involved until now
284 TargetSplitable = new G4DiffractiveSplitableHadron( *TargetNucleon );
285 TargetNucleon->Hit( TargetSplitable );
286 TargetSplitable->SetStatus( 1 ); // It takes part in the interaction
287 } else { // Target nucleon was involved before.
288 TargetSplitable = TargetNucleon->GetSplitableHadron();
289 }
290
291 G4InteractionContent* anInteraction = new G4InteractionContent( ProjectileSplitable );
292 anInteraction->SetTarget( TargetSplitable );
293 anInteraction->SetProjectileNucleon( ProjectileNucleon );
294 anInteraction->SetTargetNucleon( TargetNucleon );
295 anInteraction->SetInteractionTime( ( ProjectileNucleon->GetPosition().z() +
296 TargetNucleon->GetPosition().z() ) / betta_z );
297 anInteraction->SetStatus( 1 );
298
299 #ifdef debugFTFparticipant
300 G4cout << "Part anInteraction->GetInteractionTime() [fm] "
301 << anInteraction->GetInteractionTime()/fermi << G4endl
302 << "Splitable Pr* Tr* " << ProjectileSplitable << " "
303 << TargetSplitable << G4endl;
304 #endif
305
306 theInteractions.push_back( anInteraction );
307
308 } // End of an Interaction has happend!
309
310 #ifdef debugFTFparticipant
311 TrNuclN++;
312 #endif
313
314 } // End of while ( ( TargetNucleon = theNucleus->GetNextNucleon() ) )
315
316 #ifdef debugFTFparticipant
317 PrNuclN++;
318 #endif
319
320 } // End of while ( ( ProjectileNucleon = theProjectileNucleus->GetNextNucleon() ) )
321
322 if ( theInteractions.size() != 0 ) theProjectileNucleus->DoTranslation( theBeamPosition );
323
324 } while ( ( theInteractions.size() == 0 ) &&
325 ++loopCounter < maxNumberOfLoops ); /* Loop checking, 10.08.2015, A.Ribon */
326 if ( loopCounter >= maxNumberOfLoops ) {
327 #ifdef debugFTFparticipant
328 G4cout << "BAD situation: forced exit from the while loop!" << G4endl;
329 #endif
330 return;
331 }
332
335
336 #ifdef debugFTFparticipant
337 G4cout << G4endl << "Number of primary collisions " << theInteractions.size()
338 << "\t Bx[fm] " << impactX/fermi << "\t By[fm] " << impactY/fermi
339 << "\t B[fm] " << std::sqrt( sqr( impactX ) + sqr( impactY ) )/fermi << G4endl
340 << "FTF participant End. #######################" << G4endl << G4endl;
341 #endif
342 return;
343}
G4double B(G4double temperature)
static constexpr double twopi
Definition: G4SIunits.hh:56
static constexpr double fermi
Definition: G4SIunits.hh:83
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
#define G4UniformRand()
Definition: Randomize.hh:52
double z() const
double x() const
double y() const
G4double GetProbabilityOfInteraction(const G4double impactsquare)
G4double GetBmin2() const
G4bool SampleBinInterval() const
void SetImpactParameter(const G4double b_value)
G4double GetBmax2() const
G4double GetInteractionTime() const
void SetTargetNucleon(G4Nucleon *aNucleon)
void SetTarget(G4VSplitableHadron *aTarget)
void SetStatus(G4int aValue)
void SetInteractionTime(G4double aValue)
void SetProjectileNucleon(G4Nucleon *aNucleon)
const G4ThreeVector & GetPosition() const
Definition: G4Nucleon.hh:140
G4bool AreYouHit() const
Definition: G4Nucleon.hh:98
G4VSplitableHadron * GetSplitableHadron() const
Definition: G4Nucleon.hh:97
virtual const G4LorentzVector & Get4Momentum() const
Definition: G4Nucleon.hh:72
void Hit(G4VSplitableHadron *aHit)
Definition: G4Nucleon.hh:91
G4double GetTotalEnergy() const
G4ThreeVector GetMomentum() const
virtual G4double GetOuterRadius()=0
virtual G4Nucleon * GetNextNucleon()=0
virtual G4bool StartLoop()=0
virtual void DoTranslation(const G4ThreeVector &theShift)=0
std::pair< G4double, G4double > ChooseImpactXandY(G4double maxImpact)
Definition: G4V3DNucleus.hh:87
G4V3DNucleus * theProjectileNucleus
G4V3DNucleus * theNucleus
void SetStatus(const G4int aStatus)
const G4ThreeVector & GetPosition() const
void SetPosition(const G4ThreeVector &aPosition)
G4bool nucleon(G4int ityp)
T sqr(const T &x)
Definition: templates.hh:128

References G4Nucleon::AreYouHit(), B(), G4V3DNucleus::ChooseImpactXandY(), G4V3DNucleus::DoTranslation(), fermi, G4cout, G4endl, G4UniformRand, G4Nucleon::Get4Momentum(), GetBmax2(), GetBmin2(), G4InteractionContent::GetInteractionTime(), G4ReactionProduct::GetMomentum(), G4V3DNucleus::GetNextNucleon(), G4V3DNucleus::GetOuterRadius(), G4VSplitableHadron::GetPosition(), G4Nucleon::GetPosition(), G4FTFParameters::GetProbabilityOfInteraction(), G4Nucleon::GetSplitableHadron(), G4ReactionProduct::GetTotalEnergy(), G4Nucleon::Hit(), G4InuclParticleNames::nucleon(), SampleBinInterval(), SetImpactParameter(), G4InteractionContent::SetInteractionTime(), G4VSplitableHadron::SetPosition(), G4InteractionContent::SetProjectileNucleon(), G4VSplitableHadron::SetStatus(), G4InteractionContent::SetStatus(), G4InteractionContent::SetTarget(), G4InteractionContent::SetTargetNucleon(), ShiftInteractionTime(), SortInteractionsIncT(), sqr(), StartLoop(), G4V3DNucleus::StartLoop(), theInteractions, G4VParticipants::theNucleus, G4VParticipants::theProjectileNucleus, twopi, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by G4FTFModel::GetStrings().

◆ GetProjectileNucleus()

G4V3DNucleus * G4VParticipants::GetProjectileNucleus ( ) const
inlinevirtualinherited

Reimplemented in G4QGSParticipants.

Definition at line 78 of file G4VParticipants.hh.

79{
81}

References G4VParticipants::theProjectileNucleus.

Referenced by G4FTFModel::GetProjectileNucleus(), and G4FTFModel::Init().

◆ GetWoundedNucleus()

G4V3DNucleus * G4VParticipants::GetWoundedNucleus ( ) const
inlinevirtualinherited

Definition at line 73 of file G4VParticipants.hh.

74{
75 return theNucleus;
76}

References G4VParticipants::theNucleus.

Referenced by G4FTFModel::GetTargetNucleus(), and G4FTFModel::GetWoundedNucleus().

◆ Init()

void G4VParticipants::Init ( G4int  theZ,
G4int  theA 
)
virtualinherited

Definition at line 54 of file G4VParticipants.cc.

55{
56 if ( theNucleus == nullptr ) theNucleus = new G4Fancy3DNucleus();
57 theNucleus->Init(theA, theZ);
59}
virtual void SortNucleonsIncZ()=0
virtual void Init(G4int theA, G4int theZ, G4int numberOfLambdas=0)=0

References G4V3DNucleus::Init(), G4V3DNucleus::SortNucleonsIncZ(), and G4VParticipants::theNucleus.

Referenced by G4FTFModel::Init().

◆ InitProjectileNucleus()

void G4VParticipants::InitProjectileNucleus ( G4int  theZ,
G4int  theA,
G4int  numberOfLambdasOrAntiLambdas = 0 
)
virtualinherited

Definition at line 67 of file G4VParticipants.cc.

69{
70 if ( theProjectileNucleus == nullptr )
72
73 theProjectileNucleus->Init(theA, theZ, numberOfLambdasOrAntiLambdas);
75}
virtual void SortNucleonsDecZ()=0

References G4V3DNucleus::Init(), G4V3DNucleus::SortNucleonsDecZ(), and G4VParticipants::theProjectileNucleus.

Referenced by G4FTFModel::Init().

◆ Next()

G4bool G4FTFParticipants::Next ( )
inline

Definition at line 127 of file G4FTFParticipants.hh.

127 {
128 return ++currentInteraction < static_cast< G4int >( theInteractions.size() );
129}

References theInteractions.

Referenced by G4FTFModel::BuildStrings(), G4FTFModel::ExciteParticipants(), G4FTFModel::GetStrings(), and G4FTFModel::PutOnMassShell().

◆ operator!=() [1/2]

G4bool G4FTFParticipants::operator!= ( const G4FTFParticipants right) const
delete

◆ operator!=() [2/2]

G4bool G4VParticipants::operator!= ( const G4VParticipants right) const
deleteinherited

◆ operator=()

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

◆ operator==() [1/2]

G4bool G4FTFParticipants::operator== ( const G4FTFParticipants right) const
delete

◆ operator==() [2/2]

G4bool G4VParticipants::operator== ( const G4VParticipants right) const
deleteinherited

◆ SampleBinInterval()

G4bool G4FTFParticipants::SampleBinInterval ( ) const
inline

Definition at line 111 of file G4FTFParticipants.hh.

111 {
112 return BinInterval;
113}

References BinInterval.

Referenced by GetList().

◆ SetBminBmax()

void G4FTFParticipants::SetBminBmax ( const G4double  bmin_value,
const G4double  bmax_value 
)
inline

Definition at line 103 of file G4FTFParticipants.hh.

103 {
104 BinInterval = false;
105 if ( bmin_value < 0.0 || bmax_value < 0.0 || bmax_value < bmin_value ) return;
106 BinInterval = true;
107 Bmin2 = bmin_value * bmin_value;
108 Bmax2 = bmax_value * bmax_value;
109}

References BinInterval, Bmax2, and Bmin2.

Referenced by G4FTFModel::Init().

◆ SetImpactParameter()

void G4FTFParticipants::SetImpactParameter ( const G4double  b_value)
inline

Definition at line 95 of file G4FTFParticipants.hh.

95 {
96 Bimpact = b_value;
97}

References Bimpact.

Referenced by GetList().

◆ SetNucleus()

void G4VParticipants::SetNucleus ( G4V3DNucleus aNucleus)
virtualinherited

Definition at line 61 of file G4VParticipants.cc.

62{
63 if (theNucleus != nullptr) delete theNucleus;
64 theNucleus = aNucleus;
65}

References G4VParticipants::theNucleus.

◆ SetProjectileNucleus()

void G4VParticipants::SetProjectileNucleus ( G4V3DNucleus aNucleus)
virtualinherited

Definition at line 77 of file G4VParticipants.cc.

78{
79 if (theProjectileNucleus != nullptr) delete theProjectileNucleus;
80 theProjectileNucleus = aNucleus;
81}

References G4VParticipants::theProjectileNucleus.

Referenced by G4QGSParticipants::BuildInteractions(), and G4FTFModel::Init().

◆ ShiftInteractionTime()

void G4FTFParticipants::ShiftInteractionTime ( )

Definition at line 364 of file G4FTFParticipants.cc.

364 {
365 G4double InitialTime = theInteractions[0]->GetInteractionTime();
366 for ( unsigned int i = 1; i < theInteractions.size(); i++ ) {
367 G4double InterTime = theInteractions[i]->GetInteractionTime() - InitialTime;
368 theInteractions[i]->SetInteractionTime( InterTime );
369 G4InteractionContent* aCollision = theInteractions[i];
370 G4VSplitableHadron* projectile = aCollision->GetProjectile();
371 G4VSplitableHadron* target = aCollision->GetTarget();
372 G4ThreeVector prPosition = projectile->GetPosition();
373 prPosition.setZ( target->GetPosition().z() );
374 projectile->SetPosition( prPosition );
375 projectile->SetTimeOfCreation( InterTime );
376 target->SetTimeOfCreation( InterTime );
377 }
378 return;
379}
void setZ(double)
G4VSplitableHadron * GetProjectile() const
G4VSplitableHadron * GetTarget() const
void SetTimeOfCreation(G4double aTime)

References G4VSplitableHadron::GetPosition(), G4InteractionContent::GetProjectile(), G4InteractionContent::GetTarget(), G4VSplitableHadron::SetPosition(), G4VSplitableHadron::SetTimeOfCreation(), CLHEP::Hep3Vector::setZ(), theInteractions, and CLHEP::Hep3Vector::z().

Referenced by GetList().

◆ SortInteractionsIncT()

void G4FTFParticipants::SortInteractionsIncT ( )

Definition at line 356 of file G4FTFParticipants.cc.

356 { // on increased T
357 if ( theInteractions.size() < 2 ) return; // Avoid unnecesary work
358 std::sort( theInteractions.begin(), theInteractions.end(), G4FTFPartHelperForSortInT );
359}
bool G4FTFPartHelperForSortInT(const G4InteractionContent *Int1, const G4InteractionContent *Int2)

References G4FTFPartHelperForSortInT(), and theInteractions.

Referenced by GetList().

◆ StartLoop()

void G4FTFParticipants::StartLoop ( )
inline

Field Documentation

◆ Bimpact

G4double G4FTFParticipants::Bimpact
private

Definition at line 85 of file G4FTFParticipants.hh.

Referenced by GetImpactParameter(), and SetImpactParameter().

◆ BinInterval

G4bool G4FTFParticipants::BinInterval
private

Definition at line 86 of file G4FTFParticipants.hh.

Referenced by SampleBinInterval(), and SetBminBmax().

◆ Bmax2

G4double G4FTFParticipants::Bmax2
private

Definition at line 88 of file G4FTFParticipants.hh.

Referenced by GetBmax2(), and SetBminBmax().

◆ Bmin2

G4double G4FTFParticipants::Bmin2
private

Definition at line 87 of file G4FTFParticipants.hh.

Referenced by GetBmin2(), and SetBminBmax().

◆ currentInteraction

G4int G4FTFParticipants::currentInteraction
private

Definition at line 91 of file G4FTFParticipants.hh.

Referenced by Clean(), GetInteraction(), and StartLoop().

◆ theInteractions

std::vector< G4InteractionContent* > G4FTFParticipants::theInteractions
private

◆ theNucleus

G4V3DNucleus* G4VParticipants::theNucleus
protectedinherited

◆ theProjectileNucleus

G4V3DNucleus* G4VParticipants::theProjectileNucleus
protectedinherited

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