G4CollisionComposite Class Reference

#include <G4CollisionComposite.hh>

Inheritance diagram for G4CollisionComposite:

Data Structures
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 (G4ParticleDefinition *aP) const
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

(

)

Definition at line 44 of file G4CollisionComposite.cc.

{ 
}

G4CollisionComposite::~G4CollisionComposite ( )

virtual

Definition at line 49 of file G4CollisionComposite.cc.

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

Member Function Documentation

void G4CollisionComposite::AddComponent ( G4VCollision *  aC )

inline

Definition at line 72 of file G4CollisionComposite.hh.

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

{components.push_back(aC);}

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

virtual

Reimplemented from G4VCollision.

Reimplemented in G4CollisionNN.

Definition at line 55 of file G4CollisionComposite.cc.

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

{
  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
  {
    // 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;
}

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

virtual

Implements G4VCollision.

Definition at line 76 of file G4CollisionComposite.cc.

References G4UniformRand, and IsInCharge().

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

virtual const G4VAngularDistribution* G4CollisionComposite::GetAngularDistribution ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionPN, and G4CollisionNN.

Definition at line 76 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 78 of file G4CollisionComposite.hh.

Referenced by IsInCharge().

{ return &components;}

virtual const G4VCrossSectionSource* G4CollisionComposite::GetCrossSectionSource ( ) const

inlinevirtual

Implements G4VCollision.

Reimplemented in G4CollisionPN, and G4CollisionNN.

Definition at line 75 of file G4CollisionComposite.hh.

Referenced by CrossSection().

{ return 0; }

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

virtual

Implements G4VCollision.

Reimplemented in G4GeneralNNCollision.

Definition at line 118 of file G4CollisionComposite.cc.

References GetComponents().

Referenced by FinalState().

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

---

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

• G4CollisionComposite.hh
• G4CollisionComposite.cc