- Float_t dFiducialEta=rHeader->GetFiducialEtaMax()-rHeader->GetFiducialEtaMin();
- Double_t dMeanDist=TMath::Sqrt(2*dFiducialEta*pi/fNin);
- Double_t *xEta = xData[0];
- Double_t *xPhi = xData[1];
+ Float_t dFidEtaMax = rHeader->GetFiducialEtaMax();
+ Float_t dFidEtaMin = rHeader->GetFiducialEtaMin();
+ Float_t dFiducialEta= dFidEtaMax - dFidEtaMin;
+ Double_t *xEta = xData[0];
+ Double_t *xPhi = xData[1];
+ Int_t nEff = 0;
+ for (Int_t i=0; i<fNeff; i++) if (xEta[i]<dFidEtaMax && xEta[i]>dFidEtaMin) nEff++;
+ Double_t dMeanDist=TMath::Sqrt(2*dFiducialEta*pi/nEff);