--- /dev/null
+/**************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN. *
+ * All rights reserved. *
+ * *
+ * Author: The ALICE Off-line Project. *
+ * Contributors are mentioned in the code where appropriate. *
+ * *
+ * Permission to use, copy, modify and distribute this *
+ * software and its documentation strictly for non-commercial *
+ * purposes is hereby granted without fee, provided that the *
+ * above copyright notice appears in all copies and that both *
+ * the copyright notice and this permission notice appear in *
+ * the supporting documentation. The authors make no claims *
+ * about the suitability of this software for any purpose. It *
+ * is provided "as is" without express or implied warranty. *
+ **************************************************************/
+/* $Id$ */
+//__________________________________________________________
+//
+// Map of per strip Float_t information
+// the floats are indexed by the coordinates
+// DETECTOR # (1-3)
+// RING ID ('I' or 'O', any case)
+// SECTOR # (0-39)
+// STRIP # (0-511)
+//
+//
+// Created Mon Nov 8 12:51:51 2004 by Christian Holm Christensen
+//
+#include "AliFMDFloatMap.h" //ALIFMDFLOATMAP_H
+namespace {
+ class Printer : public AliFMDMap::ForOne
+ {
+ public:
+ Printer(const char* format)
+ : fFormat(format), fOldD(0), fOldR('-'), fOldS(1024) {}
+ Bool_t operator()(UShort_t d, Char_t r, UShort_t s, UShort_t t,
+ Float_t m)
+ {
+ if (d != fOldD) {
+ fOldD = d;
+ fOldR = '-';
+ if (d != 0) printf("\n");
+ printf("FMD%d", fOldD);
+ }
+ if (r != fOldR) {
+ fOldR = r;
+ fOldS = 1024;
+ printf("\n %s ring", (r == 'I' ? "Inner" : "Outer"));
+ }
+ if (s != fOldS) {
+ fOldS = s;
+ printf("\n Sector %2d", fOldS);
+ }
+ if (t % 4 == 0) printf("\n %3d-%3d ", t, t+3);
+ printf(fFormat, m);
+ // if (t % 4 == 3) printf("\n");
+
+ return kTRUE;
+ }
+ Bool_t operator()(UShort_t, Char_t, UShort_t, UShort_t, Int_t)
+ {
+ return kTRUE;
+ }
+ Bool_t operator()(UShort_t, Char_t, UShort_t, UShort_t, UShort_t)
+ {
+ return kTRUE;
+ }
+ Bool_t operator()(UShort_t, Char_t, UShort_t, UShort_t, Bool_t)
+ {
+ return kTRUE;
+ }
+ private:
+ Printer(const Printer& p)
+ : AliFMDMap::ForOne(p),
+ fFormat(p.fFormat),
+ fOldD(p.fOldD),
+ fOldR(p.fOldR),
+ fOldS(p.fOldS)
+ {}
+ Printer& operator=(const Printer&) { return *this; }
+ const char* fFormat;
+ UShort_t fOldD;
+ Char_t fOldR;
+ UShort_t fOldS;
+ };
+}
+//__________________________________________________________
+ClassImp(AliFMDFloatMap)
+#if 0
+ ; // This is here to keep Emacs for indenting the next line
+#endif
+
+//__________________________________________________________
+AliFMDFloatMap::AliFMDFloatMap(const AliFMDMap& other)
+ : AliFMDMap(other),
+ fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
+ fData(0)
+{
+ if (fTotal == 0) fTotal = 51200;
+ fData = new Float_t[fTotal];
+ // Copy constructor
+ if (!other.IsFloat()) return;
+ for (Int_t i = 0; i < fTotal; i++) fData[i] = other.AtAsFloat(i);
+}
+
+//__________________________________________________________
+AliFMDFloatMap::AliFMDFloatMap(const AliFMDFloatMap& other)
+ : AliFMDMap(other.fMaxDetectors,
+ other.fMaxRings,
+ other.fMaxSectors,
+ other.fMaxStrips),
+ fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
+ fData(0)
+{
+ if (fTotal == 0) fTotal = 51200;
+ fData = new Float_t[fTotal];
+ // Copy constructor
+ for (Int_t i = 0; i < fTotal; i++)
+ fData[i] = other.fData[i];
+}
+
+//__________________________________________________________
+AliFMDFloatMap::AliFMDFloatMap()
+ : AliFMDMap(),
+ fTotal(0),
+ fData(0)
+{
+ // Constructor.
+ // Parameters:
+ // None
+}
+
+//__________________________________________________________
+AliFMDFloatMap::AliFMDFloatMap(Int_t maxDet,
+ Int_t maxRing,
+ Int_t maxSec,
+ Int_t maxStr)
+ : AliFMDMap(maxDet, maxRing, maxSec, maxStr),
+ fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
+ fData(0)
+{
+ // Constructor.
+ // Parameters:
+ // maxDet Maximum number of detectors
+ // maxRing Maximum number of rings per detector
+ // maxSec Maximum number of sectors per ring
+ // maxStr Maximum number of strips per sector
+ if (fTotal == 0) fTotal = 51200;
+ fData = new Float_t[fTotal];
+ Reset(0);
+}
+
+//__________________________________________________________
+AliFMDFloatMap&
+AliFMDFloatMap::operator=(const AliFMDFloatMap& other)
+{
+ // Assignment operator
+ if(&other != this){
+ if(fMaxDetectors!= other.fMaxDetectors||
+ fMaxRings != other.fMaxRings||
+ fMaxSectors != other.fMaxSectors||
+ fMaxStrips != other.fMaxStrips){
+ // allocate new memory only if the array size is different....
+ fMaxDetectors = other.fMaxDetectors;
+ fMaxRings = other.fMaxRings;
+ fMaxSectors = other.fMaxSectors;
+ fMaxStrips = other.fMaxStrips;
+ fTotal = fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips;
+ if (fTotal == 0) fTotal = 51200;
+ if (fData) delete [] fData;
+ fData = new Float_t[fTotal];
+ }
+ for (Int_t i = 0; i < fTotal; i++) fData[i] = other.fData[i];
+ }
+ return *this;
+}
+
+
+//__________________________________________________________
+void
+AliFMDFloatMap::Reset(const Float_t& val)
+{
+ // Reset map to val
+ for (Int_t i = 0; i < fTotal; i++) fData[i] = val;
+}
+
+//__________________________________________________________
+void
+AliFMDFloatMap::Print(Option_t* option) const
+{
+ // Print contents of map
+ if (!option || option[0] == '\0') TObject::Print();
+ Printer p(option);
+ ForEach(p);
+ printf("\n");
+}
+//__________________________________________________________
+Float_t&
+AliFMDFloatMap::operator()(UShort_t det,
+ Char_t ring,
+ UShort_t sec,
+ UShort_t str)
+{
+ // Get data
+ // Parameters:
+ // det Detector #
+ // ring Ring ID
+ // sec Sector #
+ // str Strip #
+ // Returns appropriate data
+ return fData[CalcIndex(det, ring, sec, str)];
+}
+
+//__________________________________________________________
+const Float_t&
+AliFMDFloatMap::operator()(UShort_t det,
+ Char_t ring,
+ UShort_t sec,
+ UShort_t str) const
+{
+ // Get data
+ // Parameters:
+ // det Detector #
+ // ring Ring ID
+ // sec Sector #
+ // str Strip #
+ // Returns appropriate data
+ return fData[CalcIndex(det, ring, sec, str)];
+}
+
+//__________________________________________________________
+//
+// EOF
+//
+