/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/RecoTBCalo/HcalTBObjectUnpacker/src/HcalTBSlowDataUnpacker.cc

Go to the documentation of this file.

#include "RecoTBCalo/HcalTBObjectUnpacker/interface/HcalTBSlowDataUnpacker.h"
#include "TBDataFormats/HcalTBObjects/interface/HcalTBRunData.h"
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include <iostream>
#include <string>
#include <map>
#include "FWCore/Utilities/interface/Exception.h"

using namespace std;

// structure for Slow Data
struct xdaqSlowDataFormat {
  uint32_t cdfHeader[4];
  uint16_t n_doubles;
  uint16_t n_strings;
  uint16_t key_length;
  uint16_t string_value_length;
  char     start_of_data; // see below
  // char[n_doubles*key_length] doubles names
  // double[n_doubles] doubles values
  // char[n_strings*key_length] strings names
  // char[n_strings*string_value_length] strings values
  // xdaqCommonDataFormatTrailer
};

namespace hcaltb {

  void HcalTBSlowDataUnpacker::unpack(const FEDRawData&  raw,
                                      HcalTBRunData&          htbrd,
                                      HcalTBEventPosition&    htbep) {

  if (raw.size()<3*8) {
    throw cms::Exception("Missing Data") << "No data in the slow data block";
  }

    const struct xdaqSlowDataFormat *sd =
      (const struct xdaqSlowDataFormat *)(raw.data());

    map<string,double> sd_dblmap;
    map<string,string> sd_strmap;

#ifdef DEBUG
    cout << "#doubles = "   << sd->n_doubles << endl;;
    cout << "#strings = "   << sd->n_strings << endl;
    cout << "key_length = " << sd->key_length << endl;
    cout << "string_value_length = " << sd->string_value_length << endl;
#endif

    // List of doubles:

    const char   *keyptr = &sd->start_of_data;
    const double *valptr =
      (const double *)(&sd->start_of_data + sd->n_doubles*sd->key_length);

    for (int i=0; i<sd->n_doubles; i++) {
#ifdef DEBUG
      cout << keyptr << " = " << *valptr << endl;
#endif
      sd_dblmap[keyptr] = *valptr;
      keyptr += sd->key_length;
      valptr++;
    }

    // List of strings:

    keyptr = (const char *)valptr;
    const char *strptr = (keyptr + sd->n_strings*sd->key_length);

    for (int i=0; i<sd->n_strings; i++) {
#ifdef DEBUG
      cout << keyptr << " = " << strptr << endl;
#endif
      sd_strmap[keyptr] = strptr;
      keyptr += sd->key_length;
      strptr += sd->string_value_length;
    }

    // Now fill the input objects:
    htbrd.setRunData(sd_strmap["RunType"].c_str(),
                     sd_strmap["Beam.Mode"].c_str(),
                     sd_dblmap["Beam.Energy"]);

    htbep.setHFtableCoords(sd_dblmap["HFTable.X"],
                           sd_dblmap["HFTable.Y"],
                           sd_dblmap["HFTable.V"]);

    htbep.setHBHEtableCoords(sd_dblmap["Table.Eta"],
                             sd_dblmap["Table.Phi"]);
  }
}