00001 {
00002 gROOT->Reset();
00003
00004 TCanvas* c1 = new TCanvas("c1", " ");
00005 c1->SetLogy(1);
00006 c1->cd();
00007 c1->Update();
00008
00009
00010 TFile f = TFile("histograms.root");
00011 TH1D* hist = (TH1D*)f.Get("1");
00012 hist->SetTitle("Total energy deposited in astronaut");
00013 hist->GetXaxis()->SetTitle("edep (MeV)");
00014 hist->GetYaxis()->SetTitle("nb events");
00015 hist->Draw("HIST");
00016 }