#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::Event &e, 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...

virtual void	produce (edm::Event &, edm::EventSetup const &)

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

virtual	~DTDDUFileReader ()
	Destructor. More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 22 of file DTDDUFileReader.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 33 of file DTDDUFileReader.cc.

References gather_cfg::cout, 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=1;
   inputFile.read(dataPointer<uint32_t>( &runNumber_tmp ), 4);
   runNumber = runNumber_tmp;
 
   inputFile.ignore(4*(numberOfHeaderWords-1));
   
   if (skipEvents) { 
     cout<<""<<endl;
     cout<<"   Dear user, please be patient, "<<skipEvents<<" are being skipped .."<<endl;
     cout<<""<<endl;
   }
 
   produces<FEDRawDataCollection>();
 }

DTDDUFileReader::~DTDDUFileReader ( )

virtual

Destructor.

Definition at line 66 of file DTDDUFileReader.cc.

References RawFile::close(), and inputFile.

                                  {
       inputFile.close();
 }

Member Function Documentation

bool DTDDUFileReader::checkEndOfFile ( )

virtual

Definition at line 246 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 189 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::Event &	e,
		FEDRawDataCollection *&	data
	)

virtual

Generate and fill FED raw data for a full event.

Definition at line 70 of file DTDDUFileReader.cc.

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

Referenced by produce().

                                              {
   EventID eID = e.id();
   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 216 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 231 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::produce	(	edm::Event &	e,
		edm::EventSetup const &	es
	)

virtual

Implements edm::EDProducer.

Definition at line 171 of file DTDDUFileReader.cc.

References fillRawData(), edm::Event::put(), and lumiPlot::rawdata.

                                                          {
    edm::Handle<FEDRawDataCollection> rawdata;
    FEDRawDataCollection *fedcoll = 0;
    fillRawData(e,fedcoll);
    std::auto_ptr<FEDRawDataCollection> bare_product(fedcoll);
    e.put(bare_product);
  }

void DTDDUFileReader::swap ( uint64_t & word )

swapping the lsBits with the msBits

Definition at line 179 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;
 }