{
	TFile* f = new TFile("sim_digi.root");
   	TTree* t = (TTree*)f->Get("er");
   	t->AddFriend("er","reco.root");

	TCanvas *c1 = new TCanvas("c1", "Telescope Coordinate reconstruction", 1400, 1000);
	c1->Divide(3,2);

	c1->cd(1);
	gPad->SetLogy();
	TH1F *T1X = new TH1F("T1X","T1 Xreco-Xsim",100,-1,1);
	T1X->GetXaxis()->SetTitle("Xreco-Xsim[cm]");
	//T1deE->GetYaxis()->SetTitle("de");
	t->Draw("ERQTelescopeTrack_T1_DoubleSi_SD2_XY_0.fTelescopeGlobalX-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fX>>T1X","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() < 2","gof");
	T1X->Draw();

	c1->cd(2);
	gPad->SetLogy();
	TH1F *T1Y = new TH1F("T1Y","T1 Yreco-Ysim",100,-1,1);
	T1Y->GetXaxis()->SetTitle("Yreco-Ysim[cm]");
	//T1deE->GetYaxis()->SetTitle("de");
	t->Draw("ERQTelescopeTrack_T1_DoubleSi_SD2_XY_0.fTelescopeGlobalY-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fY>>T1Y","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() < 2","gof");
	T1Y->Draw();

	c1->cd(3);
	gPad->SetLogy();
	TH1F *T1Z = new TH1F("T1Z","T1 Zreco-Zsim",100,-1,1);
	T1Z->GetXaxis()->SetTitle("Zreco-Zsim[cm]");
	//T1deE->GetYaxis()->SetTitle("de");
	t->Draw("ERQTelescopeTrack_T1_DoubleSi_SD2_XY_0.fTelescopeGlobalZ-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fZ>>T1Z","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() < 2","gof");
	T1Z->Draw();

	c1->cd(4);
	gPad->SetLogy();
	TH1F *T2X = new TH1F("T2X","T2 Xreco-Xsim",100,-1,1);
	T2X->GetXaxis()->SetTitle("Xreco-Xsim[cm]");
	//T1deE->GetYaxis()->SetTitle("de");
	t->Draw("ERQTelescopeTrack_T2_DoubleSi_SD2_XY_1.fTelescopeGlobalX-ERQTelescopeSiPoint_T2_DoubleSi_SD2_XY_1_Y.fX>>T2X","ERQTelescopeSiDigi_T2_DoubleSi_SD2_XY_1_Y@.GetEntries() < 2","gof");
	T2X->Draw();

	c1->cd(5);
	gPad->SetLogy();
	TH1F *T2Y = new TH1F("T2Y","T2 Yreco-Ysim",100,-1,1);
	T2Y->GetXaxis()->SetTitle("Yreco-Ysim[cm]");
	//T1deE->GetYaxis()->SetTitle("de");
	t->Draw("ERQTelescopeTrack_T2_DoubleSi_SD2_XY_1.fTelescopeGlobalY-ERQTelescopeSiPoint_T2_DoubleSi_SD2_XY_1_Y.fY>>T2Y","ERQTelescopeSiDigi_T2_DoubleSi_SD2_XY_1_Y@.GetEntries() < 2","gof");
	T2Y->Draw();

	c1->cd(6);
	gPad->SetLogy();
	TH1F *T2Z = new TH1F("T2Z","T2 Zreco-Zsim",100,-1,1);
	T2Z->GetXaxis()->SetTitle("Zreco-Zsim[cm]");
	//T1deE->GetYaxis()->SetTitle("de");
	t->Draw("ERQTelescopeTrack_T2_DoubleSi_SD2_XY_1.fTelescopeGlobalZ-ERQTelescopeSiPoint_T2_DoubleSi_SD2_XY_1_Y.fZ>>T2Z","ERQTelescopeSiDigi_T2_DoubleSi_SD2_XY_1_Y@.GetEntries() < 2","gof");
	T2Z->Draw();



}