-
- // set hits
- segmentation->SetHit(xhit,yhit,zhit);
- // get the list of nearest neighbours
- Int_t nList, xList[10], yList[10];
- segmentation->Neighbours(ix,iy,&nList,xList,yList);
-
- qp = 0;
- for (Int_t j=0; j<nList; j++){ // loop over neighbours
- if (xList[j]!=0) { // existing neighbour
- if (j==0||j==5||qp!=0) { // built-up cluster-size
-
- // neighbour real coordinates (just for checks here)
- Float_t x,y,z;
- segmentation->GetPadC(xList[j],yList[j],x,y,z);
- // set pad (fx fy & fix fiy are the current pad coord. & Id.)
- segmentation->SetPad(xList[j],yList[j]);
- // get the chamber (i.e. current strip) response
- qp=fResponse->IntXY(segmentation);
-
- if (qp > 0.5) {
- // --- store signal information for neighbours
- newclust[0][nnew]=qp; // total charge
- newclust[1][nnew]=segmentation->Ix(); // ix-pos. of pad
- newclust[2][nnew]=segmentation->Iy(); // iy-pos. of pad
- newclust[3][nnew]=twentyNano; // time of flight
- newclust[4][nnew]=segmentation->ISector(); // sector id
- newclust[5][nnew]=(Float_t) i; // counter
- nnew++;
- } // qp > 0.5
- } // built-up cluster-size
- } // existing neighbour
- } // loop over neighbours
- } // endif hit w/o strip
+
+ // set hits
+ segmentation->SetHit(id,xhit,yhit,zhit);
+ // get the list of nearest neighbours
+ Int_t nList, xList[10], yList[10];
+ segmentation->Neighbours(id,ix,iy,&nList,xList,yList);
+
+ qp = 0;
+ for (Int_t j=0; j<nList; j++){ // loop over neighbours
+ if (xList[j]!=0) { // existing neighbour
+ if (j==0||j==5||qp!=0) { // built-up cluster-size
+
+ // neighbour real coordinates (just for checks here)
+ Float_t x,y,z;
+ segmentation->GetPadC(id,xList[j],yList[j],x,y,z);
+ // set pad (fx fy & fix fiy are the current pad coord. & Id.)
+ segmentation->SetPad(id,xList[j],yList[j]);
+ // get the chamber (i.e. current strip) response
+ qp=fResponse->IntXY(id,segmentation);
+
+ if (qp > 0.5) {
+ // --- store signal information for neighbours
+ newclust[0][nnew]=qp; // total charge
+ newclust[1][nnew]=segmentation->Ix(); // ix-pos. of pad
+ newclust[2][nnew]=segmentation->Iy(); // iy-pos. of pad
+ newclust[3][nnew]=twentyNano; // time of flight
+ newclust[4][nnew]=segmentation->ISector(); // sector id
+ newclust[5][nnew]=(Float_t) i; // counter
+ nnew++;
+ } // qp > 0.5
+ } // built-up cluster-size
+ } // existing neighbour
+ } // loop over neighbours
+ } // endif hit w/o strip