*/
//____________________________________________________________________
//
+// This class stores the sample rate (that is, how many times the
+// ATLRO's sample each VA1 channel). In principle these can be
+// controlled per half ring, but in real life it's most likely that
+// this value will be the same for all detectors. This value must be
+// retrived from DCS or the like.
//
+// IMPORTANT: The member function WriteToFile writes out the entries
+// in the format
+//
+// det,ring,id,rate
+//
+// Here, id is a number from 0 to 1, which represents the division in
+// half-rings. The mapping is as follows:
+//
+// Inner rings: Outer Rings
+// id Sectors Board id Sectors Board
+// ----+---------+------- ----+---------+-------
+// 0 | 0 - 9 | 0x10 0 | 0 - 19 | 0x11
+// 1 | 10 - 19 | 0x0 1 | 20 - 39 | 0x1
+//
+// The same mapping is used in the ReadFromFile member function
//
#include "AliFMDCalibSampleRate.h" // ALIFMDCALIBGAIN_H
-#include "AliFMDParameters.h" // ALIFMDPARAMETERS_H
-#include <AliLog.h>
+// #include "AliFMDParameters.h" // ALIFMDPARAMETERS_H
+// #include <AliLog.h>
+#include "TString.h"
+#include "AliFMDDebug.h" // Better debug macros
+#include <iostream>
//____________________________________________________________________
ClassImp(AliFMDCalibSampleRate)
: fRates(AliFMDMap::kMaxDetectors, AliFMDMap::kMaxRings, 2, 1)
// fRates(3)
{
+ // CTOR
fRates.Reset(1);
}
//____________________________________________________________________
AliFMDCalibSampleRate::AliFMDCalibSampleRate(const AliFMDCalibSampleRate& o)
: TObject(o), fRates(o.fRates)
-{}
+{
+ // Copy ctor
+}
//____________________________________________________________________
AliFMDCalibSampleRate&
AliFMDCalibSampleRate::operator=(const AliFMDCalibSampleRate& o)
{
+ // Assignment operator
fRates = o.fRates;
return (*this);
}
AliFMDCalibSampleRate::Set(UShort_t det, Char_t ring,
UShort_t sector, UShort_t, UShort_t rate)
{
- UInt_t nSec = (ring == 'I' ? 20 : 40);
+ // Set values. Strip argument is ignored
+ UInt_t nSec = (ring == 'I' ? 10 : 20);
UInt_t board = sector / nSec;
fRates(det, ring, board, 0) = rate;
+ AliFMDDebug(15, ("Setting sample rate for FMD%d%c[%2d,0] (board %d): %d",
+ det, ring, sector, board, rate));
+
}
//____________________________________________________________________
AliFMDCalibSampleRate::Rate(UShort_t det, Char_t ring,
UShort_t sec, UShort_t) const
{
- UInt_t nSec = (ring == 'I' ? 20 : 40);
- UInt_t board = sec / nSec;
- AliDebug(10, Form("Getting sample rate for FMD%d%c[%2d,0] (board %d)",
- det, ring, sec, board));
- return fRates(det, ring, board, 0);
+ // Get the sample rate
+ UInt_t nSec = (ring == 'I' ? 10 : 20);
+ UInt_t board = sec / nSec;
+ UShort_t ret = fRates(det, ring, board, 0);
+ AliFMDDebug(15, ("Getting sample rate for FMD%d%c[%2d,0] (board %d): %d",
+ det, ring, sec, board, ret));
+ return ret;
+}
+//____________________________________________________________________
+void
+AliFMDCalibSampleRate::WriteToFile(std::ostream &outFile, Bool_t* detectors)
+{
+ outFile.write("# SampleRate \n",14);
+ for(Int_t det=1;det<=3;det++) {
+ if (detectors && !detectors[det-1]) {
+ continue;
+ }
+ UShort_t FirstRing = (det == 1 ? 1 : 0);
+ for (UShort_t ir = FirstRing; ir < 2; ir++) {
+ Char_t ring = (ir == 0 ? 'O' : 'I');
+ UShort_t nsec = (ir == 0 ? 40 : 20) / 2;
+
+ for(UShort_t board = 0; board < 2; board++) {
+ UShort_t sector = board*nsec;
+ outFile << det << ','
+ << ring << ','
+ << board << ','
+ << Rate(det,ring,sector) << "\n";
+
+
+ }
+ }
+ }
+
+
+}
+//____________________________________________________________________
+void
+AliFMDCalibSampleRate::ReadFromFile(std::istream &inFile)
+{
+ TString line;
+ Bool_t readData=kFALSE;
+
+ while(line.ReadLine(inFile)) {
+ if(line.Contains("samplerate",TString::kIgnoreCase)) {
+ readData = kTRUE;
+ break;
+ }
+
+ }
+
+ UShort_t det, board;
+ Char_t ring;
+ UShort_t sampleRate;
+ Int_t thisline = inFile.tellg();
+ Char_t c[3];
+
+ while( readData ) {
+ thisline = inFile.tellg();
+ line.ReadLine(inFile);
+ if(line.Contains("# ",TString::kIgnoreCase)) {
+ readData = kFALSE;
+ continue;
+ }
+
+ inFile.seekg(thisline);
+ inFile >> det >> c[0]
+ >> ring >> c[1]
+ >> board >> c[2]
+ >> sampleRate;
+
+ UInt_t nSec = (ring == 'I' ? 20 : 40)/2;
+ UShort_t sector = board*nSec;
+ Set(det,ring,sector,0,sampleRate);
+
+
+ }
+
+ inFile.seekg(0);
+
+
}
//____________________________________________________________________