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 **************************************************************************/
18 Revision 1.5 2001/11/06 17:19:41 cblume
19 Add detailed geometry and simple simulator
21 Revision 1.4 2001/10/21 18:30:02 hristov
22 Several pointers were set to zero in the default constructors to avoid memory management problems
24 Revision 1.3 2000/10/15 23:40:01 cblume
27 Revision 1.2 2000/10/06 16:49:46 cblume
30 Revision 1.1.2.2 2000/10/04 16:34:58 cblume
31 Replace include files by forward declarations
33 Revision 1.1.2.1 2000/09/22 14:47:52 cblume
38 //////////////////////////////////////////////////////////////////////
40 // Hit compression class //
41 // Adapted from AliTPCTimeBin by Marian //
43 //////////////////////////////////////////////////////////////////////
45 #include "AliTRDcluster.h"
46 #include "AliTRDtimeBin.h"
48 ClassImp(AliTRDtimeBin)
50 //______________________________________________________
52 AliTRDtimeBin::AliTRDtimeBin() {
55 for (UInt_t i=0; i<kMaxClusterPerTimeBin; i++)
58 //______________________________________________________
60 void AliTRDtimeBin::InsertCluster(AliTRDcluster* c, UInt_t index) {
62 // Insert cluster in TimeBin cluster array.
63 // Clusters are sorted according to Y coordinate.
65 if (fN==kMaxClusterPerTimeBin) {
66 printf("AliTRDtimeBin::InsertCluster(): Too many clusters !\n");
69 if (fN==0) {fIndex[0]=index; fClusters[fN++]=c; return;}
70 Int_t i=Find(c->GetY());
71 memmove(fClusters+i+1 ,fClusters+i,(fN-i)*sizeof(AliTRDcluster*));
72 memmove(fIndex +i+1 ,fIndex +i,(fN-i)*sizeof(UInt_t));
73 fIndex[i]=index; fClusters[i]=c; fN++;
76 //______________________________________________________
78 Int_t AliTRDtimeBin::Find(Double_t y) const {
80 // Returns index of the cluster nearest in Y
82 if (y <= fClusters[0]->GetY()) return 0;
83 if (y > fClusters[fN-1]->GetY()) return fN;
84 Int_t b=0, e=fN-1, m=(b+e)/2;
85 for (; b<e; m=(b+e)/2) {
86 if (y > fClusters[m]->GetY()) b=m+1;
92 //______________________________________________________
93 AliTRDcluster *AliTRDtimeBin::operator[](Int_t i)