G4INCL::CoulombNone Class Reference

#include <G4INCLCoulombNone.hh>

Inheritance diagram for G4INCL::CoulombNone:

Public Member Functions
	CoulombNone ()
virtual	~CoulombNone ()
ParticleEntryAvatar *	bringToSurface (Particle *const p, Nucleus *const n) const
	Position the particle on the surface of the nucleus. More...
IAvatarList	bringToSurface (Cluster *const c, Nucleus *const n) const
	Position the cluster on the surface of the nucleus. More...
void	distortOut (ParticleList const &, Nucleus const *const ) const
	Modify the momenta of the outgoing particles. More...
G4double	maxImpactParameter (ParticleSpecies const &p, const G4double, Nucleus const *const n) const
	Return the maximum impact parameter for Coulomb-distorted trajectories. More...
Public Member Functions inherited from G4INCL::ICoulomb
	ICoulomb ()
virtual	~ICoulomb ()

Detailed Description

Definition at line 54 of file G4INCLCoulombNone.hh.

Constructor & Destructor Documentation

G4INCL::CoulombNone::CoulombNone ( )

inline

Definition at line 57 of file G4INCLCoulombNone.hh.

{}

virtual G4INCL::CoulombNone::~CoulombNone ( )

inlinevirtual

Definition at line 58 of file G4INCLCoulombNone.hh.

{}

Member Function Documentation

ParticleEntryAvatar * G4INCL::CoulombNone::bringToSurface ( Particle *const  p, Nucleus *const  n  ) const

virtual

Position the particle on the surface of the nucleus.

This method does not perform any distortion.

Parameters

p	incoming particle
n	distorting nucleus

Implements G4INCL::ICoulomb.

Definition at line 49 of file G4INCLCoulombNone.cc.

References G4INCL::Intersection::exists, G4INCL::IntersectionFactory::getEarlierTrajectoryIntersection(), G4INCL::Particle::getPosition(), G4INCL::Particle::getPropagationVelocity(), G4INCL::Nucleus::getUniverseRadius(), G4INCL::Intersection::position, and G4INCL::Particle::setPosition().

Referenced by G4INCL::CoulombNonRelativistic::bringToSurface().

                                                                                              {
    Intersection intersection = IntersectionFactory::getEarlierTrajectoryIntersection(p->getPosition(), p->getPropagationVelocity(), n->getUniverseRadius());
    if(intersection.exists) { // If the particle enters the nucleus
      p->setPosition(intersection.position);
      return new ParticleEntryAvatar(0.0, n, p);
    } else // If the particle does NOT enter the nucleus
      return NULL;
  }

IAvatarList G4INCL::CoulombNone::bringToSurface ( Cluster *const  c, Nucleus *const  n  ) const

virtual

Position the cluster on the surface of the nucleus.

This method does not perform any distortion.

Parameters

c	incoming cluster
n	distorting nucleus

Implements G4INCL::ICoulomb.

Definition at line 58 of file G4INCLCoulombNone.cc.

References G4INCL::Intersection::exists, G4INCL::IntersectionFactory::getEarlierTrajectoryIntersection(), G4INCL::Cluster::getParticles(), G4INCL::Nucleus::getUniverseRadius(), G4INCL::Intersection::position, and G4INCL::Intersection::time.

                                                                                    {
    // The avatar list that we will return
    IAvatarList theAvatarList;

    // Loop over the particles in the cluster
    ParticleList projectiles = c->getParticles();
    std::list<Intersection> theIntersections;
    G4double theFirstEntryTime = 1E+60; // a large time
    G4int theFirstID = 0;
    for(ParticleIter p=projectiles.begin(), e=projectiles.end(); p!=e; ++p) {
      // Check if the particle enters the nucleus
      Intersection intersection(IntersectionFactory::getEarlierTrajectoryIntersection(
            (*p)->getPosition(),
            (*p)->getPropagationVelocity(),
            n->getUniverseRadius()));
      // Store the intersections
      theIntersections.push_back(intersection);
      if(intersection.exists) {
        // Position the particle at the entry point
        (*p)->setPosition(intersection.position);

        // Keep track of the first entering particle
        if(intersection.time < theFirstEntryTime) {
          theFirstEntryTime = intersection.time;
          theFirstID = (*p)->getID();
        }
      }
    }

    std::list<Intersection>::const_iterator intIter = theIntersections.begin();
    for(ParticleIter p=projectiles.begin(), e=projectiles.end(); p!=e; ++p, ++intIter) {

      if((*intIter).exists) {
        // If the particle enters the nucleus, generate a ParticleEntryAvatar
        // for it and add it to the list of avatars that we will return
        if((*p)->getID() == theFirstID) {
          // The first particle always enters exactly at t=0 (in order to
          // avoid negative entry times due to rounding)
          theAvatarList.push_back(new ParticleEntryAvatar(0.0, n, *p));
        } else
          theAvatarList.push_back(new ParticleEntryAvatar(intIter->time - theFirstEntryTime, n, *p));
      }

   }

    return theAvatarList;
  }

void G4INCL::CoulombNone::distortOut ( ParticleList const &  , Nucleus const *  const  ) const

inlinevirtual

Modify the momenta of the outgoing particles.

This method does not perform any distortion.

Implements G4INCL::ICoulomb.

Definition at line 82 of file G4INCLCoulombNone.hh.

{}

G4double G4INCL::CoulombNone::maxImpactParameter ( ParticleSpecies const &  p, const G4double  , Nucleus const *const  n  ) const

inlinevirtual

Return the maximum impact parameter for Coulomb-distorted trajectories.

Implements G4INCL::ICoulomb.

Definition at line 86 of file G4INCLCoulombNone.hh.

References G4INCL::Composite, G4INCL::ParticleTable::getLargestNuclearRadius(), G4INCL::Nucleus::getUniverseRadius(), G4INCL::ParticleSpecies::theA, G4INCL::ParticleSpecies::theType, and G4INCL::ParticleSpecies::theZ.

                       {
      if(p.theType == Composite)
        return 2.*ParticleTable::getLargestNuclearRadius(p.theA, p.theZ)
          + n->getUniverseRadius();
      else
        return n->getUniverseRadius();
    }

---

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

• G4INCLCoulombNone.hh
• G4INCLCoulombNone.cc