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.

---

Detailed Description

Definition at line 55 of file G4INCLCoulombNonRelativistic.hh.

---

Constructor & Destructor Documentation

G4INCL::CoulombNonRelativistic::CoulombNonRelativistic

(

)

[inline]

Definition at line 57 of file G4INCLCoulombNonRelativistic.hh.

{}

virtual G4INCL::CoulombNonRelativistic::~CoulombNonRelativistic

(

)

[inline, virtual]

Definition at line 58 of file G4INCLCoulombNonRelativistic.hh.

{}

---

Member Function Documentation

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 62 of file G4INCLCoulombNonRelativistic.cc.

References G4INCL::CoulombNone::bringToSurface(), and CLHEP::detail::n.

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

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 49 of file G4INCLCoulombNonRelativistic.cc.

References G4INCL::CoulombNone::bringToSurface(), G4INCL::Particle::getZ(), and CLHEP::detail::n.

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

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 77 of file G4INCLCoulombNonRelativistic.cc.

References G4INCL::PhysicalConstants::eSquared, G4INCL::Nucleus::getDensity(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::Particle::getZ(), G4INCL::ThreeVector::mag(), G4INCL::Math::piOverTwo, and position.

                                           {

    for(ParticleIter particle=pL.begin(); particle!=pL.end(); ++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 = std::acos(cosTheta)-Math::piOverTwo + thr +
            std::acos(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);
        }
      }
    }
  }

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 123 of file G4INCLCoulombNonRelativistic.cc.

References G4INCL::Composite, G4INCL::ParticleTable::getNuclearRadius(), CLHEP::detail::n, G4INCL::ParticleSpecies::theA, G4INCL::ParticleSpecies::theType, and G4INCL::ParticleSpecies::theZ.

                                     {
    G4double theMaxImpactParameter = maxImpactParameterParticle(p, kinE, n);
    if(theMaxImpactParameter <= 0.)
      return 0.;
    if(p.theType == Composite)
      theMaxImpactParameter +=  2.*ParticleTable::getNuclearRadius(p.theA, p.theZ);
    return theMaxImpactParameter;
  }

---

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

• G4INCLCoulombNonRelativistic.hh
• G4INCLCoulombNonRelativistic.cc

---