#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, unsigned int wordCount)
	check for a 64 bits word to be a DDU trailer More...

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

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

	~DTDDUFileReader () override
	Destructor. More...

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

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

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

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 31 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 ( )

override

Destructor.

Definition at line 60 of file DTDDUFileReader.cc.

References RawFile::close(), and inputFile.

60 { inputFile.close(); }

DTDDUFileReader::inputFile

RawFile inputFile

Definition: DTDDUFileReader.h:53

RawFile::close

int close()

Close file if necessary.

Definition: RawFile.cc:52

Member Function Documentation

bool DTDDUFileReader::checkEndOfFile ( )

virtual

Definition at line 228 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 176 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 62 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 = nullptr;
         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 = nullptr;
         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.fragmentLength()<<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 202 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,
		unsigned int	wordCount
	)

check for a 64 bits word to be a DDU trailer

Definition at line 215 of file DTDDUFileReader.cc.

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

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.fragmentLength()<<endl;
     if (wordCount == candidate.fragmentLength())
       it_is = true;
   }
 
   return it_is;
 }

void DTDDUFileReader::produce	(	edm::Event &	e,
		edm::EventSetup const &	es
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 160 of file DTDDUFileReader.cc.

References fillRawData(), eostools::move(), and edm::Event::put().

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

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;
 }