Adding more bins in QA (Alis)
[u/mrichter/AliRoot.git] / MUON / AliMUONPreClusterFinder.cxx
CommitLineData
18466557 1/**************************************************************************
2* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3* *
4* Author: The ALICE Off-line Project. *
5* Contributors are mentioned in the code where appropriate. *
6* *
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**************************************************************************/
15
16// $Id$
17
18#include "AliMUONPreClusterFinder.h"
19
2103dc37 20#include "AliCodeTimer.h"
21
18466557 22#include "AliMUONCluster.h"
18466557 23#include "AliMUONPad.h"
3887d11e 24#include "AliMUONVDigit.h"
25#include "AliMUONVDigitStore.h"
e22f734a 26
27#include "AliMpArea.h"
28#include "AliMpConstants.h"
29#include "AliMpVSegmentation.h"
30
31#include "AliLog.h"
32
33#include <Riostream.h>
9e41a340 34#include <TObjArray.h>
e22f734a 35#include <TVector2.h>
18466557 36
3d1463c8 37//-----------------------------------------------------------------------------
18466557 38/// \class AliMUONPreClusterFinder
39///
40/// Implementation of AliMUONVClusterFinder
41///
42/// This class simply find adjacent pads to form clusters
43///
44/// \author Laurent Aphecetche
3d1463c8 45//-----------------------------------------------------------------------------
18466557 46
47ClassImp(AliMUONPreClusterFinder)
48
49//_____________________________________________________________________________
50AliMUONPreClusterFinder::AliMUONPreClusterFinder()
51: AliMUONVClusterFinder(),
9ec6a948 52 fClusters("AliMUONCluster"),
e22f734a 53 fPads(0x0),
54 fDetElemId(0),
861d6ce8 55 fArea(),
56 fShouldAbort(kFALSE)
18466557 57{
e22f734a 58 /// ctor
18466557 59}
60
61//_____________________________________________________________________________
62AliMUONPreClusterFinder::~AliMUONPreClusterFinder()
63{
e22f734a 64 /// dtor : note we're owner of the clusters, but not of the pads
18466557 65}
66
67//_____________________________________________________________________________
68Bool_t
69AliMUONPreClusterFinder::UsePad(const AliMUONPad& pad)
70{
71 /// Add a pad to the list of pads to be considered
72 if ( pad.DetElemId() != fDetElemId )
73 {
74 AliError(Form("Cannot add pad from DE %d to this cluster finder which is "
75 "currently dealing with DE %d",pad.DetElemId(),fDetElemId));
76 return kFALSE;
77 }
78
9e41a340 79 fPads[pad.Cathode()]->Add(new AliMUONPad(pad));
18466557 80 // FIXME: should set the ClusterId of that new pad to be -1
81 return kTRUE;
82}
83
84//_____________________________________________________________________________
85Bool_t
e22f734a 86AliMUONPreClusterFinder::Prepare(Int_t detElemId,
9e41a340 87 TObjArray* pads[2],
e22f734a 88 const AliMpArea& area)
18466557 89{
90 /// Prepare for clustering, by giving access to segmentations and digit lists
e22f734a 91
ca0b0ca1 92 fClusters.Delete();
9ec6a948 93
e22f734a 94 fPads = pads;
95 fDetElemId = detElemId;
96 fArea = area;
d592a053 97
861d6ce8 98 fShouldAbort = kFALSE;
99
18466557 100 return kTRUE;
101}
102
103//_____________________________________________________________________________
104void
105AliMUONPreClusterFinder::AddPad(AliMUONCluster& cluster, AliMUONPad* pad)
106{
107 /// Add a pad to a cluster
861d6ce8 108
2103dc37 109 if ( cluster.IsMonoCathode() && cluster.Multiplicity() > 199 )
861d6ce8 110 {
2103dc37 111 /// FIXME : we should at that point really find a better way to remove "bad" preclusters,
112 /// like e.g. computing the charge dispersion (the lower, the most probably we have noise cluster)
113 /// and/or mean charge per pad (if too close to LowestPadCharge, again that's a noise cluster...
114 /// *BUT* this needs carefull testing !
861d6ce8 115 fShouldAbort = kTRUE;
116 return;
117 }
118
ed90feae 119 AliMUONPad* addedPad = cluster.AddPad(*pad);
18466557 120
121 Int_t cathode = pad->Cathode();
9e41a340 122 TObjArray& padArray = *fPads[cathode];
d592a053 123 // WARNING: this Remove method uses the AliMUONPad::IsEqual if that method is
124 // present (otherwise just compares pointers) : so that one must be correct
125 // if implemented !
ed90feae 126 TObject* o = padArray.Remove(pad);
d592a053 127// if (!o)
128// {
129// AliFatal("Oups. Could not remove pad from pads to consider. Aborting as anyway "
130// " we'll get an infinite loop. Please check the AliMUONPad::IsEqual method"
131// " as the first suspect for failed remove");
132// }
ed90feae 133 delete o;
134
18466557 135 TIter next(&padArray);
136 AliMUONPad* testPad;
d592a053 137
18466557 138 while ( ( testPad = static_cast<AliMUONPad*>(next())))
139 {
ed90feae 140 if ( AliMUONPad::AreNeighbours(*testPad,*addedPad) )
18466557 141 {
142 AddPad(cluster,testPad);
143 }
144 }
145}
146
147//_____________________________________________________________________________
148Bool_t
149AreOverlapping(const AliMUONPad& pad, const AliMUONCluster& cluster)
150{
151 /// Whether the pad overlaps with the cluster
152
153 static Double_t precision = 1E-4; // cm
154 static TVector2 precisionAdjustment(precision,precision);//-precision,-precision);
155 for ( Int_t i = 0; i < cluster.Multiplicity(); ++i )
156 {
157 AliMUONPad* testPad = cluster.Pad(i);
158 // Note: we use negative precision numbers, meaning
159 // the area of the pads will be *increased* by these small numbers
160 // prior to check the overlap by the AreOverlapping method,
161 // so pads touching only by the corners will be considered as
162 // overlapping.
163 if ( AliMUONPad::AreOverlapping(*testPad,pad,precisionAdjustment) )
164 {
165 return kTRUE;
166 }
167 }
168 return kFALSE;
169}
170
171//_____________________________________________________________________________
e22f734a 172AliMUONPad*
173AliMUONPreClusterFinder::GetNextPad(Int_t cathode) const
174{
5bec63e7 175/// Return the next unused pad of given cathode, which is within fArea
176
e22f734a 177 TIter next(fPads[cathode]);
178
179 if ( !fArea.IsValid() )
180 {
181 return static_cast<AliMUONPad*>(next());
182 }
183 else
184 {
185 AliMUONPad* pad;
186 while ( ( pad = static_cast<AliMUONPad*>(next())) )
187 {
6e97fbb8 188 AliMpArea padArea(pad->X(), pad->Y(), pad->DX(), pad->DY());
e22f734a 189
190 if (fArea.Overlap(padArea)) return pad;
191
192 }
193 return 0x0;
194 }
195}
196
197//_____________________________________________________________________________
18466557 198AliMUONCluster*
57f3728e 199AliMUONPreClusterFinder::NewCluster()
200{
201 /// Create a new (empty) cluster
202 Int_t id = fClusters.GetLast()+1;
203 AliMUONCluster* cluster = new (fClusters[id]) AliMUONCluster;
204 cluster->SetUniqueID(id);
205 return cluster;
206}
207
208//_____________________________________________________________________________
209void
210AliMUONPreClusterFinder::RemoveCluster(AliMUONCluster* cluster)
211{
212 /// Remove a cluster
2103dc37 213 /// Note that we are *not* releasing the pads, so they won't be used further on
57f3728e 214 fClusters.Remove(cluster);
215 fClusters.Compress();
216}
217
218//_____________________________________________________________________________
219AliMUONCluster*
18466557 220AliMUONPreClusterFinder::NextCluster()
221{
222 /// Builds the next cluster, and returns it.
d592a053 223
18466557 224 // Start a new cluster
18466557 225
e22f734a 226 AliMUONPad* pad = GetNextPad(0);
18466557 227
57f3728e 228 AliMUONCluster* cluster(0x0);
229
18466557 230 if (!pad) // protection against no pad in first cathode, which might happen
231 {
232 // try other cathode
e22f734a 233 pad = GetNextPad(1);
18466557 234 if (!pad)
235 {
18466557 236 return 0x0;
237 }
57f3728e 238 else
239 {
240 cluster = NewCluster();
241 // Builds (recursively) a cluster on second cathode only
242 AddPad(*cluster,pad);
243 }
18466557 244 }
245 else
246 {
247 // Builds (recursively) a cluster on first cathode only
57f3728e 248
249 cluster = NewCluster();
18466557 250 AddPad(*cluster,pad);
251
861d6ce8 252 if ( !ShouldAbort() )
18466557 253 {
861d6ce8 254 // On the 2nd cathode, only add pads overlapping with the current cluster
9e41a340 255 TObjArray& padArray = *fPads[1];
861d6ce8 256 TIter next(&padArray);
257 AliMUONPad* testPad;
258
259 while ( ( testPad = static_cast<AliMUONPad*>(next())) && !ShouldAbort() )
18466557 260 {
861d6ce8 261 if (AreOverlapping(*testPad,*cluster) )
262 {
263 AddPad(*cluster,testPad);
264 }
18466557 265 }
266 }
267 }
268
861d6ce8 269 if ( ShouldAbort() )
270 {
fdee8d81 271 AliCodeTimerAuto(Form("Skipping a precluster in DE %d because it got too many pads",fDetElemId),0);
57f3728e 272 RemoveCluster(cluster);
2103dc37 273 return NextCluster();
861d6ce8 274 }
275
18466557 276 if ( cluster->Multiplicity() <= 1 )
277 {
278 if ( cluster->Multiplicity() == 0 )
279 {
280 // no pad is suspicious
281 AliWarning("Got an empty cluster...");
282 }
283 // else only 1 pad (not suspicious, but kind of useless, probably noise)
284 // so we remove it from our list
57f3728e 285 RemoveCluster(cluster);
18466557 286 // then proceed further
287 return NextCluster();
288 }
289
290 return cluster;
291}