1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
17 #include "AliITSClusterFinderSPD.h"
18 #include "AliITSMapA1.h"
20 #include "AliITSdigit.h"
21 #include "AliITSRawCluster.h"
22 #include "AliITSRecPoint.h"
23 #include "AliITSsegmentation.h"
24 #include "AliITSresponse.h"
29 ClassImp(AliITSClusterFinderSPD)
31 //----------------------------------------------------------
32 AliITSClusterFinderSPD::AliITSClusterFinderSPD
33 (AliITSsegmentation *seg, TClonesArray *digits, TClonesArray *recp)
39 fNclusters= fClusters->GetEntriesFast();
42 fMap=new AliITSMapA1(fSegmentation,fDigits);
46 //_____________________________________________________________________________
47 AliITSClusterFinderSPD::AliITSClusterFinderSPD()
61 //_____________________________________________________________________________
62 AliITSClusterFinderSPD::~AliITSClusterFinderSPD()
65 if (fMap) delete fMap;
69 //__________________________________________________________________________
70 AliITSClusterFinderSPD::AliITSClusterFinderSPD(const AliITSClusterFinderSPD &source){
72 if(&source == this) return;
73 this->fClusters = source.fClusters ;
74 this->fNclusters = source.fNclusters ;
75 this->fMap = source.fMap ;
76 this->fDz = source.fDz ;
77 this->fDx = source.fDx ;
78 this->fMinNCells = source.fMinNCells ;
82 //_________________________________________________________________________
83 AliITSClusterFinderSPD&
84 AliITSClusterFinderSPD::operator=(const AliITSClusterFinderSPD &source) {
85 // Assignment operator
86 if(&source == this) return *this;
87 this->fClusters = source.fClusters ;
88 this->fNclusters = source.fNclusters ;
89 this->fMap = source.fMap ;
90 this->fDz = source.fDz ;
91 this->fDx = source.fDx ;
92 this->fMinNCells = source.fMinNCells ;
96 //_____________________________________________________________________________
97 void AliITSClusterFinderSPD::SetMap()
101 if(!fMap) fMap=new AliITSMapA1(fSegmentation,fDigits);
105 //_____________________________________________________________________________
107 void AliITSClusterFinderSPD::Find1DClusters()
109 // Find one dimensional clusters, i.e.
110 // in r*phi(x) direction for each colunm in z direction
112 AliITS *iTS=(AliITS*)gAlice->GetModule("ITS");
114 // retrieve the parameters
115 Int_t fNofPixels = fSegmentation->Npz();
116 Int_t fMaxNofSamples = fSegmentation->Npx();
118 // read in digits -> do not apply threshold
119 // signal in fired pixels is always 1
123 Int_t nofFoundClusters = 0;
126 for(k=0;k<fNofPixels;k++) {
128 Int_t mmax = 10; // a size of the window for the cluster finding
130 for(it=0;it<fMaxNofSamples;it++) {
138 for(m=0;m<mmax;m++) { // find the cluster inside the window
140 if(id >= fMaxNofSamples) break; // ! no possible for the fadc
142 if(fMap->TestHit(k,id) == kUnused) { // start of the cluster
144 if(lclx == 1) xstart = id;
148 if(lclx > 0 && fMap->TestHit(k,id) == kEmpty) {
149 // end of cluster if a gap exists
157 if(lclx == 0 && ilcl == 0) it = id; // no cluster in the window,
158 // continue the "it" loop
160 if(id >= fMaxNofSamples && lclx == 0) break; // the x row finished
162 if(id < fMaxNofSamples && ilcl == 0 && lclx > 0) {
163 // cluster end is outside of the window,
164 mmax += 5; // increase mmax and repeat the cluster
169 if(id >= fMaxNofSamples && lclx > 0) { // the x row finished but
170 xstop = fMaxNofSamples - 1; // the end cluster exists
174 // --- Calculate z and x coordinates for one dimensional clusters
176 if(ilcl == 1) { // new cluster exists
180 Float_t clusterCharge = 0.;
181 Float_t zpitch = fSegmentation->Dpz(k+1);
182 Float_t clusterZ, dummyX;
184 fSegmentation->GetPadCxz(dummy,k,dummyX,clusterZ);
185 Float_t zstart = clusterZ - 0.5*zpitch;
186 Float_t zstop = clusterZ + 0.5*zpitch;
187 Float_t clusterX = 0.;
188 Int_t xstartfull = xstart;
189 Int_t xstopfull = xstop;
190 Int_t clusterSizeX = lclx;
191 Int_t clusterSizeZ = 1;
194 for(its=xstart; its<=xstop; its++) {
196 if (fMap->GetHitIndex(k,its)>=0) firedpixel=1;
197 clusterCharge += firedpixel;
198 clusterX +=its + 0.5;
200 Float_t fRphiPitch = fSegmentation->Dpx(dummy);
201 clusterX /= (clusterSizeX/fRphiPitch); // center of gravity for x
204 //printf("ClusterZ ClusterX %f %f \n",clusterZ, clusterX);
206 // Int_t nclusters = fClusters->GetEntriesFast();
207 // cout << nclusters << " clusters" << endl;
208 // cout<< "Create point"<<endl;
210 // Write the points (coordinates and some cluster information) to the
211 // AliITSRawClusterSPD object
213 AliITSRawClusterSPD *clust = new AliITSRawClusterSPD(clusterZ,clusterX,clusterCharge,clusterSizeZ,clusterSizeX,xstart,xstop,xstartfull,xstopfull,zstart,zstop,k);
214 // fClusters->Add(point);
215 iTS->AddCluster(0,clust);
216 // cout << "Cluster at Ladder: " << fLadder << ", Detector: " <<fDetector<<endl;
218 // cout<<" end of cluster finding for Z pixel "<<endl;
220 } // new cluster (ilcl=1)
221 } // X direction loop (it)
222 } // Z direction loop (k)
229 //_____________________________________________________________________________
230 void AliITSClusterFinderSPD::GroupClusters()
232 // Find two dimensional clusters, i.e. group one dimensional clusters
233 // into two dimensional ones (go both in x and z directions).
235 // get number of clusters for this module
236 Int_t nofClusters = fClusters->GetEntriesFast();
237 nofClusters -= fNclusters;
239 AliITSRawClusterSPD *clusterI;
240 AliITSRawClusterSPD *clusterJ;
242 Int_t *label=new Int_t[nofClusters];
244 for(i=0; i<nofClusters; i++) label[i] = 0;
245 for(i=0; i<nofClusters; i++) {
246 if(label[i] != 0) continue;
247 for(j=i+1; j<nofClusters; j++) {
248 if(label[j] != 0) continue;
249 clusterI = (AliITSRawClusterSPD*) fClusters->At(i);
250 clusterJ = (AliITSRawClusterSPD*) fClusters->At(j);
251 Bool_t pair = clusterI->Brother(clusterJ,fDz,fDx);
254 // if((clusterI->XStop() == clusterJ->XStart()-1)||(clusterI->XStart()==clusterJ->XStop()+1)) cout<<"!! Diagonal cluster"<<endl;
256 cout << "clusters " << i << "," << j << " before grouping" << endl;
257 clusterI->PrintInfo();
258 clusterJ->PrintInfo();
260 clusterI->Add(clusterJ);
261 // cout << "remove cluster " << j << endl;
263 fClusters->RemoveAt(j);
265 cout << "cluster " << i << " after grouping" << endl;
266 clusterI->PrintInfo();
272 fClusters->Compress();
273 // Int_t totalNofClusters = fClusters->GetEntriesFast();
274 // cout << " Nomber of clusters at the group end ="<< totalNofClusters<<endl;
282 //_____________________________________________________________________________
284 void AliITSClusterFinderSPD::TracksInCluster()
287 // Find tracks creating one cluster
289 // get number of clusters for this module
290 Int_t nofClusters = fClusters->GetEntriesFast();
291 nofClusters -= fNclusters;
293 //printf("nofClusters fNclusters %d %d\n",nofClusters,fNclusters); //commented 18-09-00
295 Int_t i, ix, iz, jx, jz, xstart, xstop, zstart, zstop, nclx, nclz;
296 // Int_t signal, track0, track1, track2;
297 const Int_t trmax = 100;
298 Int_t cltracks[trmax], itr, tracki, ii, is, js, ie, ntr, tr0, tr1, tr2;
300 for(i=0; i<nofClusters; i++) {
302 memset(cltracks,-1,sizeof(int)*trmax);
305 AliITSRawClusterSPD *clusterI = (AliITSRawClusterSPD*) fClusters->At(i);
306 //printf("clusterI %p\n",clusterI); //commented 18-09-00
308 nclx = clusterI->NclX();
309 nclz = clusterI->NclZ();
310 xstart = clusterI->XStartf();
311 xstop = clusterI->XStopf();
312 zstart = clusterI->Zend()-nclz+1;
313 zstop = clusterI->Zend();
317 for(iz=0; iz<nclz; iz++) {
319 for(ix=0; ix<nclx; ix++) {
321 ind = fMap->GetHitIndex(jz,jx);
322 //printf("ind %d\n",ind); // commented 18-09-00
323 // this part must be wrong - the index of the digit
324 // can be zero but not -1 !!! comment out this part
325 // and replace it - see below
327 if(ind == 0 && iz >= 0 && ix > 0) {
330 if(ind == 0 && iz > 0 && ix >= 0) {
333 if(ind == 0 && iz == 0 && ix == 0 && i > 0) {
337 // change the conditions above to :
338 if (ind < 0) continue;
340 AliITSdigitSPD *dig = (AliITSdigitSPD*)fMap->GetHit(jz,jx);
341 if (!dig) printf("SPD: something wrong ! dig %p\n",dig);
344 track0=dig->fTracks[0];
345 track1=dig->fTracks[1];
346 track2=dig->fTracks[2];
348 for(itr=0; itr<3; itr++) {
349 tracki = dig->fTracks[itr];
350 //printf("tracki %d\n",tracki); //commented 18-09-00
353 cltracks[ii-1] = tracki;
359 for(is=0; is<trmax; is++) {
360 if(cltracks[is]<0) continue;
361 for(js=is+1; js<trmax; js++) {
362 if(cltracks[js]<0) continue;
363 if(cltracks[js]==cltracks[is]) cltracks[js]=-5;
368 for(ie=0; ie<trmax; ie++) {
369 if(cltracks[ie] >= 0) {
371 if(ntr==1) tr0=cltracks[ie];
372 if(ntr==2) tr1=cltracks[ie];
373 if(ntr==3) tr2=cltracks[ie];
377 if(ntr == 0) ntr = 1;
379 clusterI->SetNTracks(ntr);
380 clusterI->SetTracks(tr0,tr1,tr2);
385 //_____________________________________________________________________________
387 void AliITSClusterFinderSPD::GetRecPoints()
390 AliITS *iTS=(AliITS*)gAlice->GetModule("ITS");
392 // get number of clusters for this module
393 Int_t nofClusters = fClusters->GetEntriesFast();
394 nofClusters -= fNclusters;
395 const Float_t kconv = 1.0e-4;
396 const Float_t kRMSx = 12.0*kconv; // microns -> cm ITS TDR Table 1.3
397 const Float_t kRMSz = 70.0*kconv; // microns -> cm ITS TDR Table 1.3
399 Float_t spdLength = fSegmentation->Dz();
400 Float_t spdWidth = fSegmentation->Dx();
403 Int_t track0, track1, track2;
405 for(i=0; i<nofClusters; i++) {
407 AliITSRawClusterSPD *clusterI = (AliITSRawClusterSPD*) fClusters->At(i);
408 clusterI->GetTracks(track0, track1, track2);
411 rnew.SetX((clusterI->X() - spdWidth/2)*kconv);
412 rnew.SetZ((clusterI->Z() - spdLength/2)*kconv);
415 rnew.SetSigmaX2(kRMSx*kRMSx);
416 rnew.SetSigmaZ2(kRMSz*kRMSz);
417 rnew.fTracks[0]=track0;
418 rnew.fTracks[1]=track1;
419 rnew.fTracks[2]=track2;
420 iTS->AddRecPoint(rnew);
426 //_____________________________________________________________________________
428 void AliITSClusterFinderSPD::FindRawClusters()