TH1->Divide problem displaying 1?

Hi,

I’m trying to divide two histograms with the “B” option. When the numerator and denominator are equal the resulting histogram does not show an entry at 1? I’m using 4.00/08 (I can’t use a more recent release). In the example code below, only the fist bin has an entry, bins 2-5 should have content=1, but it doesn’t get displayed.

Thanks,
Andrew

TH1F* den_hist = new TH1F(“den_hist”,“denominator”,5,0,5);
den_hist->Sumw2();
den_hist->SetBinContent(1,1.0);
den_hist->SetBinContent(2,2.0);
den_hist->SetBinContent(3,3.0);
den_hist->SetBinContent(4,4.0);
den_hist->SetBinContent(5,5.0);

TH1F* num_hist = new TH1F(“num_hist”,“numerator”,5,0,5);
num_hist->Sumw2();
num_hist->SetBinContent(1,0.9);
num_hist->SetBinContent(2,2.0);
num_hist->SetBinContent(3,3.0);
num_hist->SetBinContent(4,4.0);
num_hist->SetBinContent(5,5.0);

TH1F* efficiency_hist = new TH1F(“efficiency_hist”,“num/den”,5,0,5);
efficiency_hist->Sumw2();
efficiency_hist->SetLineWidth(3);
efficiency_hist->SetMarkerStyle(8);
efficiency_hist->Divide( num_hist, den_hist,1,1,“B” );
efficiency_hist->Draw();

You should use version 4.04 or newer

Rene