1 #ifndef ALIMUONTRACKERDATA_H
2 #define ALIMUONTRACKERDATA_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
10 /// \class AliMUONTrackerData
11 /// \brief Implementation of AliMUONVTrackerData
13 // Author Laurent Aphecetche, Subatech
15 #ifndef ALIMUONVTRACKERDATA_H
16 # include "AliMUONVTrackerData.h"
19 class AliMUONSparseHisto;
20 class AliMUONVCalibParam;
22 class AliMpDetElement;
24 class AliMUONTrackerData : public AliMUONVTrackerData
27 AliMUONTrackerData(const char* name="", const char* title="",
29 Bool_t issingleevent=kFALSE);
31 AliMUONTrackerData(const char* name, const char* title,
32 const AliMUONVStore& manuValues);
34 AliMUONTrackerData(const char* name, const char* title,
35 const AliMUONVStore& deValues, Int_t val);
37 virtual ~AliMUONTrackerData();
39 Bool_t Add(const AliMUONTrackerData& data);
41 virtual Bool_t Add(const AliMUONVStore& channelValues, TArrayI* nofEventsPerDDL=0x0);
43 virtual Bool_t Replace(const AliMUONVStore& channelValues);
45 virtual Double_t BusPatch(Int_t busPatchId, Int_t dim=0) const;
47 virtual Double_t Chamber(Int_t chamberId, Int_t dim=0) const;
49 virtual Double_t Channel(Int_t detElemId, Int_t manuId, Int_t manuChannel,
52 virtual void Clear(Option_t* opt="");
54 virtual Double_t Count(Int_t detElemId, Int_t manuId, Int_t manuChannel) const;
56 virtual Double_t DetectionElement(Int_t detElemId, Int_t dim=0) const;
58 virtual TString DimensionName(Int_t dim) const;
60 virtual TString ExternalDimensionName(Int_t dim) const;
62 virtual Bool_t HasChamber(Int_t chamberId) const;
64 virtual Bool_t HasBusPatch(Int_t busPatchId) const;
66 virtual Bool_t HasDetectionElement(Int_t detElemId) const;
68 virtual Bool_t HasManu(Int_t detElemId, Int_t manuId) const;
70 virtual Bool_t HasPCB(Int_t detElemId, Int_t pcbIndex) const;
72 /// Whether we can be run
73 virtual Bool_t IsSingleEvent() const { return fIsSingleEvent; }
75 virtual Double_t Manu(Int_t detElemId, Int_t manuId, Int_t dim=0) const;
77 /// Returns the number of dimensions (i.e. the number of values) each element has
78 virtual Int_t NumberOfDimensions() const;
80 /// The number of values we are inputting
81 virtual Int_t ExternalDimension() const { return fExternalDimension; }
83 /// Convert from internal to external dimension
84 virtual Int_t InternalToExternal(Int_t dim) const { return dim/2; }
86 /// Returns the number of events we have seen so far
87 virtual Int_t NumberOfEvents(Int_t ddlNumber) const;
89 virtual Double_t PCB(Int_t detElemId, Int_t pcbIndex, Int_t dim=0) const;
93 /// Print, with option, all objects whose name matches wildcard
94 virtual void Print(Option_t* wildcard, Option_t* opt) const;
96 virtual void SetDimensionName(Int_t index, const char* value);
98 Bool_t CanHistogram() const { return kTRUE; }
100 void MakeHistogramForDimension(Int_t index, Bool_t value, Double_t xmin=0.0, Double_t xmax=4096.0);
102 virtual void HistogramRange(Double_t& xmin, Double_t& xmax) const { xmin = fXmin; xmax = fXmax; }
104 AliMUONSparseHisto* GetManuSparseHisto(Int_t detElemId, Int_t manuId,
107 AliMUONSparseHisto* GetManuSparseHisto(Int_t detElemId, Int_t manuId,
110 AliMUONSparseHisto* GetChannelSparseHisto(Int_t detElemId, Int_t manuId,
111 Int_t manuChannel, Int_t dim=0);
113 virtual AliMUONSparseHisto* GetChannelSparseHisto(Int_t detElemId, Int_t manuId,
114 Int_t manuChannel, Int_t dim=0) const;
116 /// Disable storing values at the channel level
117 virtual void DisableChannelLevel();
119 /// Whether we store values at the channel level or not
120 virtual Bool_t IsChannelLevelEnabled() const { return fIsChannelLevelEnabled; }
122 /// Disable storing values at the manu level
123 virtual void DisableManuLevel();
125 /// Whether we store values at the manu level or not
126 virtual Bool_t IsManuLevelEnabled() const { return fIsManuLevelEnabled; }
128 /// Whether we store values at the bus patch level or not
129 virtual Bool_t IsBusPatchLevelEnabled() const { return fIsBustPatchLevelEnabled; }
131 /// Whether we store values at the PCB level or not
132 virtual Bool_t IsPCBLevelEnabled() const { return fIsPCBLevelEnabled; }
134 /// To allow merging of different objects
135 virtual Long64_t Merge(TCollection* list);
137 Bool_t ExportAsASCIIOccupancyFile(const char* filename, Int_t runNumber) const;
141 void FillHisto(Int_t detElemId, Int_t manuId, Int_t manuChannel,
142 Int_t dim, Double_t value);
144 AliMUONVCalibParam* BusPatchParam(Int_t busPatch, Bool_t create=kFALSE) const;
146 AliMUONVCalibParam* CreateBusPatchParam(Int_t busPatch) const;
148 AliMUONVCalibParam* ChamberParam(Int_t chamberId, Bool_t create=kFALSE) const;
150 AliMUONVCalibParam* CreateChamberParam(Int_t chamberId) const;
152 AliMUONVCalibParam* ChannelParam(Int_t detElemId, Int_t manuId,
153 const AliMUONVCalibParam* external=0x0) const;
155 AliMUONVCalibParam* DetectionElementParam(Int_t detElemId, Bool_t create=kFALSE) const;
157 AliMUONVCalibParam* CreateDetectionElementParam(Int_t detElemId) const;
159 AliMUONVCalibParam* ManuParam(Int_t detElemId, Int_t manuId, Bool_t create=kFALSE) const;
161 AliMUONVCalibParam* CreateManuParam(Int_t detElemInd, Int_t manuId) const;
163 AliMUONVCalibParam* PCBParam(Int_t detElemId, Int_t pcbIndex, Bool_t create=kFALSE) const;
165 AliMUONVCalibParam* CreatePCBParam(Int_t detElemId, Int_t pcbIndex) const;
167 /// Index of the dimension containing the number of time an item was hit
168 virtual Int_t IndexOfNumberDimension() const { return fDimension - 1; }
170 /// Index of the dimension containing the occupancy number
171 virtual Int_t IndexOfOccupancyDimension() const { return fDimension - 2; }
173 /// Whether we have histograms for a given dimension, or not
174 virtual Bool_t IsHistogrammed(Int_t dim) const { return ( fHistogramming[dim] > 0 ); }
176 Int_t DdlIdFromBusPatchId(Int_t buspatchid) const;
177 Int_t DdlIdFromDetElemId(Int_t detelemid) const;
178 Int_t DdlIdFromChamberId(Int_t chamberid) const;
181 AliMUONTrackerData(const AliMUONTrackerData& rhs);
183 AliMUONTrackerData& operator=(const AliMUONTrackerData& rhs);
185 AliMUONVCalibParam* CreateDouble(const AliMUONVCalibParam& param, Int_t detElemId, Int_t manuId) const;
187 Int_t GetParts(AliMUONVCalibParam* external,
188 AliMUONVCalibParam*& chamber,
189 AliMUONVCalibParam*& de,
190 AliMUONVCalibParam*& busPatch,
191 AliMUONVCalibParam*& pcb,
192 AliMUONVCalibParam*& manu,
193 AliMUONVCalibParam*& channel,
194 AliMpDetElement*& mpde);
196 /// Convert from external to internal index
197 Int_t External2Internal(Int_t index) const;
199 void SetInternalDimensionName(Int_t index, const char* value);
201 void SetExternalDimensionName(Int_t index, const char* value);
203 Double_t Value(const AliMUONVCalibParam& param, Int_t i, Int_t dim, Int_t ddlId) const;
205 /// The number of values we actually *store* for each item
206 Int_t Dimension() const { return fDimension; }
208 Bool_t InternalAdd(const AliMUONVStore& store, TArrayI* nevents, Bool_t replace);
210 void GetDEManu(const AliMUONVCalibParam& param,
211 Int_t& detElemId, Int_t& manuId) const;
213 void AddCalibParams(const AliMUONVCalibParam& src, AliMUONVCalibParam& dest) const;
215 void Add2D(const AliMUONVStore& src, AliMUONVStore& dest) const;
217 void Add1D(const AliMUONVStore& src, AliMUONVStore& dest) const;
221 Bool_t UpdateNumberOfEvents(TArrayI* nevents);
225 Bool_t fIsSingleEvent; ///< whether we can deal with more than one event
226 AliMUONVStore* fChannelValues; ///< the channel store
227 AliMUONVStore* fManuValues; ///< the manu store
228 AliMUONVStore* fBusPatchValues; ///< the bus patch store
229 AliMUONVStore* fDEValues; ///< the detection element store
230 AliMUONVStore* fChamberValues; ///< the chamber store
231 AliMUONVStore* fPCBValues; ///< the pcb store
232 Int_t fDimension; ///< the dimension of the data
233 Int_t fNevents; ///< the number of events treated
234 TObjArray* fDimensionNames; ///< the names of the (internal) dimensions
235 TObjArray* fExternalDimensionNames; ///< the names of the external (i.e. original) dimensions
236 Int_t fExternalDimension; ///< number of interface values per item
237 /// whether we should histogram the dimension(s)
238 Int_t* fHistogramming; //[fExternalDimension] whether we should histogram the dimension(s)
239 AliMUONVStore* fHistos; ///< the lowest histograms we have
240 Double_t fXmin; ///< min x value for histograms
241 Double_t fXmax; ///< max x value for histograms
242 static const Int_t fgkExtraDimension; ///< to hold extra information
243 static const Int_t fgkVirtualExtraDimension; ///< to give access to information not stored, but computed on the fly
244 Bool_t fIsChannelLevelEnabled; ///< whether we allow storing of channel (fChannelValues) values
245 Bool_t fIsManuLevelEnabled; ///< whether we allow storing of manu (fManuValues) values
246 Bool_t fIsBustPatchLevelEnabled; ///< whether we allow storing of bus patches (fBusPatchValues) values
247 Bool_t fIsPCBLevelEnabled; ///< whether we allow storing of PCB values (fPCBValues)
248 Int_t fNofDDLs; ///< nof of DDLs we're dealing with
249 /// the number of events treated (per DDL)
250 Int_t* fNofEventsPerDDL; //[fNofDDLs] the number of events treated (per DDL)
252 ClassDef(AliMUONTrackerData,8) // Implementation of AliMUONVTrackerData