ROOT examples

make_binned_and_fitted_to_root.C

#include "TFile.h"
#include "TGraphErrors.h"
#include "TCanvas.h"
#include "TAxis.h"
//g++ `root-config --libs --cflags --ldflags` make_graph.C -o make_graph
//.x ./make_binned_and_fitted_to_root.C("binned_and_fitted_channel_signal_1_channel_trigger_2_pulse_length_3.csv", "graphs.root", "title", "x_name", "y_name")
#include "Riostream.h"
using namespace std;

void make_binned_and_fitted_to_root(TString csv_file, TString root_file, TString Graph_Name="", TString x_name="", TString y_name=""){
	ifstream in;
	in.open(csv_file);
	TFile *f = new TFile(root_file, "UPDATE");
	TGraphErrors *g1 = new TGraphErrors();
	TCanvas *c1 = new TCanvas(Graph_Name);
	Double_t dac, A, mu, sigma;
	Int_t point = 0;

	while(1){
		in>>dac>>A>>mu>>sigma;
		point++;
		if (!in.good()){
			break;
		}
		g1->Set(point);
		g1->SetPoint(point, dac, mu);
		g1->SetPointError(point, 50, sigma);
	}
	c1->cd();
	g1->SetTitle(Graph_Name);
	g1->GetXaxis()->SetTitle(x_name);
	g1->GetYaxis()->SetTitle(y_name);
	g1->Draw("ap");
	f->cd();
	c1->Write();
	f->Close();
	delete f;
	delete c1;
	delete g1;
}

int main(int argc, char **argv){
	switch(argc){
		case 3:
			make_binned_and_fitted_to_root(argv[1], argv[2]);
			break;
		case 4: 
			make_binned_and_fitted_to_root(argv[1], argv[2], argv[3]);
			break;
		case 5:
			make_binned_and_fitted_to_root(argv[1], argv[2], argv[3],argv[4]);
			break;
		case 6:
			make_binned_and_fitted_to_root(argv[1], argv[2], argv[3],argv[4],argv[5]);
			break;
		default:
			return 1;
	}
	return 0;
}

make_graph.C

#include "TFile.h"
#include "TGraph.h"
#include "TCanvas.h"
#include "TAxis.h"
//g++ `root-config --libs --cflags --ldflags` make_graph.C -o make_graph
#include "Riostream.h"
using namespace std;

void make_graph(TString csv_file, TString root_file, TString Graph_Name="", TString x_name="", TString y_name=""){
	ifstream in;
	in.open(csv_file);
	TFile *f = new TFile(root_file, "UPDATE");
	TGraph *g1 = new TGraph();
	TCanvas *c1 = new TCanvas(Graph_Name);
	Double_t x, y;
	Int_t point = 0;
	
	while(1){
		in>>x>>y;
		point++;
		if (!in.good()){
			break;
		}
		g1->Set(point);
		g1->SetPoint(point, x, y);
	}
	c1->cd();
	g1->SetTitle(Graph_Name);
	g1->GetXaxis()->SetTitle(x_name);
	g1->GetYaxis()->SetTitle(y_name);
	g1->Draw("ap");
	f->cd();
	c1->Write();
	f->Close();
}

int main(int argc, char **argv){
	switch(argc){
		case 3:
			make_graph(argv[1], argv[2]);
			break;
		case 4: 
			make_graph(argv[1], argv[2], argv[3]);
			break;
		case 5:
			make_graph(argv[1], argv[2], argv[3],argv[4]);
			break;
		case 6:
			make_graph(argv[1], argv[2], argv[3],argv[4],argv[5]);
			break;
		default:
			return 1;
	}
	return 0;
}