#include <ScatteringPhysics.h>

Inheritance diagram for ScatteringPhysics:

Collaboration diagram for ScatteringPhysics:

[legend]

Public Member Functions
	ScatteringPhysics (const std::string &name, ElementBase *element, std::string &mat, bool enableRutherford, double lowEnergyThr)

	~ScatteringPhysics ()

virtual void	apply (PartBunchBase< double, 3 > *bunch, const std::pair< Vector_t, double > &boundingSphere)

virtual const std::string	getType () const

virtual void	print (Inform &os)

virtual bool	stillActive ()

virtual double	getTime ()

virtual std::string	getName ()

virtual size_t	getParticlesInMat ()

virtual unsigned	getRediffused ()

virtual unsigned int	getNumEntered ()

void	computeInteraction (PartBunchBase< double, 3 > *bunch)

virtual bool	computeEnergyLoss (PartBunchBase< double, 3 > *bunch, Vector_t &P, const double deltat, bool includeFluctuations=true) const

Public Member Functions inherited from ParticleMatterInteractionHandler
	ParticleMatterInteractionHandler (std::string name, ElementBase *elref)

virtual	~ParticleMatterInteractionHandler ()

void	setFlagAllParticlesIn (bool p)

bool	getFlagAllParticlesIn () const

void	updateElement (ElementBase *newref)

ElementBase *	getElement ()

Private Member Functions
void	configureMaterialParameters ()
	The material of the collimator. More...

void	computeCoulombScattering (Vector_t &R, Vector_t &P, double dt)

void	applyRotation (Vector_t &P, Vector_t &R, double xplane, double thetacou)

void	applyRandomRotation (Vector_t &P, double theta0)

void	copyFromBunch (PartBunchBase< double, 3 > *bunch, const std::pair< Vector_t, double > &boundingSphere)

void	addBackToBunch (PartBunchBase< double, 3 > *bunch)

void	deleteParticleFromLocalVector ()

void	calcStat (double Eng)

void	gatherStatistics ()

void	push ()

void	resetTimeStep ()

void	setTimeStepForLeavingParticles ()

Private Attributes
double	T_m

double	dT_m

double	mass_m

double	charge_m

gsl_rng *	rGen_m

std::string	material_m

std::unique_ptr< InsideTester >	hitTester_m

ElementBase::ElementType	collshape_m

std::string	collshapeStr_m

double	Z_m

double	A_m

double	rho_m

double	X0_m

double	I_m

double	A1_c

double	A2_c

double	A3_c

double	A4_c

double	A5_c

double	B1_c

double	B2_c

double	B3_c

double	B4_c

double	B5_c

unsigned int	bunchToMatStat_m

unsigned int	stoppedPartStat_m

unsigned int	rediffusedStat_m

unsigned int	totalPartsInMat_m

double	Eavg_m

double	Emax_m

double	Emin_m

std::vector< PART >	locParts_m

std::unique_ptr< LossDataSink >	lossDs_m

bool	enableRutherford_m

double	lowEnergyThr_m

IpplTimings::TimerRef	DegraderApplyTimer_m

IpplTimings::TimerRef	DegraderLoopTimer_m

IpplTimings::TimerRef	DegraderDestroyTimer_m

Additional Inherited Members
Protected Attributes inherited from ParticleMatterInteractionHandler
ElementBase *	element_ref_m

bool	allParticleInMat_m
	if all particles are in matter stay inside the particle matter interaction More...

const std::string	name_m

Detailed Description

Definition at line 66 of file ScatteringPhysics.h.

Constructor & Destructor Documentation

◆ ScatteringPhysics()

ScatteringPhysics::ScatteringPhysics	(	const std::string &	name,
		ElementBase *	element,
		std::string &	mat,
		bool	enableRutherford,
		double	lowEnergyThr
	)

Definition at line 87 of file ScatteringPhysics.cpp.

References Options::asciidump, ElementBase::CCOLLIMATOR, collshape_m, configureMaterialParameters(), ElementBase::DEGRADER, DegraderApplyTimer_m, DegraderDestroyTimer_m, DegraderLoopTimer_m, ParticleMatterInteractionHandler::element_ref_m, ElementBase::FLEXIBLECOLLIMATOR, getName(), IpplTimings::getTimer(), ElementBase::getType(), hitTester_m, lossDs_m, IpplInfo::myNode(), rGen_m, and Options::seed.

Here is the call graph for this function:

◆ ~ScatteringPhysics()

ScatteringPhysics::~ScatteringPhysics ( )

Definition at line 163 of file ScatteringPhysics.cpp.

References locParts_m, lossDs_m, and rGen_m.

Member Function Documentation

◆ addBackToBunch()

void ScatteringPhysics::addBackToBunch ( PartBunchBase< double, 3 > * bunch )

private

Definition at line 485 of file ScatteringPhysics.cpp.

Referenced by apply().

Here is the call graph for this function:

◆ apply()

void ScatteringPhysics::apply	(	PartBunchBase< double, 3 > *	bunch,
		const std::pair< Vector_t, double > &	boundingSphere
	)

virtual

Implements ParticleMatterInteractionHandler.

Definition at line 193 of file ScatteringPhysics.cpp.

References addBackToBunch(), ParticleMatterInteractionHandler::allParticleInMat_m, ALPHA, bunchToMatStat_m, charge_m, computeInteraction(), copyFromBunch(), DegraderApplyTimer_m, DegraderLoopTimer_m, DEUTERON, dT_m, Eavg_m, Emax_m, Emin_m, gatherStatistics(), PartBunchBase< T, Dim >::getdT(), PartBunchBase< T, Dim >::getM(), PartBunchBase< T, Dim >::getPType(), PartBunchBase< T, Dim >::getPTypeString(), PartBunchBase< T, Dim >::getQ(), PartBunchBase< T, Dim >::getT(), H2P, HMINUS, mass_m, MUON, PROTON, push(), rediffusedStat_m, resetTimeStep(), setTimeStepForLeavingParticles(), IpplTimings::startTimer(), stoppedPartStat_m, IpplTimings::stopTimer(), T_m, and totalPartsInMat_m.

Here is the call graph for this function:

◆ applyRandomRotation()

void ScatteringPhysics::applyRandomRotation	(	Vector_t &	P,
		double	theta0
	)

private

Definition at line 416 of file ScatteringPhysics.cpp.

References abs(), atan(), cos(), cross(), dot(), euclidean_norm(), rGen_m, sin(), sqrt(), Physics::two_pi, Attrib::Distribution::W, and X.

Referenced by computeCoulombScattering().

Here is the call graph for this function:

◆ applyRotation()

void ScatteringPhysics::applyRotation	(	Vector_t &	P,
		Vector_t &	R,
		double	xplane,
		double	thetacou
	)

private

Definition at line 400 of file ScatteringPhysics.cpp.

References atan2(), cos(), fmod(), Attrib::Distribution::R, sin(), and Physics::two_pi.

Referenced by computeCoulombScattering().

Here is the call graph for this function:

◆ calcStat()

void ScatteringPhysics::calcStat ( double Eng )

inlineprivate

Definition at line 183 of file ScatteringPhysics.h.

References Eavg_m, Emax_m, and Emin_m.

◆ computeCoulombScattering()

void ScatteringPhysics::computeCoulombScattering	(	Vector_t &	R,
		Vector_t &	P,
		double	dt
	)

private

Coulomb Scattering: Including Multiple Coulomb Scattering and large angle Rutherford Scattering. Using the distribution given in Classical Electrodynamics, by J. D. Jackson.

Definition at line 441 of file ScatteringPhysics.cpp.

References abs(), applyRandomRotation(), applyRotation(), Physics::c, charge_m, cos(), enableRutherford_m, euclidean_norm(), log(), mass_m, Physics::pi, pow(), Attrib::Distribution::R, rGen_m, CoordinateSystemTrafo::rotateFrom(), CoordinateSystemTrafo::rotateTo(), sin(), sqrt(), CoordinateSystemTrafo::transformFrom(), Physics::two_pi, and X0_m.

Referenced by computeInteraction().

Here is the call graph for this function:

◆ computeEnergyLoss()

bool ScatteringPhysics::computeEnergyLoss	(	PartBunchBase< double, 3 > *	bunch,
		Vector_t &	P,
		const double	deltat,
		bool	includeFluctuations = `true`
	)		const

virtual

Energy Loss: using the Bethe-Bloch equation. In low-energy region use Andersen-Ziegler fitting (only for protons and alpha) Energy straggling: For relatively thick absorbers such that the number of collisions is large, the energy loss distribution is shown to be Gaussian in form. See Particle Physics Booklet, chapter 'Passage of particles through matter' or Review of Particle Physics, DOI: 10.1103/PhysRevD.86.010001, page 329 ff

Implements ParticleMatterInteractionHandler.

Definition at line 315 of file ScatteringPhysics.cpp.

References A1_c, A2_c, A3_c, A4_c, A5_c, A_m, ALPHA, Physics::amu, Physics::Avo, B1_c, B2_c, B3_c, B4_c, B5_c, Physics::c, charge_m, Physics::e, endl(), euclidean_norm(), Util::getGamma(), PartBunchBase< T, Dim >::getPType(), I_m, INFOMSG, level4(), log(), Physics::m_e, Physics::m_p, mass_m, Physics::pi, pow(), PROTON, Physics::r_e, rGen_m, rho_m, sqrt(), Attrib::Legacy::Distribution::T, and Z_m.

Referenced by computeInteraction().

Here is the call graph for this function:

◆ computeInteraction()

void ScatteringPhysics::computeInteraction ( PartBunchBase< double, 3 > * bunch )

Definition at line 268 of file ScatteringPhysics.cpp.

References computeCoulombScattering(), computeEnergyLoss(), deleteParticleFromLocalVector(), hitTester_m, locParts_m, lossDs_m, Attrib::Distribution::R, stoppedPartStat_m, and T_m.

Referenced by apply().

Here is the call graph for this function:

◆ configureMaterialParameters()

void ScatteringPhysics::configureMaterialParameters ( )

private

The material of the collimator.

Definition at line 173 of file ScatteringPhysics.cpp.

References Physics::Material::A1, A1_c, Physics::Material::A2, A2_c, Physics::Material::A3, A3_c, Physics::Material::A4, A4_c, Physics::Material::A5, A5_c, A_m, Physics::Material::B1, B1_c, Physics::Material::B2, B2_c, Physics::Material::B3, B3_c, Physics::Material::B4, B4_c, Physics::Material::B5, B5_c, Physics::Material::getMaterial(), I_m, material_m, rho_m, X0_m, and Z_m.

Referenced by ScatteringPhysics().

Here is the call graph for this function:

◆ copyFromBunch()

void ScatteringPhysics::copyFromBunch	(	PartBunchBase< double, 3 > *	bunch,
		const std::pair< Vector_t, double > &	boundingSphere
	)

private

Definition at line 529 of file ScatteringPhysics.cpp.

Referenced by apply().

Here is the call graph for this function:

◆ deleteParticleFromLocalVector()

void ScatteringPhysics::deleteParticleFromLocalVector ( )

private

Definition at line 630 of file ScatteringPhysics.cpp.

References Eavg_m, Emax_m, Emin_m, and locParts_m.

Referenced by addBackToBunch(), and computeInteraction().

◆ gatherStatistics()

void ScatteringPhysics::gatherStatistics ( )

private

Definition at line 703 of file ScatteringPhysics.cpp.

References allreduce(), bunchToMatStat_m, locParts_m, rediffusedStat_m, stoppedPartStat_m, and totalPartsInMat_m.

Referenced by apply().

Here is the call graph for this function:

◆ getName()

std::string ScatteringPhysics::getName ( )

inlinevirtual

Implements ParticleMatterInteractionHandler.

Definition at line 197 of file ScatteringPhysics.h.

References ParticleMatterInteractionHandler::element_ref_m, ElementBase::getName(), and ParticleMatterInteractionHandler::name_m.

Referenced by ScatteringPhysics().

Here is the call graph for this function:

◆ getNumEntered()

unsigned int ScatteringPhysics::getNumEntered ( )

inlinevirtual

Implements ParticleMatterInteractionHandler.

Definition at line 212 of file ScatteringPhysics.h.

References bunchToMatStat_m.

◆ getParticlesInMat()

size_t ScatteringPhysics::getParticlesInMat ( )

inlinevirtual

Implements ParticleMatterInteractionHandler.

Definition at line 202 of file ScatteringPhysics.h.

References totalPartsInMat_m.

◆ getRediffused()

unsigned int ScatteringPhysics::getRediffused ( )

inlinevirtual

Implements ParticleMatterInteractionHandler.

Definition at line 207 of file ScatteringPhysics.h.

References rediffusedStat_m.

◆ getTime()

double ScatteringPhysics::getTime ( )

inlinevirtual

Implements ParticleMatterInteractionHandler.

Definition at line 192 of file ScatteringPhysics.h.

References T_m.

◆ getType()

const std::string ScatteringPhysics::getType ( ) const

inlinevirtual

Implements ParticleMatterInteractionHandler.

Definition at line 217 of file ScatteringPhysics.h.

◆ print()

void ScatteringPhysics::print ( Inform & os )

virtual

Implements ParticleMatterInteractionHandler.

Definition at line 595 of file ScatteringPhysics.cpp.

References bunchToMatStat_m, ParticleMatterInteractionHandler::element_ref_m, endl(), Inform::flags(), ElementBase::getName(), level2(), material_m, ParticleMatterInteractionHandler::name_m, rediffusedStat_m, stoppedPartStat_m, OPALTimer::Timer::time(), Util::toStringWithThousandSep(), and totalPartsInMat_m.

Here is the call graph for this function:

◆ push()

void ScatteringPhysics::push ( )

private

Definition at line 655 of file ScatteringPhysics.cpp.

References Physics::c, dT_m, Util::getGamma(), locParts_m, and Attrib::Distribution::R.

Referenced by apply().

Here is the call graph for this function:

◆ resetTimeStep()

void ScatteringPhysics::resetTimeStep ( )

private

Definition at line 695 of file ScatteringPhysics.cpp.

References dT_m, and locParts_m.

Referenced by apply().

◆ setTimeStepForLeavingParticles()

void ScatteringPhysics::setTimeStepForLeavingParticles ( )

private

Definition at line 671 of file ScatteringPhysics.cpp.

References Physics::c, dT_m, ParticleMatterInteractionHandler::element_ref_m, ElementBase::getElementLength(), Util::getGamma(), locParts_m, PAssert_GE, PAssert_LE, and Attrib::Distribution::R.

Referenced by apply().

Here is the call graph for this function:

◆ stillActive()

bool ScatteringPhysics::stillActive ( )

virtual

Implements ParticleMatterInteractionHandler.

Definition at line 620 of file ScatteringPhysics.cpp.

References totalPartsInMat_m.

Member Data Documentation

◆ A1_c

double ScatteringPhysics::A1_c

private

Definition at line 145 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ A2_c

double ScatteringPhysics::A2_c

private

Definition at line 146 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ A3_c

double ScatteringPhysics::A3_c

private

Definition at line 147 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ A4_c

double ScatteringPhysics::A4_c

private

Definition at line 148 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ A5_c

double ScatteringPhysics::A5_c

private

Definition at line 149 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ A_m

double ScatteringPhysics::A_m

private

Definition at line 135 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ B1_c

double ScatteringPhysics::B1_c

private

Definition at line 150 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ B2_c

double ScatteringPhysics::B2_c

private

Definition at line 151 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ B3_c

double ScatteringPhysics::B3_c

private

Definition at line 152 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ B4_c

double ScatteringPhysics::B4_c

private

Definition at line 153 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ B5_c

double ScatteringPhysics::B5_c

private

Definition at line 154 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ bunchToMatStat_m

unsigned int ScatteringPhysics::bunchToMatStat_m

private

Definition at line 157 of file ScatteringPhysics.h.

Referenced by apply(), copyFromBunch(), gatherStatistics(), getNumEntered(), and print().

◆ charge_m

double ScatteringPhysics::charge_m

private

Definition at line 125 of file ScatteringPhysics.h.

Referenced by apply(), computeCoulombScattering(), and computeEnergyLoss().

◆ collshape_m

ElementBase::ElementType ScatteringPhysics::collshape_m

private

Definition at line 130 of file ScatteringPhysics.h.

Referenced by ScatteringPhysics().

◆ collshapeStr_m

std::string ScatteringPhysics::collshapeStr_m

private

Definition at line 131 of file ScatteringPhysics.h.

◆ DegraderApplyTimer_m

IpplTimings::TimerRef ScatteringPhysics::DegraderApplyTimer_m

private

Definition at line 177 of file ScatteringPhysics.h.

Referenced by apply(), and ScatteringPhysics().

◆ DegraderDestroyTimer_m

IpplTimings::TimerRef ScatteringPhysics::DegraderDestroyTimer_m

private

Definition at line 179 of file ScatteringPhysics.h.

Referenced by copyFromBunch(), and ScatteringPhysics().

◆ DegraderLoopTimer_m

IpplTimings::TimerRef ScatteringPhysics::DegraderLoopTimer_m

private

Definition at line 178 of file ScatteringPhysics.h.

Referenced by apply(), and ScatteringPhysics().

◆ dT_m

double ScatteringPhysics::dT_m

private

Definition at line 122 of file ScatteringPhysics.h.

Referenced by addBackToBunch(), apply(), push(), resetTimeStep(), and setTimeStepForLeavingParticles().

◆ Eavg_m

double ScatteringPhysics::Eavg_m

private

Definition at line 166 of file ScatteringPhysics.h.

Referenced by apply(), calcStat(), and deleteParticleFromLocalVector().

◆ Emax_m

double ScatteringPhysics::Emax_m

private

Definition at line 167 of file ScatteringPhysics.h.

Referenced by apply(), calcStat(), and deleteParticleFromLocalVector().

◆ Emin_m

double ScatteringPhysics::Emin_m

private

Definition at line 168 of file ScatteringPhysics.h.

Referenced by apply(), calcStat(), and deleteParticleFromLocalVector().

◆ enableRutherford_m

bool ScatteringPhysics::enableRutherford_m

private

Definition at line 174 of file ScatteringPhysics.h.

Referenced by computeCoulombScattering().

◆ hitTester_m

std::unique_ptr<InsideTester> ScatteringPhysics::hitTester_m

private

Definition at line 129 of file ScatteringPhysics.h.

Referenced by computeInteraction(), copyFromBunch(), and ScatteringPhysics().

◆ I_m

double ScatteringPhysics::I_m

private

Definition at line 138 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ locParts_m

std::vector<PART> ScatteringPhysics::locParts_m

private

Definition at line 170 of file ScatteringPhysics.h.

Referenced by addBackToBunch(), computeInteraction(), copyFromBunch(), deleteParticleFromLocalVector(), gatherStatistics(), push(), resetTimeStep(), setTimeStepForLeavingParticles(), and ~ScatteringPhysics().

◆ lossDs_m

std::unique_ptr<LossDataSink> ScatteringPhysics::lossDs_m

private

Definition at line 172 of file ScatteringPhysics.h.

Referenced by computeInteraction(), ScatteringPhysics(), and ~ScatteringPhysics().

◆ lowEnergyThr_m

double ScatteringPhysics::lowEnergyThr_m

private

Definition at line 175 of file ScatteringPhysics.h.

◆ mass_m

double ScatteringPhysics::mass_m

private

Definition at line 124 of file ScatteringPhysics.h.

Referenced by apply(), computeCoulombScattering(), and computeEnergyLoss().

◆ material_m

std::string ScatteringPhysics::material_m

private

Definition at line 128 of file ScatteringPhysics.h.

Referenced by configureMaterialParameters(), and print().

◆ rediffusedStat_m

unsigned int ScatteringPhysics::rediffusedStat_m

private

Definition at line 161 of file ScatteringPhysics.h.

Referenced by addBackToBunch(), apply(), gatherStatistics(), getRediffused(), and print().

◆ rGen_m

gsl_rng* ScatteringPhysics::rGen_m

private

Definition at line 127 of file ScatteringPhysics.h.

Referenced by applyRandomRotation(), computeCoulombScattering(), computeEnergyLoss(), ScatteringPhysics(), and ~ScatteringPhysics().

◆ rho_m

double ScatteringPhysics::rho_m

private

Definition at line 136 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

◆ stoppedPartStat_m

unsigned int ScatteringPhysics::stoppedPartStat_m

private

Definition at line 159 of file ScatteringPhysics.h.

Referenced by apply(), computeInteraction(), gatherStatistics(), and print().

◆ T_m

double ScatteringPhysics::T_m

private

Definition at line 121 of file ScatteringPhysics.h.

Referenced by apply(), computeInteraction(), and getTime().

◆ totalPartsInMat_m

unsigned int ScatteringPhysics::totalPartsInMat_m

private

Definition at line 163 of file ScatteringPhysics.h.

Referenced by apply(), gatherStatistics(), getParticlesInMat(), print(), and stillActive().

◆ X0_m

double ScatteringPhysics::X0_m

private

Definition at line 137 of file ScatteringPhysics.h.

Referenced by computeCoulombScattering(), and configureMaterialParameters().

◆ Z_m

double ScatteringPhysics::Z_m

private

Definition at line 134 of file ScatteringPhysics.h.

Referenced by computeEnergyLoss(), and configureMaterialParameters().

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

/Users/gsell/src/OPAL/src/src/Classic/Solvers/ScatteringPhysics.h
/Users/gsell/src/OPAL/src/src/Classic/Solvers/ScatteringPhysics.cpp

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ScatteringPhysics()

◆ ~ScatteringPhysics()

Member Function Documentation

◆ addBackToBunch()

◆ apply()

◆ applyRandomRotation()

◆ applyRotation()

◆ calcStat()

◆ computeCoulombScattering()

◆ computeEnergyLoss()

◆ computeInteraction()

◆ configureMaterialParameters()

◆ copyFromBunch()

◆ deleteParticleFromLocalVector()

◆ gatherStatistics()

◆ getName()

◆ getNumEntered()

◆ getParticlesInMat()

◆ getRediffused()

◆ getTime()

◆ getType()

◆ print()

◆ push()

◆ resetTimeStep()

◆ setTimeStepForLeavingParticles()

◆ stillActive()

Member Data Documentation

◆ A1_c

◆ A2_c

◆ A3_c

◆ A4_c

◆ A5_c

◆ A_m

◆ B1_c

◆ B2_c

◆ B3_c

◆ B4_c

◆ B5_c

◆ bunchToMatStat_m

◆ charge_m

◆ collshape_m

◆ collshapeStr_m

◆ DegraderApplyTimer_m

◆ DegraderDestroyTimer_m

◆ DegraderLoopTimer_m

◆ dT_m

◆ Eavg_m

◆ Emax_m

◆ Emin_m

◆ enableRutherford_m

◆ hitTester_m

◆ I_m

◆ locParts_m

◆ lossDs_m

◆ lowEnergyThr_m

◆ mass_m

◆ material_m

◆ rediffusedStat_m

◆ rGen_m

◆ rho_m

◆ stoppedPartStat_m

◆ T_m

◆ totalPartsInMat_m

◆ X0_m

◆ Z_m