#include "TLine.h"
#include "TPad.h"
#include "TCutG.h"

void showT1Sim()
{

	TFile *fc = new TFile("cut.root", "READ");
	TCutG *cHe = (TCutG*)fc->Get("cut3He");

	TFile *fcPos = new TFile("cutsPos.root", "READ");
	TCutG *c3He = (TCutG*)fcPos->Get("c3HePos");
	TCutG *c6He = (TCutG*)fcPos->Get("c6HePos");

	TFile *fr = new TFile("sim_digi.root", "READ");
	TTree *tr = (TTree*)fr->Get("er");



	TCanvas *c3 = new TCanvas("c3", "Telescope info", 1200, 800);
	c3->Divide(3,2);

	c3->cd(1);
//	tr->Draw("(MCTrack.fEnergy - MCTrack.fMass)*1000", "MCTrack.fPdgCode == 1000020060");
	tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",
			"ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 && ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70","");
	cHe->Draw("same");

	c3->cd(2);
//	tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",
//			"ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 && ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70","");
//	cHe->Draw("same");
	tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "cut3He", "col");

	c3->cd(3);
	tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "", "col");
	c3He->SetLineWidth(3);
	c3He->SetLineColor(kRed);
	c3He->Draw("same");
	c6He->Draw("same");
	c6He->SetLineWidth(3);
	c6He->SetLineColor(kGreen);


	c3->cd(4);
	tr->Draw("(MCTrack.fEnergy - MCTrack.fMass)*1000", "MCTrack.fPdgCode == 1000010010");
	tr->Draw("(MCTrack.fEnergy - MCTrack.fMass)*1000", "MCTrack.fPdgCode == 2212");
	tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",
			"ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 && ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70","");
	cHe->Draw("same");
	tr->SetMarkerColor(kRed);
	tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",
			"ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 "
			"&& ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70 "
			"&& cut3He","same");


	c3->cd(5);
	tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000","","col");


	c3->cd(6);
	tr->Draw("ERQTelescopeSiPoint_DoubleSi_SD2_XY_1_X.fY:ERQTelescopeSiPoint_DoubleSi_SD2_XY_1_X.fX", "", "col");

	return;

	TCanvas *c1 = new TCanvas("c1", "Telescope info", 800, 800);
//	c1->Divide(3,2);

	std::cout << tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "", "col") << std::endl;
	c3He->Draw("same");
//	std::cout << tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "c3HePos", "col") << std::endl;
//	std::cout << tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "c6HePos", "col") << std::endl;

}