5 TDirectory *
dir = (TDirectory*)gDirectory->FindObjectAny(
"TAGH_doubles");
8 TH1F *doubles =(TH1F*)gDirectory->Get(
"BeforeMergingDoubles/BM2_Energy");
9 TH1F *
total =(TH1F*)gDirectory->Get(
"BeforeMergingDoubles/BM1_Energy");
11 TH1F *f_doubles = (TH1F*)doubles->Clone();
13 f_doubles->Divide(doubles, total);
16 TCanvas *
c1 =
new TCanvas(
"c1",
"TAGH double-hit fraction",150,10,990,660);
23 TCanvas*
c1 = gPad->GetCanvas();
25 gStyle->SetOptStat(
"");
26 f_doubles->SetTitle(
"TAGH double-hit fraction vs. energy");
27 f_doubles->SetTitleSize(0.045,
"XY");
28 f_doubles->GetXaxis()->SetTitle(
"TAGH energy [GeV]");
29 f_doubles->GetYaxis()->SetTitle(
"double-hit fraction");
30 f_doubles->SetAxisRange(0.0,1.0,
"Y");