MUON/AliMUONTrackerData.h

   1 #ifndef ALIMUONTRACKERDATA_H
   2 #define ALIMUONTRACKERDATA_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
   9 /// \ingroup calib
  10 /// \class AliMUONTrackerData
  11 /// \brief Implementation of AliMUONVTrackerData
  12 ///
  13 // Author Laurent Aphecetche, Subatech
  14
  15 #ifndef ALIMUONVTRACKERDATA_H
  16 #  include "AliMUONVTrackerData.h"
  17 #endif
  18
  19 class AliMUONSparseHisto;
  20 class AliMUONVCalibParam;
  21 class AliMUONVStore;
  22 class AliMpDetElement;
  23 class AliMUONRejectList;
  24
  25 class AliMUONTrackerData : public AliMUONVTrackerData
  26 {
  27 public:
  28   AliMUONTrackerData(const char* name="", const char* title="",
  29                      Int_t dimension=0,
  30                      Bool_t issingleevent=kFALSE);
  31
  32   AliMUONTrackerData(const char* name, const char* title,
  33                      const AliMUONVStore& manuValues);
  34
  35   AliMUONTrackerData(const char* name, const char* title,
  36                      const AliMUONVStore& deOrBpValues, Int_t val);
  37
  38   AliMUONTrackerData(const char* name, const char* title,
  39                      const AliMUONRejectList& rejectList);
  40
  41   virtual ~AliMUONTrackerData();
  42
  43   Bool_t Add(const AliMUONTrackerData& data);
  44
  45   virtual Bool_t Add(const AliMUONVStore& channelValues, TArrayI* nofEventsPerDDL=0x0);
  46
  47   virtual Bool_t Replace(const AliMUONVStore& channelValues);
  48
  49   virtual Double_t BusPatch(Int_t busPatchId, Int_t dim=0) const;
  50
  51   virtual Double_t Chamber(Int_t chamberId, Int_t dim=0) const;
  52
  53   virtual Double_t Channel(Int_t detElemId, Int_t manuId, Int_t manuChannel,
  54                            Int_t dim=0) const;
  55
  56   virtual void Clear(Option_t* opt="");
  57
  58   virtual Double_t Count(Int_t detElemId, Int_t manuId, Int_t manuChannel) const;
  59
  60   virtual Double_t DetectionElement(Int_t detElemId, Int_t dim=0) const;
  61
  62   virtual TString DimensionName(Int_t dim) const;
  63
  64   virtual TString ExternalDimensionName(Int_t dim) const;
  65
  66   virtual Bool_t HasChamber(Int_t chamberId) const;
  67
  68   virtual Bool_t HasBusPatch(Int_t busPatchId) const;
  69
  70   virtual Bool_t HasDetectionElement(Int_t detElemId) const;
  71
  72   virtual Bool_t HasManu(Int_t detElemId, Int_t manuId) const;
  73
  74   virtual Bool_t HasPCB(Int_t detElemId, Int_t pcbIndex) const;
  75
  76   /// Whether we can be run
  77   virtual Bool_t IsSingleEvent() const { return fIsSingleEvent; }
  78
  79   virtual Double_t Manu(Int_t detElemId, Int_t manuId, Int_t dim=0) const;
  80
  81   /// Returns the number of dimensions (i.e. the number of values) each element has
  82   virtual Int_t NumberOfDimensions() const;
  83
  84   /// The number of values we are inputting
  85   virtual Int_t ExternalDimension() const { return fExternalDimension; }
  86
  87   /// Convert from internal to external dimension
  88   virtual Int_t InternalToExternal(Int_t dim) const { return dim/2; }
  89
  90   /// Returns the number of events we have seen so far
  91   virtual Int_t NumberOfEvents(Int_t ddlNumber) const;
  92
  93   virtual Double_t PCB(Int_t detElemId, Int_t pcbIndex, Int_t dim=0) const;
  94
  95   using TObject::Print;
  96
  97   /// Print, with option, all objects whose name matches wildcard
  98   virtual void Print(Option_t* wildcard, Option_t* opt) const;
  99
 100   virtual void SetDimensionName(Int_t index, const char* value);
 101
 102   Bool_t CanHistogram() const { return kTRUE; }
 103
 104   void MakeHistogramForDimension(Int_t index, Bool_t value, Double_t xmin=0.0, Double_t xmax=4096.0);
 105
 106   virtual void HistogramRange(Double_t& xmin, Double_t& xmax) const { xmin = fXmin; xmax = fXmax; }
 107
 108   AliMUONSparseHisto* GetManuSparseHisto(Int_t detElemId, Int_t manuId,
 109                                          Int_t dim=0);
 110
 111   AliMUONSparseHisto* GetManuSparseHisto(Int_t detElemId, Int_t manuId,
 112                                          Int_t dim=0) const;
 113
 114   AliMUONSparseHisto* GetChannelSparseHisto(Int_t detElemId, Int_t manuId,
 115                                             Int_t manuChannel, Int_t dim=0);
 116
 117   virtual AliMUONSparseHisto* GetChannelSparseHisto(Int_t detElemId, Int_t manuId,
 118                                                     Int_t manuChannel, Int_t dim=0) const;
 119
 120         /// Disable storing values at the channel level
 121         virtual void DisableChannelLevel();
 122
 123         /// Whether we store values at the channel level or not
 124         virtual Bool_t IsChannelLevelEnabled() const { return fIsChannelLevelEnabled; }
 125
 126   /// Disable storing values at the manu level
 127         virtual void DisableManuLevel();
 128
 129         /// Whether we store values at the manu level or not
 130         virtual Bool_t IsManuLevelEnabled() const { return fIsManuLevelEnabled; }
 131
 132   /// Whether we store values at the bus patch level or not
 133   virtual Bool_t IsBusPatchLevelEnabled() const { return fIsBustPatchLevelEnabled; }
 134
 135   /// Whether we store values at the PCB level or not
 136   virtual Bool_t IsPCBLevelEnabled() const { return fIsPCBLevelEnabled; }
 137
 138   /// To allow merging of different objects
 139   virtual Long64_t Merge(TCollection* list);
 140
 141   Bool_t ExportAsASCIIOccupancyFile(const char* filename, Int_t runNumber) const;
 142
 143 private:
 144
 145   void BuildFromDEStore(const AliMUONVStore& deStore);
 146
 147   void BuildFromBPStore(const AliMUONVStore& bpStore);
 148
 149   void FillHisto(Int_t detElemId, Int_t manuId, Int_t manuChannel,
 150                  Int_t dim, Double_t value);
 151
 152   AliMUONVCalibParam* BusPatchParam(Int_t busPatch, Bool_t create=kFALSE) const;
 153
 154   AliMUONVCalibParam* CreateBusPatchParam(Int_t busPatch) const;
 155
 156   AliMUONVCalibParam* ChamberParam(Int_t chamberId, Bool_t create=kFALSE) const;
 157
 158   AliMUONVCalibParam* CreateChamberParam(Int_t chamberId) const;
 159
 160   AliMUONVCalibParam* ChannelParam(Int_t detElemId, Int_t manuId,
 161                                    const AliMUONVCalibParam* external=0x0) const;
 162
 163   AliMUONVCalibParam* DetectionElementParam(Int_t detElemId, Bool_t create=kFALSE) const;
 164
 165   AliMUONVCalibParam* CreateDetectionElementParam(Int_t detElemId) const;
 166
 167   AliMUONVCalibParam* ManuParam(Int_t detElemId, Int_t manuId, Bool_t create=kFALSE) const;
 168
 169   AliMUONVCalibParam* CreateManuParam(Int_t detElemInd, Int_t manuId) const;
 170
 171   AliMUONVCalibParam* PCBParam(Int_t detElemId, Int_t pcbIndex, Bool_t create=kFALSE) const;
 172
 173   AliMUONVCalibParam* CreatePCBParam(Int_t detElemId, Int_t pcbIndex) const;
 174
 175   /// Index of the dimension containing the number of time an item was hit
 176   virtual Int_t IndexOfNumberDimension() const { return fDimension - 1; }
 177
 178   /// Index of the dimension containing the occupancy number
 179   virtual Int_t IndexOfOccupancyDimension() const { return fDimension - 2; }
 180
 181   /// Whether we have histograms for a given dimension, or not
 182   virtual Bool_t IsHistogrammed(Int_t dim) const { return ( fHistogramming[dim] > 0 ); }
 183
 184   Int_t DdlIdFromBusPatchId(Int_t buspatchid) const;
 185   Int_t DdlIdFromDetElemId(Int_t detelemid) const;
 186   Int_t DdlIdFromChamberId(Int_t chamberid) const;
 187
 188   /// Not implemented
 189   AliMUONTrackerData(const AliMUONTrackerData& rhs);
 190   /// Not implemented
 191   AliMUONTrackerData& operator=(const AliMUONTrackerData& rhs);
 192
 193   AliMUONVCalibParam* CreateDouble(const AliMUONVCalibParam& param, Int_t detElemId, Int_t manuId) const;
 194
 195   Int_t GetParts(AliMUONVCalibParam* external,
 196                  AliMUONVCalibParam*& chamber,
 197                  AliMUONVCalibParam*& de,
 198                  AliMUONVCalibParam*& busPatch,
 199                  AliMUONVCalibParam*& pcb,
 200                  AliMUONVCalibParam*& manu,
 201                  AliMUONVCalibParam*& channel,
 202                  AliMpDetElement*& mpde);
 203
 204   /// Convert from external to internal index
 205   Int_t External2Internal(Int_t index) const;
 206
 207   void SetInternalDimensionName(Int_t index, const char* value);
 208
 209   void SetExternalDimensionName(Int_t index, const char* value);
 210
 211   Double_t Value(const AliMUONVCalibParam& param, Int_t i, Int_t dim, Int_t ddlId) const;
 212
 213   /// The number of values we actually *store* for each item
 214   Int_t Dimension() const { return fDimension; }
 215
 216   Bool_t InternalAdd(const AliMUONVStore& store, TArrayI* nevents, Bool_t replace);
 217
 218   void GetDEManu(const AliMUONVCalibParam& param,
 219                   Int_t& detElemId, Int_t& manuId) const;
 220
 221   void AddCalibParams(const AliMUONVCalibParam& src, AliMUONVCalibParam& dest) const;
 222
 223   void Add2D(const AliMUONVStore& src, AliMUONVStore& dest) const;
 224
 225   void Add1D(const AliMUONVStore& src, AliMUONVStore& dest) const;
 226
 227   void AssertStores();
 228
 229   Bool_t UpdateNumberOfEvents(TArrayI* nevents);
 230
 231   static void DispatchValue(AliMUONVCalibParam& param, Int_t index, Double_t y, Double_t ey, Int_t nchannels);
 232
 233 private:
 234
 235   Bool_t fIsSingleEvent; ///< whether we can deal with more than one event
 236   AliMUONVStore* fChannelValues; ///< the channel store
 237   AliMUONVStore* fManuValues; ///< the manu store
 238   AliMUONVStore* fBusPatchValues; ///< the bus patch store
 239   AliMUONVStore* fDEValues; ///< the detection element store
 240   AliMUONVStore* fChamberValues; ///< the chamber store
 241   AliMUONVStore* fPCBValues; ///< the pcb store
 242   Int_t fDimension; ///< the dimension of the data
 243   Int_t fNevents; ///< the number of events treated
 244   TObjArray* fDimensionNames; ///< the names of the (internal) dimensions
 245   TObjArray* fExternalDimensionNames; ///< the names of the external (i.e. original) dimensions
 246   Int_t fExternalDimension; ///< number of interface values per item
 247   /// whether we should histogram the dimension(s)
 248   Int_t* fHistogramming; //[fExternalDimension] whether we should histogram the dimension(s)
 249   AliMUONVStore* fHistos; ///< the lowest histograms we have
 250   Double_t fXmin; ///< min x value for histograms
 251   Double_t fXmax; ///< max x value for histograms
 252   static const Int_t fgkExtraDimension; ///< to hold extra information
 253   static const Int_t fgkVirtualExtraDimension; ///< to give access to information not stored, but computed on the fly
 254   Bool_t fIsChannelLevelEnabled; ///< whether we allow storing of channel (fChannelValues) values
 255   Bool_t fIsManuLevelEnabled; ///< whether we allow storing of manu (fManuValues) values
 256   Bool_t fIsBustPatchLevelEnabled; ///< whether we allow storing of bus patches (fBusPatchValues) values
 257   Bool_t fIsPCBLevelEnabled; ///< whether we allow storing of PCB values (fPCBValues)
 258   Int_t fNofDDLs; ///< nof of DDLs we're dealing with
 259   /// the number of events treated (per DDL)
 260   Int_t* fNofEventsPerDDL; //[fNofDDLs] the number of events treated (per DDL)
 261
 262   ClassDef(AliMUONTrackerData,8) // Implementation of AliMUONVTrackerData
 263 };
 264
 265 #endif