-// clusterPads[1] = col + 0.5
-// + (clusterSignal[2] - clusterSignal[0])
-// / clusterCharge;
- clusterPads[1] = 0.5
- + (clusterSignal[2] - clusterSignal[0])
- / clusterCharge;
+ for (Int_t i=0;i<5;i++) padSignal[i]=0;
+ padSignal[2] = TMath::Abs(digits->GetDataUnchecked(row,col,time)); // central pad
+ padSignal[1] = TMath::Abs(digits->GetDataUnchecked(row,col-1,time)); // left pad
+ padSignal[3] = TMath::Abs(digits->GetDataUnchecked(row,col+1,time)); // right pad
+ if (col>2 &&TMath::Abs(digits->GetDataUnchecked(row,col-2,time)<padSignal[1])){
+ padSignal[0] = TMath::Abs(digits->GetDataUnchecked(row,col-2,time));
+ }
+ if (col<nColMax-3 &&TMath::Abs(digits->GetDataUnchecked(row,col+2,time)<padSignal[3])){
+ padSignal[4] = TMath::Abs(digits->GetDataUnchecked(row,col+2,time));
+ }
+ clusterPads[1] = GetCOG(padSignal);
+ Double_t check = fPar->LUTposition(iplan,clusterSignal[0]
+ ,clusterSignal[1]
+ ,clusterSignal[2]);
+ // Float_t diff = clusterPads[1] - check;