00001 {
00002 TFile* f = new TFile("out.root");
00003 TH1* h1 = (TH1*)f->Get("GeSpectrum/hEnergy");
00004 TH1* h2 = (TH1*)f->Get("GeSpectrum/hEnergyFarOOT");
00005 TH1* h3 = (TH1*)f->Get("GeSpectrum/hEnergyOOT");
00006 TH1* h4 = ((TH2*)f->Get("GeSpectrum/hTimeEnergy"))->ProjectionY();
00007 h1->SetStats(0);
00008 h2->SetStats(0);
00009 h3->SetStats(0);
00010 h4->SetStats(0);
00011 h1->GetXaxis()->SetRangeUser(100., 15000.);
00012 h2->GetXaxis()->SetRangeUser(100., 15000.);
00013 h3->GetXaxis()->SetRangeUser(100., 15000.);
00014 h4->GetXaxis()->SetRangeUser(100., 15000.);
00015 h1->SetTitle("Ge Si16P");
00016 h2->SetTitle("Ge Si16P");
00017 h3->SetTitle("Ge Si16P");
00018 h4->SetTitle("Ge Si16P");
00019 h1->GetXaxis()->SetTitle("Energy (ADC)"); h1->GetYaxis()->SetTitle("Counts");
00020 h2->GetXaxis()->SetTitle("Energy (ADC)"); h2->GetYaxis()->SetTitle("Counts");
00021 h3->GetXaxis()->SetTitle("Energy (ADC)"); h3->GetYaxis()->SetTitle("Counts");
00022 h4->GetXaxis()->SetTitle("Energy (ADC)"); h4->GetYaxis()->SetTitle("Counts");
00023
00024 TLegend *leg = new TLegend(0.7,0.7,0.9,0.9);
00025 leg->AddEntry(h1, "All");
00026 leg->AddEntry(h2, "Out of time (> +/- 5 us)");
00027 leg->AddEntry(h3, "Out of time (> +/-500 ns)");
00028 leg->AddEntry(h4, "In Time (< +/-500 ns)");
00029
00030 h1->SetLineColor(kBlack);
00031 h2->SetLineColor(kBlue);
00032 h3->SetLineColor(kRed);
00033 h4->SetLineColor(kGreen);
00034
00035 new TCanvas;
00036 h1->DrawClone();
00037 h2->DrawClone("SAME");
00038 h3->DrawClone("SAME");
00039 h4->DrawClone("SAME");
00040 leg->DrawClone("SAME");
00041 (new TCanvas)->SetLogy();
00042 h1->GetYaxis()->SetRangeUser(1., 10000.);
00043 h2->GetYaxis()->SetRangeUser(1., 10000.);
00044 h3->GetYaxis()->SetRangeUser(1., 10000.);
00045 h4->GetYaxis()->SetRangeUser(1., 10000.);
00046 h1->Draw();
00047 h2->Draw("SAME");
00048 h3->Draw("SAME");
00049 h4->Draw("SAME");
00050 leg->Draw("SAME");
00051 }
