G4INCL::CoulombNonRelativistic Class Reference

#include <G4INCLCoulombNonRelativistic.hh>

Inheritance diagram for G4INCL::CoulombNonRelativistic:

Public Member Functions
	CoulombNonRelativistic ()
virtual	~CoulombNonRelativistic ()
ParticleEntryAvatar *	bringToSurface (Particle *const p, Nucleus *const n) const
	Modify the momentum of the particle and position it on the surface of the nucleus.
IAvatarList	bringToSurface (Cluster *const c, Nucleus *const n) const
	Modify the momentum of the incoming cluster and position it on the surface of the nucleus.
void	distortOut (ParticleList const &pL, Nucleus const *const n) const
	Modify the momenta of the outgoing particles.
G4double	maxImpactParameter (ParticleSpecies const &p, const G4double kinE, 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 56 of file G4INCLCoulombNonRelativistic.hh.

Constructor & Destructor Documentation

CoulombNonRelativistic()

G4INCL::CoulombNonRelativistic::CoulombNonRelativistic ( )

inline

Definition at line 58 of file G4INCLCoulombNonRelativistic.hh.

{}

~CoulombNonRelativistic()

virtual G4INCL::CoulombNonRelativistic::~CoulombNonRelativistic ( )

inlinevirtual

Definition at line 59 of file G4INCLCoulombNonRelativistic.hh.

{}

Member Function Documentation

bringToSurface() [1/2]

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

virtual

Modify the momentum of the incoming cluster and position it on the surface of the nucleus.

This method performs non-relativistic distortion. The momenta of the particles that compose the cluster are also distorted.

Parameters

c	incoming cluster
n	distorting nucleus

Implements G4INCL::ICoulomb.

Definition at line 63 of file G4INCLCoulombNonRelativistic.cc.

                                                                                               {
    // Neutral clusters?!
// assert(c->getZ()>0);
 
    // Perform the actual Coulomb deviation
    const G4bool success = coulombDeviation(c, n);
    if(!success) {
      return IAvatarList();
    }
 
    // Rely on the CoulombNone slave to compute the straight-line intersection
    // and actually bring the particle to the surface of the nucleus
    return theCoulombNoneSlave.bringToSurface(c,n);
  }

bringToSurface() [2/2]

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

virtual

Modify the momentum of the particle and position it on the surface of the nucleus.

This method performs non-relativistic distortion.

Parameters

p	incoming particle
n	distorting nucleus

Implements G4INCL::ICoulomb.

Definition at line 50 of file G4INCLCoulombNonRelativistic.cc.

                                                                                                         {
    // No distortion for neutral particles
    if(p->getZ()!=0) {
      const G4bool success = coulombDeviation(p, n);
      if(!success) // transparent
        return NULL;
    }
 
    // Rely on the CoulombNone slave to compute the straight-line intersection
    // and actually bring the particle to the surface of the nucleus
    return theCoulombNoneSlave.bringToSurface(p,n);
  }

distortOut()

void G4INCL::CoulombNonRelativistic::distortOut ( ParticleList const & pL, Nucleus const *const n ) const

virtual

Modify the momenta of the outgoing particles.

This method performs non-relativistic distortion.

Parameters

pL	list of outgoing particles
n	distorting nucleus

Implements G4INCL::ICoulomb.

Definition at line 78 of file G4INCLCoulombNonRelativistic.cc.

                                           {
 
    for(ParticleIter particle=pL.begin(), e=pL.end(); particle!=e; ++particle) {
 
      const G4int Z = (*particle)->getZ();
      if(Z == 0) continue;
 
      const G4double tcos=1.-0.000001;
 
      const G4double et1 = PhysicalConstants::eSquared * nucleus->getZ();
      const G4double transmissionRadius =
        nucleus->getDensity()->getTransmissionRadius(*particle);
 
      const ThreeVector position = (*particle)->getPosition();
      ThreeVector momentum = (*particle)->getMomentum();
      const G4double r = position.mag();
      const G4double p = momentum.mag();
      const G4double cosTheta = position.dot(momentum)/(r*p);
      if(cosTheta < 0.999999) {
        const G4double sinTheta = std::sqrt(1.-cosTheta*cosTheta);
        const G4double eta = et1 * Z / (*particle)->getKineticEnergy();
        if(eta > transmissionRadius-0.0001) {
          // If below the Coulomb barrier, radial emission:
          momentum = position * (p/r);
          (*particle)->setMomentum(momentum);
        } else {
          const G4double b0 = 0.5 * (eta + std::sqrt(eta*eta +
                4. * std::pow(transmissionRadius*sinTheta,2)
                * (1.-eta/transmissionRadius)));
          const G4double bInf = std::sqrt(b0*(b0-eta));
          const G4double thr = std::atan(eta/(2.*bInf));
          G4double uTemp = (1.-b0/transmissionRadius) * std::sin(thr) +
            b0/transmissionRadius;
          if(uTemp>tcos) uTemp=tcos;
          const G4double thd = Math::arcCos(cosTheta)-Math::piOverTwo + thr +
            Math::arcCos(uTemp);
          const G4double c1 = std::sin(thd)*cosTheta/sinTheta + std::cos(thd);
          const G4double c2 = -p*std::sin(thd)/(r*sinTheta);
          const ThreeVector newMomentum = momentum*c1 + position*c2;
          (*particle)->setMomentum(newMomentum);
        }
      }
    }
  }

maxImpactParameter()

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

virtual

Return the maximum impact parameter for Coulomb-distorted trajectories.

Implements G4INCL::ICoulomb.

Definition at line 124 of file G4INCLCoulombNonRelativistic.cc.

                                                                                   {
    const G4double theMinimumDistance = minimumDistance(p, kinE, n);
    G4double rMax = n->getUniverseRadius();
    if(p.theType == Composite)
      rMax +=  2.*ParticleTable::getLargestNuclearRadius(p.theA, p.theZ);
    const G4double theMaxImpactParameterSquared = rMax*(rMax-theMinimumDistance);
    if(theMaxImpactParameterSquared<=0.)
      return 0.;
    const G4double theMaxImpactParameter = std::sqrt(theMaxImpactParameterSquared);
    return theMaxImpactParameter;
  }

---

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

• G4INCLCoulombNonRelativistic.hh
• G4INCLCoulombNonRelativistic.cc