MatacqDataFormatter Class Reference

#include <MatacqDataFormatter.h>

Public Member Functions
void	interpretRawData (const FEDRawData &data, EcalMatacqDigiCollection &matacqDigiCollection)
void	interpretRawData (const MatacqRawEvent &data, EcalMatacqDigiCollection &matacqDigiCollection)
	MatacqDataFormatter ()

Private Member Functions
void	printData (std::ostream &out, const MatacqRawEvent &event) const

Detailed Description

Class to interpret ECAL MATACQ raw data and produce the MATACQ digis. This class is used by the MatacqProducer module.

Author

: Ph. Gras (CEA/Saclay)

Definition at line 16 of file MatacqDataFormatter.h.

Constructor & Destructor Documentation

MatacqDataFormatter::MatacqDataFormatter ( )

inline

Definition at line 18 of file MatacqDataFormatter.h.

{};

Member Function Documentation

void MatacqDataFormatter::interpretRawData	(	const FEDRawData &	data,
		EcalMatacqDigiCollection &	matacqDigiCollection
	)

Callback method for decoding raw data

Parameters

data	raw data
matacqDigiCollection	[out] digi collection object to fill with the decoded data

Definition at line 18 of file MatacqDataFormatter.cc.

References gather_cfg::cout, FEDRawData::data(), mps_fire::i, SiStripPI::max, and FEDRawData::size().

Referenced by MatacqProducer::addMatacqData().

                                                                                                                 {
#if MATACQ_DEBUG
  cout << "****************************************************************\n";
  cout << "********************** MATACQ decoder **************************\n";
  cout << "****************************************************************\n";
  cout << "FEDRawData: \n";
  char oldPad = cout.fill('0');
  for (int i = 0; i < max(100, (int)data.size()); ++i) {
    cout << hex << setw(2) << (int)(data.data()[i]) << ((i + 1) % 8 ? " " : "\n");
  }
  cout.fill(oldPad);
  cout << "======================================================================\n";
#endif  //MATACQ_DEBUG defined
  interpretRawData(MatacqRawEvent(data.data(), data.size()), matacqDigiCollection);
}

void MatacqDataFormatter::interpretRawData	(	const MatacqRawEvent &	data,
		EcalMatacqDigiCollection &	matacqDigiCollection
	)

Callback method for decoding raw data

Parameters

data	raw data
matacqDigiCollection	[out] digi collection object to fill with the decoded data

Definition at line 34 of file MatacqDataFormatter.cc.

References EcalMatacqDigi::attenuation_dB(), edm::SortedCollection< T, SORT >::back(), EcalMatacqDigi::bxId(), filterCSVwithJSON::copy, gather_cfg::cout, EcalMatacqDigi::delayA(), relativeConstraints::empty, EcalMatacqDigi::emtcDelay(), EcalMatacqDigi::emtcPhase(), MatacqRawEvent::getAttenuation_dB(), MatacqRawEvent::getBxId(), MatacqRawEvent::getChannelData(), MatacqRawEvent::getDelayA(), MatacqRawEvent::getEmtcDelay(), MatacqRawEvent::getEmtcPhase(), MatacqRawEvent::getEventId(), MatacqRawEvent::getFreqGHz(), MatacqRawEvent::getLaserPower(), MatacqRawEvent::getMatacqDataFormatVersion(), MatacqRawEvent::getOrbitId(), MatacqRawEvent::getPostTrig(), MatacqRawEvent::getTimeStamp(), MatacqRawEvent::getTriggerType(), MatacqRawEvent::getTrigRec(), MatacqRawEvent::getTTrigPs(), MatacqRawEvent::getVernier(), EcalMatacqDigi::l1a(), EcalMatacqDigi::laserPower(), SiStripPI::max, hgc_digi::nSamples, EcalMatacqDigi::orbitId(), EcalMatacqDigi::postTrig(), edm::SortedCollection< T, SORT >::push_back(), submitPVValidationJobs::t, EcalMatacqDigi::timeStamp(), EcalMatacqDigi::triggerType(), EcalMatacqDigi::trigRec(), EcalMatacqDigi::vernier(), and relval_steps::version.

                                                                                           {
#if MATACQ_DEBUG
  printData(cout, matacq);
#endif  //MATACQ_DEBUG defined

  const double ns = 1.e-9;   //ns->s
  const double ps = 1.e-12;  //ps->s
  double ts = ns / matacq.getFreqGHz();
  double tTrig = matacq.getTTrigPs() < .5 * numeric_limits<int>::max() ? ps * matacq.getTTrigPs() : 999.;
  int version = matacq.getMatacqDataFormatVersion();

  vector<int16_t> samples;
  //FIXME: the interpretRawData method should fill an EcalMatacqDigiCollection
  //instead of an EcalMatacqDigi because Matacq channels are several.
  //In the meamtime copy only the first channel appearing in data:
  const vector<MatacqRawEvent::ChannelData>& chData = matacq.getChannelData();
  for (unsigned iCh = 0; iCh < chData.size(); ++iCh) {
    //copy time samples into a vector:
    samples.resize(chData[iCh].nSamples);
    copy(chData[iCh].samples, chData[iCh].samples + chData[iCh].nSamples, samples.begin());
    int chId = chData[iCh].chId;
    vector<int16_t> empty;
    EcalMatacqDigi matacqDigi(empty, chId, ts, version, tTrig);
#if (defined(ECAL_MATACQ_DIGI_VERS) && (ECAL_MATACQ_DIGI_VERS >= 2))
    matacqDigi.bxId(matacq.getBxId());
    matacqDigi.l1a(matacq.getEventId());
    matacqDigi.triggerType(matacq.getTriggerType());
    matacqDigi.orbitId(matacq.getOrbitId());
    matacqDigi.trigRec(matacq.getTrigRec());
    matacqDigi.postTrig(matacq.getPostTrig());
    matacqDigi.vernier(matacq.getVernier());
    matacqDigi.delayA(matacq.getDelayA());
    matacqDigi.emtcDelay(matacq.getEmtcDelay());
    matacqDigi.emtcPhase(matacq.getEmtcPhase());
    matacqDigi.attenuation_dB(matacq.getAttenuation_dB());
    matacqDigi.laserPower(matacq.getLaserPower());
    timeval t;
    matacq.getTimeStamp(t);
    matacqDigi.timeStamp(t);
#endif  //matacq digi version >=2
    matacqDigiCollection.push_back(matacqDigi);
    matacqDigiCollection.back().swap(samples);  //swap is more efficient than a copy
  }
}

void MatacqDataFormatter::printData ( std::ostream &  out, const MatacqRawEvent &  event  ) const

private

Definition at line 80 of file MatacqDataFormatter.cc.

References gpuClustering::adc, gather_cfg::cout, TauDecayModes::dec, MatacqRawEvent::getBxId(), MatacqRawEvent::getChannelCount(), MatacqRawEvent::getChannelData(), MatacqRawEvent::getDccErrors(), MatacqRawEvent::getDccLen(), MatacqRawEvent::getEventId(), MatacqRawEvent::getFedId(), MatacqRawEvent::getFov(), MatacqRawEvent::getFreqGHz(), MatacqRawEvent::getMatacqDataFormatVersion(), MatacqRawEvent::getRunNum(), MatacqRawEvent::getStatus(), MatacqRawEvent::getTimeStamp(), MatacqRawEvent::getTriggerType(), and mps_fire::i.

                                                                                    {
  cout << "FED id: " << hex << "0x" << matacq.getFedId() << dec << "\n";
  cout << "Event id (lv1): " << hex << "0x" << matacq.getEventId() << dec << "\n";
  cout << "FOV: " << hex << "0x" << matacq.getFov() << dec << "\n";
  cout << "BX id: " << hex << "0x" << matacq.getBxId() << dec << "\n";
  cout << "Trigger type: " << hex << "0x" << matacq.getTriggerType() << dec << "\n";
  cout << "DCC Length: " << matacq.getDccLen() << "\n";
  cout << "Run number: " << matacq.getRunNum() << "\n";
  cout << "Field 'DCC errors': " << hex << "0x" << matacq.getDccErrors() << dec << "\n";

  if (matacq.getStatus()) {
    cout << "Error in matacq data. Errot code: " << hex << "0x" << matacq.getStatus() << dec << "\n";
  }

  cout << "MATACQ data format version: " << matacq.getMatacqDataFormatVersion() << "\n";
  cout << "Sampling frequency: " << matacq.getFreqGHz() << " GHz\n";
  cout << "MATACQ channel count: " << matacq.getChannelCount() << "\n";
  time_t timeStamp = matacq.getTimeStamp();
  cout << "Data acquired on : " << ctime(&timeStamp);

  const vector<MatacqRawEvent::ChannelData>& channels = matacq.getChannelData();
  for (unsigned i = 0; i < channels.size(); ++i) {
    cout << "-------------------------------------- Channel " << channels[i].chId << ": " << setw(4)
         << channels[i].nSamples << " samples --------------------------------------\n";

    for (int iSample = 0; iSample < channels[i].nSamples; ++iSample) {
      MatacqRawEvent::int16le_t adc = (channels[i].samples)[iSample];
      cout << setw(4) << adc << ((iSample % 20 == 19) ? "\n" : " ");
    }
  }
  cout << "=================================================="
          "==================================================\n\n";
}