Implemented a new version of cluster (with its store and iterator):
[u/mrichter/AliRoot.git] / MUON / AliMUONRawCluster.h
1 #ifndef ALIMUONRAWCLUSTER_H
2 #define ALIMUONRAWCLUSTER_H
3
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 /* $Id$ */
8 // Revision of includes 07/05/2004
9
10 /// \ingroup base
11 /// \class AliMUONRawCluster
12 /// \brief MUON raw cluster
13 ///
14 /// Class for the MUON RecPoint
15 /// It contains the properties of the physics cluters found in the tracking chambers
16 /// RawCluster contains also the information from the both cathode of the chambers.
17
18
19 #include <TObject.h>
20 #include <TMath.h> // because of inline funtion GetRadius
21 #include <TArrayF.h>
22 class AliMUONRawCluster : public TObject {
23
24 public:
25    AliMUONRawCluster();
26    virtual ~AliMUONRawCluster();
27    
28    /// Return radius
29    Float_t      GetRadius(Int_t i) {return TMath::Sqrt(fX[i]*fX[i]+fY[i]*fY[i]);}
30    /// Return true as the function Compare() is implemented
31    Bool_t       IsSortable() const {return kTRUE;}
32    Int_t        Compare(const TObject *obj) const;
33    Int_t        PhysicsContribution() const;
34    virtual void Print(Option_t* opt="") const;
35    static Int_t BinarySearch(Float_t r, TArrayF ccord, Int_t from, Int_t upto);
36    static void  SortMin(Int_t *idx,Float_t *xdarray, Float_t *xarray, Float_t *yarray, Float_t *qarray,Int_t ntr);
37    void         DumpIndex();
38
39    Int_t        AddCharge(Int_t i, Float_t Q);
40    Int_t        AddX(Int_t i, Float_t X);
41    Int_t        AddY(Int_t i, Float_t Y);
42    Int_t        AddZ(Int_t i, Float_t Z);
43
44    /// Return the ID number of the detection element (slat) on which the cluster is found
45    Int_t DetElemId() const { return fDetElemId; }
46    
47    Float_t        GetCharge(Int_t i=0) const;
48    Float_t      GetX(Int_t i=0) const;
49    Float_t      GetY(Int_t i=0) const;
50    Float_t      GetZ(Int_t i=0) const;
51    Int_t        GetTrack(Int_t i=0) const;
52    Float_t        GetPeakSignal(Int_t i=0) const;
53    Int_t        GetMultiplicity(Int_t i=0) const;
54    Int_t        GetClusterType() const;
55    Int_t        GetGhost() const;
56    Int_t        GetNcluster(Int_t i=0) const;
57    Float_t      GetChi2(Int_t i=0) const;
58    Int_t        GetIndex(Int_t i, Int_t j) const;
59    Int_t        GetOffset(Int_t i, Int_t j) const;
60    Float_t      GetContrib(Int_t i, Int_t j) const;
61    Int_t        GetPhysics(Int_t i) const;
62    Int_t        GetDetElemId() const ; 
63    Float_t      GetError(Int_t iXY) const;
64    Float_t      GetErrX() const;
65    Float_t      GetErrY() const;
66
67    Int_t        SetCharge(Int_t i, Float_t Q);
68    Int_t        SetX(Int_t i, Float_t X);
69    Int_t        SetY(Int_t i, Float_t Y);
70    Int_t        SetZ(Int_t i, Float_t Z);
71    void         SetDetElemId(Int_t Id); 
72    Int_t        SetTrack(Int_t i, Int_t track);
73    Int_t        SetPeakSignal(Int_t i, Float_t peaksignal);
74    Int_t        SetMultiplicity(Int_t i, Int_t mul);
75    Int_t        SetClusterType(Int_t type);
76    Int_t        SetGhost(Int_t ghost);
77    Int_t        SetNcluster(Int_t i, Int_t ncluster);
78    Int_t        SetChi2(Int_t i, Float_t chi2);
79    void         SetIndex(Int_t i, Int_t j, Int_t index);
80    void         SetOffset(Int_t i, Int_t j, Int_t offset);
81    void         SetContrib(Int_t i, Int_t j, Float_t contrib);
82    void         SetPhysics(Int_t i, Int_t physics);
83    void         SetError(Int_t iXY, Float_t err);
84    void         SetErrX(Float_t err);
85    void         SetErrY(Float_t err);
86
87 private:
88    Int_t       fIndexMap[50][2];  ///< Indices of digits
89    Int_t       fOffsetMap[50][2]; ///< Emmanuel special
90    Float_t     fContMap[50][2];   ///< Contribution from digit
91    Int_t     fPhysicsMap[50];   ///< Distinguish signal and background contr.
92   
93    Float_t     fQ[2]  ;           ///< Q of cluster (in ADC counts)     
94    Float_t     fX[2]  ;           ///< X of cluster
95    Float_t     fY[2]  ;           ///< Y of cluster
96    Float_t     fZ[2]  ;           ///< Z of cluster
97    Int_t       fTracks[3];        ///< Labels of overlapped tracks
98    Float_t     fPeakSignal[2];    ///< Peak signal 
99    Int_t       fMultiplicity[2];  ///< Cluster multiplicity
100    Int_t       fClusterType;      ///< Cluster type
101    Int_t       fGhost;            ///< Ghost info
102                                   // 0 if not a ghost or ghost problem solved
103                                   // >0 if ghost problem remains because
104                                   // 1 both (true and ghost) satify 
105                                   //   charge chi2 compatibility
106                                   // 2 none give satisfactory chi2
107    Int_t       fNcluster[2];      ///< Number of clusters
108    Float_t     fChi2[2];          ///< Chi**2 of fit
109    Int_t       fDetElemId;        ///< ID number of the detection element (slat) on which the cluster is found. 
110    Float_t     fErrXY[2];         ///< coordinate errors
111    ClassDef(AliMUONRawCluster,2)  //Cluster class for MUON
112 };
113
114 // inline functions
115
116 /// Return Indices of digits
117 inline  Int_t  AliMUONRawCluster::GetIndex(Int_t i, Int_t j) const
118 { return fIndexMap[i][j]; }
119
120 /// Return Emmanuel special offset map
121 inline  Int_t  AliMUONRawCluster::GetOffset(Int_t i, Int_t j) const
122 { return fOffsetMap[i][j]; }
123
124 /// Return Contribution from digit
125 inline  Float_t  AliMUONRawCluster::GetContrib(Int_t i, Int_t j) const
126 { return fContMap[i][j]; }
127
128 /// Return Distinguish signal and background contr.
129 inline  Int_t  AliMUONRawCluster::GetPhysics(Int_t i) const
130 { return fPhysicsMap[i]; }
131
132 /// Set Indices of digits
133 inline  void  AliMUONRawCluster::SetIndex(Int_t i, Int_t j, Int_t index)
134 { fIndexMap[i][j] = index; }
135
136 /// Set Emmanuel special offset map
137 inline  void  AliMUONRawCluster::SetOffset(Int_t i, Int_t j, Int_t offset)
138 { fOffsetMap[i][j] = offset; }
139
140 /// Set Contribution from digit
141 inline  void  AliMUONRawCluster::SetContrib(Int_t i, Int_t j, Float_t contrib)
142 { fContMap[i][j] = contrib; }
143
144 /// Set Distinguish signal and background contr.
145 inline  void  AliMUONRawCluster::SetPhysics(Int_t i, Int_t physics)
146 { fPhysicsMap[i] = physics; }
147
148 /// Set ID number of the detection element (slat) on which the cluster is found.
149 inline void AliMUONRawCluster::SetDetElemId(Int_t Id)
150 { fDetElemId = Id; }
151
152 /// Return ID number of the detection element (slat) on which the cluster is found.
153 inline Int_t AliMUONRawCluster::GetDetElemId() const
154 { return fDetElemId;}
155
156 /// Set coordinate errors
157 inline void AliMUONRawCluster::SetError(Int_t iXY, Float_t err)
158 { fErrXY[iXY] = err; }
159
160 /// Set x coordinate error
161 inline void AliMUONRawCluster::SetErrX(Float_t err)
162 { SetError(0, err); }
163
164 /// Set y coordinate error
165 inline void AliMUONRawCluster::SetErrY(Float_t err)
166 { SetError(1, err); }
167
168 /// Return coordinate errors
169 inline Float_t AliMUONRawCluster::GetError(Int_t iXY) const
170 { return fErrXY[iXY]; }
171
172 /// Return x coordinate error
173 inline Float_t AliMUONRawCluster::GetErrX() const
174 { return GetError(0); }
175
176 /// Return y coordinate error
177 inline Float_t AliMUONRawCluster::GetErrY() const
178 { return GetError(1); }
179
180 #endif
181
182
183
184
185
186