#include <G4CollisionComposite.hh>

Inheritance diagram for G4CollisionComposite:

Data Structures
struct	Register

struct	Resolve

Public Member Functions
void	AddComponent (G4VCollision *aC)

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

void	establish_G4MT_TLS_G4VCollision ()

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

	G4CollisionComposite ()

virtual const G4VAngularDistribution *	GetAngularDistribution () const

virtual const G4CollisionVector *	GetComponents () const

virtual const G4VCrossSectionSource *	GetCrossSectionSource () const

virtual G4String	GetName () const =0

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

G4bool	operator!= (const G4VCollision &right) const

G4bool	operator== (const G4VCollision &right) const

virtual void	Print () const

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

virtual	~G4CollisionComposite ()

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

G4int	GetNumberOfPartons (const G4ParticleDefinition *aP) const

Private Member Functions
void	BufferCrossSection (const G4ParticleDefinition aP, const G4ParticleDefinition bP)

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

	G4CollisionComposite (const G4CollisionComposite &right)

const G4CollisionComposite &	operator= (const G4CollisionComposite &right)

Private Attributes
G4Mutex	bufferMutex

G4CollisionVector	components

std::vector< G4CrossSectionBuffer >	theBuffer

Static Private Attributes
static const G4int	nPoints = 32

static const G4double	theT []

Detailed Description

Definition at line 57 of file G4CollisionComposite.hh.

Constructor & Destructor Documentation

◆ G4CollisionComposite() [1/2]

G4CollisionComposite::G4CollisionComposite ( )

Definition at line 44 of file G4CollisionComposite.cc.

{ 
  G4MUTEXINIT( bufferMutex );
}

References bufferMutex, and G4MUTEXINIT.

◆ ~G4CollisionComposite()

G4CollisionComposite::~G4CollisionComposite ( )

virtual

Definition at line 50 of file G4CollisionComposite.cc.

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

References bufferMutex, components, and G4MUTEXDESTROY.

◆ G4CollisionComposite() [2/2]

G4CollisionComposite::G4CollisionComposite ( const G4CollisionComposite & right )

private

Member Function Documentation

◆ AddComponent()

void G4CollisionComposite::AddComponent ( G4VCollision * aC )

inline

Definition at line 72 of file G4CollisionComposite.hh.

72{components.push_back(aC);}

References components.

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

◆ BufferCrossSection()

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

private

Definition at line 141 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();
}

References components, G4ParticleDefinition::GetPDGMass(), GeV, IsInCharge(), nPoints, G4CrossSectionBuffer::push_back(), theBuffer, and theT.

Referenced by CrossSection().

◆ BufferedCrossSection()

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

private

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

References G4KineticTrack::GetDefinition(), and theBuffer.

Referenced by CrossSection().

◆ CrossSection()

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

virtual

Reimplemented from G4VCollision.

Reimplemented in G4CollisionMesonBaryonToResonance, and 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;
}

References BufferCrossSection(), BufferedCrossSection(), bufferMutex, G4VCrossSectionSource::CrossSection(), GetCrossSectionSource(), and G4KineticTrack::GetDefinition().

Referenced by G4CollisionMesonBaryonToResonance::CrossSection().

◆ establish_G4MT_TLS_G4VCollision()

void G4VCollision::establish_G4MT_TLS_G4VCollision ( )

inherited

Definition at line 143 of file G4VCollision.cc.

143{ }

Referenced by G4VScatteringCollision::establish_G4MT_TLS_G4VScatteringCollision().

◆ 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;
}

References components, G4UniformRand, and IsInCharge().

◆ 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;}

References 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().

◆ GetListOfColliders()

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

protectedpure virtualinherited

◆ GetName()

virtual G4String G4VCollision::GetName ( ) const

pure virtualinherited

Referenced by G4VCollision::Print().

◆ GetNumberOfPartons()

G4int G4VCollision::GetNumberOfPartons ( const G4ParticleDefinition * aP ) const

inlineprotectedinherited

Definition at line 68 of file G4VCollision.hh.

  {
    G4int result = 0;
    for(G4int i=0; i<6; i++) 
    {
      result += aP->GetQuarkContent(i+1);
      result += aP->GetAntiQuarkContent(i+1);
    }
    return result;
  }

References G4ParticleDefinition::GetAntiQuarkContent(), and G4ParticleDefinition::GetQuarkContent().

Referenced by G4CollisionMesonBaryonElastic::IsInCharge().

◆ 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;
}

References GetComponents(), and IsInCharge().

Referenced by BufferCrossSection(), FinalState(), and IsInCharge().

◆ operator!=()

G4bool G4VCollision::operator!= ( const G4VCollision & right ) const

inherited

Definition at line 48 of file G4VCollision.cc.

{
    return (this != (G4VCollision *) &right);
}

◆ operator=()

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

private

◆ operator==()

G4bool G4VCollision::operator== ( const G4VCollision & right ) const

inherited

Definition at line 42 of file G4VCollision.cc.

{
    return (this == (G4VCollision *) &right);
}

◆ Print() [1/2]

void G4VCollision::Print ( ) const

virtualinherited

Definition at line 70 of file G4VCollision.cc.

{
    G4String name = GetName();
 
    G4cout << "---- " << name << "---- Cross section" << G4endl;
 
    const G4VCrossSectionSource* xSource = GetCrossSectionSource();
    if (xSource) xSource->Print();
 
    G4int nComponents = 0;
    const G4CollisionVector* components = GetComponents();
    if (components)
    {
        nComponents = components->size();
    }
    G4cout << "---- " << name << "---- has " << nComponents << " components" <<G4endl;
    G4int i = 0;
    G4CollisionVector::const_iterator iter;
    if (components)
    {
        for (iter = components->begin(); iter != components->end(); ++iter)
        {
            G4cout << "---- " << name << " ---- Component " << i << G4endl;
            ((*iter))->Print();
            i++;
        }
    }
 
}

References G4cout, G4endl, G4VCollision::GetComponents(), G4VCollision::GetCrossSectionSource(), G4VCollision::GetName(), G4InuclParticleNames::name(), G4VCollision::Print(), and G4VCrossSectionSource::Print().

Referenced by G4VCollision::Print().

◆ Print() [2/2]

void G4VCollision::Print	(	const G4KineticTrack &	trk1,
		const G4KineticTrack &	trk2
	)		const

virtualinherited

Definition at line 101 of file G4VCollision.cc.

{
    G4String name = GetName();
 
    if (IsInCharge(trk1,trk2))
    {
        G4cout << "---- " << name << "is in charge ---- " << G4endl;
    }
    else
    {
        G4cout << "---- " << name << "is not in charge ---- " << G4endl;
    }
 
    G4cout << "---- " << name << "---- Cross section" << G4endl;
 
    const G4VCrossSectionSource* xSource = GetCrossSectionSource();
    if (xSource) xSource->Print();
    G4cout << "Cross section = " << CrossSection(trk1,trk2) << G4endl;
 
    G4int nComponents = 0;
    const G4CollisionVector* components = GetComponents();
    if (components)
    {
        nComponents = components->size();
    }
    G4cout << "---- " << name << "has " << nComponents << " components" <<G4endl;
 
    G4int i = 0;
    G4CollisionVector::const_iterator iter;
    if (components)
    {
        for (iter = components->begin(); iter != components->end(); ++iter)
        {
            G4cout << "Component " << i << G4endl;
            ((*iter))->Print();
            i++;
        }
    }
}

References G4VCollision::CrossSection(), G4cout, G4endl, G4VCollision::GetComponents(), G4VCollision::GetCrossSectionSource(), G4VCollision::GetName(), G4VCollision::IsInCharge(), G4InuclParticleNames::name(), G4VCollision::Print(), and G4VCrossSectionSource::Print().

Field Documentation

◆ bufferMutex

G4Mutex G4CollisionComposite::bufferMutex

mutableprivate

Definition at line 126 of file G4CollisionComposite.hh.

Referenced by CrossSection(), G4CollisionComposite(), and ~G4CollisionComposite().

◆ components

G4CollisionVector G4CollisionComposite::components

private

Definition at line 120 of file G4CollisionComposite.hh.

Referenced by AddComponent(), BufferCrossSection(), FinalState(), GetComponents(), and ~G4CollisionComposite().

◆ nPoints

const G4int G4CollisionComposite::nPoints = 32

staticprivate

Definition at line 123 of file G4CollisionComposite.hh.

Referenced by BufferCrossSection().

◆ theBuffer

std::vector<G4CrossSectionBuffer> G4CollisionComposite::theBuffer

private

Definition at line 121 of file G4CollisionComposite.hh.

Referenced by BufferCrossSection(), and BufferedCrossSection().

◆ 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 124 of file G4CollisionComposite.hh.

Referenced by BufferCrossSection().

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

source/processes/hadronic/models/im_r_matrix/include/G4CollisionComposite.hh
source/processes/hadronic/models/im_r_matrix/src/G4CollisionComposite.cc

Data Structures

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4CollisionComposite() [1/2]

◆ ~G4CollisionComposite()

◆ G4CollisionComposite() [2/2]

Member Function Documentation

◆ AddComponent()

◆ BufferCrossSection()

◆ BufferedCrossSection()

◆ CrossSection()

◆ establish_G4MT_TLS_G4VCollision()

◆ FinalState()

◆ GetAngularDistribution()

◆ GetComponents()

◆ GetCrossSectionSource()

◆ GetListOfColliders()

◆ GetName()

◆ GetNumberOfPartons()

◆ IsInCharge()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ Print() [1/2]

◆ Print() [2/2]

Field Documentation

◆ bufferMutex

◆ components

◆ nPoints

◆ theBuffer

◆ theT