+//____________________________________________________________________
+Int_t
+AliFMDRawReader::NewDDL(AliAltroRawStreamV3& input, UShort_t& det)
+{
+ // Process a new DDL. Sets the internal data members fZeroSuppress,
+ // fSampleRate, and fNoiseFactor based on information in the RCU trailer.
+ //
+ // Parameters:
+ // input Input stream
+ // det On return, the detector number
+ //
+ // Return value:
+ // negative value in case of problems, the DDL number otherwise
+
+ // Get the DDL number
+ UInt_t ddl = input.GetDDLNumber();
+ AliFMDDebug(2, ("DDL number %d", ddl));
+
+ // Get zero suppression flag
+ fZeroSuppress[ddl] = input.GetZeroSupp();
+ AliFMDDebug(3, ("RCU @ DDL %d zero suppression: %s",
+ ddl, (fZeroSuppress[ddl] ? "yes" : "no")));
+
+ // WARNING: We store the noise factor in the 2nd baseline
+ // filters excluded post samples, since we'll never use that
+ // mode.
+ fNoiseFactor[ddl] = input.GetNPostsamples();
+ AliFMDDebug(3, ("RCU @ DDL %d noise factor: %d", ddl,fNoiseFactor[ddl]));
+
+ // WARNING: We store the sample rate in the number of pre-trigger
+ // samples, since we'll never use that mode.
+ fSampleRate[ddl] = input.GetNPretriggerSamples();
+ AliFMDDebug(3, ("RCU @ DDL %d sample rate: %d", ddl,fSampleRate[ddl]));
+
+ // Get Errors seen
+ Int_t nChAddrMismatch = input.GetNChAddrMismatch();
+ Int_t nChLenMismatch = input.GetNChLengthMismatch();
+ if (nChAddrMismatch != 0)
+ AliWarning(Form("Got %d channels with address mis-matches for 0x%03x",
+ nChAddrMismatch, ddl));
+ if (nChLenMismatch != 0)
+ AliWarning(Form("Got %d channels with length mis-matches for 0x%03x",
+ nChLenMismatch, ddl));
+
+ // Map DDL number to the detector number
+ AliFMDParameters* pars = AliFMDParameters::Instance();
+ AliFMDAltroMapping* map = pars->GetAltroMap();
+ if (map->DDL2Detector(ddl) < 0) return -1;
+ det = map->DDL2Detector(ddl);
+
+ if (AliLog::GetDebugLevel("FMD", 0) > 5)
+ input.PrintRCUTrailer();
+ return ddl;
+}
+
+//____________________________________________________________________
+Int_t
+AliFMDRawReader::NewChannel(AliAltroRawStreamV3& input, UShort_t det,
+ Char_t& ring, UShort_t& sec, Short_t& strbase)
+{
+ // Processs a new channel. Sets the internal data members
+ // fMinStrip, fMaxStrip, and fPreSamp.
+ //
+ // Parameter:
+ // input Input stream
+ // ring On return, the ring identifier
+ // sec On return, the sector number
+ // strbase On return, the strip base
+ //
+ // Return value
+ // negative value in case of problems, hardware address otherwise
+
+ // Get the hardware address, and map that to detector coordinates
+ UShort_t board, chip, channel;
+ Int_t ddl = input.GetDDLNumber();
+ Int_t hwaddr = input.GetHWAddress();
+ if (input.IsChannelBad()) {
+ AliError(Form("Channel 0x%03x is marked as bad!", hwaddr));
+ }
+
+ AliFMDParameters* pars = AliFMDParameters::Instance();
+ AliFMDAltroMapping* map = pars->GetAltroMap();
+ // Map to hardware stuff
+ map->ChannelAddress(hwaddr, board, chip, channel);
+ // Then try to map to detector address
+ if (!map->Channel2StripBase(board, chip, channel, ring, sec, strbase)) {
+ AliError(Form("Failed to get detector id from DDL %d, "
+ "hardware address 0x%03x", ddl, hwaddr));
+ return -1;
+ }
+ AliFMDDebug(3, ("Board: 0x%02x, Altro: 0x%x, Channel: 0x%x",
+ board, chip, channel));
+
+ // Get the 'conditions'
+ fMinStrip = pars->GetMinStrip(det, ring, sec, strbase);
+ fMaxStrip = pars->GetMaxStrip(det, ring, sec, strbase);
+ fPreSamp = pars->GetPreSamples(det, ring, sec, strbase);
+ if (fSampleRate[ddl] == 0)
+ fSampleRate[ddl] = pars->GetSampleRate(det, ring, sec, strbase);
+
+ return hwaddr;
+}
+
+//____________________________________________________________________
+Int_t
+AliFMDRawReader::NewSample(AliAltroRawStreamV3& input,
+ Int_t i, UShort_t t, UShort_t sec,
+ UShort_t strbase, Short_t& str, UShort_t& samp)
+{
+ // Process a new timebin
+ //
+ // Parameters:
+ // input Input stream
+ // i Index into bunch data
+ // t Time
+ // strbase Base of strip numbers for this channel
+ // str On return, the strip number
+ // samp On return, the sample number
+ //
+ // Return value
+ // negative value in case of problems, ADC value otherwise
+ if (t < fPreSamp) return -1;
+
+ Int_t ddl = input.GetDDLNumber();
+ Int_t hwa = input.GetHWAddress();
+ const UShort_t* data = input.GetSignals();
+ Short_t adc = data[i];
+ AliFMDDebug(10, ("0x%04x/0x%03x/%04d %4d", ddl, hwa, t, adc));
+
+ AliFMDParameters* pars = AliFMDParameters::Instance();
+ AliFMDAltroMapping* map = pars->GetAltroMap();
+
+ samp = 0;
+ Short_t stroff = 0;
+ map->Timebin2Strip(sec, t, fPreSamp, fSampleRate[ddl], stroff, samp);
+ str = strbase + stroff;
+
+ AliFMDDebug(20, ("0x%04x/0x%03x/%04d=%4d maps to strip %3d sample %d "
+ "(pre: %d, min: %d, max: %d, rate: %d)",
+ ddl, hwa, t, adc, str, samp, fPreSamp,
+ fMinStrip, fMaxStrip, fSampleRate[ddl]));
+ if (str < 0) {
+ AliFMDDebug(10, ("Got presamples at timebin %d", i));
+ return -1;
+ }
+
+ // VA1 Local strip number
+ Short_t lstrip = (t - fPreSamp) / fSampleRate[ddl] + fMinStrip;
+
+ AliFMDDebug(15, ("Checking if strip %d (%d) in range [%d,%d]",
+ lstrip, str, fMinStrip, fMaxStrip));
+ if (lstrip < fMinStrip || lstrip > fMaxStrip) {
+ AliFMDDebug(10, ("Strip %03d-%d (%d,%d) from t=%d out of range (%3d->%3d)",
+ str, samp, lstrip, stroff, t, fMinStrip, fMaxStrip));
+ adc = -1;
+ }
+ // Possibly do pedestal subtraction of signal
+ if (adc > 1023)
+ AliWarning(Form("ADC value out of range: %4d", adc));
+ return adc;
+}
+
+//____________________________________________________________________
+Bool_t
+AliFMDRawReader::NextSample(UShort_t& det, Char_t& rng, UShort_t& sec,
+ UShort_t& str, UShort_t& sam, UShort_t& rat,
+ Short_t& adc, Bool_t& zs, UShort_t& fac)
+{
+ // Scan current event for next signal. It returns kFALSE when
+ // there's no more data in the event.
+ static AliAltroRawStreamV3 stream(fReader); // = 0;
+ static Int_t ddl = -1;
+ static UShort_t tdet = 0;
+ static Char_t trng = '\0';
+ static UShort_t tsec = 0;
+ static Short_t tstr = 0;
+ static Short_t bstr = -1;
+ static UShort_t tsam = 0;
+ static UInt_t trate = 0;
+ static Int_t hwaddr = -1;
+ static UShort_t start = 0;
+ static UShort_t length = 0;
+ static Short_t t = -1;
+ static Int_t i = 0;
+ // First entry!
+ if (stream.GetDDLNumber() < 0) {
+ fReader->Select("FMD");
+
+ // Reset variables
+ ddl = -1;
+ trate = 0;
+ tdet = 0;
+ trng = '\0';
+ tsec = 0;
+ tstr = 0;
+ tsam = -1;
+ hwaddr = -1;
+ }
+
+ do {
+ if (t < start - length + 1) {
+ if (!stream.NextBunch()) {
+ if (!stream.NextChannel()) {
+ if (!stream.NextDDL()) {
+ stream.Reset();
+ return kFALSE;
+ }
+ ddl = NewDDL(stream, tdet);
+ }
+ hwaddr = NewChannel(stream, tdet, trng, tsec, bstr);
+ }
+ start = stream.GetStartTimeBin();
+ length = stream.GetBunchLength();
+ t = start;
+ i = 0;
+ }
+ Int_t tadc = NewSample(stream, i, t, tsec, bstr, tstr, tsam);
+ if (tadc >= 0) {
+ det = tdet;
+ rng = trng;
+ sec = tsec;
+ str = tstr;
+ sam = tsam;
+ adc = tadc;
+ rat = fSampleRate[ddl];
+ zs = fZeroSuppress[ddl];
+ fac = fNoiseFactor[ddl];
+ t--;
+ i++;
+ break;
+ }
+ } while (true);
+
+ return kTRUE;
+}
+
+
+//____________________________________________________________________
+Bool_t
+AliFMDRawReader::NextSignal(UShort_t& det, Char_t& rng,
+ UShort_t& sec, UShort_t& str,
+ Short_t& adc, Bool_t& zs,
+ UShort_t& fac)
+{
+
+ do {
+ UShort_t samp, rate;
+ if (!NextSample(det, rng, sec, str, samp, rate, adc, zs, fac))
+ return kFALSE;
+
+ Bool_t take = kFALSE;
+ switch (rate) {
+ case 1: take = kTRUE; break;
+ case 2: if (samp == 1) take = kTRUE; break;
+ case 3: if (samp == 1) take = kTRUE; break;
+ case 4: if (samp == 2) take = kTRUE; break;
+ default: if (samp == rate-2) take = kTRUE; break;
+ }
+ if (!take) continue;
+ break;
+ } while (true);
+ return kTRUE;
+}
+
+//____________________________________________________________________
+Bool_t
+AliFMDRawReader::SelectSample(UShort_t samp, UShort_t rate)
+{
+ // Check if the passed sample is the one we need
+ Bool_t take = kFALSE;
+ switch (rate) {
+ case 1: take = kTRUE; break;
+ case 2: if (samp == 1) take = kTRUE; break;
+ case 3: if (samp == 1) take = kTRUE; break;
+ case 4: if (samp == 2) take = kTRUE; break;
+ default: if (samp == rate-2) take = kTRUE; break;
+ }
+
+ return take;
+}
+
+//____________________________________________________________________
+Bool_t
+AliFMDRawReader::ReadAdcs(TClonesArray* array)
+{
+ // Read ADC values from raw input into passed TClonesArray of AliFMDDigit
+ // objects.
+ AliFMDDebug(2, ("Reading ADC values into a TClonesArray"));
+
+ // Read raw data into the digits array, using AliFMDAltroReader.
+ if (!array) {
+ AliError("No TClonesArray passed");
+ return kFALSE;
+ }
+ const UShort_t kUShortMax = (1 << 16) - 1;
+ fSeen.Reset(kUShortMax);
+
+ AliAltroRawStreamV3 input(fReader);
+ input.Reset();
+ input.SelectRawData("FMD");
+
+ // Loop over input RORCs
+ while (input.NextDDL()) {
+ UShort_t det = 0;
+ Int_t ddl = NewDDL(input, det);
+ if (ddl < 0) break;
+
+ while (input.NextChannel()) {
+ // Get the hardware address, and map that to detector coordinates
+ Char_t ring;
+ UShort_t sec;
+ Short_t strbase;
+ Int_t hwaddr = NewChannel(input, det, ring, sec, strbase);
+ if (hwaddr < 0) break;
+
+ // Loop over bunches
+ while (input.NextBunch()) {
+ // Get Lenght of bunch, and pointer to the data
+ const UShort_t* data = input.GetSignals();
+ UShort_t start = input.GetStartTimeBin();
+ UShort_t length = input.GetBunchLength();
+
+ // Loop over the data and store it.
+ for (Int_t i = 0; i < length; i++) {
+ // Time
+ Short_t str;
+ UShort_t samp;
+ Int_t t = start - i;
+ Int_t adc = NewSample(input, i, t, sec, strbase, str, samp);
+ if (adc < 0) continue;
+ UShort_t counts = adc;
+
+ AliFMDDebug(10, ("FMD%d%c[%02d,%03d]-%d: %4d",
+ det, ring, sec, str, samp, counts));
+ // Check the cache of indicies
+ Int_t idx = fSeen(det, ring, sec, str);
+ AliFMDDigit* digit = 0;
+ if (idx == kUShortMax) {
+ // We haven't seen this strip yet.
+ fSeen(det, ring, sec, str) = idx = array->GetEntriesFast();
+ AliFMDDebug(7,("making digit for FMD%d%c[%2d,%3d]-%d "
+ "from timebin %4d",
+ det, ring, sec, str, samp, t));
+ digit = new ((*array)[idx]) AliFMDDigit(det, ring, sec, str);
+ digit->SetDefaultCounts(fSampleRate[ddl]);
+ }
+ else
+ digit = static_cast<AliFMDDigit*>(array->At(idx));
+ AliFMDDebug(10, ("Setting FMD%d%c[%2d,%3d]-%d "
+ "from timebin %4d=%4d (%4d)",
+ det, ring, sec, str, samp, t, counts, data[i]));
+ digit->SetCount(samp, counts);
+ } // for (i)
+ } // while (bunch)
+ } // while (channel)
+ } // while (ddl)
+ return kTRUE;
+}
+
+//____________________________________________________________________
+Bool_t AliFMDRawReader::ReadSODevent(AliFMDCalibSampleRate* sampleRate,
+ AliFMDCalibStripRange* stripRange,
+ TArrayS &pulseSize,
+ TArrayS &pulseLength,
+ Bool_t* detectors)
+{
+
+ AliFMDDebug(0, ("Start of SOD/EOD"));
+
+ UInt_t shift_clk[18];
+ UInt_t sample_clk[18];
+ UInt_t strip_low[18];
+ UInt_t strip_high[18];
+ UInt_t pulse_size[18];
+ UInt_t pulse_length[18];
+ AliFMDParameters* param = AliFMDParameters::Instance();
+ AliFMDAltroMapping* map = param->GetAltroMap();
+
+ AliAltroRawStreamV3 streamer(fReader);
+ streamer.Reset();
+ streamer.SelectRawData("FMD");
+ //fReader->GetDDLID();
+ while (streamer.NextDDL()) {
+ Int_t ddl = streamer.GetDDLNumber();
+ Int_t detID = fReader->GetDetectorID();
+ if (detectors) detectors[map->DDL2Detector(ddl)-1] = kTRUE;
+ AliFMDDebug(0, (" From reader: DDL number is %d , det ID is %d",ddl,detID));
+
+ ULong_t nPayloadWords = streamer.GetRCUPayloadSizeInSOD();
+ UChar_t* payloadData = streamer.GetRCUPayloadInSOD();
+ UInt_t* payloadWords = reinterpret_cast<UInt_t*>(payloadData);
+ //UInt_t* payloadWords = streamer.GetRCUPayloadInSOD();
+
+ //std::cout<<nPayloadWords<<" "<<ddl<<std::endl;
+ for (ULong_t i = 1; i <= nPayloadWords ; i++, payloadWords++) {
+ UInt_t payloadWord = *payloadWords; // Get32bitWord(i);
+
+ //std::cout<<i<<Form(" word: 0x%x",payloadWord)<<std::endl;
+ // address is only 24 bit
+ UInt_t address = (0xffffff & payloadWord);
+ UInt_t type = ((address >> 21) & 0xf);
+ UInt_t error = ((address >> 20) & 0x1);
+ UInt_t bcast = ((address >> 18) & 0x1);
+ UInt_t bc_not_altro = ((address >> 17) & 0x1);
+ UInt_t board = ((address >> 12) & 0x1f);
+ UInt_t instruction = 0;
+ UInt_t chip = 0;
+ UInt_t channel = 0;
+ if(bc_not_altro)
+ instruction = address & 0xfff;
+ else {
+ chip = ((address >> 9) & 0x7);
+ channel = ((address >> 5) & 0x5);
+ instruction = (address & 0x1f);