#include <DTDDUFileReader.h>

Inheritance diagram for DTDDUFileReader:

Public Member Functions
virtual bool	checkEndOfFile ()

uint64_t	dmaUnpack (bool &isData, int &nread)
	pre-unpack the data if read via DMA More...

	DTDDUFileReader (const edm::ParameterSet &pset)
	Constructor. More...

virtual int	fillRawData (edm::EventID &eID, edm::Timestamp &tstamp, FEDRawDataCollection *&data)
	Generate and fill FED raw data for a full event. More...

bool	isHeader (uint64_t word, bool dataTag)
	check for a 64 bits word to be a DDU header More...

bool	isTrailer (uint64_t word, bool dataTag, int wordCount)
	check for a 64 bits word to be a DDU trailer More...

void	swap (uint64_t &word)
	swapping the lsBits with the msBits More...

virtual	~DTDDUFileReader ()
	Destructor. More...

Public Member Functions inherited from DaqBaseReader
	DaqBaseReader ()

virtual void	setRunNumber (edm::RunNumber_t runNumber)
	set the run number More...

virtual	~DaqBaseReader ()

Private Attributes
int	dduID

edm::EventNumber_t	eventNumber

RawFile	inputFile

int	numberOfHeaderWords

bool	readFromDMA

edm::RunNumber_t	runNumber

int	skipEvents

Static Private Attributes
static const int	dduWordLength = 8

Detailed Description

Read DT ROS8 raw data files

Date:: 2010/02/03 16:58:24

Revision:: 1.11

Author: M. Zanetti - INFN Padova

Definition at line 20 of file DTDDUFileReader.h.

Constructor & Destructor Documentation

DTDDUFileReader::DTDDUFileReader ( const edm::ParameterSet & pset )

Constructor.

Definition at line 32 of file DTDDUFileReader.cc.

References gather_cfg::cout, edm::hlt::Exception, RawFile::fail(), lut2db_cfg::filename, edm::ParameterSet::getUntrackedParameter(), RawFile::ignore(), inputFile, numberOfHeaderWords, RawFile::open(), RawFile::read(), readFromDMA, runNumber, and skipEvents.

                                                             : 
   runNumber(1), eventNumber(1) {
 
   const string & filename = pset.getUntrackedParameter<string>("fileName");
 
   readFromDMA = pset.getUntrackedParameter<bool>("readFromDMA",true);
   numberOfHeaderWords = pset.getUntrackedParameter<int>("numberOfHeaderWords",10);
   skipEvents = pset.getUntrackedParameter<int>("skipEvents",0);  
 
   inputFile.open(filename.c_str());
   if( inputFile.fail() ) {
     throw cms::Exception("InputFileMissing") 
       << "[DTDDUFileReader]: the input file: " << filename <<" is not present";
   } else {
     cout << "DTDDUFileReader: DaqSource file '" << filename << "' was succesfully opened" << endl;
   }
 
   uint32_t runNumber_tmp;
   inputFile.read(dataPointer<uint32_t>( &runNumber_tmp ), 4);
   runNumber = runNumber_tmp;
 
   inputFile.ignore(4*(numberOfHeaderWords-1));
   
   if (skipEvents) { 
     cout<<""<<endl;
     cout<<"   Dear user, pleas be patient, "<<skipEvents<<" are being skipped .."<<endl;
     cout<<""<<endl;
   }
 
 }

DTDDUFileReader::~DTDDUFileReader ( )

virtual

Destructor.

Definition at line 64 of file DTDDUFileReader.cc.

References RawFile::close(), and inputFile.

                                  {
       inputFile.close();
 }

Member Function Documentation

bool DTDDUFileReader::checkEndOfFile ( )

virtual

Definition at line 235 of file DTDDUFileReader.cc.

References RawFile::eof(), and inputFile.

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

uint64_t DTDDUFileReader::dmaUnpack	(	bool &	isData,
		int &	nread
	)

pre-unpack the data if read via DMA

Definition at line 178 of file DTDDUFileReader.cc.

References data, inputFile, and RawFile::read().

Referenced by fillRawData().

                                                               {
   
   uint64_t dduWord = 0;
 
   uint32_t td[4];
   // read 4 32-bits word from the file;
   nread = inputFile.read(dataPointer<uint32_t>( &td[0] ), 4);
   nread += inputFile.read(dataPointer<uint32_t>( &td[1] ), 4);
   nread += inputFile.read(dataPointer<uint32_t>( &td[2] ), 4);
   nread += inputFile.read(dataPointer<uint32_t>( &td[3] ), 4);
 
   uint32_t data[2] = {0, 0};
   // adjust 4 32-bits words  into 2 32-bits words
   data[0] |= td[3] & 0x3ffff;
   data[0] |= (td[2] << 18 ) & 0xfffc0000;
   data[1] |= (td[2] >> 14 ) & 0x0f;
   data[1] |= (td[1] << 4 ) & 0x3ffff0;
   data[1] |= (td[0] << 22 ) & 0xffc00000;
 
   isData = ( td[0] >> 10 ) & 0x01;
 
   // push_back to a 64 word
   dduWord = (uint64_t(data[1]) << 32) | data[0];
 
   return dduWord;
 }

int DTDDUFileReader::fillRawData	(	edm::EventID &	eID,
		edm::Timestamp &	tstamp,
		FEDRawDataCollection *&	data
	)

virtual

Generate and fill FED raw data for a full event.

Implements DaqBaseReader.

Definition at line 68 of file DTDDUFileReader.cc.

References filterCSVwithJSON::copy, gather_cfg::cout, FEDRawData::data(), data, dduID, dduWordLength, dmaUnpack(), eventNumber, FEDRawDataCollection::FEDData(), inputFile, isHeader(), isTrailer(), FEDNumbering::MINDTFEDID, RawFile::read(), readFromDMA, FEDRawData::resize(), runNumber, skipEvents, and swap().

                                              {
   data = new FEDRawDataCollection();
 
   vector<uint64_t> eventData;
   size_t estimatedEventDimension = 102400; // dimensione hardcoded
   eventData.reserve(estimatedEventDimension); 
   uint64_t word = 0;
 
   bool haederTag = false;
   bool dataTag = true;
   bool headerAlreadyFound = false;
   
   int wordCount = 0;
   
   // getting the data word by word from the file
   // do it until you get the DDU trailer
   while ( !isTrailer(word, dataTag, wordCount) ) {
     //while ( !isTrailer(word) ) { 
     
     if (readFromDMA) {
       int nread;
       word = dmaUnpack(dataTag,nread);
       if ( nread<=0 ) {
     cout<<"[DTDDUFileReader]: ERROR! no more words and failed to get the trailer"<<endl;
     delete data; data=0;
     return false;
       }
     }
     
     else {
       int nread = inputFile.read(dataPointer<uint64_t>( &word ), dduWordLength);
       dataTag = false;
       if ( nread<=0 ) {
     cout<<"[DTDDUFileReader]: ERROR! failed to get the trailer"<<endl;
     delete data; data=0;
     return false;
       }
     }
     
     // get the DDU header
     if (!headerAlreadyFound) 
       if ( isHeader(word,dataTag)) {
     headerAlreadyFound=true; 
     haederTag=true;
       }
     
     // from now on fill the eventData with the ROS data
     if (haederTag) {
       
       if (readFromDMA) {
     // swapping only the 8 byte words
     if (dataTag) {
       swap(word);
     } // WARNING also the ddu status words have been swapped!
     // Control the correct interpretation in DDUUnpacker
       }
       
       eventData.push_back(word);
       wordCount++;
     }
     
   } 
   
   //     FEDTrailer candidate(reinterpret_cast<const unsigned char*>(&word));
   //     cout<<"EventSize: pushed back "<<eventData.size()
   //    <<";  From trailer "<<candidate.lenght()<<endl;
   
   // next event reading will start with meaningless trailer+header from DTLocalDAQ
   // those will be skipped automatically when seeking for the DDU header
   //if (eventData.size() > estimatedEventDimension) throw 2;
   
   // Eventually skipping events
   if ((int)eventNumber >= skipEvents) {
 
     // Setting the Event ID
     eID = EventID( runNumber, 1U, eventNumber);
     
     // eventDataSize = (Number Of Words)* (Word Size)
     int eventDataSize = eventData.size()*dduWordLength;
     
 
     if ( dduID<770 || dduID > 775 ) {
       cout<<"[DTDDUFileReader]: ERROR. DDU ID out of range. DDU id="<<dduID<<endl;
       // The FED ID is always the first in the DT range
       dduID = FEDNumbering::MINDTFEDID;
     } 
     FEDRawData& fedRawData = data->FEDData( dduID );
     fedRawData.resize(eventDataSize);
     
     copy(reinterpret_cast<unsigned char*>(&eventData[0]),
      reinterpret_cast<unsigned char*>(&eventData[0]) + eventDataSize, fedRawData.data());
 
   }
 
   return true;
     
 }

bool DTDDUFileReader::isHeader	(	uint64_t	word,
		bool	dataTag
	)

check for a 64 bits word to be a DDU header

Definition at line 205 of file DTDDUFileReader.cc.

References FEDHeader::check(), dduID, eventNumber, and FEDHeader::sourceID().

Referenced by fillRawData().

                                                           {
 
   bool it_is = false;
   FEDHeader candidate(reinterpret_cast<const unsigned char*>(&word));
   if ( candidate.check() ) {
     // if ( candidate.check() && !dataTag) {
     it_is = true;
     dduID = candidate.sourceID();
     eventNumber++;
   }
  
   return it_is;
 }

bool DTDDUFileReader::isTrailer	(	uint64_t	word,
		bool	dataTag,
		int	wordCount
	)

check for a 64 bits word to be a DDU trailer

Definition at line 220 of file DTDDUFileReader.cc.

References FEDTrailer::check(), and FEDTrailer::lenght().

Referenced by fillRawData().

                                                                           {
 
   bool it_is = false;
   FEDTrailer candidate(reinterpret_cast<const unsigned char*>(&word));
   if ( candidate.check() ) {
     //  if ( candidate.check() && !dataTag) {
     //cout<<"[DTDDUFileReader] "<<wordCount<<" - "<<candidate.lenght()<<endl;
     if ( wordCount == candidate.lenght())
       it_is = true;
   }
  
   return it_is;
 }

void DTDDUFileReader::swap ( uint64_t & word )

swapping the lsBits with the msBits

Definition at line 168 of file DTDDUFileReader.cc.

References twoNibble64::lsBits, and twoNibble64::msBits.

Referenced by fillRawData().

                                           {
   
   twoNibble64* newWorld = reinterpret_cast<twoNibble64*>(&word);
 
   uint32_t msBits_tmp = newWorld->msBits;
   newWorld->msBits = newWorld->lsBits;
   newWorld->lsBits = msBits_tmp;
 }