G4CollisionComposite Class Reference

#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

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::G4CollisionComposite

(

)

Definition at line 45 of file G4CollisionComposite.cc.

{ 
  G4MUTEXINIT( bufferMutex );
}

G4CollisionComposite::~G4CollisionComposite ( )

virtual

Definition at line 51 of file G4CollisionComposite.cc.

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

Member Function Documentation

void G4CollisionComposite::AddComponent ( G4VCollision *  aC )

inline

Definition at line 73 of file G4CollisionComposite.hh.

{components.push_back(aC);}

Here is the caller graph for this function:

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

virtual

Reimplemented from G4VCollision.

Reimplemented in 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:

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:

virtual const G4VAngularDistribution* G4CollisionComposite::GetAngularDistribution ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionPN, and G4CollisionNN.

Definition at line 77 of file G4CollisionComposite.hh.

{ return 0; }

virtual const G4CollisionVector* G4CollisionComposite::GetComponents ( ) const

inlinevirtual

Reimplemented from G4VCollision.

Reimplemented in G4CollisionNN, G4CollisionNNToDeltaDelta1700, G4CollisionNNToDeltaDelta1600, G4CollisionNNToDeltaDelta1620, G4CollisionNNToNDelta1620, G4CollisionNNToNDelta1700, G4CollisionNNToNDelta1900, G4CollisionNNToNDelta1905, G4CollisionNNToNDelta1910, G4CollisionNNToNDelta1920, G4CollisionNNToNDelta1930, G4CollisionNNToNDelta1950, G4CollisionNNToDeltaDelta1900, G4CollisionNNToDeltaDelta1905, G4CollisionNNToDeltaDelta1910, G4CollisionNNToDeltaDelta1920, G4CollisionNNToDeltaDelta1930, G4CollisionNNToDeltaDelta1950, G4CollisionNNToNDelta, and G4CollisionNNToNDelta1600.

Definition at line 79 of file G4CollisionComposite.hh.

{ return &components;}

Here is the caller graph for this function:

virtual const G4VCrossSectionSource* G4CollisionComposite::GetCrossSectionSource ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionPN, and G4CollisionNN.

Definition at line 76 of file G4CollisionComposite.hh.

{ return 0; }

Here is the caller graph for this function:

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:

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/im_r_matrix/include/G4CollisionComposite.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/im_r_matrix/src/G4CollisionComposite.cc