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 **************************************************************************/
20 ///////////////////////////////////////////////////////////////////////////////
22 // This class contains the points for the ALICE event display //
26 <img src="gif/AliRICHpointsClass.gif">
31 ///////////////////////////////////////////////////////////////////////////////
32 #include "AliRICHdisplay.h"
33 #include "AliRICHpoints.h"
39 const Int_t MAX_Nipx=400, MAX_Nipy=800;
41 ClassImp(AliRICHpoints)
43 //_____________________________________________________________________________
44 AliRICHpoints::AliRICHpoints()
47 // Default constructor
54 //_____________________________________________________________________________
55 AliRICHpoints::AliRICHpoints(Int_t npoints)
59 // Standard constructor
66 //_____________________________________________________________________________
67 AliRICHpoints::~AliRICHpoints()
77 //_____________________________________________________________________________
78 void AliRICHpoints::DumpHit()
81 // Dump hit corresponding to this point
83 AliRICHhit *hit = GetHit();
87 //_____________________________________________________________________________
88 void AliRICHpoints::DumpDigit()
91 // Dump digit corresponding to this point
93 AliRICHdigit *digit = GetDigit();
94 if (digit) digit->Dump();
97 //_____________________________________________________________________________
98 void AliRICHpoints::InspectHit()
101 // Inspect hit corresponding to this point
103 AliRICHhit *hit = GetHit();
104 if (hit) hit->Inspect();
107 //_____________________________________________________________________________
108 void AliRICHpoints::InspectDigit()
111 // Inspect digit corresponding to this point
113 AliRICHdigit *digit = GetDigit();
114 if (digit) digit->Inspect();
117 //_____________________________________________________________________________
118 Int_t AliRICHpoints::GetTrackIndex()
121 // Dump digit corresponding to this point
123 printf("GetTrackIndex - fTrackIndex %d \n",fTrackIndex);
128 //_____________________________________________________________________________
129 AliRICHhit *AliRICHpoints::GetHit() const
132 // Returns pointer to hit index in AliRun::fParticles
134 AliRICH *RICH = (AliRICH*)gAlice->GetDetector("RICH");
135 gAlice->TreeH()->GetEvent(fTrackIndex);
136 TClonesArray *RICHhits = RICH->Hits();
137 Int_t nhits = RICHhits->GetEntriesFast();
138 if (fHitIndex < 0 || fHitIndex >= nhits) return 0;
139 return (AliRICHhit*)RICHhits->UncheckedAt(fHitIndex);
142 //_____________________________________________________________________________
143 AliRICHdigit *AliRICHpoints::GetDigit() const
146 // Returns pointer to digit index in AliRun::fParticles
149 AliRICHdisplay *display=(AliRICHdisplay*)gAlice->Display();
150 Int_t chamber=display->GetChamber();
151 Int_t cathode=display->GetCathode();
153 AliRICH *RICH = (AliRICH*)gAlice->GetDetector("RICH");
154 TClonesArray *RICHdigits = RICH->DigitsAddress(chamber-1);
155 gAlice->TreeD()->GetEvent(cathode);
156 Int_t ndigits = RICHdigits->GetEntriesFast();
157 if (fDigitIndex < 0 || fDigitIndex >= ndigits) return 0;
158 return (AliRICHdigit*)RICHdigits->UncheckedAt(fDigitIndex);
160 //_____________________________________________________________________________
162 const AliRICHdigit *dig;
164 Bin() {dig=0; idx=-1;}
167 struct PreCluster : public AliRICHreccluster {
168 const AliRICHdigit* summit;
172 PreCluster() : AliRICHreccluster() {cut=npeaks=0;}
174 //_____________________________________________________________________________
176 static void FindCluster(AliRICHchamber *iChamber, AliRICHsegmentation *segmentation, int i, int j, Bin bins[MAX_Nipx][MAX_Nipy], PreCluster &c)
184 printf("I'm in FindCluster \n");
187 Int_t q=b.dig->fSignal;
189 printf("FindCluster - i j q %d %d %d\n",i,j,q);
195 if (b.idx >= 0 && b.idx != c.idx) {
200 if (q > TMath::Abs(c.summit->fSignal)) c.summit=b.dig;
202 Int_t npx = segmentation->Npx();
203 Int_t npy = segmentation->Npy();
205 segmentation->GetPadCxy(i-npx, j-npy, x,y);
206 printf("FindCluster - x y %f %f \n",x,y);
213 b.dig = 0; b.idx = c.idx;
215 if (bins[i-1][j].dig) FindCluster(iChamber,segmentation,i-1,j,bins,c);
216 if (bins[i][j-1].dig) FindCluster(iChamber,segmentation,i,j-1,bins,c);
217 if (bins[i+1][j].dig) FindCluster(iChamber,segmentation,i+1,j,bins,c);
218 if (bins[i][j+1].dig) FindCluster(iChamber,segmentation,i,j+1,bins,c);
222 //_____________________________________________________________________________
224 void AliRICHpoints::GetCenterOfGravity()
227 // simple RICH cluster finder from digits -- finds neighbours and
228 // calculates center of gravity for the cluster
230 const Int_t MAX_nipx=400, MAX_nipy=800;
231 printf("\n Hallo world");
232 AliRICHdisplay *display=(AliRICHdisplay*)gAlice->Display();
233 Int_t chamber=display->GetChamber();
234 Int_t cathode=display->GetCathode();
236 AliRICH *RICH = (AliRICH*)gAlice->GetDetector("RICH");
237 AliRICHchamber *iChamber;
238 AliRICHsegmentation *segmentation;
239 iChamber =&(RICH->Chamber(chamber-1));
240 segmentation=iChamber->GetSegmentationModel(cathode);
241 Int_t npx = segmentation->Npx();
242 Int_t npy = segmentation->Npy();
243 Float_t zpos=iChamber->ZPosition();
245 TClonesArray *RICHdigits = RICH->DigitsAddress(chamber-1);
246 gAlice->TreeD()->GetEvent(cathode);
247 Int_t ndigits = RICHdigits->GetEntriesFast();
248 if (fDigitIndex < 0 || fDigitIndex >= ndigits) return;
251 dig=(AliRICHdigit*)RICHdigits->UncheckedAt(fDigitIndex);
252 Int_t ipx=dig->fPadX;
253 Int_t ipy=dig->fPadY;
254 Bin bins[MAX_nipx][MAX_nipy];
255 bins[ipx+npx][ipy+npy].dig=dig;
259 for (ndig=0; ndig<ndigits; ndig++) {
260 dig = (AliRICHdigit*)RICHdigits->UncheckedAt(ndig);
261 int i=dig->fPadX, j=dig->fPadY;
262 bins[i+npx][j+npy].dig=dig;
265 PreCluster c; c.summit=bins[ipx+npx][ipy+npy].dig; c.idx=ncls;
266 FindCluster(iChamber,segmentation,ipx+npx, ipy+npy, bins, c);
268 printf("GetCenterOfGravity -- more than one peak");
272 printf("GetCenterOfGravity - c.fX c.fY c.fQ %f %f %d \n",c.fX,c.fY,c.fQ);
274 c.fTracks[0]=c.summit->fTracks[0];
275 c.fTracks[1]=c.summit->fTracks[1];
276 c.fTracks[2]=c.summit->fTracks[2];
278 AliRICHpoints *points = 0;
279 points = new AliRICHpoints(1);
280 points->SetMarkerColor(kYellow);
281 points->SetMarkerStyle(5);
282 points->SetMarkerSize(1.);
283 points->SetPoint(0,c.fX,c.fY,zpos);
286 printf("GetCenterOfGravity -- ncls %d \n",ncls);