Class to process Pede result file(s) More...

#include <PedeResult.h>

Classes
struct	parameterData
	Struct to hold data for a parameter. More...

Public Member Functions
	PedeResult ()
	Default constructor, use read(...) to init the object.

	PedeResult (const std::string &filename)
	Constructor which loads given file.

void	read (std::string filename="millepede.res")
	Reads the result file and inits the object.

void	readEigenFile (std::string filename="millepede.eve")
	Reads file with eigen-vector/numbers.

bool	isValid ()
	Was the object initialized properly from a result file?

void	dump ()
	Dump the content to std::cout.

int	getNoParameters () const
	Get number of parameters in result (for looping over)

int	getNoDeterminedParameters () const
	Get number of determined parameters (with correction)

int	getParameterIndex (int parameterLabel)
	Get index of parameter with given label.

unsigned int	getParameterLabel (unsigned int parameterIndex)
	Get label of parameter at index.

double	getParameterCorrection (unsigned int parameterIndex)
	Get determined correction of parameter at index.

double	getParameterError (unsigned int parameterIndex)
	Get correction error of parameter at index.

double	getParameterPresigma (unsigned int parameterIndex)
	Get presigma of parameter at index.

bool	isParameterFixed (unsigned int parameterIndex)
	Is parameter at given index fixed?

bool	isParameterDetermined (unsigned int parameterIndex)
	Is parameter at given index determined?

int	getNoEigenPairs () const
	Get the number of eigenvectors(numbers) in eigen file.

double	getEigenNumber (unsigned int eigenPairIndex)
	Get eigennumber at given index.

double	getEigenVectorElement (unsigned int eigenPairIndex, unsigned int parameterIndex)
	Get eigenvector at given index.

Private Member Functions
bool	outOfBounds (int iparam)
	Is index out of valid range?

Private Attributes
std::vector< parameterData >	data
	Vector with all the parameter data.

std::vector< double >	eigenNumbers
	Vector of loaded eigennumbers.

bool	valid
	Flag to check if data wa loaded.

std::map< int, int >	labelIndices
	Map to link parameter labels and their indices in result.

int	noDeterminedParams {0}
	Number of parameters actually determined.

Detailed Description

Class to process Pede result file(s)

Definition at line 22 of file PedeResult.h.

Constructor & Destructor Documentation

◆ PedeResult() [1/2]

PedeResult ( )

inline

Default constructor, use read(...) to init the object.

Definition at line 25 of file PedeResult.h.

25: data(), eigenNumbers(), valid(false), labelIndices() {}

Belle2::alignment::PedeResult::valid

bool valid

Flag to check if data wa loaded.

Definition: PedeResult.h:105

Belle2::alignment::PedeResult::labelIndices

std::map< int, int > labelIndices

Map to link parameter labels and their indices in result.

Definition: PedeResult.h:107

Belle2::alignment::PedeResult::eigenNumbers

std::vector< double > eigenNumbers

Vector of loaded eigennumbers.

Definition: PedeResult.h:103

Belle2::alignment::PedeResult::data

std::vector< parameterData > data

Vector with all the parameter data.

Definition: PedeResult.h:101

◆ PedeResult() [2/2]

PedeResult ( const std::string & filename )

explicit

Constructor which loads given file.

Parameters

filename Name of the result file (millepede.res)

Definition at line 126 of file PedeResult.cc.

                                                : data(), eigenNumbers(), valid(false), labelIndices()
    {
      read(filename);
    }

Member Function Documentation

◆ dump()

void dump ( )

Dump the content to std::cout.

Definition at line 176 of file PedeResult.cc.

    {
      for (auto pdata : data) {
        std::cout << pdata.index << " : " << pdata.label << " corr: " << pdata.correction << " err: " << pdata.error << " presigma: " <<
                  pdata.presigma << endl;
      }
    }

◆ getEigenNumber()

double getEigenNumber ( unsigned int eigenPairIndex )

Get eigennumber at given index.

Definition at line 161 of file PedeResult.cc.

    {
      if (eigenPairIndex >= eigenNumbers.size()) return 0.;
      return eigenNumbers[eigenPairIndex];
    }

◆ getEigenVectorElement()

double getEigenVectorElement	(	unsigned int	eigenPairIndex,
		unsigned int	parameterIndex
	)

Get eigenvector at given index.

Definition at line 166 of file PedeResult.cc.

    {
      if (outOfBounds(parameterIndex))
        return 0.;
      parameterData& element = data[parameterIndex];
      if (eigenPairIndex >= element.eigenweights.size())
        return 0.;
 
      return data[parameterIndex].eigenweights[eigenPairIndex];
    }

◆ getNoDeterminedParameters()

int getNoDeterminedParameters ( ) const

inline

Get number of determined parameters (with correction)

Definition at line 40 of file PedeResult.h.

40{return noDeterminedParams;}

Belle2::alignment::PedeResult::noDeterminedParams

int noDeterminedParams

Number of parameters actually determined.

Definition: PedeResult.h:109

◆ getNoEigenPairs()

int getNoEigenPairs ( ) const

inline

Get the number of eigenvectors(numbers) in eigen file.

Definition at line 56 of file PedeResult.h.

56{return eigenNumbers.size();}

◆ getNoParameters()

int getNoParameters ( ) const

inline

Get number of parameters in result (for looping over)

Definition at line 38 of file PedeResult.h.

38{return data.size();}

◆ getParameterCorrection()

double getParameterCorrection ( unsigned int parameterIndex )

Get determined correction of parameter at index.

Definition at line 136 of file PedeResult.cc.

    {
      if (outOfBounds(parameterIndex)) return 0.;
      return data[parameterIndex].correction;
    }

◆ getParameterError()

double getParameterError ( unsigned int parameterIndex )

Get correction error of parameter at index.

Definition at line 141 of file PedeResult.cc.

    {
      if (outOfBounds(parameterIndex)) return 0.;
      return data[parameterIndex].error;
    }

◆ getParameterIndex()

int getParameterIndex ( int parameterLabel )

Get index of parameter with given label.

Definition at line 27 of file PedeResult.cc.

    {
      auto index = labelIndices.find(parameterLabel);
      if (index != labelIndices.end())
        return index->second;
      else
        return -1;
    }

◆ getParameterLabel()

unsigned int getParameterLabel ( unsigned int parameterIndex )

Get label of parameter at index.

Definition at line 131 of file PedeResult.cc.

    {
      if (outOfBounds(parameterIndex)) return 0;
      return data[parameterIndex].label;
    }

◆ getParameterPresigma()

double getParameterPresigma ( unsigned int parameterIndex )

Get presigma of parameter at index.

Definition at line 146 of file PedeResult.cc.

    {
      if (outOfBounds(parameterIndex)) return 0.;
      return data[parameterIndex].presigma;
    }

◆ isParameterDetermined()

bool isParameterDetermined ( unsigned int parameterIndex )

Is parameter at given index determined?

Definition at line 156 of file PedeResult.cc.

    {
      if (outOfBounds(parameterIndex)) return false;
      return !isParameterFixed(parameterIndex) && data[parameterIndex].correction != 0. && data[parameterIndex].error != 0.;
    }

◆ isParameterFixed()

bool isParameterFixed ( unsigned int parameterIndex )

Is parameter at given index fixed?

Definition at line 151 of file PedeResult.cc.

    {
      if (outOfBounds(parameterIndex)) return false;
      return data[parameterIndex].presigma < 0;
    }

◆ isValid()

bool isValid ( )

inline

Was the object initialized properly from a result file?

Definition at line 34 of file PedeResult.h.

34{return valid;}

◆ outOfBounds()

bool outOfBounds ( int iparam )

inlineprivate

Is index out of valid range?

Definition at line 64 of file PedeResult.h.

64{if (iparam < 0) return true; if ((unsigned int)iparam >= data.size()) return true; return false;}

◆ read()

void read ( std::string filename = "millepede.res" )

Reads the result file and inits the object.

Definition at line 95 of file PedeResult.cc.

    {
      string line;
      ifstream res(filename);
      // skip comment line
      getline(res, line);
 
      noDeterminedParams = 0;
      data.clear();
      labelIndices.clear();
      int index(-1);
 
      while (getline(res, line)) {
        ++index;
        int label = -1;
        double param = 0.;
        double presigma = 0.;
        double differ = 0.;
        double error = 0.;
        stringstream ss;
        ss << line;
        ss >> label >> param >> presigma >> differ >> error;
        data.push_back(parameterData(index, label, param, error, presigma));
        labelIndices.insert(make_pair(label, index));
        if (isParameterDetermined(index))
          ++noDeterminedParams;
      }
 
      if (data.size() > 0)
        valid = true;
    }

◆ readEigenFile()

void readEigenFile ( std::string filename = "millepede.eve" )

Reads file with eigen-vector/numbers.

Definition at line 35 of file PedeResult.cc.

    {
      if (!valid)
        return;
 
      string line;
      ifstream file(filename);
 
      if (!file.is_open())
        return;
 
      int eigenPairIndex = -1;
 
      for (auto& pData : data) {
        pData.eigenweights.clear();
        pData.eigenweights.reserve(20);
      }
 
      while (getline(file, line)) {
        cout << line << endl;
        stringstream ss;
        ss << line;
        string tmp1, tmp2, tmp3;
        int eigenPairNumber;
        double eigenvalue;
        ss >> tmp1 >> eigenPairNumber >> tmp2 >> tmp3 >> eigenvalue;
        cout << tmp1 << " " << tmp2 << " " << tmp3 << " " << eigenPairNumber << " " << eigenvalue << endl;
        if (tmp1 == "Eigenvector")
          ++eigenPairIndex;
        else
          continue;
 
        // Read elements
        vector<pair<int, double>> elements;
        while (getline(file, line) && line.length() > 2) {
          int lab1 = 0;
          int lab2 = 0;
          int lab3 = 0;
          double w1, w2, w3;
          stringstream triplet;
          triplet << line;
          triplet >> lab1 >> w1 >> lab2 >> w2 >> lab3 >> w3;
          if (lab1) elements.push_back({lab1, w1});
          if (lab2) elements.push_back({lab2, w2});
          if (lab3) elements.push_back({lab3, w3});
          //cout << lab1 << w1 << lab2 << w2 << lab3 << w3 << endl;
        }
        for (auto& pData : data) {
          pData.eigenweights.push_back(0.);
        }
        for (auto el : elements) {
          int index = getParameterIndex(el.first);
          if (index >= 0)
            data[index].eigenweights[eigenPairIndex] = el.second;
        }
        eigenNumbers.push_back(eigenvalue);
      }
 
 
    }

Member Data Documentation

◆ data

std::vector<parameterData> data

private

Vector with all the parameter data.

Definition at line 101 of file PedeResult.h.

◆ eigenNumbers

std::vector<double> eigenNumbers

private

Vector of loaded eigennumbers.

Definition at line 103 of file PedeResult.h.

◆ labelIndices

std::map<int, int> labelIndices

private

Map to link parameter labels and their indices in result.

Definition at line 107 of file PedeResult.h.

◆ noDeterminedParams

int noDeterminedParams {0}

private

Number of parameters actually determined.

Definition at line 109 of file PedeResult.h.

◆ valid

bool valid

private

Flag to check if data wa loaded.

Definition at line 105 of file PedeResult.h.

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

alignment/include/PedeResult.h
alignment/src/PedeResult.cc

Classes

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PedeResult() [1/2]

◆ PedeResult() [2/2]

Member Function Documentation

◆ dump()

◆ getEigenNumber()

◆ getEigenVectorElement()

◆ getNoDeterminedParameters()

◆ getNoEigenPairs()

◆ getNoParameters()

◆ getParameterCorrection()

◆ getParameterError()

◆ getParameterIndex()

◆ getParameterLabel()

◆ getParameterPresigma()

◆ isParameterDetermined()

◆ isParameterFixed()

◆ isValid()

◆ outOfBounds()

◆ read()

◆ readEigenFile()

Member Data Documentation

◆ data

◆ eigenNumbers

◆ labelIndices

◆ noDeterminedParams

◆ valid