+ //printf("AliEMCALRecPoint::EvalDispersion() : Dispersion %f \n",fDispersion);
+}
+
+//____________________________________________________________________________
+void AliEMCALRecPoint::EvalDistanceToBadChannels(AliCaloCalibPedestal* caloped)
+{
+ //For each EMC rec. point set the distance to the nearest bad channel.
+ //AliInfo(Form("%d bad channel(s) found.\n", caloped->GetDeadTowerCount()));
+ //It is done in cell units and not in global or local position as before (Sept 2010)
+
+ if(!caloped->GetDeadTowerCount()) return;
+
+ //Get channels map of the supermodule where the cluster is.
+ TH2D* hMap = caloped->GetDeadMap(fSuperModuleNumber);
+
+ Int_t dRrow, dReta;
+ Float_t minDist = 10000.;
+ Float_t dist = 0.;
+ Int_t nSupMod, nModule;
+ Int_t nIphi, nIeta;
+ Int_t iphi, ieta;
+ fDigitIndMax = GetMaximalEnergyIndex();
+ fGeomPtr->GetCellIndex(fAbsIdList[fDigitIndMax], nSupMod,nModule,nIphi,nIeta);
+ fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
+
+ //Loop on tower status map
+ for(Int_t irow = 0; irow < AliEMCALGeoParams::fgkEMCALRows; irow++){
+ for(Int_t icol = 0; icol < AliEMCALGeoParams::fgkEMCALCols; icol++){
+ //Check if tower is bad.
+ if(hMap->GetBinContent(icol,irow)==AliCaloCalibPedestal::kAlive) continue;
+ //printf("AliEMCALRecPoint::EvalDistanceToBadChannels() - Bad channel in SM %d, col %d, row %d\n",iSM,icol, irow);
+
+ dRrow=TMath::Abs(irow-iphi);
+ dReta=TMath::Abs(icol-ieta);
+ dist=TMath::Sqrt(dRrow*dRrow+dReta*dReta);
+ if(dist < minDist) minDist = dist;
+
+ }
+ }
+
+ //In case the cluster is shared by 2 SuperModules, need to check the map of the second Super Module
+ if (fSharedCluster) {
+ TH2D* hMap2 = 0;
+ Int_t nSupMod2 = -1;
+
+ //The only possible combinations are (0,1), (2,3) ... (10,11)
+ if(fSuperModuleNumber%2) nSupMod2 = fSuperModuleNumber-1;
+ else nSupMod2 = fSuperModuleNumber+1;
+ hMap2 = caloped->GetDeadMap(nSupMod2);
+
+ //Loop on tower status map of second super module
+ for(Int_t irow = 0; irow < AliEMCALGeoParams::fgkEMCALRows; irow++){
+ for(Int_t icol = 0; icol < AliEMCALGeoParams::fgkEMCALCols; icol++){
+ //Check if tower is bad.
+ if(hMap2->GetBinContent(icol,irow)==AliCaloCalibPedestal::kAlive) continue;
+ //printf("AliEMCALRecPoint::EvalDistanceToBadChannels() - Bad channel in SM %d, col %d, row %d\n",iSM,icol, irow);
+ dRrow=TMath::Abs(irow-iphi);
+
+ if(fSuperModuleNumber%2) {
+ dReta=TMath::Abs(icol-(AliEMCALGeoParams::fgkEMCALCols+ieta));
+ }
+ else {
+ dReta=TMath::Abs(AliEMCALGeoParams::fgkEMCALCols+icol-ieta);
+ }
+
+ dist=TMath::Sqrt(dRrow*dRrow+dReta*dReta);
+ if(dist < minDist) minDist = dist;
+
+ }
+ }
+
+ }// shared cluster in 2 SuperModules
+
+ fDistToBadTower = minDist;
+ //printf("AliEMCALRecPoint::EvalDistanceToBadChannel() - Distance to Bad is %f cm, shared cluster? %d \n",fDistToBadTower,fSharedCluster);