+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDReconstructor::IsDigSurvive(Int_t *pUserCut, AliHMPIDDigit *pDig)const
+{
+//Check if the current digit survive to a riapllied sigma cut
+//Arguments: pDig pointer to the current digit
+// Returns: kTRUE if charge > mean+n*sigma
+ Int_t iCh = pDig->Ch();
+ Int_t iDaqSigCut =(Int_t)fDaqSig->At(iCh)->GetUniqueID();
+ if(pUserCut[iCh]<=iDaqSigCut) return kTRUE;
+ TMatrixF *pM = (TMatrixF*)fDaqSig->At(pDig->Ch());
+ Float_t sig = (*pM)(pDig->PadChX(),pDig->PadChY());
+ if(pDig->Q()>pUserCut[iCh]*sig) return kTRUE;
+ else return kFALSE;
+}