+
+ if(mult==0){ //too soft cluster
+ fMulDigit =0 ;
+ fAmp = 0.; //Recalculate total energy
+ }
+
+ //Remove non-connected cells
+ Int_t index[mult] ;
+ Bool_t used[mult] ;
+ for(Int_t i=0; i<mult; i++) used[i]=0 ;
+ Int_t inClu=0 ;
+ Double_t eMax=0. ;
+ //find maximum
+ for(Int_t iDigit=0; iDigit<mult; iDigit++) {
+ if(eMax<tempoE[iDigit]){
+ eMax=tempoE[iDigit];
+ index[0]=iDigit ;
+ inClu=1 ;
+ }
+ }
+ if(mult>0)
+ used[index[0]]=kTRUE ; //mark as used
+ for(Int_t i=0; i<inClu; i++){
+ AliPHOSDigit * digit = (AliPHOSDigit *) digits->At(tempo[index[i]]) ;
+ for(Int_t iDigit=0 ;iDigit<mult; iDigit++){
+ if(used[iDigit]) //already used
+ continue ;
+ AliPHOSDigit * digitN = (AliPHOSDigit *) digits->At(tempo[iDigit]) ;
+ if(AreNeighbours(digit,digitN)){
+ index[inClu]= iDigit ;
+ inClu++ ;
+ used[iDigit]=kTRUE ;
+ }
+ }
+ }
+
+ fMulDigit = inClu ;