{
	gROOT->Reset();
  
	delete gROOT->GetListOfFiles()->FindObject("200keV-nopol.root");
	delete gROOT->GetListOfFiles()->FindObject("200keV-x-tarpol.root");
	delete gROOT->GetListOfFiles()->FindObject("200keV-y-tarpol.root");
	delete gROOT->GetListOfFiles()->FindObject("200keV-z-tarpol.root");
	delete gROOT->GetListOfFiles()->FindObject("200keV-x-beamtarpol.root");
	delete gROOT->GetListOfFiles()->FindObject("200keV-y-beamtarpol.root");
	delete gROOT->GetListOfFiles()->FindObject("200keV-z-beamtarpol.root");

	TFile f = TFile("200keV-nopol.root");  

	TFile fx = TFile("200keV-x-tarpol.root");  
	TFile fx2 = TFile("200keV-x-beamtarpol.root");  

	TFile fy = TFile("200keV-y-tarpol.root");  
	TFile fy2 = TFile("200keV-y-beamtarpol.root");  

	TFile fz = TFile("200keV-z-tarpol.root");  
	TFile fz2 = TFile("200keV-z-beamtarpol.root");  





	TCanvas *c1 = new TCanvas("c1", " ");
	c1->Divide(4,2);
	
	// T(0,0,0)-B(0,0,0)

	c1->cd(1);
	gPad->SetLogy();

	//Gammas
	
	TH1D *dE_gam=(TH1D*)f.Get("7");
	dE_gam->SetLineColor(kMagenta);
	dE_gam->Draw("HIST");
	dE_gam->GetXaxis()->SetTitle("Energy(MeV)(gammas));
	dE_gam->GetYaxis()->SetTitle("Events");

	//Electrons
	
	TH1D *dE_ele=(TH1D*)f.Get("15");
	dE_ele->SetLineColor(kSpring);
	dE_ele->Draw("same");
	dE_ele->GetXaxis()->SetTitle("Energy(MeV)(electrons));
	dE_ele->GetYaxis()->SetTitle("Events");

	//Positirons
	
	TH1D *dE_pos=(TH1D*)f.Get("23");
	dE_pos->SetLineColor(kOrange);
	dE_pos->Draw("same");
	dE_pos->GetXaxis()->SetTitle("Energy(MeV)(positrons));
	dE_pos->GetYaxis()->SetTitle("Events");

	THStack *st = new THStack();
	st->Add(dE_gam);
	st->Add(dE_ele);
	st->Add(dE_pos);

	st->Draw("no stack");



	

	// T(1,0,0)-B(0,0,0)
	c1->cd(2);
	gPad->SetLogy();

	//Gammas
	
	TH1D *dE_gam_x=(TH1D*)fx.Get("7");
	dE_gam_x->SetLineColor(kMagenta);
	dE_gam_x->Draw("HIST");
	dE_gam_x->GetXaxis()->SetTitle("Energy(MeV)(gammas));
	dE_gam_x->GetYaxis()->SetTitle("Events");
	
	//Electrons
	
	TH1D *dE_ele_x=(TH1D*)fx.Get("15");
	dE_ele_x->SetLineColor(kSpring);
	dE_ele_x->Draw("same");
	dE_ele_x->GetXaxis()->SetTitle("Energy(MeV)(electrons));
	dE_ele_x->GetYaxis()->SetTitle("Events");
	
	//Positirons
	
	TH1D *dE_pos_x=(TH1D*)fx.Get("23");
	dE_pos_x->SetLineColor(kOrange);
	dE_pos_x->Draw("same");
	dE_pos_x->GetXaxis()->SetTitle("Energy(MeV)(positrons));
	dE_pos_x->GetYaxis()->SetTitle("Events");


	
	// T(0,1,0)-B(0,0,0)

	c1->cd(3);
	gPad->SetLogy();

	//Gammas
	
	TH1D *dE_gam_y=(TH1D*)fy.Get("7");
	dE_gam_y->SetLineColor(kMagenta);
	dE_gam_y->Draw("HIST");
	dE_gam_y->GetXaxis()->SetTitle("Energy(MeV)(gammas));
	dE_gam_y->GetYaxis()->SetTitle("Events");

	//Electrons
	
	TH1D *dE_ele_y=(TH1D*)fy.Get("15");
	dE_ele_y->SetLineColor(kSpring);
	dE_ele_y->Draw("same");
	dE_ele_y->GetXaxis()->SetTitle("Energy(MeV)(electrons));
	dE_ele_y->GetYaxis()->SetTitle("Events");

	//Positirons
	
	TH1D *dE_pos-y=(TH1D*)fy.Get("23");
	dE_pos_y->SetLineColor(kOrange);
	dE_pos_y->Draw("same");
	dE_pos_y->GetXaxis()->SetTitle("Energy(MeV)(positrons));
	dE_pos_y->GetYaxis()->SetTitle("Events");


	// T(0,0,1)-B(0,0,0)

	c1->cd(4);
	gPad->SetLogy();

	//Gammas
	
	TH1D *dE_gam_z=(TH1D*)fz.Get("7");
	dE_gam_z->SetLineColor(kMagenta);
	dE_gam_z->Draw("HIST");
	dE_gam_z->GetXaxis()->SetTitle("Energy(MeV)(gammas));
	dE_gam_z->GetYaxis()->SetTitle("Events");

	//Electrons
	
	TH1D *dE_ele-z=(TH1D*)fz.Get("15");
	dE_ele_z->SetLineColor(kSpring);
	dE_ele_z->Draw("same");
	dE_ele_z->GetXaxis()->SetTitle("Energy(MeV)(electrons));
	dE_ele_z->GetYaxis()->SetTitle("Events");

	//Positirons
	
	TH1D *dE_pos_z=(TH1D*)fz.Get("23");
	dE_pos_z->SetLineColor(kOrange);
	dE_pos_z->Draw("same");
	dE_pos_z->GetXaxis()->SetTitle("Energy(MeV)(positrons));
	dE_pos_z->GetYaxis()->SetTitle("Events");


	// T(1,0,0)-B(1,0,0)

	c1->cd(5);
	gPad->SetLogy();

	//Gammas
	
	TH1D *dE_gam_x2=(TH1D*)fx2.Get("7");
	dE_gam_x2->SetLineColor(kMagenta);
	dE_gam_x2->Draw("HIST");
	dE_gam_x2->GetXaxis()->SetTitle("Energy(MeV)(gammas));
	dE_gam_x2->GetYaxis()->SetTitle("Events");
	
	//Electrons
	
	TH1D *dE_ele_x2=(TH1D*)fx2.Get("15");
	dE_ele_x2->SetLineColor(kSpring);
	dE_ele_x2->Draw("same");
	dE_ele_x2->GetXaxis()->SetTitle("Energy(MeV)(electrons));
	dE_ele_x2->GetYaxis()->SetTitle("Events");

	//Positirons
	
	TH1D *dE_pos_x2=(TH1D*)fx2.Get("23");
	dE_pos_x2->SetLineColor(kOrange);
	dE_pos_x2->Draw("same");
	dE_pos_x2->GetXaxis()->SetTitle("Energy(MeV)(positrons));
	dE_pos_x2->GetYaxis()->SetTitle("Events");

	// T(0,1,0)-B(0,1,0)

	c1->cd(6);
	gPad->SetLogy();

	//Gammas
	
	TH1D *dE_gam_y2=(TH1D*)fy2.Get("7");
	dE_gam_y2->SetLineColor(kMagenta);
	dE_gam_y2->Draw("HIST");
	dE_gam_y2->GetXaxis()->SetTitle("Energy(MeV)(gammas));
	dE_gam_y2->GetYaxis()->SetTitle("Events");
	
	//Electrons
	
	TH1D *dE_ele_y2=(TH1D*)fy2.Get("15");
	dE_ele_y2->SetLineColor(kSpring);
	dE_ele_y2->Draw("same");
	dE_ele_y2->GetXaxis()->SetTitle("Energy(MeV)(electrons));
	dE_ele_y2->GetYaxis()->SetTitle("Events");

	//Positirons
	
	TH1D *dE_pos_y2=(TH1D*)fy2.Get("23");
	dE_pos_y2->SetLineColor(kOrange);
	dE_pos_y2->Draw("same");
	dE_pos_y2->GetXaxis()->SetTitle("Energy(MeV)(positrons));
	dE_pos_y2->GetYaxis()->SetTitle("Events");

	// T(0,0,1)-B(0,0,1)

	c1->cd(7);
	gPad->SetLogy();

	//Gammas
	
	TH1D *dE_gam_z2=(TH1D*)fz2.Get("7");
	dE_gam_z2->SetLineColor(kMagenta);
	dE_gam_z2->Draw("HIST");
	dE_gam_z2->GetXaxis()->SetTitle("Energy(MeV)(gammas));
	dE_gam_z2->GetYaxis()->SetTitle("Events");
	
	//Electrons
	
	TH1D *dE_ele_z2=(TH1D*)fz2.Get("15");
	dE_ele_z2->SetLineColor(kSpring);
	dE_ele_z2->Draw("same");
	dE_ele_z2->GetXaxis()->SetTitle("Energy(MeV)(electrons));
	dE_ele_z2->GetYaxis()->SetTitle("Events");

	//Positirons
	
	TH1D *dE_pos_z2=(TH1D*)fz2.Get("23");
	dE_pos_z2->SetLineColor(kOrange);
	dE_pos_z2->Draw("same");
	dE_pos_z2->GetXaxis()->SetTitle("Energy(MeV)(positrons));
	dE_pos_z2->GetYaxis()->SetTitle("Events");

	//c1->Update();
	//c1->cd();
	//c1->SaveAs("T_000-B_000.pdf");
	
}