// divides generated by measured to get the correction
//
-// if (!fhCorr) {
-// fhCorr = (TH1*)fhGene->Clone("correction");
-// fhCorr->SetTitle(Form("%s correction",GetTitle()));
-// fhCorr->Reset();
-// }
-
if (!fhMeas || !fhGene || !fhCorr) {
AliDebug(AliLog::kError, "measured or generated histograms not available");
return;
if (!fhCorr)
return;
- for (Int_t x=1; x<=fhCorr->GetNbinsX(); ++x)
- for (Int_t y=1; y<=fhCorr->GetNbinsY(); ++y)
- for (Int_t z=1; z<=fhCorr->GetNbinsZ(); ++z)
+ for (Int_t x=0; x<=fhCorr->GetNbinsX()+1; ++x)
+ for (Int_t y=0; y<=fhCorr->GetNbinsY()+1; ++y)
+ for (Int_t z=0; z<=fhCorr->GetNbinsZ()+1; ++z)
+ {
fhCorr->SetBinContent(x, y, z, 1);
+ fhCorr->SetBinError(x, y, z, 0);
+ }
}
//____________________________________________________________________
fhMeas->Scale(factor);
fhGene->Scale(factor);
}
+
+//____________________________________________________________________
+void AliCorrectionMatrix::ResetErrorsOnCorrections()
+{
+ // set the errors on the correction matrix to 0
+
+ if (!fhCorr)
+ return;
+
+ for (Int_t x=0; x<=fhCorr->GetNbinsX()+1; ++x)
+ for (Int_t y=0; y<=fhCorr->GetNbinsY()+1; ++y)
+ for (Int_t z=0; z<=fhCorr->GetNbinsZ()+1; ++z)
+ fhCorr->SetBinError(x, y, z, 0);
+}