#include <TH2.h>
#include <iostream>
#include <TCanvas.h>
#include <TChain.h>
using namespace std;


void winkelvergleich(){
	
	
	TString infile = "/net/user/eleser/L5_6062_Erzgebirge_Ext_Cut.root";
	TString Treename = "aTree";
	TString Zenith = "MCPrimary_Zenith";
	TString Energy = "MPEFit_Zenith";
	
	double mpefitzen, primaryzen;
	
	
	TChain* inChain = new TChain(Treename);
	inChain->Add(infile);
	inChain->SetBranchStatus("*",1);
	inChain->SetBranchAddress(mpefitzenith,&mpefitzen);
	inChain->SetBranchAddress(primaryzenith,&primaryzen);
	
		
	TCanvas *can = new TCanvas("can","",600,400);
	TH2D *hzenith = new TH2D("hzenith2D","2D", 50, 1.5, 3.2 , 50, 1.5, 3.2);  //"h2p","2D", n, MPEFitMinwinkel, MPEFitMaxwinkel, m, Primaryzenmin, primaryzenmax
	//*histo2D->SetName("Box");                                             //Name der Fehlerbox im Fenster
	 //histo2D->SetTitle("2D-Histo"); */                                      //Name des Titels
	gStyle->SetPalette(1);                                              //Setzt die Farbpalette fest, 2 ist grau
	
	for(int i = 0 ; i< inChain->GetEntries(); ++i){
		inChain->GetEntry(i);
		hzenith->Fill(mpefitzen, primaryzen);							//MPEFItzen = MPEFitzen*(360/(2*Pi));								//soll aus radiant werten grad werte machen
	}															//Primaryzen = Primaryzen*(360/(2*Pi)
	
	hzenith->GetXaxis()->SetTitle("MPEFit Zenith");
	hzenith->GetYaxis()->SetTitle("Primary Zenith");
	//histo2D->GetZaxis()->SetTitle("Counts");
	hzenith->Draw("COLZ");                                                  //COLZ ist die Grafik wie das Histogramm spŠter aussehen soll
}