+ // ---- equalize radial clusters distribution ----------------------
+
+ for (Int_t iPlane=0; iPlane<fNPlanes; iPlane++) {
+ for (Int_t iBin=1; iBin<=fHistClusterRadialPosition[iPlane]->GetNbinsX(); iBin++) {
+ Double_t rMin = fHistClusterRadialPosition[iPlane]->GetBinLowEdge(iBin); // in cm
+ Double_t rMax = fHistClusterRadialPosition[iPlane]->GetBinWidth(iBin) + rMin; // in cm
+ Double_t area = 100.*TMath::Pi()*(rMax*rMax - rMin*rMin); // in mm^2
+ Double_t density = fHistClusterRadialPosition[iPlane]->GetBinContent(iBin)/area;
+ fHistClusterRadialPosition[iPlane]->SetBinContent(iBin, density);
+ fHistClusterRadialPosition[iPlane]->SetBinError(iBin, fHistClusterRadialPosition[iPlane]->GetBinError(iBin)/area);
+ }
+ fHistClusterRadialPosition[iPlane] -> SetBinContent(1, fEv); // "scaler" bin
+ }
+
+ // -----------------------------------------------------------------
+