1 #ifndef ALIMUONCLUSTER_H
2 #define ALIMUONCLUSTER_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
10 /// \class AliMUONCluster
11 /// \brief A group of adjacent pads
13 // Author Laurent Aphecetche
19 # include "TVector2.h"
25 # include "AliMpArea.h"
27 #ifndef ALI_MP_DIRECTION_H
28 # include "AliMpDirection.h"
30 #ifndef ALI_MP_ENCODE_PAIR_H
31 # include "AliMpEncodePair.h"
34 #ifndef ROOT_TObjArray
35 # include "TObjArray.h"
40 class AliMUONCluster : public TObject
44 AliMUONCluster(const AliMUONCluster& rhs);
45 AliMUONCluster& operator=(const AliMUONCluster& rhs);
47 virtual ~AliMUONCluster();
49 Bool_t Contains(const AliMUONPad& pad) const;
51 TString AsString() const;
53 static Bool_t AreOverlapping(const AliMUONCluster& c1, const AliMUONCluster& c2);
55 AliMUONPad* AddPad(const AliMUONPad& pad);
57 /// Area that contains all the pads (whatever the cathode)
58 AliMpArea Area() const;
60 /// Area that contains all the pads of a given cathode
61 AliMpArea Area(Int_t cathode) const;
63 Float_t Charge() const;
64 Float_t Charge(Int_t cathode) const;
66 /// Return the cathode's charges asymmetry
67 Float_t ChargeAsymmetry() const;
69 /// Return chi2 of the RawCharge fit (if any)
70 Float_t Chi2() const { return fChi2; }
72 virtual void Copy(TObject& obj) const;
74 /// Return false for pre-cluster
75 Bool_t HasPosition() const { return fHasPosition; }
77 /// Whether we have at least one saturated pad in a given cathode
78 Bool_t IsSaturated(Int_t cathode) const { return fIsSaturated[cathode]; }
80 /// Whether we have one saturated pad on *each* cathode
81 Bool_t IsSaturated() const { return IsSaturated(0) && IsSaturated(1); }
83 /// Return the max charge on the chathod
84 Int_t MaxChargeCathode() const { return Charge(0) > Charge(1) ? 0:1; }
86 /// Return the max raw charge on the chathod
87 Int_t MaxRawChargeCathode() const { return RawCharge(0) > RawCharge(1) ? 0:1; }
89 /// Return the biggest pad dimensions for a given cathode
90 TVector2 MaxPadDimensions(Int_t cathode, Int_t statusMask, Bool_t matchMask) const;
92 /// Return the biggest pad dimensions
93 TVector2 MaxPadDimensions(Int_t statusMask, Bool_t matchMask) const;
95 /// Return the smallest pad dimensions for a given cathode
96 TVector2 MinPadDimensions(Int_t cathode, Int_t statusMask, Bool_t matchMask) const;
98 /// Return the smallest pad dimensions
99 TVector2 MinPadDimensions(Int_t statusMask, Bool_t matchMask) const;
101 Int_t Multiplicity() const;
102 Int_t Multiplicity(Int_t cathode) const;
104 /// Compute number of pads in X and Y direction for a given cathode.
105 Long_t NofPads(Int_t cathode, Int_t statusMask, Bool_t matchMask) const;
107 /// Number of pads in (X,Y) direction, whatever the cathode.
108 Long_t NofPads(Int_t statusMask, Bool_t matchMask=kTRUE) const;
110 /// Return true as the function Compare is implemented
111 Bool_t IsSortable() const { return kTRUE; }
113 // Bool_t IsEqual(const TObject* obj) const;
115 virtual Int_t Compare(const TObject* obj) const;
117 AliMUONPad* Pad(Int_t index) const;
119 virtual void Paint(Option_t* opt="");
121 /// Return (x,y) of that cluster
122 TVector2 Position() const { return fPosition; }
123 /// Return errors on (x,y)
124 TVector2 PositionError() const { return fPositionError; }
126 virtual void Print(Option_t* opt="") const;
128 /// By default, return the average of both cathode RawCharges.
129 Float_t RawCharge() const;
131 /// Returns the RawCharge on the given cathode.
132 Float_t RawCharge(Int_t cathode) const;
134 /// Return the cathode's raw charges asymmetry
135 Float_t RawChargeAsymmetry() const;
137 void RemovePad(AliMUONPad* pad);
139 /// Set cathode (re)computed charges
140 void SetCharge(Float_t chargeCath0, Float_t chargeCath1)
141 { fHasCharge = kTRUE; fCharge[0]=chargeCath0; fCharge[1]=chargeCath1; }
143 /// Set chi2 of the RawCharge fit
144 void SetChi2(Float_t chi2) { fChi2 = chi2; }
146 /// Set (x,y) of that cluster and errors
147 void SetPosition(const TVector2& pos, const TVector2& errorOnPos)
148 { fHasPosition = kTRUE; fPosition = pos; fPositionError = errorOnPos; }
150 Int_t Cathode() const;
152 void AddCluster(const AliMUONCluster& cluster);
154 void Clear(Option_t* opt="");
160 TObjArray fPads; ///< AliMUONPad(s) composing this cluster
161 Bool_t fHasPosition; ///< false for pre-cluster (i.e. not yet computed)
162 TVector2 fPosition; ///< (x,y) of that cluster (only valid if fHasPosition is kTRUE)
163 TVector2 fPositionError; ///< errors on (x,y)
164 Int_t fMultiplicity[2]; ///< number of pads in each cathode
165 Float_t fRawCharge[2]; ///< cathode RawCharges
166 Bool_t fHasCharge; ///< false if SetCharge has not been called
167 Float_t fCharge[2]; ///< cathode (re)computed charges
168 Float_t fChi2; ///< chi2 of the RawCharge fit (if any)
169 Bool_t fIsSaturated[2]; ///< saturation status of cathodes
171 ClassDef(AliMUONCluster,3) // A cluster of AliMUONPad