Set of optical functions corresponding to one scoring plane along LHC, including splines for interpolation performance. More...

#include <LHCInterpolatedOpticalFunctionsSet.h>

Inheritance diagram for LHCInterpolatedOpticalFunctionsSet:

Classes
struct	Kinematics
	proton kinematics description More...

Public Member Functions
void	initializeSplines ()
	builds splines from m_*_values fields More...

	LHCInterpolatedOpticalFunctionsSet ()=default

	LHCInterpolatedOpticalFunctionsSet (const LHCOpticalFunctionsSet &src)

const std::vector< std::shared_ptr< const TSpline3 > > &	splines () const

void	transport (const Kinematics &input, Kinematics &output, bool calculateAngles=false) const
	transports proton according to the splines More...

	~LHCInterpolatedOpticalFunctionsSet ()=default

Public Member Functions inherited from LHCOpticalFunctionsSet
const std::vector< std::vector< double > > &	getFcnValues () const

double	getScoringPlaneZ () const
	returns the position of the scoring plane (LHC/TOTEM convention) More...

const std::vector< double > &	getXiValues () const

	LHCOpticalFunctionsSet ()=default

	LHCOpticalFunctionsSet (const std::string &fileName, const std::string &directoryName, double z)
	fills m_*_values fields from a ROOT file More...

	~LHCOpticalFunctionsSet ()=default

Protected Attributes
friend	CTPPSInterpolatedOpticalFunctionsESSource

friend	CTPPSModifiedOpticalFunctionsESSource

std::vector< std::shared_ptr< const TSpline3 > >	m_splines

Protected Attributes inherited from LHCOpticalFunctionsSet
std::vector< std::vector< double > >	m_fcn_values
	length unit cm More...

std::vector< double >	m_xi_values

double	m_z
	position of the scoring plane, in LHC/TOTEM convention, cm More...

Additional Inherited Members
Public Types inherited from LHCOpticalFunctionsSet
enum	{ evx, eLx, e14, exd, evpx, eLpx, e24, expd, e32, evy, eLy, eyd, e42, evpy, eLpy, eypd, nFunctions }
	indices for m_fcn_values and m_splines data members More...

Detailed Description

Set of optical functions corresponding to one scoring plane along LHC, including splines for interpolation performance.

Definition at line 14 of file LHCInterpolatedOpticalFunctionsSet.h.

Constructor & Destructor Documentation

LHCInterpolatedOpticalFunctionsSet::LHCInterpolatedOpticalFunctionsSet ( )

default

LHCInterpolatedOpticalFunctionsSet::LHCInterpolatedOpticalFunctionsSet ( const LHCOpticalFunctionsSet & src )

inline

Definition at line 19 of file LHCInterpolatedOpticalFunctionsSet.h.

References ~LHCInterpolatedOpticalFunctionsSet().

19 : LHCOpticalFunctionsSet(src) {}

LHCOpticalFunctionsSet::LHCOpticalFunctionsSet

LHCOpticalFunctionsSet()=default

LHCInterpolatedOpticalFunctionsSet::~LHCInterpolatedOpticalFunctionsSet ( )

default

Referenced by LHCInterpolatedOpticalFunctionsSet().

Member Function Documentation

void LHCInterpolatedOpticalFunctionsSet::initializeSplines ( )

builds splines from m_*_values fields

Definition at line 7 of file LHCInterpolatedOpticalFunctionsSet.cc.

References mps_fire::i, LHCOpticalFunctionsSet::m_fcn_values, m_splines, and LHCOpticalFunctionsSet::m_xi_values.

Referenced by CTPPSInterpolatedOpticalFunctionsESSource::produce(), and splines().

 {
   const unsigned int num_xi_vals = m_xi_values.size();
 
   m_splines.resize(m_fcn_values.size());
   for (unsigned int i = 0; i < m_fcn_values.size(); ++i)
     m_splines[i] = std::make_shared<TSpline3>("", m_xi_values.data(), m_fcn_values[i].data(), num_xi_vals);
 }

const std::vector<std::shared_ptr<const TSpline3> >& LHCInterpolatedOpticalFunctionsSet::splines ( ) const

inline

Definition at line 23 of file LHCInterpolatedOpticalFunctionsSet.h.

References initializeSplines(), and m_splines.

Referenced by ProtonReconstructionAlgorithm::init().

23 { return m_splines; }

LHCInterpolatedOpticalFunctionsSet::m_splines

std::vector< std::shared_ptr< const TSpline3 > > m_splines

Definition: LHCInterpolatedOpticalFunctionsSet.h:45

void LHCInterpolatedOpticalFunctionsSet::transport	(	const Kinematics &	input,
		LHCInterpolatedOpticalFunctionsSet::Kinematics &	output,
		bool	calculateAngles = `false`
	)		const

transports proton according to the splines

Definition at line 18 of file LHCInterpolatedOpticalFunctionsSet.cc.

Referenced by CTPPSOpticsPlotter::analyze(), and ProtonReconstructionAlgorithm::init().

 {
   const double xi = input.xi;
 
   output.x = m_splines[exd]->Eval(xi) + m_splines[evx]->Eval(xi) * input.x
     + m_splines[eLx]->Eval(xi) * input.th_x + m_splines[e14]->Eval(xi) * input.th_y;
 
   output.th_x = (!calculateAngles) ? 0. : m_splines[expd]->Eval(xi) + m_splines[evpx]->Eval(xi) * input.x
     + m_splines[eLpx]->Eval(xi) * input.th_x + m_splines[e24]->Eval(xi) * input.th_y;
 
   output.y = m_splines[eyd]->Eval(xi) + m_splines[evy]->Eval(xi) * input.y
     + m_splines[eLy]->Eval(xi) * input.th_y + m_splines[e32]->Eval(xi) * input.th_x;
 
   output.th_y = (!calculateAngles) ? 0. : m_splines[eypd]->Eval(xi) + m_splines[evpy]->Eval(xi) * input.y
     + m_splines[eLpy]->Eval(xi) * input.th_y + m_splines[e42]->Eval(xi) * input.th_x;
 
   output.xi = input.xi;
 }