#ifndef ROOT_TObject
# include <TObject.h>
#endif
+class TFile;
+
//____________________________________________________________________
-//
-// Base class for caches of per-strip information.
-// This is used to index a strip.
-// Data stored depends on derived class.
-//
+/** @class AliFMDMap
+ @brief Base class for caches of per-strip information.
+ @ingroup FMD_data
+ This is used to index a strip. Data stored depends on derived
+ class. */
class AliFMDMap : public TObject
{
public:
enum {
+ /** Default maximum detector number */
kMaxDetectors = 3,
+ /** Default maximum number of rings */
kMaxRings = 2,
+ /** Default maximum number of sectors */
kMaxSectors = 40,
+ /** Default maximum number of strips */
kMaxStrips = 512
};
- AliFMDMap(size_t maxDet = kMaxDetectors,
- size_t maxRing= kMaxRings,
- size_t maxSec = kMaxSectors,
- size_t maxStr = kMaxStrips);
+ /** Constructor
+ @param maxDet Maximum allowed detector number
+ @param maxRing Maximum number of rings
+ @param maxSec Maximum number of sectors
+ @param maxStr Maximum number of strips
+ @return */
+ AliFMDMap(UShort_t maxDet = kMaxDetectors,
+ UShort_t maxRing= kMaxRings,
+ UShort_t maxSec = kMaxSectors,
+ UShort_t maxStr = kMaxStrips);
+ /** Destructor */
virtual ~AliFMDMap() {}
- size_t MaxDetectors() const { return fMaxDetectors; }
- size_t MaxRings() const { return fMaxRings; }
- size_t MaxSectors() const { return fMaxSectors; }
- size_t MaxStrips() const { return fMaxStrips; }
- Int_t CheckIndex(size_t det, Char_t ring, size_t sec, size_t str) const;
+ /** @return Maximum detector number */
+ UShort_t MaxDetectors() const { return fMaxDetectors; }
+ /** @return Maximum number of rings */
+ UShort_t MaxRings() const { return fMaxRings; }
+ /** @return Maximum number of sectors */
+ UShort_t MaxSectors() const { return fMaxSectors; }
+ /** @return Maximum number of strip */
+ UShort_t MaxStrips() const { return fMaxStrips; }
+ /** Calculate, check, and return index for strip. If the index is
+ invalid, -1 is returned
+ @param det Detector number
+ @param ring Ring identifier
+ @param sec Sector number
+ @param str Strip number
+ @return Unique index, or -1 in case of errors */
+ Int_t CheckIndex(UShort_t det, Char_t ring,
+ UShort_t sec, UShort_t str) const;
+ /** Check if we need UShort_t hack
+ @param file File this object was read from */
+ void CheckNeedUShort(TFile* file);
+ enum {
+ /** In case of version 2 of this class, this bit should be set. */
+ kNeedUShort = 14
+ };
protected:
- size_t CalcIndex(size_t det, Char_t ring, size_t sec, size_t str) const;
- size_t fMaxDetectors; // Maximum # of detectors
- size_t fMaxRings; // Maximum # of rings
- size_t fMaxSectors; // Maximum # of sectors
- size_t fMaxStrips; // Maximum # of strips
- ClassDef(AliFMDMap, 1) // Cache of per strip information
-};
+ /** Calculate index and return
+ @param det Detector number
+ @param ring Ring identifier
+ @param sec Sector number
+ @param str Strip number
+ @return Index (not checked) */
+ Int_t CalcIndex(UShort_t det, Char_t ring,
+ UShort_t sec, UShort_t str) const;
+ UShort_t fMaxDetectors; // Maximum # of detectors
+ UShort_t fMaxRings; // Maximum # of rings
+ UShort_t fMaxSectors; // Maximum # of sectors
+ UShort_t fMaxStrips; // Maximum # of strips
-#ifdef MAY_USE_TEMPLATES
-#ifndef __VECTOR__
-# include <vector>
-#endif
-//____________________________________________________________________
-//
-// Class template for classes that cache per strip information.
-// Access to the data goes via
-//
-// Type& AliFMDMap<Type>::operator()(size_t detector,
-// Char_t ring,
-// size_t sector,
-// size_t strip);
-//
-// (as well as a const version of this member function).
-// The elements can be reset to the default value by calling
-// AliFMDMap<Type>::Clear(). This resets the values to `Type()'.
-//
-template <typename Type>
-class AliFMDMap : public TObject
-{
-public:
- AliFMDMap(size_t maxDet=3, size_t maxRing=2, size_t maxSec=40,
- size_t maxStr=512);
- virtual ~AliFMDMap() {}
- void Clear(const Type& val=Type());
- Type& operator()(size_t det, Char_t ring, size_t sec, size_t str);
- const Type& operator()(size_t det, Char_t ring, size_t sec, size_t str)const;
-private:
- typedef std::vector<Type> ValueVector; // Type of container
- ValueVector fValues; // Contained values
- size_t fMaxDetectors; // Maximum # of detectors
- size_t fMaxRings; // Maximum # of rings
- size_t fMaxSectors; // Maximum # of sectors
- size_t fMaxStrips; // Maximum # of strips
-
- size_t CalcIndex(size_t det, Char_t ring, size_t sec, size_t str) const;
- ClassDef(AliFMDMap, 0); // Map of FMD index's to values
+ ClassDef(AliFMDMap, 3) // Cache of per strip information
};
-
-//____________________________________________________________________
-template <typename Type>
-inline
-AliFMDMap<Type>::AliFMDMap(size_t maxDet,
- size_t maxRing,
- size_t maxSec,
- size_t maxStr)
- : fValues(maxDet * maxRing * maxSec * maxStr),
- fMaxDetectors(maxDet),
- fMaxRings(maxRing),
- fMaxSectors(maxSec),
- fMaxStrips(maxStr)
-{
- // Construct a map
- //
- // Parameters:
- // maxDet Maximum # of detectors
- // maxRinf Maximum # of rings
- // maxSec Maximum # of sectors
- // maxStr Maximum # of strips
-}
-
-
-//____________________________________________________________________
-template <typename Type>
-inline size_t
-AliFMDMap<Type>::CalcIndex(size_t det, Char_t ring, size_t sec, size_t str) const
-{
- // Calculate index into storage from arguments.
- //
- // Parameters:
- // det Detector #
- // ring Ring ID
- // sec Sector #
- // str Strip #
- //
- // Returns appropriate index into storage
- //
- size_t ringi = (ring == 'I' || ring == 'i' ? 0 : 1);
- size_t idx =
- (det + fMaxDetectors * (ringi + fMaxRings * (sec + fMaxSectors * str)));
- if (idx >= fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips) {
- Fatal("CalcIndex", "Index (%d,'%c',%d,%d) out of bounds, "
- "in particular the %s index",
- det, ring, sec, str,
- (det >= fMaxDetectors ? "Detector" :
- (ringi >= fMaxRings ? "Ring" :
- (sec >= fMaxSectors ? "Sector" : "Strip"))));
- return 0;
- }
- return idx;
-}
-
-//____________________________________________________________________
-template <typename Type>
-inline void
-AliFMDMap<Type>::Clear(const Type& val)
-{
- // Resets stored values to the default value for that type
- for (size_t i = 0; i < fValues.size(); ++i) fValues[i] = val;
-}
-
-//____________________________________________________________________
-template <typename Type>
-inline Type&
-AliFMDMap<Type>::operator()(size_t det, Char_t ring, size_t sec, size_t str)
-{
- // Parameters:
- // det Detector #
- // ring Ring ID
- // sec Sector #
- // str Strip #
- //
- // Returns data[det][ring][sec][str]
- return fValues[CalcIndex(det, ring, sec, str)];
-}
-
-//____________________________________________________________________
-template <typename Type>
-inline const Type&
-AliFMDMap<Type>::operator()(size_t det,
- Char_t ring,
- size_t sec,
- size_t str) const
-{
- // Parameters:
- // det Detector #
- // ring Ring ID
- // sec Sector #
- // str Strip #
- //
- // Returns data[det][ring][sec][str]
- return fValues[CalcIndex(det, ring, sec, str)];
-}
-
-
-//____________________________________________________________________
-//
-// Some specialisations
-//
-typedef AliFMDMap<UShort_t> AliFMDAdcMap;
-typedef AliFMDMap<std::pair<Float_t, UShort_t> > AliFMDEdepMap;
-
-#endif
#endif
//____________________________________________________________________
//