RawDataCollectorByLabel.cc

Go to the documentation of this file.

#include "EventFilter/RawDataCollector/src/RawDataCollectorByLabel.h"
#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h" 
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Provenance/interface/ProcessHistory.h" 
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <iostream>

using namespace edm;

RawDataCollectorByLabel::RawDataCollectorByLabel(const edm::ParameterSet& pset) {

  inputTags_ = pset.getParameter<std::vector<InputTag> >("RawCollectionList");
  verbose_ = pset.getUntrackedParameter<int>("verbose",0);

  inputTokens_.reserve(inputTags_.size());
  for(tag_iterator_t inputTag = inputTags_.begin(); inputTag != inputTags_.end(); ++inputTag ) {
    inputTokens_.push_back(consumes<FEDRawDataCollection>(*inputTag));
  }
  produces<FEDRawDataCollection>();
}

RawDataCollectorByLabel::~RawDataCollectorByLabel(){

}


void RawDataCollectorByLabel::produce(Event & e, const EventSetup& c){

 std::vector< Handle<FEDRawDataCollection> > rawData;
 rawData.reserve(inputTokens_.size());
 for(tok_iterator_t inputTok = inputTokens_.begin(); inputTok != inputTokens_.end(); ++inputTok ) {
   Handle<FEDRawDataCollection> input;
   if (e.getByToken(*inputTok,input)){
     rawData.push_back(input);
   }
   //else{     //skipping the inputtag requested. but this is a normal operation to bare data & MC. silent warning   }
 }

 std::auto_ptr<FEDRawDataCollection> producedData(new FEDRawDataCollection);

 for (unsigned int i=0; i< rawData.size(); ++i ) { 

   const FEDRawDataCollection *rdc=rawData[i].product();

   if ( verbose_ > 0 ) {
     std::cout << "\nRAW collection #" << i+1 << std::endl;
     std::cout << "branch name = " << rawData[i].provenance()->branchName() << std::endl;
     std::cout << "process index = " << rawData[i].provenance()->productID().processIndex() << std::endl;
   }

   for ( int j=0; j< FEDNumbering::MAXFEDID; ++j ) {
     const FEDRawData & fedData = rdc->FEDData(j);
     size_t size=fedData.size();

     if ( size > 0 ) {
       // this fed has data -- lets copy it
       if(verbose_ > 1) std::cout << "Copying data from FED #" << j << std::endl;
       FEDRawData & fedDataProd = producedData->FEDData(j);
       if ( fedDataProd.size() != 0 ) {
     if(verbose_ > 1) {
       std::cout << " More than one FEDRawDataCollection with data in FED ";
       std::cout << j << " Skipping the 2nd\n";
     }
     continue;
       } 
       fedDataProd.resize(size);
       unsigned char *dataProd=fedDataProd.data();
       const unsigned char *data=fedData.data();
       for ( unsigned int k=0; k<size; ++k ) {
     dataProd[k]=data[k];
       }
     }
   }
 }

 // Insert the new product in the event  
 e.put(producedData);  

}