/*
 * PndLmdAlignQA.h
 *
 *  Created on: Jul 6, 2015
 *      Author: Roman Klasen, roklasen@uni-mainz.de or klasen@kph.uni-mainz.de
 *
 *      helper class that plots stuff, use one instance for desired plot, then delete and create new.
 */

#ifndef LMD_LMDSENSORALIGNMENT_PNDLMDALIGNQA_H_
#define LMD_LMDSENSORALIGNMENT_PNDLMDALIGNQA_H_

#include <sstream>
#include <string>
#include <vector>

#include <PndLmdAlignManager.h>
#include <matrix.h>

struct histParams{
	bool printCMPXinPathName=true;
	double scaleFactor=1.0;
	std::string fileName, title, xtitle, ytitle, path;
	int bins=30;
	double xMin=-1;
	double xMax=-1;
	int vectorIndex=0;
};

class PndLmdAlignQA {

private:

	std::vector<std::string> _inputFiles;
	std::string outputPath, pdfOutPath, binaryMatPath, LMDMatPath;
	double infoMomentum;
	bool infoAbsolute, infoRelative, byPlane, _inCentimeters, _enableHelperMatrix;
	bool alignOptionBool;
	int pairsRequired;
	PndLmdAlignManager manager;

	//contains number of pairs on overlap area
	std::map<int, int> matrixInfo;

	void createHist(std::vector<std::vector<double> > &vec, histParams &parameters);

	void histPixelDistances(int sen1, int sen2, bool aligned=true);
	double calculateOverlappingArea(int id1, int id2, bool aligned=true);

	int noOfPairs(int id1, int id2);

	PndLmdDim *dimension;

public:

	PndLmdAlignQA();
	virtual ~PndLmdAlignQA();
	void init();


	// ========== QA functions

	void checkCyclicMatrices(bool inCentimeters=true);
	void checkCombinedMatrices(bool inCentimeters=true);
	void plotCMvsPXmatrices();
	void plotMatrixresiduals(bool inCentimeters=true);
	void plotPXvsCMmatricesResiduals();
	void histogramPairDistances();

	// ========== end of QA


	// ========== helper functions

	void calculatePixelDistancesFrontToBack();
	void readMatrixInfo();
	void checkIOpaths();
	bool checkForMatrixFiles();
	void calculateOverlapingAreas();
	void compareMatrices(int param);

	// ========== end of helper functions



	// ========== settings for this QA task

	void setInCentimeters(bool inCentimeters) {	this->_inCentimeters = inCentimeters; }
	void setCorrection(bool correction) { _enableHelperMatrix = correction;	}
	void setLmdMatPath(const std::string& path) { LMDMatPath = path; }
	void setBinaryMatPath(const std::string& path) { binaryMatPath = path; }
	void setPdfOutPath(const std::string& path) { pdfOutPath = path; }
	void setAlignedGeometry(bool aligned) {	alignOptionBool = aligned; }
	void setPairsRequired(int number) {	pairsRequired = number;	}
	void setOutputPath(std::string path){ outputPath = path; }
	void setInfoAbsolute(bool info) { infoAbsolute = info; }
	void setInfoMomentum(double info) { infoMomentum = info; }
	void setInfoRelative(bool info) { infoRelative = info; }

	// ========== end of settings
};

#endif /* LMD_LMDSENSORALIGNMENT_PNDLMDALIGNQA_H_ */