#include <TString.h>
#include <TSystem.h>
#include "ROOT/RVec.hxx"
#include "ROOT/RDataFrame.hxx"
#include <string>
#include "TTreeReaderValue.h"
void plot_multi()
{
	ROOT::EnableImplicitMT(4);
	ROOT::RDataFrame d("tree", "./data.root");

	auto Si = d.Filter("std::string s(ParticleName.begin(), ParticleName.end()); return s == \"proton\";").Histo1D({"Si","Si",4000,0,400},"particleEnergy");
	
	TCanvas *c1 = new TCanvas;
	
	c1->SetTopMargin(0.1);
	c1->SetBottomMargin(0.13);
	c1->SetLogy();
	Si->SetLineColor(2);
	Si->SetLineWidth(3);
	Si->SetStats(0);

	Si->SetTitle(" ");
	Si->GetXaxis()->SetTitle("Proton Energy [MeV]");
	Si->GetYaxis()->SetTitle("counts");

	Si->GetYaxis()->SetLabelFont(12);
	Si->GetYaxis()->SetLabelSize(0.05);
	Si->GetYaxis()->SetTitleSize(0.05);
	Si->GetYaxis()->SetTitleOffset(0.8);
	Si->GetYaxis()->SetTitleFont(22);
//	Si->GetYaxis()->SetRangeUser(1,1e7);

	Si->GetXaxis()->SetLabelFont(12);
	Si->GetXaxis()->SetLabelSize(0.05);
	Si->GetXaxis()->SetTitleSize(0.05);
	Si->GetXaxis()->SetTitleOffset(1.2);
	Si->GetXaxis()->SetTitleFont(22);

	Si->GetXaxis()->CenterTitle();
	Si->GetYaxis()->CenterTitle();
	Si->DrawCopy();
}