Detailed Description

Dedicated optimised inversion of range-energy function for KVedaLoss.

We generate a TGraph of particle energy as a function of range from the VEDALOSS tables, and then use a TSpline3 interpolation in order to calculate rapidly the particle energy for any range value.

See also: KVedaLoss

Definition at line 24 of file KVedaLossInverseRangeFunction.h.

#include <KVedaLossInverseRangeFunction.h>

Inheritance diagram for KVedaLossInverseRangeFunction:

[legend]

Public Member Functions
	KVedaLossInverseRangeFunction ()

	KVedaLossInverseRangeFunction (TF1 *range_func, Int_t A, Double_t riso, Int_t ninter=50)

virtual	~KVedaLossInverseRangeFunction ()

Double_t	GetEnergyPerNucleon (Double_t range, Double_t riso)

Private Attributes
std::unique_ptr< TSpline3 >	fInterpol
	interpolation of inverse range-energy curve More...

Constructor & Destructor Documentation

◆ KVedaLossInverseRangeFunction() [1/2]

KVedaLossInverseRangeFunction::KVedaLossInverseRangeFunction ( )

inline

Definition at line 33 of file KVedaLossInverseRangeFunction.h.

◆ KVedaLossInverseRangeFunction() [2/2]

KVedaLossInverseRangeFunction::KVedaLossInverseRangeFunction	(	TF1 *	range_func,
		Int_t	A,
		Double_t	riso,
		Int_t	ninter = `50`
	)

Initialise from VEDALOSS range function object

Parameters

A	is mass of ion used in call to KVedaLossMaterial::GetRangeFunction
riso	is value corresponding to previous call to KVedaLossMaterial::GetRangeFunction
ninter	number of points in graph used for interpolation

We make a graph from 0.01 MeV/nucleon to EMAX.

We divide the calculated range by the current value of riso, thereby effectively ensuring that we always calculate range values corresponding to riso=1

Definition at line 27 of file KVedaLossInverseRangeFunction.cpp.

◆ ~KVedaLossInverseRangeFunction()

virtual KVedaLossInverseRangeFunction::~KVedaLossInverseRangeFunction ( )

inlinevirtual

Definition at line 42 of file KVedaLossInverseRangeFunction.h.

Member Function Documentation

◆ GetEnergyPerNucleon()

Double_t KVedaLossInverseRangeFunction::GetEnergyPerNucleon	(	Double_t	range,
		Double_t	riso
	)

Parameters

range	of ion in \(g/cm^2\)
riso	takes into account any change in ion mass and/or material mass

Returns: corresponding energy in MeV/nucleon

Definition at line 78 of file KVedaLossInverseRangeFunction.cpp.

Member Data Documentation

◆ fInterpol

std::unique_ptr<TSpline3> KVedaLossInverseRangeFunction::fInterpol

private

interpolation of inverse range-energy curve

Definition at line 29 of file KVedaLossInverseRangeFunction.h.