Class representing transport parameterisation as four Equation. More...

#include <Parameterization.h>

Collaboration diagram for AFP::Parameterization:

Classes
class	Equation
	Local class for storing parameterization equation. More...

Public Member Functions
	Parameterization (std::string_view filename)
	Takes filename of parametrization file as argument. More...

double	x (double x0, double y0, double z0, double sx0, double sy0, double E) const
	Evaluates value of horizontal position equation. More...

double	y (double x0, double y0, double z0, double sx0, double sy0, double E) const
	Evaluates value of vertical position equation. More...

double	sx (double x0, double y0, double z0, double sx0, double sy0, double E) const
	Evaluates value of horizontal slope equation. More...

double	sy (double x0, double y0, double z0, double sx0, double sy0, double E) const
	Evaluates value of vertical slope equation. More...

double	parametrizationPosition () const
	Returns position for which parametrization was performed. More...

double	energy () const
	Returns nominal energy for which parametrization was performed. More...

const Equation *	getEquation (int iEqu) const
	Returns const raw pointer to equation of number iEqu. More...

bool	yIsUsed () const
	Returns m_yIsUsed. More...

Private Attributes
std::array< double, 2 >	m_slopeInitIP
	Crossing angles in rad. More...

std::array< double, 3 >	m_positionInitIP
	Initial vertex postition. More...

double	m_energy
	Nominal beam energy. More...

double	m_parametrizationPosition
	Position for which the parameterisation was performed. More...

std::vector< std::unique_ptr< Equation > >	m_equation
	Vector of equations: 0 - x, 1 - y, 2 - sx, 3 - sy. More...

bool	m_yIsUsed
	True if y measurements are used in reconstruction. More...

Static Private Attributes
static constexpr int	m_numberOfPolynomials = 8
	Number of polynomials in equation. More...

static constexpr int	m_numberOfEquations = 4
	Number of equations - each for one measured value. More...

Detailed Description

Class representing transport parameterisation as four Equation.

Definition at line 23 of file Parameterization.h.

Constructor & Destructor Documentation

◆ Parameterization()

AFP::Parameterization::Parameterization ( std::string_view filename )

Takes filename of parametrization file as argument.

Loads parametrization file.

Definition at line 9 of file Parameterization.cxx.

                                                             : m_slopeInitIP{}, m_positionInitIP{} {
     std::cout << "Reading parameterization from file: \"" << filename << "\"" << '\n';
  
     std::ifstream file;
     file.open(filename.data());
     if(!file.is_open()) {
       std::stringstream errorMsg;
       errorMsg << "Parameterization::load(): \tCould not open parameterization file: \"" << filename << "\"";
       throw std::runtime_error(errorMsg.str().data());
     }
  
     while (file.peek() == '#') {
       char ch;
       std::string name;
       file >> ch >> name;
  
       if (name == "m_xPositionInitIP")  { double value; file >> value; std::cout << '\t' << name << ": " << value << '\n'; m_positionInitIP[0] = value; }
       if (name == "m_yPositionInitIP")  { double value; file >> value; std::cout << '\t' << name << ": " << value << '\n'; m_positionInitIP[1] = value; }
       if (name == "m_zPositionInitIP")  { double value; file >> value; std::cout << '\t' << name << ": " << value << '\n'; m_positionInitIP[2] = value; }
       if (name == "m_xSlopeInitIP")     { double value; file >> value; std::cout << '\t' << name << ": " << value << '\n'; m_slopeInitIP[0] = value; }
       if (name == "m_ySlopeInitIP")     { double value; file >> value; std::cout << '\t' << name << ": " << value << '\n'; m_slopeInitIP[1] = value; }
       if (name == "m_detectorPosition") { double value; file >> value; std::cout << '\t' << name << ": " << value << '\n'; m_parametrizationPosition = value; }
       if (name == "m_nominalEnergy")    { double value; file >> value; std::cout << '\t' << name << ": " << value << '\n'; m_energy = value; }
  
       file.ignore(1000, '\n');
     }
  
     std::vector<std::unique_ptr<TF1>> poly;
  
     m_yIsUsed = false;
  
     for (int iEqu = 0; iEqu < m_numberOfEquations; iEqu++) {
       file.ignore(1000, '\n');
       std::vector<double> coef;
  
       for (int iPoly = 0; iPoly < m_numberOfPolynomials; iPoly++) {
         int N = 0;
         file >> N;
         double value;
  
         for (int i = 0; i <= N; i++) {
           if (file.eof()) {
             std::stringstream errorMsg;
             errorMsg << "Parameterization::load(): \tFile \"" << filename
                      << "\" does not contain a correct parameterization";
             throw std::runtime_error(errorMsg.str().data());
           }
  
           file >> value;
           coef.push_back(value);
         }
  
         file.ignore(1000, '\n');
  
         TString name;
         name.Form("pol%d", N);
         auto temp = std::make_unique<TF1>(name, name);
         for (int iCoef = 0; iCoef <= N; ++iCoef) {
           temp->SetParameter(iCoef, coef.at(iCoef));
         }
  
         // For y-dependent equation check if polynomials are non-zero order
         if (iEqu == 1 or iEqu == 3) {
           m_yIsUsed = m_yIsUsed or not (N == 0 and coef[0] == 0);
         }
  
         poly.push_back(std::move(temp));
         coef.clear();
       }
  
       m_equation.push_back(std::make_unique<Equation>(m_energy, std::move(poly)));
       poly.clear();
     }
  
     std::cout << "Finished reading parameterization: \"" << filename << "\"\n\n";
   }

Member Function Documentation

◆ energy()

double AFP::Parameterization::energy ( ) const

inline

Returns nominal energy for which parametrization was performed.

Definition at line 66 of file Parameterization.h.

66 { return m_energy; }

◆ getEquation()

const Equation* AFP::Parameterization::getEquation ( int iEqu ) const

inline

Returns const raw pointer to equation of number iEqu.

Definition at line 69 of file Parameterization.h.

69 { return m_equation.at(iEqu).get(); }

◆ parametrizationPosition()

double AFP::Parameterization::parametrizationPosition ( ) const

inline

Returns position for which parametrization was performed.

Definition at line 63 of file Parameterization.h.

63 { return m_parametrizationPosition; }

◆ sx()

double AFP::Parameterization::sx	(	double	x0,
		double	y0,
		double	z0,
		double	sx0,
		double	sy0,
		double	E
	)		const

Evaluates value of horizontal slope equation.

Definition at line 94 of file Parameterization.cxx.

                                                                                                      {
     return getEquation(2)->evaluate(x0 + m_positionInitIP.at(0), y0 + m_positionInitIP.at(1), z0 + m_positionInitIP.at(2), sx0 + m_slopeInitIP.at(0), sy0 + m_slopeInitIP.at(1), E);
   }

◆ sy()

double AFP::Parameterization::sy	(	double	x0,
		double	y0,
		double	z0,
		double	sx0,
		double	sy0,
		double	E
	)		const

Evaluates value of vertical slope equation.

Definition at line 98 of file Parameterization.cxx.

                                                                                                      {
     return getEquation(3)->evaluate(x0 + m_positionInitIP.at(0), y0 + m_positionInitIP.at(1), z0 + m_positionInitIP.at(2), sx0 + m_slopeInitIP.at(0), sy0 + m_slopeInitIP.at(1), E);
   }

◆ x()

double AFP::Parameterization::x	(	double	x0,
		double	y0,
		double	z0,
		double	sx0,
		double	sy0,
		double	E
	)		const

Evaluates value of horizontal position equation.

Definition at line 86 of file Parameterization.cxx.

                                                                                                     {
     return getEquation(0)->evaluate(x0 + m_positionInitIP.at(0), y0 + m_positionInitIP.at(1), z0 + m_positionInitIP.at(2), sx0 + m_slopeInitIP.at(0), sy0 + m_slopeInitIP.at(1), E);
   }

◆ y()

double AFP::Parameterization::y	(	double	x0,
		double	y0,
		double	z0,
		double	sx0,
		double	sy0,
		double	E
	)		const

Evaluates value of vertical position equation.

Definition at line 90 of file Parameterization.cxx.

                                                                                                     {
     return getEquation(1)->evaluate(x0 + m_positionInitIP.at(0), y0 + m_positionInitIP.at(1), z0 + m_positionInitIP.at(2), sx0 + m_slopeInitIP.at(0), sy0 + m_slopeInitIP.at(1), E);
   }

◆ yIsUsed()

bool AFP::Parameterization::yIsUsed ( ) const

inline

Returns m_yIsUsed.

Definition at line 72 of file Parameterization.h.

72 { return m_yIsUsed; }

Member Data Documentation

◆ m_energy

double AFP::Parameterization::m_energy

private

Nominal beam energy.

Definition at line 88 of file Parameterization.h.

◆ m_equation

std::vector<std::unique_ptr<Equation> > AFP::Parameterization::m_equation

private

Vector of equations: 0 - x, 1 - y, 2 - sx, 3 - sy.

Definition at line 94 of file Parameterization.h.

◆ m_numberOfEquations

constexpr int AFP::Parameterization::m_numberOfEquations = 4

staticconstexprprivate

Number of equations - each for one measured value.

Definition at line 79 of file Parameterization.h.

◆ m_numberOfPolynomials

constexpr int AFP::Parameterization::m_numberOfPolynomials = 8

staticconstexprprivate

Number of polynomials in equation.

Definition at line 76 of file Parameterization.h.

◆ m_parametrizationPosition

double AFP::Parameterization::m_parametrizationPosition

private

Position for which the parameterisation was performed.

Definition at line 91 of file Parameterization.h.

◆ m_positionInitIP

std::array<double, 3> AFP::Parameterization::m_positionInitIP

private

Initial vertex postition.

Definition at line 85 of file Parameterization.h.

◆ m_slopeInitIP

std::array<double, 2> AFP::Parameterization::m_slopeInitIP

private

Crossing angles in rad.

Definition at line 82 of file Parameterization.h.

◆ m_yIsUsed

bool AFP::Parameterization::m_yIsUsed

private

True if y measurements are used in reconstruction.

Definition at line 97 of file Parameterization.h.

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

Classes

Public Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Parameterization()

Member Function Documentation

◆ energy()

◆ getEquation()

◆ parametrizationPosition()

◆ sx()

◆ sy()

◆ x()

◆ y()

◆ yIsUsed()

Member Data Documentation

◆ m_energy

◆ m_equation

◆ m_numberOfEquations

◆ m_numberOfPolynomials

◆ m_parametrizationPosition

◆ m_positionInitIP

◆ m_slopeInitIP

◆ m_yIsUsed