// Copyright 2019-2020 CERN and copyright holders of ALICE O2.

// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.

// All rights not expressly granted are reserved.

//

// This software is distributed under the terms of the GNU General Public

// License v3 (GPL Version 3), copied verbatim in the file "COPYING".

//

// In applying this license CERN does not waive the privileges and immunities

// granted to it by virtue of its status as an Intergovernmental Organization

// or submit itself to any jurisdiction.

/// @file Mille.h

/// @author Gero Flucke, ruben.shahoyan@cern.ch, michael.lettrich@cern.ch

/// @since October 2006

/// @brief Class to write C binary file.

/**

* Class to write a C binary (cf. below) file of a given name and to fill it

* with information used as input to **pede**.

* Use its member functions \c mille(), \c special(), \c kill() and \c end()

* as you would use the fortran \ref mille.f90 "MILLE"

* and its entry points \c MILLSP, \c KILLE and \c ENDLE.

*

* For debugging purposes constructor flags enable switching to text output and/or

* to write also derivatives and labels which are ==0.

* But note that **pede** will not be able to read text output and has not been tested with

* derivatives/labels ==0.

*/

#ifndef MILLE_H

#define MILLE_H

#include <fstream>

#include <TArrayI.h>

#include <TArrayF.h>

namespace o2

{

namespace align

{

class Mille

{

public:

Mille(const char* outFileName, bool asBinary = true, bool writeZero = false);

~Mille();

void mille(int NLC, const float* derLc, int NGL, const float* derGl,

const int* label, float rMeas, float sigma);

void special(int nSpecial, const float* floatings, const int* integers);

void kill();

int end();

private:

void newSet();

bool checkBufferSize(int nLocal, int nGlobal);

std::ofstream myOutFile; ///< C-binary for output

bool myAsBinary; ///< if false output as text

bool myWriteZero; ///< if true also write out derivatives/labels ==0

/// buffer size for ints and floats

int myBufferSize; ///< buffer size for ints and floats

TArrayI myBufferInt; ///< to collect labels etc.

TArrayF myBufferFloat; ///< to collect derivatives etc.

int myBufferPos; ///< position in buffer

bool myHasSpecial; ///< if true, special(..) already called for this record

/// largest label allowed: 2^31 - 1

enum { myMaxLabel = (0xFFFFFFFF - (1 << 31)) };

};

} // namespace align

} // namespace o2

#endif