]>
Commit | Line | Data |
---|---|---|
a9e2aefa | 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$ */ | |
30178c30 | 8 | // Revision of includes 07/05/2004 |
a9e2aefa | 9 | |
692de412 | 10 | /// \ingroup base |
11 | /// \class AliMUONRawCluster | |
12 | /// \brief MUON raw cluster | |
13 | /// | |
14 | /// Class for the MUON RecPoint | |
585d648a | 15 | /// It contains the properties of the physics cluters found in the tracking chambers |
692de412 | 16 | /// RawCluster contains also the information from the both cathode of the chambers. |
3afdba21 | 17 | |
18 | ||
2060b217 | 19 | #include "AliMUONVCluster.h" |
ecfa008b | 20 | #include <TMath.h> // because of inline funtion GetRadius |
f480df09 | 21 | #include <TArrayF.h> |
2060b217 | 22 | |
23 | class AliMUONRawCluster : public AliMUONVCluster { | |
3afdba21 | 24 | |
a9e2aefa | 25 | public: |
3afdba21 | 26 | AliMUONRawCluster(); |
71a2d3aa | 27 | virtual ~AliMUONRawCluster(); |
28 | ||
2060b217 | 29 | /// Clear method (used by TClonesArray) |
30 | virtual void Clear(Option_t* = "") {} | |
31 | ||
96ebe67e | 32 | /// Create a copy of the current cluster |
1467f4ba | 33 | virtual AliMUONRawCluster* Clone(const char* = "") const {return new AliMUONRawCluster(*this);} |
96ebe67e | 34 | |
35 | /// Set coordinates (cm) | |
2060b217 | 36 | virtual void SetXYZ(Double_t x, Double_t y, Double_t z) {fX[0] = x; fY[0] = y; fZ[0] = z;} |
37 | /// Return coordinate X (cm) | |
38 | virtual Double_t GetX() const {return fX[0];} | |
39 | /// Return coordinate Y (cm) | |
40 | virtual Double_t GetY() const {return fY[0];} | |
41 | /// Return coordinate Z (cm) | |
42 | virtual Double_t GetZ() const {return fZ[0];} | |
43 | ||
44 | /// Set resolution (cm) on coordinates (X,Y) | |
45 | virtual void SetErrXY(Double_t errX, Double_t errY) {fErrXY[0] = errX; fErrXY[1] = errY;} | |
46 | /// Return resolution (cm) on coordinate X | |
47 | virtual Double_t GetErrX() const {return fErrXY[0];} | |
48 | /// Return resolution**2 (cm**2) on coordinate X | |
49 | virtual Double_t GetErrX2() const {return fErrXY[0] * fErrXY[0];} | |
50 | /// Return resolution (cm) on coordinate Y | |
51 | virtual Double_t GetErrY() const {return fErrXY[1];} | |
52 | /// Return resolution**2 (cm**2) on coordinate Y | |
53 | virtual Double_t GetErrY2() const {return fErrXY[1] * fErrXY[1];} | |
54 | ||
55 | /// Set the cluster charge | |
56 | virtual void SetCharge(Double_t q) {fQ[0] = q;} | |
57 | /// Set the cluster charge | |
58 | virtual Double_t GetCharge() const {return fQ[0];} | |
59 | ||
60 | /// Return chamber Id | |
61 | virtual Int_t GetChamberId() const {return fDetElemId/100 - 1;} | |
62 | /// Set detection element Id | |
63 | void SetDetElemId(Int_t id) {fDetElemId = id;} | |
64 | /// Return detection element Id | |
65 | virtual Int_t GetDetElemId() const {return fDetElemId;} | |
66 | ||
67 | virtual void SetDigitsId(Int_t nDigits, const UInt_t *digitsId); | |
68 | /// Add a digit Id to the array of associated digits | |
69 | virtual void AddDigitId(UInt_t id) {fIndexMap[fMultiplicity[0]++][0] = id;} | |
70 | ||
71 | /// Return number of associated digits | |
72 | virtual Int_t GetNDigits() const {return fMultiplicity[0];} | |
73 | /// Return Id of digits i | |
74 | virtual UInt_t GetDigitId(Int_t i) const {return (i < fMultiplicity[0] && i < 50) ? (UInt_t)fIndexMap[i][0] : 0;} | |
75 | ||
76 | /// Set chi2 of cluster | |
77 | virtual void SetChi2( Double_t chi2) {fChi2[0] = chi2;} | |
78 | /// Return chi2 of cluster | |
79 | virtual Double_t GetChi2() const {return fChi2[0];} | |
80 | ||
2e2d0c44 | 81 | /// Set the corresponding MC track number |
82 | virtual void SetMCLabel(Int_t label) {SetTrack(0, label);} | |
83 | /// Return the corresponding MC track number | |
84 | virtual Int_t GetMCLabel() const {return GetTrack(0);} | |
85 | ||
86 | /// Return radius | |
9ee1d6ff | 87 | Float_t GetRadius(Int_t i) const {return TMath::Sqrt(fX[i]*fX[i]+fY[i]*fY[i]);} |
71a2d3aa | 88 | /// Return true as the function Compare() is implemented |
3afdba21 | 89 | Bool_t IsSortable() const {return kTRUE;} |
90 | Int_t Compare(const TObject *obj) const; | |
91 | Int_t PhysicsContribution() const; | |
6464217e | 92 | virtual void Print(Option_t* opt="") const; |
3afdba21 | 93 | static Int_t BinarySearch(Float_t r, TArrayF ccord, Int_t from, Int_t upto); |
94 | static void SortMin(Int_t *idx,Float_t *xdarray, Float_t *xarray, Float_t *yarray, Float_t *qarray,Int_t ntr); | |
95 | void DumpIndex(); | |
a9e2aefa | 96 | |
a31c3507 | 97 | Int_t AddCharge(Int_t i, Float_t Q); |
3afdba21 | 98 | Int_t AddX(Int_t i, Float_t X); |
99 | Int_t AddY(Int_t i, Float_t Y); | |
100 | Int_t AddZ(Int_t i, Float_t Z); | |
101 | ||
2060b217 | 102 | Float_t GetCharge(Int_t i) const; |
103 | Float_t GetX(Int_t i) const; | |
104 | Float_t GetY(Int_t i) const; | |
105 | Float_t GetZ(Int_t i) const; | |
6464217e | 106 | Int_t GetTrack(Int_t i=0) const; |
2060b217 | 107 | Float_t GetPeakSignal(Int_t i=0) const; |
6464217e | 108 | Int_t GetMultiplicity(Int_t i=0) const; |
9e993f2a | 109 | Int_t GetClusterType() const; |
3b5272e3 | 110 | Int_t GetGhost() const; |
6464217e | 111 | Int_t GetNcluster(Int_t i=0) const; |
2060b217 | 112 | Float_t GetChi2(Int_t i) const; |
0164904a | 113 | Int_t GetIndex(Int_t i, Int_t j) const; |
114 | Int_t GetOffset(Int_t i, Int_t j) const; | |
115 | Float_t GetContrib(Int_t i, Int_t j) const; | |
116 | Int_t GetPhysics(Int_t i) const; | |
3afdba21 | 117 | |
a31c3507 | 118 | Int_t SetCharge(Int_t i, Float_t Q); |
3afdba21 | 119 | Int_t SetX(Int_t i, Float_t X); |
120 | Int_t SetY(Int_t i, Float_t Y); | |
121 | Int_t SetZ(Int_t i, Float_t Z); | |
9e993f2a | 122 | Int_t SetTrack(Int_t i, Int_t track); |
a31c3507 | 123 | Int_t SetPeakSignal(Int_t i, Float_t peaksignal); |
9e993f2a | 124 | Int_t SetMultiplicity(Int_t i, Int_t mul); |
125 | Int_t SetClusterType(Int_t type); | |
3b5272e3 | 126 | Int_t SetGhost(Int_t ghost); |
127 | Int_t SetNcluster(Int_t i, Int_t ncluster); | |
128 | Int_t SetChi2(Int_t i, Float_t chi2); | |
0164904a | 129 | void SetIndex(Int_t i, Int_t j, Int_t index); |
130 | void SetOffset(Int_t i, Int_t j, Int_t offset); | |
131 | void SetContrib(Int_t i, Int_t j, Float_t contrib); | |
132 | void SetPhysics(Int_t i, Int_t physics); | |
3afdba21 | 133 | |
0164904a | 134 | private: |
18b6b8c7 | 135 | Int_t fIndexMap[50][2]; ///< Indices of digits |
136 | Int_t fOffsetMap[50][2]; ///< Emmanuel special | |
137 | Float_t fContMap[50][2]; ///< Contribution from digit | |
2060b217 | 138 | Int_t fPhysicsMap[50]; ///< Distinguish signal and background contr. |
3b5272e3 | 139 | |
a31c3507 | 140 | Float_t fQ[2] ; ///< Q of cluster (in ADC counts) |
18b6b8c7 | 141 | Float_t fX[2] ; ///< X of cluster |
142 | Float_t fY[2] ; ///< Y of cluster | |
143 | Float_t fZ[2] ; ///< Z of cluster | |
144 | Int_t fTracks[3]; ///< Labels of overlapped tracks | |
a31c3507 | 145 | Float_t fPeakSignal[2]; ///< Peak signal |
18b6b8c7 | 146 | Int_t fMultiplicity[2]; ///< Cluster multiplicity |
147 | Int_t fClusterType; ///< Cluster type | |
148 | Int_t fGhost; ///< Ghost info | |
149 | // 0 if not a ghost or ghost problem solved | |
3b5272e3 | 150 | // >0 if ghost problem remains because |
151 | // 1 both (true and ghost) satify | |
152 | // charge chi2 compatibility | |
153 | // 2 none give satisfactory chi2 | |
18b6b8c7 | 154 | Int_t fNcluster[2]; ///< Number of clusters |
155 | Float_t fChi2[2]; ///< Chi**2 of fit | |
156 | Int_t fDetElemId; ///< ID number of the detection element (slat) on which the cluster is found. | |
157 | Float_t fErrXY[2]; ///< coordinate errors | |
2060b217 | 158 | |
159 | ClassDef(AliMUONRawCluster,3) //Cluster class for MUON | |
a9e2aefa | 160 | }; |
0164904a | 161 | |
162 | // inline functions | |
163 | ||
71a2d3aa | 164 | /// Return Indices of digits |
0164904a | 165 | inline Int_t AliMUONRawCluster::GetIndex(Int_t i, Int_t j) const |
166 | { return fIndexMap[i][j]; } | |
167 | ||
71a2d3aa | 168 | /// Return Emmanuel special offset map |
0164904a | 169 | inline Int_t AliMUONRawCluster::GetOffset(Int_t i, Int_t j) const |
170 | { return fOffsetMap[i][j]; } | |
171 | ||
71a2d3aa | 172 | /// Return Contribution from digit |
0164904a | 173 | inline Float_t AliMUONRawCluster::GetContrib(Int_t i, Int_t j) const |
174 | { return fContMap[i][j]; } | |
175 | ||
71a2d3aa | 176 | /// Return Distinguish signal and background contr. |
0164904a | 177 | inline Int_t AliMUONRawCluster::GetPhysics(Int_t i) const |
178 | { return fPhysicsMap[i]; } | |
179 | ||
71a2d3aa | 180 | /// Set Indices of digits |
0164904a | 181 | inline void AliMUONRawCluster::SetIndex(Int_t i, Int_t j, Int_t index) |
182 | { fIndexMap[i][j] = index; } | |
183 | ||
71a2d3aa | 184 | /// Set Emmanuel special offset map |
0164904a | 185 | inline void AliMUONRawCluster::SetOffset(Int_t i, Int_t j, Int_t offset) |
186 | { fOffsetMap[i][j] = offset; } | |
187 | ||
71a2d3aa | 188 | /// Set Contribution from digit |
0164904a | 189 | inline void AliMUONRawCluster::SetContrib(Int_t i, Int_t j, Float_t contrib) |
190 | { fContMap[i][j] = contrib; } | |
191 | ||
71a2d3aa | 192 | /// Set Distinguish signal and background contr. |
0164904a | 193 | inline void AliMUONRawCluster::SetPhysics(Int_t i, Int_t physics) |
194 | { fPhysicsMap[i] = physics; } | |
195 | ||
87e6e2dd | 196 | |
a9e2aefa | 197 | #endif |
198 |