EcalFEDWithCRCErrorProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    EcalFEDWithCRCErrorProducer
// Class:      EcalFEDWithCRCErrorProducer
// 
//
// Original Author:  Giovanni FRANZONI
//         Created:  Tue Jan 22 13:55:00 CET 2008
// $Id: EcalFEDWithCRCErrorProducer.cc,v 1.5 2010/10/05 13:39:30 vlimant Exp $
//
//


// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Framework/interface/EDProducer.h"

#include <string>
#include <iostream>
#include <vector>
#include <iomanip>

#include <DataFormats/FEDRawData/interface/FEDRawData.h>
#include <DataFormats/FEDRawData/interface/FEDRawDataCollection.h>
#include <DataFormats/FEDRawData/interface/FEDNumbering.h>
//
// class declaration
//

class EcalFEDWithCRCErrorProducer : public edm::EDProducer 
{
public:
  explicit EcalFEDWithCRCErrorProducer(const edm::ParameterSet&);
  ~EcalFEDWithCRCErrorProducer();
  
private:
  virtual void produce(edm::Event&, const edm::EventSetup&);
  
  // ----------member data ---------------------------
  
  edm::InputTag     DataLabel_;
  std::vector<int> fedUnpackList_;
  bool writeAllEcalFEDs_;
};



//
// constructors and destructor
//
EcalFEDWithCRCErrorProducer::EcalFEDWithCRCErrorProducer(const edm::ParameterSet& iConfig)
{
  //now do what ever initialization is needed

  DataLabel_     = iConfig.getParameter<edm::InputTag>("InputLabel");
  fedUnpackList_ = iConfig.getUntrackedParameter< std::vector<int> >("FEDs", std::vector<int>());
  writeAllEcalFEDs_ = iConfig.getUntrackedParameter<bool >("writeAllEcalFED", false);
  if (fedUnpackList_.empty()) 
    for (int i=FEDNumbering::MINECALFEDID; i<=FEDNumbering::MAXECALFEDID; i++)
      fedUnpackList_.push_back(i);

  produces< FEDRawDataCollection >();
}


EcalFEDWithCRCErrorProducer::~EcalFEDWithCRCErrorProducer()
{
 
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called on each new Event  ------------
void
EcalFEDWithCRCErrorProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  using namespace edm;

  edm::Handle<FEDRawDataCollection> rawdata;  
  iEvent.getByLabel(DataLabel_,rawdata);

  std::auto_ptr<FEDRawDataCollection> producedData(new FEDRawDataCollection);
  // get fed raw data and SM id

  // loop over FEDS
  for (std::vector<int>::const_iterator i=fedUnpackList_.begin(); i!=fedUnpackList_.end(); i++) 
    {

      // get fed raw data and SM id
      const FEDRawData & fedData = rawdata->FEDData(*i);
      int length = fedData.size()/sizeof(uint64_t);

      //    LogDebug("EcalRawToDigi") << "raw data length: " << length ;
      //if data size is not null interpret data
      if ( length >= 1 )
    {
            uint64_t * pData = (uint64_t *)(fedData.data());
        //When crc error is found return true
        uint64_t * fedTrailer = pData + (length - 1);
        bool crcError = (*fedTrailer >> 2 ) & 0x1; 
        // this fed has data -- lets copy it
        if (writeAllEcalFEDs_ || crcError)
          {
        FEDRawData & fedDataProd = producedData->FEDData(*i);
        if ( fedDataProd.size() != 0 ) 
          {
            //                std::cout << " More than one FEDRawDataCollection with data in FED ";
            //                std::cout << j << " Skipping the 2nd\n";
            continue;
          }
        fedDataProd.resize(fedData.size());
        unsigned char *dataProd=fedDataProd.data();
        const unsigned char *data=fedData.data();
        for ( unsigned int k=0; k<fedData.size(); ++k ) 
          {
            dataProd[k]=data[k];
          }
          }
    }
    }
  
  iEvent.put(producedData);  
}

//define this as a plug-in
DEFINE_FWK_MODULE(EcalFEDWithCRCErrorProducer);