#include <G4CollisionComposite.hh>

Inheritance diagram for G4CollisionComposite:

Classes
struct	Register

struct	Resolve

Public Member Functions
	G4CollisionComposite ()

virtual	~G4CollisionComposite ()

virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

virtual G4KineticTrackVector *	FinalState (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

virtual G4bool	IsInCharge (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

void	AddComponent (G4VCollision *aC)

virtual const G4VCrossSectionSource *	GetCrossSectionSource () const

virtual const G4VAngularDistribution *	GetAngularDistribution () const

virtual const G4CollisionVector *	GetComponents () const

Public Member Functions inherited from G4VCollision
	G4VCollision ()

void	establish_G4MT_TLS_G4VCollision ()

	G4VCollision (void s1, void s2, void s3, void s4, void s5, void s6, void *s7)

virtual	~G4VCollision ()

G4bool	operator== (const G4VCollision &right) const

G4bool	operator!= (const G4VCollision &right) const

virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

virtual G4KineticTrackVector *	FinalState (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0

virtual G4bool	IsInCharge (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0

virtual G4String	GetName () const =0

virtual void	Print () const

virtual void	Print (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

Additional Inherited Members
Protected Member Functions inherited from G4VCollision
G4int	GetNumberOfPartons (const G4ParticleDefinition *aP) const

virtual const G4CollisionVector *	GetComponents () const

virtual const G4VCrossSectionSource *	GetCrossSectionSource () const =0

virtual const G4VAngularDistribution *	GetAngularDistribution () const =0

virtual const std::vector< G4String > &	GetListOfColliders (G4int whichOne) const =0

Detailed Description

Definition at line 57 of file G4CollisionComposite.hh.

Constructor & Destructor Documentation

◆ G4CollisionComposite()

G4CollisionComposite::G4CollisionComposite ( )

Definition at line 44 of file G4CollisionComposite.cc.

{ 
  G4MUTEXINIT( bufferMutex );
}

◆ ~G4CollisionComposite()

G4CollisionComposite::~G4CollisionComposite ( )

virtual

Definition at line 50 of file G4CollisionComposite.cc.

{
  G4MUTEXDESTROY(bufferMutex);
  std::for_each(components.begin(), components.end(), G4Delete());
}

Member Function Documentation

◆ AddComponent()

void G4CollisionComposite::AddComponent ( G4VCollision * aC )

inline

Definition at line 72 of file G4CollisionComposite.hh.

72{components.push_back(aC);}

Referenced by G4CollisionMesonBaryon::G4CollisionMesonBaryon(), G4CollisionMesonBaryonToResonance::G4CollisionMesonBaryonToResonance(), G4CollisionComposite::Register::operator()(), and G4CollisionComposite::Resolve::operator()().

◆ CrossSection()

G4double G4CollisionComposite::CrossSection	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)		const

virtual

Reimplemented from G4VCollision.

Reimplemented in G4CollisionNN.

Definition at line 57 of file G4CollisionComposite.cc.

{
  G4double crossSect = 0.;
  const G4VCrossSectionSource* xSource = GetCrossSectionSource();
  if (xSource != 0)
  // There is a total cross section for this Collision
  {
    crossSect = xSource->CrossSection(trk1,trk2);
  }
  else
  {
    G4AutoLock l(&bufferMutex);
    // waiting for mutable to enable buffering.
    const_cast<G4CollisionComposite *>(this)->BufferCrossSection(trk1.GetDefinition(), trk2.GetDefinition());
//    G4cerr << "Buffer filled, reying with sqrts = "<< (trk1.Get4Momentum()+trk2.Get4Momentum()).mag() <<G4endl;
    crossSect = BufferedCrossSection(trk1,trk2);
  }
  return crossSect;
}

◆ FinalState()

G4KineticTrackVector * G4CollisionComposite::FinalState	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)		const

virtual

Implements G4VCollision.

Definition at line 79 of file G4CollisionComposite.cc.

{
  std::vector<G4double> cxCache;
  G4double partialCxSum = 0.0;
 
  size_t i;
  for (i=0; i<components.size(); i++) 
  {
    G4double partialCx;
//    cout << "comp" << i << " " << components[i]()->GetName();
    if (components[i]->IsInCharge(trk1,trk2)) 
    {
      partialCx = components[i]->CrossSection(trk1,trk2);
    } 
    else 
    {
      partialCx = 0.0;
    }
//    cout << "   cx=" << partialCx << endl;
    partialCxSum += partialCx;
    cxCache.push_back(partialCx);
  }
 
  G4double random = G4UniformRand()*partialCxSum;
  G4double running = 0;
  for (i=0; i<cxCache.size(); i++) 
  {
    running += cxCache[i];
    if (running > random) 
    {
      return components[i]->FinalState(trk1, trk2);
    }
  }
//  G4cerr <<"in charge = "<<IsInCharge(trk1, trk2)<<G4endl;
//  G4cerr <<"Cross-section = "<<CrossSection(trk1, trk2)/millibarn<<" "<<running<<" "<<cxCache.size()<<G4endl;
//  G4cerr <<"Names = "<<trk1.GetDefinition()->GetParticleName()<<", "<<trk2.GetDefinition()->GetParticleName()<<G4endl;
//  throw G4HadronicException(__FILE__, __LINE__, "G4CollisionComposite: no final state found!");
  return NULL;
}

◆ GetAngularDistribution()

virtual const G4VAngularDistribution * G4CollisionComposite::GetAngularDistribution ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionNN, and G4CollisionPN.

Definition at line 76 of file G4CollisionComposite.hh.

76{ return 0; }

◆ GetComponents()

virtual const G4CollisionVector * G4CollisionComposite::GetComponents ( ) const

inlinevirtual

Reimplemented from G4VCollision.

Definition at line 78 of file G4CollisionComposite.hh.

78{ return &components;}

Referenced by IsInCharge().

◆ GetCrossSectionSource()

virtual const G4VCrossSectionSource * G4CollisionComposite::GetCrossSectionSource ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionNN, and G4CollisionPN.

Definition at line 75 of file G4CollisionComposite.hh.

75{ return 0; }

Referenced by CrossSection().

◆ IsInCharge()

G4bool G4CollisionComposite::IsInCharge	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)		const

virtual

Implements G4VCollision.

Reimplemented in G4GeneralNNCollision.

Definition at line 121 of file G4CollisionComposite.cc.

{
  G4bool isInCharge = false;
 
  // The composite is in charge if any of its components is in charge
 
  const G4CollisionVector* comps = GetComponents();
  if (comps)
    {
      G4CollisionVector::const_iterator iter;
      for (iter = comps->begin(); iter != comps->end(); ++iter)
    {
     if ( ((*iter))->IsInCharge(trk1,trk2) ) isInCharge = true;
    }
    }
 
  return isInCharge;
}

Referenced by FinalState(), and IsInCharge().

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

Classes

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4CollisionComposite()

◆ ~G4CollisionComposite()

Member Function Documentation

◆ AddComponent()

◆ CrossSection()

◆ FinalState()

◆ GetAngularDistribution()

◆ GetComponents()

◆ GetCrossSectionSource()

◆ IsInCharge()