#include <G4CollisionComposite.hh>

Inheritance diagram for G4CollisionComposite:

Collaboration 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 G4String	GetName () const =0

virtual void	Print () const

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

Private Member Functions
	G4CollisionComposite (const G4CollisionComposite &right)

const G4CollisionComposite &	operator= (const G4CollisionComposite &right)

void	BufferCrossSection (const G4ParticleDefinition aP, const G4ParticleDefinition bP)

G4double	BufferedCrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

Private Attributes
G4CollisionVector	components

std::vector< G4CrossSectionBuffer >	theBuffer

G4Mutex	bufferMutex

Static Private Attributes
static const G4int	nPoints = 32

static const G4double	theT []

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

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

Detailed Description

Definition at line 58 of file G4CollisionComposite.hh.

Constructor & Destructor Documentation

◆ G4CollisionComposite() [1/2]

G4CollisionComposite::G4CollisionComposite ( )

Definition at line 45 of file G4CollisionComposite.cc.

 { 
   G4MUTEXINIT( bufferMutex );
 }

Here is the caller graph for this function:

◆ ~G4CollisionComposite()

G4CollisionComposite::~G4CollisionComposite ( )

virtual

Definition at line 51 of file G4CollisionComposite.cc.

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

◆ G4CollisionComposite() [2/2]

G4CollisionComposite::G4CollisionComposite ( const G4CollisionComposite & right )

private

Member Function Documentation

◆ AddComponent()

void G4CollisionComposite::AddComponent ( G4VCollision * aC )

inline

Definition at line 73 of file G4CollisionComposite.hh.

73 {components.push_back(aC);}

G4CollisionComposite::components

G4CollisionVector components

Definition: G4CollisionComposite.hh:121

Here is the caller graph for this function:

◆ BufferCrossSection()

void G4CollisionComposite::BufferCrossSection	(	const G4ParticleDefinition *	aP,
		const G4ParticleDefinition *	bP
	)

private

Definition at line 143 of file G4CollisionComposite.cc.

 {
    // check if already buffered
    size_t i;
    for(i=0; i<theBuffer.size(); i++)
    {
      if(theBuffer[i].InCharge(aP, bP)) return;
    }
 //   G4cerr << "Buffering for "<<aP->GetParticleName()<<" "<<bP->GetParticleName()<<G4endl;
    
    // buffer the new one.
    G4CrossSectionBuffer aNewBuff(aP, bP);
    size_t maxE=nPoints;
    for(size_t tt=0; tt<maxE; tt++)
    {
      G4double aT = theT[tt]*GeV;
      G4double crossSect = 0;
      // The total cross-section is summed over all the component channels
      
      //A.R. 28-Sep-2012 Fix reproducibility problem
      //                 Assign the kinetic energy to the lightest of the
      //                 two particles, instead to the first one always.
      G4double atime = 0;
      G4double btime = 0;
      G4ThreeVector aPosition(0,0,0);
      G4ThreeVector bPosition(0,0,0);
      G4double aM = aP->GetPDGMass();
      G4double bM = bP->GetPDGMass();
      G4double aE = aM;
      G4double bE = bM;
      G4ThreeVector aMom(0,0,0);
      G4ThreeVector bMom(0,0,0);
      if ( aM <= bM ) {
       aE += aT;
       aMom = G4ThreeVector(0,0,std::sqrt(aE*aE-aM*aM));
      } else {
       bE += aT;
       bMom = G4ThreeVector(0,0,std::sqrt(bE*bE-bM*bM));
      }
      G4LorentzVector a4Momentum(aE, aMom);
      G4LorentzVector b4Momentum(bE, bMom);
      G4KineticTrack a(aP, atime, aPosition, a4Momentum);
      G4KineticTrack b(bP, btime, bPosition, b4Momentum);
      
      for (i=0; i<components.size(); i++)
      {
        if(components[i]->IsInCharge(a,b))
        {
      crossSect += components[i]->CrossSection(a,b);
        }
      }
      G4double sqrts = (a4Momentum+b4Momentum).mag();
      aNewBuff.push_back(sqrts, crossSect);
    }
    theBuffer.push_back(aNewBuff);
 //   theBuffer.back().Print();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ BufferedCrossSection()

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

private

Definition at line 203 of file G4CollisionComposite.cc.

 {
    for(size_t i=0; i<theBuffer.size(); i++)
    {
      if(theBuffer[i].InCharge(trk1.GetDefinition(), trk2.GetDefinition())) 
      {
        return theBuffer[i].CrossSection(trk1, trk2);
      }
    }
    throw G4HadronicException(__FILE__, __LINE__, "G4CollisionComposite::BufferedCrossSection - Blitz !!");
    return 0;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ CrossSection()

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

virtual

Reimplemented from G4VCollision.

Reimplemented in G4CollisionMesonBaryonToResonance, and G4CollisionNN.

Definition at line 58 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;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ FinalState()

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

virtual

Implements G4VCollision.

Definition at line 80 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;
 }

Here is the call graph for this function:

◆ GetAngularDistribution()

virtual const G4VAngularDistribution* G4CollisionComposite::GetAngularDistribution ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionPN, and G4CollisionNN.

Definition at line 77 of file G4CollisionComposite.hh.

77 { return 0; }

◆ GetComponents()

virtual const G4CollisionVector* G4CollisionComposite::GetComponents ( ) const

inlinevirtual

Reimplemented from G4VCollision.

Definition at line 79 of file G4CollisionComposite.hh.

79 { return &components;}

G4CollisionComposite::components

G4CollisionVector components

Definition: G4CollisionComposite.hh:121

Here is the caller graph for this function:

◆ GetCrossSectionSource()

virtual const G4VCrossSectionSource* G4CollisionComposite::GetCrossSectionSource ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionPN, and G4CollisionNN.

Definition at line 76 of file G4CollisionComposite.hh.

76 { return 0; }

Here is the caller graph for this function:

◆ IsInCharge()

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

virtual

Implements G4VCollision.

Reimplemented in G4GeneralNNCollision.

Definition at line 122 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;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ operator=()

const G4CollisionComposite& G4CollisionComposite::operator= ( const G4CollisionComposite & right )

private

Here is the caller graph for this function:

Member Data Documentation

◆ bufferMutex

G4Mutex G4CollisionComposite::bufferMutex

mutableprivate

Definition at line 127 of file G4CollisionComposite.hh.

◆ components

G4CollisionVector G4CollisionComposite::components

private

Definition at line 121 of file G4CollisionComposite.hh.

◆ nPoints

const G4int G4CollisionComposite::nPoints = 32

staticprivate

Definition at line 124 of file G4CollisionComposite.hh.

◆ theBuffer

std::vector<G4CrossSectionBuffer> G4CollisionComposite::theBuffer

private

Definition at line 122 of file G4CollisionComposite.hh.

◆ theT

const G4double G4CollisionComposite::theT

staticprivate

Initial value:

=

{.01, .03, .05, .1, .15, .2, .3, .4, .5, .6, .7, .8, .9, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.5, 3.0, 3.5, 4.0, 5.0, 6.0, 8.0, 10., 15, 20, 50, 100}

Definition at line 125 of file G4CollisionComposite.hh.

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

Classes

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4CollisionComposite() [1/2]

◆ ~G4CollisionComposite()

◆ G4CollisionComposite() [2/2]

Member Function Documentation

◆ AddComponent()

◆ BufferCrossSection()

◆ BufferedCrossSection()

◆ CrossSection()

◆ FinalState()

◆ GetAngularDistribution()

◆ GetComponents()

◆ GetCrossSectionSource()

◆ IsInCharge()

◆ operator=()

Member Data Documentation

◆ bufferMutex

◆ components

◆ nPoints

◆ theBuffer

◆ theT