00001 {
00002 gROOT->Reset();
00003
00004
00005
00006 TFile f = TFile("histograms.root");
00007 TCanvas* c1 = new TCanvas("c1", " ");
00008
00009 TH1D* hist1 = (TH1D*)f.Get("1");
00010 c1->SetLogy(1);
00011 c1->cd();
00012 c1->Update();
00013
00014 hist1->SetTitle("Total energy deposited in astronaut");
00015 hist1->GetXaxis()->SetTitle("edep (MeV)");
00016 hist1->GetYaxis()->SetTitle("nb events");
00017 hist1->Draw("HIST");
00018
00019 TH1D* hist2 = (TH1D*)f.Get("2");
00020 c1->SetLogy(0);
00021 c1->cd();
00022 c1->Update();
00023
00024 hist2->SetTitle("edep longitudinal profile");
00025 hist2->GetXaxis()->SetTitle("zlocal (mm)");
00026 hist2->GetYaxis()->SetTitle("Edep (MeV/mm)");
00027 hist2->Draw("HIST");
00028 }
