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


void EuWLisa(){


TString infile = "/net/user/eleser/L5_6062_Erzgebirge_Ext.root";
TString Treename = "aTree";
TString Zenith = "MPEFit_Zenith";
TString Energy = "MPEFitMuE_Energy";

double zen, en, Minen, Maxen;


TChain* inChain = new TChain(Treename);
inChain->Add(infile);
inChain->SetBranchStatus("*",1);
inChain->SetBranchAddress(Energy,&en);
inChain->SetBranchAddress(Zenith,&zen);


for(int i = 0 ; i< inChain->GetEntries(); ++i){
    inChain->GetEntry(i);
    Minen=TMath::Min(en, Minen);
    Maxen=TMath::Max(en, Maxen);
    
}

cout<<Minen<<"   "<<Maxen<<endl;


TCanvas *can = new TCanvas("can","",600,400);
TH2D *histo2D = new TH2D("histo2","2D", 50, 1.5, 3.2 , 50, 0, 10000);                                              //"h2p","2D", n, Mintime, Maxtime, m, zenmin, zenmax
/*histo2->SetName("Box");                                             //Name der Fehlerbox im Fenster
histo2->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);
    histo2D->Fill(zen, en);
}
histo2D->GetXaxis()->SetTitle("Zenithwinkel");
histo2D->GetYaxis()->SetTitle("Energie");
//histo2D->GetZaxis()->SetTitle("Counts");
histo2D->Draw("COLZ");                                                  //COLZ ist die Grafik wie das Histogramm später aussehen soll
}