Fit-AB.cpp

// this file is distributed under 
// MIT license
#include <iostream>
#include <sstream>
#include <thread>
#include <map>
#include <gnuplot_wrap.h>
#include <IO/gethist.h>
#include <JPetLargeBarrelExtensions/PetDict.h>
#include <JPetLargeBarrelExtensions/BarrelExtensions.h>
#include <Calc/SyncProcedures.h>
using namespace std;
using namespace GnuplotWrap;
using namespace MathTemplates;
int main(int argc, char **argv) {
	if(argc<2){
		cerr<<"Usage: "<<argv[0]<<" <filename> <filename> ..."<<endl;
		return -1;
	}
	vector<string> root_filenames;
	for(int i=1;i<argc;i++)
		root_filenames.push_back(string(argv[i]));
	Plotter::Instance().SetOutput(".","AB");
	auto map=make_JPetMap<SyncAB_results>();
	for(size_t layer=1;layer<=map->layersCount();layer++){
		hist<double> position,height,chisq;
		for(size_t slot=1;slot<=map->layerSize(layer);slot++){
			const auto name=LayerSlot(layer,slot);
			const auto shist=ReadHist(root_filenames,name);
			auto&item=map->item({.layer=layer,.slot=slot});
			item=Sync::Fit4SyncAB(shist,"SyncAB "+name,1);
			if(item.valid()){
			    height<<point<value<double>>(slot,item.height);
			    position<<point<value<double>>(slot,item.peak);
			    chisq<<point<value<double>>(slot,item.chi_sq);
			}
		}
		Plot<double>().Hist(height,"Peak height")
		<<"set key on"<<"set title 'Layer "+to_string(layer)+"'";
		Plot<double>().Hist(position,"Position")
		<<"set key on"<<"set title 'Layer "+to_string(layer)+"'";
		Plot<double>().Hist(chisq,"Chi^2")
		<<"set key on"<<"set title 'Layer "+to_string(layer)+"'";
	}
	cout<<(*map);
	return 0;
}