/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/EventFilter/CastorRawToDigi/src/CastorCtdcUnpacker.cc

Go to the documentation of this file.

#include "EventFilter/CastorRawToDigi/interface/CastorCtdcUnpacker.h"
#include "EventFilter/CastorRawToDigi/interface/CastorCTDCHeader.h"
#include "EventFilter/CastorRawToDigi/interface/CastorCORData.h"
#include "DataFormats/HcalDetId/interface/HcalOtherDetId.h"
#include "DataFormats/HcalDigi/interface/HcalQIESample.h"
#include "DataFormats/HcalDigi/interface/CastorDataFrame.h"
#include "DataFormats/HcalDigi/interface/HcalTriggerPrimitiveSample.h"
#include <iostream>
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace std;

CastorCtdcUnpacker::CastorCtdcUnpacker(int sourceIdOffset, int beg, int end) : sourceIdOffset_(sourceIdOffset)
{
 if ( beg >= 0 && beg <= CastorDataFrame::MAXSAMPLES -1 ) {
        startSample_ = beg;
 } else {
        startSample_ = 0;
 }
 if ( end >= 0 && end <= CastorDataFrame::MAXSAMPLES -1 && end >= beg ) {
        endSample_ = end;
 } else {
        endSample_ = CastorDataFrame::MAXSAMPLES -1;
 }
}


void CastorCtdcUnpacker::unpack(const FEDRawData& raw, const CastorElectronicsMap& emap,
                          CastorRawCollections& colls, HcalUnpackerReport& report) {

  if (raw.size()<16) {
    edm::LogWarning("Invalid Data") << "Empty/invalid DCC data, size = " << raw.size();
    return;
  }

  // get the CTDC header
  const CastorCTDCHeader* ctdcHeader=(const CastorCTDCHeader*)(raw.data());
  int ctdcid=ctdcHeader->getSourceId()-sourceIdOffset_;

  // space for unpacked data from one COR
  std::vector<unsigned short> precdata(CastorCORData::CHANNELS_PER_SPIGOT*CastorCORData::MAXIMUM_SAMPLES_PER_CHANNEL);
  std::vector<unsigned short> trigdata(CastorCORData::CHANNELS_PER_SPIGOT*CastorCORData::MAXIMUM_SAMPLES_PER_CHANNEL);
  std::vector<unsigned char>  preclen(CastorCORData::CHANNELS_PER_SPIGOT);
  std::vector<unsigned char>  triglen(CastorCORData::CHANNELS_PER_SPIGOT);
  
  // walk through the COR data...
  CastorCORData cor;
  for (int spigot=0; spigot<CastorCTDCHeader::SPIGOT_COUNT; spigot++) {
    if (!ctdcHeader->getSpigotPresent(spigot)) continue;

    int retval=ctdcHeader->getSpigotData(spigot,cor,raw.size());
    if (retval!=0) {
      if (retval==-1) {
        edm::LogWarning("Invalid Data") << "Invalid COR data (data beyond payload size) observed on spigot " << spigot << " of CTDC with source id " << ctdcHeader->getSourceId();
        report.countSpigotFormatError();
      }
      continue;
    }
    // check
    if (!cor.check()) {
      edm::LogWarning("Invalid Data") << "Invalid COR data observed on spigot " << spigot << " of CTDC with source id " << ctdcHeader->getSourceId();
      report.countSpigotFormatError();
      continue;
    }
    if (cor.isHistogramEvent()) {
      edm::LogWarning("Invalid Data") << "Histogram data passed to non-histogram unpacker on spigot " << spigot << " of CTDC with source id " << ctdcHeader->getSourceId();
      continue;

    }
    // calculate "real" number of presamples
    int nps=cor.getNPS()-startSample_;

    // new: do not use get pointers .. instead make use of CastorCORData::unpack   

     cor.unpack(&(preclen[0]),&(precdata[0]),
                            &(triglen[0]),&(trigdata[0]));
 
    
   int ichan;
   for ( ichan = 0; ichan<CastorCORData::CHANNELS_PER_SPIGOT; ichan++) {
           if ( preclen[ichan] == 0 || preclen[ichan] & 0xc0 ) continue;
                int fiber = ichan/3;
                int fiberchan = ichan%3;
        // lookup the right channel
            CastorElectronicsId partialEid(fiberchan,fiber+1,spigot,ctdcid);
            // does this partial id exist?
            CastorElectronicsId eid;
            HcalGenericDetId did;
            bool found;
            found = emap.lookup(partialEid,eid,did);

                if (found) {                    
                        CastorDataFrame digi = CastorDataFrame(HcalCastorDetId(did));
                        // set parameters - presamples
                        digi.setPresamples(nps);

                        int ntaken = 0;
                        for ( int sample = startSample_; sample <= endSample_; sample++ ) {
                                        digi.setSample(ntaken,precdata[ichan*CastorCORData::MAXIMUM_SAMPLES_PER_CHANNEL+sample]);
                                        ntaken++;
                        }
                        digi.setSize(ntaken);
                        colls.castorCont->push_back(digi);                      
                } else {
                        report.countUnmappedDigi();
                        if (unknownIds_.find(partialEid)==unknownIds_.end()) {
                          edm::LogWarning("CASTOR") << "CastorCtdcUnpacker: No match found for electronics partialEid :" << partialEid;
                          unknownIds_.insert(partialEid);
                    }
                }
    }

  }
}