#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. More...

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. More...

void	distortOut (ParticleList const &pL, Nucleus const *const n) const
	Modify the momenta of the outgoing particles. More...

G4double	maxImpactParameter (ParticleSpecies const &p, const G4double kinE, 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 55 of file G4INCLCoulombNonRelativistic.hh.

Constructor & Destructor Documentation

G4INCL::CoulombNonRelativistic::CoulombNonRelativistic ( )

inline

Definition at line 57 of file G4INCLCoulombNonRelativistic.hh.

57 {}

virtual G4INCL::CoulombNonRelativistic::~CoulombNonRelativistic ( )

inlinevirtual

Definition at line 58 of file G4INCLCoulombNonRelativistic.hh.

58 {}

Member Function Documentation

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(), and G4INCL::Particle::getZ().

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

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().

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

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 plottest35::c1, G4INCL::ThreeVector::dot(), G4INCL::PhysicalConstants::eSquared, G4INCL::Nucleus::getDensity(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::Particle::getZ(), G4INCL::ThreeVector::mag(), and G4INCL::Math::piOverTwo.

                                            {
 
     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 = 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::getLargestNuclearRadius(), G4INCL::Nucleus::getUniverseRadius(), G4INCL::ParticleSpecies::theA, G4INCL::ParticleSpecies::theType, and G4INCL::ParticleSpecies::theZ.

                                                                                    {
     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:

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation