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.
ParticleEntryAvatar *	bringToSurfaceAbar (Particle *const p, Nucleus *const n) const
	Position the particle on the surface of the nucleus. ONLY FOR ANTIDEUTERON !!! This method does not perform any distortion.
IAvatarList	bringToSurface (Cluster *const c, Nucleus *const n) const
	Position the cluster on the surface of the nucleus.
void	distortOut (ParticleList const &, Nucleus const *const) const
	Modify the momenta of the outgoing particles.
G4double	maxImpactParameter (ParticleSpecies const &p, const G4double, Nucleus const *const n) const
	Return the maximum impact parameter for Coulomb-distorted trajectories.
Public Member Functions inherited from G4INCL::ICoulomb
	ICoulomb ()
virtual	~ICoulomb ()

Detailed Description

Definition at line 55 of file G4INCLCoulombNone.hh.

Constructor & Destructor Documentation

CoulombNone()

G4INCL::CoulombNone::CoulombNone ( )

inline

Definition at line 58 of file G4INCLCoulombNone.hh.

{}

~CoulombNone()

virtual G4INCL::CoulombNone::~CoulombNone ( )

inlinevirtual

Definition at line 59 of file G4INCLCoulombNone.hh.

{}

Member Function Documentation

bringToSurface() [1/2]

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 68 of file G4INCLCoulombNone.cc.

                                                                                    {
    // The avatar list that we will return
    IAvatarList theAvatarList;
 
    // Loop over the particles in the cluster
    ParticleList const &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 = (G4int)(*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;
  }

bringToSurface() [2/2]

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 50 of file G4INCLCoulombNone.cc.

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

bringToSurfaceAbar()

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

virtual

Position the particle on the surface of the nucleus. ONLY FOR ANTIDEUTERON !!! This method does not perform any distortion.

Parameters

p	incoming particle
n	distorting nucleus

Implements G4INCL::ICoulomb.

Definition at line 59 of file G4INCLCoulombNone.cc.

                                                                                                  {
    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.001, n, p);
    } else // If the particle does NOT enter the nucleus
      return NULL;
  }

distortOut()

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 92 of file G4INCLCoulombNone.hh.

{}

maxImpactParameter()

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 96 of file G4INCLCoulombNone.hh.

                       {
      if(p.theType == Composite || p.theType == antiComposite)
        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