00001
00002
00003
00004 {
00005
00006 TFile *f;
00007 TFile *fa;
00008 TFile *fb;
00009 TFile *fc;
00010
00011 f = TFile::Open("histograms_Htogammagamma_Hm350.00_Nevt10000.root");
00012
00013 TH1F *this = (TH1F*) f->Get("Hmass");
00014
00015 this->Rebin(4);
00016 this->SetTitle("Input Higgs mass 350 GeV/c^{2}");
00017 this->GetXaxis()->SetTitle("#gamma #gamma Invariant Mass (Gev/c^{2})");
00018 this->GetXaxis()->SetRangeUser(0,2000);
00019
00020 fa = TFile::Open("histograms_Htogammagamma_Hm500.00_Nevt10000.root");
00021
00022 TH1F *thisa = (TH1F*) fa->Get("Hmass");
00023
00024 thisa->Rebin(4);
00025 thisa->SetTitle("Input Higgs mass 500 GeV/c^{2}");
00026 thisa->GetXaxis()->SetTitle("#gamma #gamma Invariant Mass (Gev/c^{2})");
00027 thisa->GetXaxis()->SetRangeUser(0,2000);
00028
00029 fb = TFile::Open("histograms_Htogammagamma_Hm650.00_Nevt10000.root");
00030
00031 TH1F *thisb = (TH1F*) fb->Get("Hmass");
00032
00033 thisb->Rebin(4);
00034 thisb->SetTitle("Input Higgs mass 650 GeV/c^{2}");
00035 thisb->GetXaxis()->SetTitle("#gamma #gamma Invariant Mass (Gev/c^{2})");
00036 thisb->GetXaxis()->SetRangeUser(0,2000);
00037
00038 fc = TFile::Open("histograms_Htogammagamma_Hm800.00_Nevt10000.root");
00039
00040 TH1F *thisc = (TH1F*) fc->Get("Hmass");
00041
00042 thisc->Rebin(4);
00043 thisc->SetTitle("Input Higgs mass 800 GeV/c^{2}");
00044 thisc->GetXaxis()->SetTitle("#gamma #gamma Invariant Mass (Gev/c^{2})");
00045 thisc->GetXaxis()->SetRangeUser(0,2000);
00046
00047 TCanvas *c1 = new TCanvas ("c1","Higgs plots");
00048 c1->Divide(2,2);
00049 c1->cd(1);
00050 this->Draw();
00051 c1->cd(2);
00052 thisa->Draw();
00053 c1->cd(3);
00054 thisb->Draw();
00055 c1->cd(4);
00056 thisc->Draw();
00057
00058 }
