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

#include <MilleBinary.h>

Public Member Functions
void	addData (double aMeas, double aErr, unsigned int numLocal, unsigned int indLocal, 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 68 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 41 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 58 of file MilleBinary.cc.

References binaryFile.

                             {
     binaryFile.close();
   }

Member Function Documentation

void gbl::MilleBinary::addData	(	double	aMeas,
		double	aErr,
		unsigned int	numLocal,
		unsigned int *	indLocal,
		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]	numLocal	Number of local labels/derivatives
[in]	indLocal	Array of labels of local parameters
[in]	derLocal	Array of derivatives for local parameters
[in]	labGlobal	List of labels of global parameters
[in]	derGlobal	List of derivatives for global parameters

Definition at line 72 of file MilleBinary.cc.

References doubleBuffer, doublePrecision, floatBuffer, mps_fire::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 < numLocal; ++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 < numLocal; ++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 113 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);
   }