{
	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(1,3);

	c1->cd(1);
	gPad->SetLogy();
	t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.P()-(ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPy**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPz**2)**0.5");
	
	c1->cd(2);
	gPad->SetLogy();
	//t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.Px()-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 ");
	t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.P()-(ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPy**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPz**2)**0.5","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0.fEdep > 0.001 && ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPDG==1000020030");

	c1->cd(3);
	gPad->SetLogy();
	//t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.Px()-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 ");
	t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.P()-(ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPy**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPz**2)**0.5","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0.fEdep > 0.001 && Sum$(ERQTelescopeCsIPoint_T1_CsI_1_0.fPx_out)==0 && ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPDG==1000020030");

}