DTROS8FileReader.cc

Go to the documentation of this file.

#include <IORawData/DTCommissioning/plugins/DTROS8FileReader.h>
#include <IORawData/DTCommissioning/plugins/DTFileReaderHelpers.h>

#include <DataFormats/FEDRawData/interface/FEDHeader.h>
#include <DataFormats/FEDRawData/interface/FEDTrailer.h>
#include <DataFormats/FEDRawData/interface/FEDNumbering.h>

#include "DataFormats/Provenance/interface/EventID.h"
#include <DataFormats/Provenance/interface/Timestamp.h>
#include <DataFormats/FEDRawData/interface/FEDRawData.h>
#include <DataFormats/FEDRawData/interface/FEDRawDataCollection.h>

#include <FWCore/ParameterSet/interface/ParameterSet.h>
#include <FWCore/Utilities/interface/Exception.h>

#include <string>
#include <iosfwd>
#include <iostream>
#include <algorithm>
   
using namespace std;
using namespace edm;


DTROS8FileReader::DTROS8FileReader(const edm::ParameterSet& pset) : 
  runNum(1), eventNum(0) {
      
  //const string & filename = pset.getParameter<string>("fileName");  // Don't work, it produces an error saying that the DTNewROS8FileReader plugin doesn't exist!!!
 
  const string & filename = pset.getUntrackedParameter<string>("fileName");

  inputFile.open(filename.c_str());
  if( inputFile.fail() ) {
    throw cms::Exception("InputFileMissing") 
      << "DTROS8FileReader: the input file: " << filename <<" is not present";
  }

 produces<FEDRawDataCollection>();
}


DTROS8FileReader::~DTROS8FileReader(){
  inputFile.close();
}


int DTROS8FileReader::fillRawData(Event& e,
//                Timestamp& tstamp, 
                  FEDRawDataCollection*& data){
  EventID eID = e.id();
  data = new FEDRawDataCollection();

  try {
    

    if( checkEndOfFile() ) throw 1; 
    

    // Get the total number of words from the 1st word in the payload
    int numberOfWords=0;
    int nread = 0;
    nread = inputFile.read(dataPointer<int>( &numberOfWords ), ros8WordLenght);
    if ( nread<=0 )  throw 1;


    // Get the event data (all words but the 1st)
    int* eventData = new int[numberOfWords];
    nread = inputFile.read(dataPointer<int>( eventData + 1 ), (numberOfWords-1) * ros8WordLenght );
    if ( nread<=0 ) throw 1;
    

    // Check that the event data size corresponds to the 1st word datum 
    if ( eventData[numberOfWords-1] != numberOfWords ) {
      cout << "[DTROS8FileReader]: word counter mismatch exception: "
       << numberOfWords << " " << eventData[numberOfWords-1] << endl;
      throw 99;
    }

    // The header added by the local DAQ occupies 8 words, starting from the 2nd 
    int* head = eventData + 1;
    
    /* 
      Header word 0: run number
      Header word 1: spill number
      Header word 2: event number
      Header word 3: reserved
      Header word 4: ROS data offset
      Header word 5: PU data offset
      Header word 6: reserved
      Header word 7: reserved
    */

    // WARNING: the event number is reset at a new spill
    eID = EventID( head[0], 1U, head[1]*head[2]);

    // The pointer to the ROS payload (the 1st word being the ROS words counter)
    int* rosData = eventData + head[4];

    // The ROS payload size
    int eventDataSize = *rosData * ros8WordLenght;
    // It has to be a multiple of 8 bytes. if not, adjust the size of the FED payload
    int adjustment = (eventDataSize/4)%2 == 1 ? 4 : 0;   
    
    //if ( (eventDataSize/4)%2 ) adjustment = 4;


    // The FED ID is always the first in the DT range
    FEDRawData& fedRawData = data->FEDData( FEDNumbering::MINDTFEDID );
    fedRawData.resize(eventDataSize+adjustment);
    
    // I pass only the ROS data to the Event
    copy(reinterpret_cast<unsigned char*>(rosData), 
     reinterpret_cast<unsigned char*>(rosData) + eventDataSize, fedRawData.data());

    // needed to get rid of memory leaks (?)
    delete[] eventData;

    return true;
  }

  catch( int i ) {

    if ( i == 1 ){
      cout << "[DTROS8FileReader]: END OF FILE REACHED. "
           << "No information read for the requested event" << endl;
      delete data; data=0;
      return false;
    }    
    else {
      cout << "[DTROS8FileReader]: PROBLEM WITH EVENT INFORMATION ON THE FILE. "
           << "EVENT DATA READING FAILED  code= " << i << endl;
      delete data; data=0;
      return false;
    }

  }

}

void DTROS8FileReader::produce(Event&e, EventSetup const&es){
   edm::Handle<FEDRawDataCollection> rawdata;
   FEDRawDataCollection *fedcoll = 0;
   fillRawData(e,fedcoll);
   std::auto_ptr<FEDRawDataCollection> bare_product(fedcoll);
   e.put(bare_product);
 }

bool DTROS8FileReader::checkEndOfFile(){

  bool retval=false;
  if ( inputFile.eof() ) retval=true;
  return retval;

}