Millepede-II (binary) record. More...

#include <MilleBinary.h>

Public Member Functions
void	addData (double aMeas, double aPrec, const std::vector< unsigned int > &indLocal, const std::vector< double > &derLocal, const std::vector< int > &labGlobal, const std::vector< double > &derGlobal)
	Add data block to (end of) record. More...

	MilleBinary (const std::string fileName="milleBinaryISN.dat", bool doublePrec=false, unsigned int aSize=2000)
	Create binary file. More...

void	writeRecord ()
	Write record to file. More...

virtual	~MilleBinary ()

Private Attributes
std::ofstream	binaryFile
	Binary File. More...

std::vector< double >	doubleBuffer
	Double buffer. More...

bool	doublePrecision
	Flag for storage in as double values. More...

std::vector< float >	floatBuffer
	Float buffer. More...

std::vector< int >	intBuffer
	Integer buffer. More...

Detailed Description

Millepede-II (binary) record.

Containing information for local (track) and global fit.

The data blocks are collected in two arrays, a real array (containing float or double values) and integer array, of same length. A positive record length indicate float and a negative one double values. The content of the record is:

*         real array              integer array
*     0   0.0                     error count (this record)
*     1   RMEAS, measured value   0                            -+
*     2   local derivative        index of local derivative     |
*     3   local derivative        index of local derivative     |
*     4    ...                                                  | block
*         SIGMA, error (>0)       0                             |
*         global derivative       label of global derivative    |
*         global derivative       label of global derivative   -+
*         RMEAS, measured value   0
*         local derivative        index of local derivative
*         local derivative        index of local derivative
*         ...
*         SIGMA, error            0
*         global derivative       label of global derivative
*         global derivative       label of global derivative
*         ...
*         global derivative       label of global derivative
*

Definition at line 46 of file MilleBinary.h.

Constructor & Destructor Documentation

gbl::MilleBinary::MilleBinary	(	const std::string	fileName = `"milleBinaryISN.dat"`,
		bool	doublePrec = `false`,
		unsigned int	aSize = `2000`
	)

Create binary file.

Parameters

[in]	fileName	File name
[in]	doublePrec	Flag for storage as double values
[in]	aSize	Buffer size

Definition at line 19 of file MilleBinary.cc.

References doubleBuffer, doublePrecision, floatBuffer, and intBuffer.

                             :
         binaryFile(fileName.c_str(), std::ios::binary | std::ios::out), intBuffer(), floatBuffer(), doubleBuffer(), doublePrecision(
                 doublePrec) {
     intBuffer.reserve(aSize);
     intBuffer.push_back(0); // first word is error counter
     if (doublePrecision) {
         doubleBuffer.reserve(aSize);
         doubleBuffer.push_back(0.);
 
     } else {
         floatBuffer.reserve(aSize);
         floatBuffer.push_back(0.);
     }
 }

gbl::MilleBinary::~MilleBinary ( )

virtual

Definition at line 35 of file MilleBinary.cc.

References binaryFile.

                           {
     binaryFile.close();
 }

Member Function Documentation

void gbl::MilleBinary::addData	(	double	aMeas,
		double	aErr,
		const std::vector< unsigned int > &	indLocal,
		const std::vector< double > &	derLocal,
		const std::vector< int > &	labGlobal,
		const std::vector< double > &	derGlobal
	)

Add data block to (end of) record.

Parameters

[in]	aMeas	Value
[in]	aErr	Error
[in]	indLocal	List of labels of local parameters
[in]	derLocal	List of derivatives for local parameters
[in]	labGlobal	List of labels of global parameters
[in]	derGlobal	List of derivatives for global parameters

Definition at line 48 of file MilleBinary.cc.

References doubleBuffer, doublePrecision, floatBuffer, i, and intBuffer.

Referenced by gbl::GblTrajectory::milleOut().

                                             {
 
     if (doublePrecision) {
         // double values
         intBuffer.push_back(0);
         doubleBuffer.push_back(aMeas);
         for (unsigned int i = 0; i < indLocal.size(); ++i) {
             intBuffer.push_back(indLocal[i]);
             doubleBuffer.push_back(derLocal[i]);
         }
         intBuffer.push_back(0);
         doubleBuffer.push_back(aErr);
         for (unsigned int i = 0; i < labGlobal.size(); ++i) {
             if (derGlobal[i]) {
                 intBuffer.push_back(labGlobal[i]);
                 doubleBuffer.push_back(derGlobal[i]);
             }
         }
     } else {
         // float values
         intBuffer.push_back(0);
         floatBuffer.push_back(aMeas);
         for (unsigned int i = 0; i < indLocal.size(); ++i) {
             intBuffer.push_back(indLocal[i]);
             floatBuffer.push_back(derLocal[i]);
         }
         intBuffer.push_back(0);
         floatBuffer.push_back(aErr);
         for (unsigned int i = 0; i < labGlobal.size(); ++i) {
             if (derGlobal[i]) {
                 intBuffer.push_back(labGlobal[i]);
                 floatBuffer.push_back(derGlobal[i]);
             }
         }
     }
 }

void gbl::MilleBinary::writeRecord ( )

Write record to file.

Definition at line 89 of file MilleBinary.cc.

References binaryFile, doubleBuffer, doublePrecision, floatBuffer, and intBuffer.

Referenced by gbl::GblTrajectory::milleOut().

                               {
 
     const int recordLength =
             (doublePrecision) ? -intBuffer.size() * 2 : intBuffer.size() * 2;
     binaryFile.write(reinterpret_cast<const char*>(&recordLength),
             sizeof(recordLength));
     if (doublePrecision)
         binaryFile.write(reinterpret_cast<char*>(&doubleBuffer[0]),
                 doubleBuffer.size() * sizeof(doubleBuffer[0]));
     else
         binaryFile.write(reinterpret_cast<char*>(&floatBuffer[0]),
                 floatBuffer.size() * sizeof(floatBuffer[0]));
     binaryFile.write(reinterpret_cast<char*>(&intBuffer[0]),
             intBuffer.size() * sizeof(intBuffer[0]));
 // start with new record
     intBuffer.resize(1);
     if (doublePrecision)
         doubleBuffer.resize(1);
     else
         floatBuffer.resize(1);
 }