Inheritance diagram for MatacqProducer:

Classes
struct	MatacqEventId

class	PosEstimator

struct	stats_t

Public Types
enum	calibTrigType_t { laserType = 4, ledType = 5, tpType = 6, pedType = 7 }

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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

Public Member Functions
	MatacqProducer (const edm::ParameterSet &params)

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

	~MatacqProducer () override

Public Member Functions inherited from edm::one::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDProducerBase () 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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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

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

std::vector< ESResolverIndex > 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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef FILE *	FILE_t

typedef off_t	filepos_t

Private Member Functions
void	addMatacqData (edm::Event &event)

int	getCalibTriggerType (edm::Event &ev) const

bool	getMatacqEvent (uint32_t runNumber, int32_t orbitId, bool fileChange)

bool	getMatacqFile (uint32_t runNumber, uint32_t orbitId, bool *fileChange=nullptr)

uint32_t	getOrbitId (edm::Event &ev) const

bool	getOrbitRange (uint32_t &firstOrb, uint32_t &lastOrb)

uint32_t	getRunNumber (edm::Event &ev) const

void	loadOrbitOffset ()

bool	mcheck (const std::string &name)

void	mclose ()

bool	meof ()

bool	misOpened ()

bool	mopen (const std::string &name)

bool	mread (char buf, size_t n, const char mess=nullptr, bool peek=false)

bool	mrewind ()

bool	mseek (filepos_t offset, int whence=SEEK_SET, const char *mess=nullptr)

bool	msize (filepos_t &s)

bool	mtell (filepos_t &pos)

void	newRun (int prevRun, int newRun)

Static Private Member Functions
static std::string	runSubDir (uint32_t runNumber)

Private Attributes
std::vector< unsigned char >	data_

std::string	digiInstanceName_

bool	disabled_

bool	doOrbitOffset_

int	eventSkipCounter_

int	fastRetrievalThresh_

std::vector< std::string >	fileNames_

MatacqDataFormatter	formatter_

bool	ignoreTriggerType_

FILE_t	inFile_

std::string	inFileName_

edm::InputTag	inputRawCollection_

edm::EDGetTokenT< FEDRawDataCollection >	inputRawCollectionToken_

int32_t	lastOrb_

std::ofstream	logFile_

std::string	logFileName_

bool	logTiming_

MatacqRawEvent	matacq_

bool	mergeRaw_

int	onErrorDisablingEvtCnt_

uint32_t	openedFileRunNumber_

std::map< uint32_t, uint32_t >	orbitOffset_

std::string	orbitOffsetFile_

PosEstimator	posEstim_

bool	produceDigis_

bool	produceRaw_

std::string	rawInstanceName_

uint32_t	runNumber_

timeval	startTime_

struct MatacqProducer::stats_t	stats_

std::ofstream	timeLog_

std::string	timeLogFile_

timeval	timer_

bool	timing_

int	verbosity_

Static Private Attributes
static const int	bufferSize = 30000

static const int	matacqFedId_ = 655

static const int	orbitTolerance_ = 80

static const stats_t	stats_init = {0, 0, 0}

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

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

Definition at line 37 of file MatacqProducer.cc.

Member Typedef Documentation

◆ FILE_t

typedef FILE* MatacqProducer::FILE_t

private

Definition at line 47 of file MatacqProducer.cc.

◆ filepos_t

typedef off_t MatacqProducer::filepos_t

private

Definition at line 46 of file MatacqProducer.cc.

Member Enumeration Documentation

◆ calibTrigType_t

enum MatacqProducer::calibTrigType_t

Enumerator
laserType
ledType
tpType
pedType

Definition at line 39 of file MatacqProducer.cc.

39 { laserType = 4, ledType = 5, tpType = 6, pedType = 7 };

MatacqProducer::ledType

Definition: MatacqProducer.cc:39

MatacqProducer::pedType

Definition: MatacqProducer.cc:39

MatacqProducer::tpType

Definition: MatacqProducer.cc:39

MatacqProducer::laserType

Definition: MatacqProducer.cc:39

Constructor & Destructor Documentation

◆ MatacqProducer()

MatacqProducer::MatacqProducer ( const edm::ParameterSet & params )

explicit

Constructor

Parameters

params seletive readout parameters

Definition at line 372 of file MatacqProducer.cc.

References gather_cfg::cout, digiInstanceName_, doOrbitOffset_, Exception, inputRawCollectionToken_, loadOrbitOffset(), logFile_, logFileName_, logTiming_, now(), orbitOffsetFile_, MillePedeFileConverter_cfg::out, submitPVValidationJobs::params, posEstim_, produceDigis_, produceRaw_, rawInstanceName_, startTime_, timeLog_, timeLogFile_, timer_, MatacqProducer::PosEstimator::verbosity(), and verbosity_.

     : fileNames_(params.getParameter<std::vector<std::string> >("fileNames")),
       digiInstanceName_(params.getParameter<string>("digiInstanceName")),
       rawInstanceName_(params.getParameter<string>("rawInstanceName")),
       timing_(params.getUntrackedParameter<bool>("timing", false)),
       disabled_(params.getParameter<bool>("disabled")),
       verbosity_(params.getUntrackedParameter<int>("verbosity", 0)),
       produceDigis_(params.getParameter<bool>("produceDigis")),
       produceRaw_(params.getParameter<bool>("produceRaw")),
       inputRawCollection_(params.getParameter<edm::InputTag>("inputRawCollection")),
       mergeRaw_(params.getParameter<bool>("mergeRaw")),
       ignoreTriggerType_(params.getParameter<bool>("ignoreTriggerType")),
       matacq_(nullptr, 0),
       inFile_(nullptr),
       data_(bufferSize),
       openedFileRunNumber_(0),
       lastOrb_(0),
       fastRetrievalThresh_(0),
       orbitOffsetFile_(params.getUntrackedParameter<std::string>("orbitOffsetFile", "")),
       inFileName_(""),
       stats_(stats_init),
       logFileName_(params.getUntrackedParameter<std::string>("logFileName", "matacqProducer.log")),
       eventSkipCounter_(0),
       onErrorDisablingEvtCnt_(params.getParameter<int>("onErrorDisablingEvtCnt")),
       timeLogFile_(params.getUntrackedParameter<std::string>("timeLogFile", "")),
       runNumber_(0) {
   if (verbosity_ >= 4)
     cout << "[Matacq " << now() << "] in MatacqProducer ctor" << endl;
 
   gettimeofday(&timer_, nullptr);
 
   if (!timeLogFile_.empty()) {
     timeLog_.open(timeLogFile_.c_str());
     if (timeLog_.fail()) {
       cout << "[LaserSorter " << now() << "] "
            << "Failed to open file " << timeLogFile_ << " to log timing.\n";
       logTiming_ = false;
     } else {
       logTiming_ = true;
     }
   }
 
   posEstim_.verbosity(verbosity_);
 
   logFile_.open(logFileName_.c_str(), ios::app | ios::out);
 
   if (logFile_.bad()) {
     throw cms::Exception("FileOpen") << "Failed to open file " << logFileName_ << " for logging.\n";
   }
 
   inputRawCollectionToken_ = consumes<FEDRawDataCollection>(params.getParameter<InputTag>("inputRawCollection"));
 
   if (produceDigis_) {
     if (verbosity_ > 0)
       cout << "[Matacq " << now()
            << "] registering new "
               "EcalMatacqDigiCollection product with instance name '"
            << digiInstanceName_ << "'\n";
     produces<EcalMatacqDigiCollection>(digiInstanceName_);
   }
 
   if (produceRaw_) {
     if (verbosity_ > 0)
       cout << "[Matacq " << now()
            << "] registering new FEDRawDataCollection "
               "product with instance name '"
            << rawInstanceName_ << "'\n";
     produces<FEDRawDataCollection>(rawInstanceName_);
   }
 
   startTime_.tv_sec = startTime_.tv_usec = 0;
   if (!orbitOffsetFile_.empty()) {
     doOrbitOffset_ = true;
     loadOrbitOffset();
   } else {
     doOrbitOffset_ = false;
   }
   if (verbosity_ >= 4)
     cout << "[Matacq " << now() << "] exiting MatacqProducer ctor" << endl;
 }

◆ ~MatacqProducer()

MatacqProducer::~MatacqProducer ( )

override

Destructor

Definition at line 1131 of file MatacqProducer.cc.

References gather_cfg::cout, logTiming_, mclose(), now(), startTime_, and submitPVValidationJobs::t.

                                 {
   mclose();
   timeval t;
   gettimeofday(&t, nullptr);
   if (logTiming_ && startTime_.tv_sec != 0) {
     //not using logger, to allow timing with different logging options
     cout << "[Matacq " << now()
          << "] Time elapsed between first event and "
             "destruction of MatacqProducer: "
          << ((t.tv_sec - startTime_.tv_sec) * 1. + (t.tv_usec - startTime_.tv_usec) * 1.e-6) << "s\n";
   }
 }

Member Function Documentation

◆ addMatacqData()

void MatacqProducer::addMatacqData ( edm::Event & event )

private

Add matacq digi to the event

Parameters

event	the event
digiInstanceName_	name to give to the matacq digi instance

Definition at line 484 of file MatacqProducer.cc.

References bufferSize, filterCSVwithJSON::copy, gather_cfg::cout, data_, DigiToRawDM_cff::digiColl, digiInstanceName_, doOrbitOffset_, eventSkipCounter_, FEDRawDataCollection::FEDData(), formatter_, getCalibTriggerType(), MatacqRawEvent::getDccLen(), getMatacqEvent(), getMatacqFile(), getOrbitId(), MatacqRawEvent::getOrbitId(), MatacqRawEvent::getParsedLen(), getRunNumber(), ignoreTriggerType_, inputRawCollectionToken_, MatacqDataFormatter::interpretRawData(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, laserType, matacq_, matacqFedId_, mergeRaw_, eostools::move(), MatacqProducer::stats_t::nLaserEventsWithMatacq, MatacqProducer::stats_t::nNonLaserEventsWithMatacq, now(), orbitOffset_, produceDigis_, produceRaw_, rawInstanceName_, convertSQLiteXML::runNumber, FEDRawData::size(), stats_, and verbosity_.

Referenced by produce().

                                                   {
   edm::Handle<FEDRawDataCollection> sourceColl;
   event.getByToken(inputRawCollectionToken_, sourceColl);
 
   std::unique_ptr<FEDRawDataCollection> rawColl;
   if (produceRaw_) {
     if (mergeRaw_) {
       rawColl = std::make_unique<FEDRawDataCollection>(*sourceColl);
     } else {
       rawColl = std::make_unique<FEDRawDataCollection>();
     }
   }
 
   auto digiColl = std::make_unique<EcalMatacqDigiCollection>();
 
   if (eventSkipCounter_ == 0) {
     if (sourceColl->FEDData(matacqFedId_).size() > 4 && !produceRaw_) {
       //input raw data collection already contains matacqData
       formatter_.interpretRawData(sourceColl->FEDData(matacqFedId_), *digiColl);
     } else {
       bool isLaserEvent = (getCalibTriggerType(event) == laserType);
 
       //      cout << "---> " << (ignoreTriggerType_?"yes":"no") << " " << getCalibTriggerType(event) << endl;
 
       if (isLaserEvent || ignoreTriggerType_) {
         const uint32_t runNumber = getRunNumber(event);
         const uint32_t orbitId = getOrbitId(event);
 
         LogInfo("Matacq") << "Run " << runNumber << "\t Orbit " << orbitId << "\n";
 
         bool fileChange;
         if (doOrbitOffset_) {
           map<uint32_t, uint32_t>::iterator it = orbitOffset_.find(runNumber);
           if (it == orbitOffset_.end()) {
             LogWarning("Matacq") << "Orbit offset not found for run " << runNumber
                                  << ". No orbit correction will be applied.";
           }
         }
 
         if (getMatacqFile(runNumber, orbitId, &fileChange)) {
           //matacq file retrieval succeeded
           LogInfo("Matacq") << "Matacq data file found for "
                             << "run " << runNumber << " orbit " << orbitId;
           if (getMatacqEvent(runNumber, orbitId, fileChange)) {
             if (produceDigis_) {
               formatter_.interpretRawData(matacq_, *digiColl);
             }
             if (produceRaw_) {
               uint32_t dataLen64 = matacq_.getParsedLen();
               if (dataLen64 > bufferSize * 8 || matacq_.getDccLen() != dataLen64) {
                 LogWarning("Matacq") << " Error in Matacq event fragment length! "
                                      << "DCC len: " << matacq_.getDccLen()
                                      << "*8 Bytes, Parsed len: " << matacq_.getParsedLen() << "*8 Bytes.  "
                                      << "Matacq data will not be included for this event.\n";
               } else {
                 rawColl->FEDData(matacqFedId_).resize(dataLen64 * 8);
                 copy(data_.begin(), data_.begin() + dataLen64 * 8, rawColl->FEDData(matacqFedId_).data());
               }
             }
             LogInfo("Matacq") << "Associating matacq data with orbit id " << matacq_.getOrbitId()
                               << " to dcc event with orbit id " << orbitId << std::endl;
             if (isLaserEvent) {
               ++stats_.nLaserEventsWithMatacq;
             } else {
               ++stats_.nNonLaserEventsWithMatacq;
             }
           } else {
             if (isLaserEvent) {
               LogWarning("Matacq") << "No matacq data found for laser event "
                                    << "of run " << runNumber << " orbit " << orbitId;
             }
           }
         } else {
           LogWarning("Matacq") << "No matacq file found for event " << event.id();
         }
       }
     }
     if (eventSkipCounter_ > 0) {  //error occured for this events
       //                       and some events will be skipped following
       //                       to this error.
       LogInfo("Matacq") << " [" << now() << "] " << eventSkipCounter_
                         << " next events will be skipped, following to an "
                         << "error on the last processed event, "
                         << "which is expected to be persistant.";
     }
   } else {
     --eventSkipCounter_;
   }
 
   if (produceRaw_) {
     if (verbosity_ > 1)
       cout << "[Matacq " << now() << "] "
            << "Adding FEDRawDataCollection collection "
            << " to event.\n";
     event.put(std::move(rawColl), rawInstanceName_);
   }
 
   if (produceDigis_) {
     if (verbosity_ > 1)
       cout << "[Matacq " << now() << "] "
            << "Adding EcalMatacqDigiCollection collection "
            << " to event.\n";
     event.put(std::move(digiColl), digiInstanceName_);
   }
 }

◆ getCalibTriggerType()

int MatacqProducer::getCalibTriggerType ( edm::Event & ev ) const

private

Definition at line 1002 of file MatacqProducer.cc.

References data, makeMEIFBenchmarkPlots::ev, Exception, FEDRawDataCollection::FEDData(), EcalPhiSymFlatTableProducers_cfi::id, inputRawCollectionToken_, FEDNumbering::inRange(), edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, and edm_modernize_messagelogger::stat.

Referenced by addMatacqData().

                                                           {
   edm::Handle<FEDRawDataCollection> rawdata;
   ev.getByToken(inputRawCollectionToken_, rawdata);
   if (!(rawdata.isValid())) {
     throw cms::Exception("NotFound") << "No FED raw data collection found. ECAL raw data are "
                                         "required to retrieve the trigger type";
   }
 
   Majority<int> stat;
   for (int id = 601; id <= 654; ++id) {
     if (!FEDNumbering::inRange(id))
       continue;
     const FEDRawData& data = rawdata->FEDData(id);
     const int detailedTrigger32 = 5;
     if (data.size() >= 4 * (detailedTrigger32 + 1)) {
       const unsigned char* pTType = data.data() + detailedTrigger32 * 4;
       int tType = pTType[1] & 0x7;
       stat.add(tType);
     }
   }
   double p;
   int tType = stat.result(&p);
   if (p < 0) {
     //throw cms::Exception("NotFound") << "No ECAL DCC data found\n";
     LogWarning("NotFound") << "No ECAL DCC data found\n";
     tType = -1;
   }
   if (p < .8) {
     //throw cms::Exception("EventCorruption") << "Inconsitency in detailed trigger type indicated in ECAL DCC data headers\n";
     LogWarning("EventCorruption") << "Inconsitency in detailed trigger type indicated in ECAL DCC data headers\n";
     tType = -1;
   }
   return tType;
 }

◆ getMatacqEvent()

bool MatacqProducer::getMatacqEvent	(	uint32_t	runNumber,
		int32_t	orbitId,
		bool	fileChange
	)

private

Definition at line 636 of file MatacqProducer.cc.

References funct::abs(), gather_cfg::cout, data_, MatacqProducer::PosEstimator::eventLength(), Exception, runEdmFileComparison::failed, fastRetrievalThresh_, newFWLiteAna::found, MatacqRawEvent::getDccLen(), MatacqRawEvent::getOrbitId(), MatacqRawEvent::getRunNum(), inFileName_, createfilelist::int, MatacqProducer::PosEstimator::invalid(), lastOrb_, matacq_, mread(), mrewind(), mseek(), msize(), mtell(), now(), HLT_IsoTrack_cff::offset, orbitTolerance_, MatacqProducer::PosEstimator::pos(), posEstim_, writedatasetfile::run, convertSQLiteXML::runNumber, and verbosity_.

Referenced by addMatacqData().

                                                                                         {
   filepos_t startPos;
   if (!mtell(startPos))
     return false;
 
   int32_t startOrb = -1;
   const size_t headerSize = 8 * 8;
   if (mread((char*)&data_[0], headerSize, "Reading matacq header", true)) {
     startOrb = MatacqRawEvent::getOrbitId(&data_[0], headerSize);
     if (startOrb < 0)
       startOrb = 0;
   } else {
     if (verbosity_ > 2) {
       cout << "[Matacq " << now()
            << "] Failed to read matacq header. Moved to start of "
               " the file.\n";
     }
     mrewind();
     if (mread((char*)&data_[0], headerSize, "Reading matacq header", true)) {
       startPos = 0;
       startOrb = MatacqRawEvent::getOrbitId(&data_[0], headerSize);
     } else {
       if (verbosity_ > 2)
         cout << "[Matacq " << now() << "] Looks like matacq file is empty"
              << "\n";
       return false;
     }
   }
 
   if (verbosity_ > 2)
     cout << "[Matacq " << now() << "] Last read orbit: " << lastOrb_ << " looking for orbit " << orbitId
          << ". Current file position: " << startPos << " Orbit at current position: " << startOrb << "\n";
 
   //  f.clear();
   bool didCoarseMove = false;
 
   //FIXME: case where posEtim_.invalid() is false
   if (!posEstim_.invalid() && (abs(lastOrb_ - orbitId) > fastRetrievalThresh_)) {
     filepos_t pos = posEstim_.pos(orbitId);
 
     //    struct stat st;
     filepos_t fsize = -1;
     //    if(0==stat(inFileName_.c_str(), &st)){
     if (msize(fsize)) {
       //      const int64_t fsize = st.st_size;
       if (0 != posEstim_.eventLength() && pos > fsize) {
         //estimated position is beyong end of file
         //-> move to beginning of last event:
         int64_t evtSize = posEstim_.eventLength() * sizeof(uint64_t);
         pos = ((int64_t)fsize / evtSize - 1) * evtSize;
         if (verbosity_ > 2) {
           cout << "[Matacq " << now()
                << "] Estimated position was beyond end of file. "
                   "Changed to "
                << pos << "\n";
         }
       }
     } else {
       LogWarning("Matacq") << "Failed to access file " << inFileName_ << ".";
     }
     if (pos >= 0) {
       if (verbosity_ > 2)
         cout << "[Matacq " << now() << "] jumping to estimated position " << pos << "\n";
       mseek(pos, SEEK_SET, "Jumping to estimated event position");
       if (mread((char*)&data_[0], headerSize, "Reading matacq header", true)) {
         didCoarseMove = true;
       } else {
         //estimated position might have been beyond the end of the file,
         //try, with original position:
         didCoarseMove = false;
         if (!mread((char*)&data_[0], headerSize, "Reading event header", true)) {
           return false;
         }
       }
     } else {
       if (verbosity_)
         cout << "[Matacq " << now()
              << "] Event orbit outside of orbit range "
                 "of matacq data file events\n";
       return false;
     }
   }
 
   int32_t orb = MatacqRawEvent::getOrbitId(&data_[0], headerSize);
 
   if (didCoarseMove) {
     //autoadjustement of threshold for coarse move:
     if (abs(orb - orbitId) > fastRetrievalThresh_) {
       if (verbosity_ > 2)
         cout << "[Matacq " << now() << "] Fast retrieval threshold increased from " << fastRetrievalThresh_;
       fastRetrievalThresh_ = 2 * abs(orb - orbitId);
       if (verbosity_ > 2)
         cout << " to " << fastRetrievalThresh_ << "\n";
     }
 
     //if coarse move did not improve situation, rolls back:
     if (startOrb > 0 && (abs(orb - orbitId) > abs(startOrb - orbitId))) {
       if (verbosity_ > 2)
         cout << "[Matacq " << now() << "] Estimation (-> orbit " << orb
              << ") "
                 "was worst than original position (-> orbit "
              << startOrb << "). Restoring position (" << startPos << ").\n";
       mseek(startPos, SEEK_SET);
       mread((char*)&data_[0], headerSize, "Reading event header", true);
       orb = MatacqRawEvent::getOrbitId(&data_[0], headerSize);
     }
   }
 
   bool searchBackward = (orb > orbitId) ? true : false;
   //BEWARE: len must be signed, because we are using latter in the code (-len)
   //expression
   int len = (int)MatacqRawEvent::getDccLen(&data_[0], headerSize);
 
   if (len == 0) {
     cout << "[Matacq " << now() << "] read DCC length is null! Cancels matacq event search "
          << " and move matacq file pointer to beginning of the file. "
          << "(" << __FILE__ << ":" << __LINE__ << ")."
          << "\n";
     //rewind(f);
     mrewind();
     return false;
   }
 
   enum state_t { searching, found, failed } state = searching;
 
   while (state == searching) {
     orb = MatacqRawEvent::getOrbitId(&data_[0], headerSize);
     len = (int)MatacqRawEvent::getDccLen(&data_[0], headerSize);
     uint32_t run = MatacqRawEvent::getRunNum(&data_[0], headerSize);
     if (verbosity_ > 3) {
       filepos_t pos = -1;
       mtell(pos);
       cout << "[Matacq " << now() << "] Header read at file position " << pos << ":  orbit = " << orb
            << " len = " << len << "x8 Byte"
            << " run = " << run << "\n";
     }
     if ((abs(orb - orbitId) < orbitTolerance_) && (runNumber == 0 || runNumber == run)) {
       state = found;
       lastOrb_ = orb;
       //reads the rest of the event:
       if ((int)data_.size() < len * 8) {
         throw cms::Exception("Matacq") << "Buffer overflow";
       }
       if (verbosity_ > 2)
         cout << "[Matacq " << now()
              << "] Event found. Reading "
                 " matacq event."
              << "\n";
       if (!mread((char*)&data_[0], len * 8, "Reading matacq event")) {
         if (verbosity_ > 2)
           cout << "[Matacq " << now() << "] Failed to read matacq event."
                << "\n";
         state = failed;
       }
       matacq_ = MatacqRawEvent((unsigned char*)&data_[0], len * 8);
     } else {
       if ((searchBackward && (orb < orbitId)) || (!searchBackward && (orb > orbitId))) {  //search ended
         lastOrb_ = orb;
         state = failed;
         if (verbosity_ > 2)
           cout << "[Matacq " << now() << "] No matacq data found for run " << run << ", orbit ID " << orbitId << "."
                << "\n";
       } else {
         off_t offset = (searchBackward ? -len : len) * 8;
         lastOrb_ = orb;
         if (verbosity_ > 3) {
           cout << "[Matacq " << now() << "] In matacq file, moving " << abs(offset) << " byte "
                << (offset > 0 ? "forward" : "backward") << ".\n";
         }
 
         if (mseek(offset, SEEK_CUR, (searchBackward ? "Moving to previous event" : "Moving to next event")) &&
             mread((char*)&data_[0], headerSize, "Reading event header", true)) {
         } else {
           if (!searchBackward)
             mseek(-len * 8, SEEK_CUR, "Moving to start of last complete event");
           state = failed;
         }
       }
     }
   }
 
   if (state == found) {
     filepos_t pos = -1;
     filepos_t fsize = -1;
     mtell(pos);
     msize(fsize);
     if (pos == fsize - 1) {  //last byte.
       if (verbosity_ > 2) {
         cout << "[Matacq " << now()
              << "] Event found was at the end of the file. Moving "
                 "stream position to beginning of this event."
              << "\n";
       }
       mseek(-(int)len * 8 - 1, SEEK_CUR, "Moving to beginning of last matacq event");
     }
   }
   return (state == found);
 }

◆ getMatacqFile()

bool MatacqProducer::getMatacqFile	(	uint32_t	runNumber,
		uint32_t	orbitId,
		bool *	fileChange = `nullptr`
	)

private

Retrieve the file containing a given matacq event

Parameters

runNumber	Number of the run the matacq event is looking from
orbitId	Id of the orbit of the matacq event
fileChange	if not null pointer, set to true if the file changed.

Returns: true if file retrieval succeeded, false otherwise. found.

Definition at line 835 of file MatacqProducer.cc.

References gather_cfg::cout, eventSkipCounter_, fileNames_, alignmentValidation::fname, newFWLiteAna::found, g, getOrbitRange(), RooUtil::FileUtil::glob(), mps_fire::i, MatacqProducer::PosEstimator::init(), lastOrb_, misOpened(), mopen(), now(), onErrorDisablingEvtCnt_, openedFileRunNumber_, posEstim_, convertSQLiteXML::runNumber, runSubDir(), and verbosity_.

Referenced by addMatacqData().

                                                                                          {
   if (openedFileRunNumber_ != 0 && openedFileRunNumber_ == runNumber) {
     uint32_t firstOrb, lastOrb;
     bool goodRange = getOrbitRange(firstOrb, lastOrb);
     //    if(orbitId < firstOrb || orbitId > lastOrb) continue;
     if (goodRange && firstOrb <= orbitId && orbitId <= lastOrb) {
       if (fileChange != nullptr)
         *fileChange = false;
       return misOpened();
     }
   }
 
   if (fileNames_.empty())
     return false;
 
   const string runNumberFormat = "%08d{,_*}";
   string sRunNumber = fmt::sprintf(runNumberFormat, runNumber);
   //cout << "Run number string: " << sRunNumber << "\n";
   bool found = false;
   string fname;
   uint32_t maxOrb = 0;
   //we make two iterations to handle the case where the event is procesed
   //before the matacq data are available. In such case we would have
   //orbitId > maxOrb (maxOrb: orbit of last written matacq event)
   //  for(int itry = 0; itry < 2 && (orbitId > maxOrb); ++itry){
   for (int itry = 0; itry < 1 && (orbitId > maxOrb); ++itry) {
     if (itry > 0) {
       int n_sec = 1;
       std::cout << "[Matacq " << now() << "] Event orbit id (" << orbitId
                 << ") goes "
                    "beyound the range of available one. Waiting for "
                 << n_sec
                 << " seconds in case "
                    "it was not written yet to disk.";
       sleep(n_sec);
     }
 
     for (unsigned i = 0; i < fileNames_.size() && !found; ++i) {
       fname = fileNames_[i];
       boost::algorithm::replace_all(fname, "%run_subdir%", runSubDir(runNumber));
       boost::algorithm::replace_all(fname, "%run_number%", sRunNumber);
 
       glob_t g;
       int rc = glob(fname.c_str(), GLOB_BRACE, nullptr, &g);
       if (rc) {
         if (verbosity_ > 1) {
           switch (rc) {
             case GLOB_NOSPACE:
               std::cout << "[Matacq " << now()
                         << "] Running out of memory while calling glob function to look for matacq file paths\n";
               break;
             case GLOB_ABORTED:
               std::cout << "[Matacq " << now()
                         << "] Read error while calling glob function to look for matacq file paths\n";
               break;
             case GLOB_NOMATCH:
               //ok. No message to report.
               break;
           }
           continue;
         }
       }  //rc
       for (unsigned iglob = 0; iglob < g.gl_pathc; ++iglob) {
         char* thePath = g.gl_pathv[iglob];
         //FIXME: add sanity check on the path
         static std::atomic<int> nOpenErrors{0};
         const int maxOpenErrors = 50;
         if (!mopen(thePath) && nOpenErrors < maxOpenErrors) {
           std::cout << "[Matacq " << now() << "] Failed to open file " << thePath;
           ++nOpenErrors;
           if (nOpenErrors == maxOpenErrors) {
             std::cout << nOpenErrors << "This is the " << maxOpenErrors
                       << "th occurence of this error. Report of this error is now disabled.\n";
           } else {
             std::cout << "\n";
           }
         }
         uint32_t firstOrb;
         uint32_t lastOrb;
         bool goodRange = getOrbitRange(firstOrb, lastOrb);
         std::cout << "Get orbit range " << (goodRange ? "succeeded" : "failed") << ". Range: " << firstOrb << "..."
                   << lastOrb << "\n";
         if (goodRange && lastOrb > maxOrb)
           maxOrb = lastOrb;
         if (goodRange && firstOrb <= orbitId && orbitId <= lastOrb) {
           found = true;
           //continue;
           fname = thePath;
           if (verbosity_ > 1)
             std::cout << "[Matacq " << now() << "] Switching to file " << fname << "\n";
           break;
         }
       }  //next iglob
       globfree(&g);
     }  //next filenames
   }  //next itry
 
   if (found) {
     LogInfo("Matacq") << "Uses matacq data file: '" << fname << "'\n";
   } else {
     if (verbosity_ >= 0)
       cout << "[Matacq " << now()
            << "] no matacq file found "
               "for run "
            << runNumber << ", orbit " << orbitId << "\n";
     eventSkipCounter_ = onErrorDisablingEvtCnt_;
     openedFileRunNumber_ = 0;
     if (fileChange != nullptr)
       *fileChange = false;
     return false;
   }
 
   if (found) {
     openedFileRunNumber_ = runNumber;
     lastOrb_ = 0;
     posEstim_.init(this);
     if (fileChange != nullptr)
       *fileChange = true;
     return true;
   } else {
     return false;
   }
 }

◆ getOrbitId()

uint32_t MatacqProducer::getOrbitId ( edm::Event & ev ) const

private

Definition at line 961 of file MatacqProducer.cc.

References funct::abs(), data, makeMEIFBenchmarkPlots::ev, Exception, FEDRawDataCollection::FEDData(), EcalPhiSymFlatTableProducers_cfi::id, inputRawCollectionToken_, FEDNumbering::inRange(), edm::HandleBase::isValid(), and orbitTolerance_.

Referenced by addMatacqData().

                                                       {
   //on CVS HEAD (June 4, 08), class Event has a method orbitNumber()
   //we could use here. The code would be shorten to:
   //return ev.orbitNumber();
   //we have to deal with what we have in current CMSSW releases:
   edm::Handle<FEDRawDataCollection> rawdata;
   ev.getByToken(inputRawCollectionToken_, rawdata);
   if (!(rawdata.isValid())) {
     throw cms::Exception("NotFound") << "No FED raw data collection found. ECAL raw data are "
                                         "required to retrieve the orbit ID";
   }
 
   int orbit = 0;
   for (int id = 601; id <= 654; ++id) {
     if (!FEDNumbering::inRange(id))
       continue;
     const FEDRawData& data = rawdata->FEDData(id);
     const int orbitIdOffset64 = 3;
     if (data.size() >= 8 * (orbitIdOffset64 + 1)) {  //orbit id is in 4th 64-bit word
       const unsigned char* pOrbit = data.data() + orbitIdOffset64 * 8;
       int thisOrbit = pOrbit[0] | (pOrbit[1] << 8) | (pOrbit[2] << 16) | (pOrbit[3] << 24);
       if (orbit != 0 && thisOrbit != 0 && abs(orbit - thisOrbit) > orbitTolerance_) {
         //throw cms::Exception("EventCorruption")
         //  << "Orbit ID inconsitency in DCC headers";
         LogWarning("EventCorruption") << "Orbit ID inconsitency in DCC headers";
         orbit = 0;
         break;
       }
       if (thisOrbit != 0)
         orbit = thisOrbit;
     }
   }
 
   if (orbit == 0) {
     //    throw cms::Exception("NotFound")
     //  << "Failed to retrieve orbit ID of event "<< ev.id();
     LogWarning("NotFound") << "Failed to retrieve orbit ID of event " << ev.id();
   }
   return orbit;
 }

◆ getOrbitRange()

bool MatacqProducer::getOrbitRange	(	uint32_t &	firstOrb,
		uint32_t &	lastOrb
	)

private

Definition at line 1441 of file MatacqProducer.cc.

References MatacqRawEvent::getDccLen(), MatacqRawEvent::getOrbitId(), RecoTauValidation_cfi::header, createfilelist::int, mread(), mseek(), msize(), mtell(), nEvts, and createJobs::tmp.

Referenced by getMatacqFile().

                                                                         {
   filepos_t pos = -1;
   filepos_t fsize = -1;
   mtell(pos);
   msize(fsize);
   const unsigned headerSize = 8 * 8;
   unsigned char header[headerSize];
   //FIXME: Don't we need here to rewind?
   mseek(0);
   if (!mread((char*)header, headerSize, nullptr, false))
     return false;
   firstOrb = MatacqRawEvent::getOrbitId(header, headerSize);
   int len = (int)MatacqRawEvent::getDccLen(header, headerSize);
   //number of complete events. If last event is partially written,
   //it won't be included in the count.
   unsigned nEvts = fsize / (len * 8);
   //Position of last complete event:
   filepos_t lastEvtPos = (filepos_t)(nEvts - 1) * len * 8;
   //  std::cout << "Move to position : " << lastEvtPos
   //        << "(" << (nEvts - 1) << "*" << len << "*" << 64 << ")"
   //<< "\n";
   mseek(lastEvtPos);
   filepos_t tmp;
   mtell(tmp);
   //std::cout << "New position, sizeof(tmp): " << tmp << "," << sizeof(tmp) << "\n";
   mread((char*)header, headerSize, nullptr, false);
   lastOrb = MatacqRawEvent::getOrbitId(header, headerSize);
 
   //restore file position:
   mseek(pos);
 
   return true;
 }

◆ getRunNumber()

uint32_t MatacqProducer::getRunNumber ( edm::Event & ev ) const

private

Definition at line 959 of file MatacqProducer.cc.

References makeMEIFBenchmarkPlots::ev.

Referenced by addMatacqData(), and produce().

959 { return ev.run(); }

makeMEIFBenchmarkPlots.ev

ev

Definition: makeMEIFBenchmarkPlots.py:55

◆ loadOrbitOffset()

void MatacqProducer::loadOrbitOffset ( )

private

Loading orbit correction table from file.

See also: orbitOffsetFile_

Definition at line 1144 of file MatacqProducer.cc.

References visDQMUpload::buf, gather_cfg::cout, Exception, f, WZElectronSkims53X_cff::max, now(), orbitOffset_, orbitOffsetFile_, writedatasetfile::run, and alignCSCRings::s.

Referenced by MatacqProducer().

                                      {
   std::ifstream f(orbitOffsetFile_.c_str());
   if (f.bad()) {
     throw cms::Exception("Matacq") << "Failed to open orbit ID correction file '" << orbitOffsetFile_ << "'\n";
   }
 
   cout << "[Matacq " << now() << "] "
        << "Offset to substract to Matacq events Orbit ID: \n"
        << "#Run Number\t Offset\n";
 
   string s;
   stringstream buf;
   while (f.eof()) {
     getline(f, s);
     if (s[0] == '#') {  //comment
       //skip line:
       f.ignore(numeric_limits<streamsize>::max(), '\n');
       continue;
     }
     buf.str("");
     buf << s;
     int run;
     int orbit;
     buf >> run;
     buf >> orbit;
     if (buf.bad()) {
       throw cms::Exception("Matacq") << "Syntax error in Orbit offset file '" << orbitOffsetFile_ << "'";
     }
     cout << run << "\t" << orbit << "\n";
     orbitOffset_.insert(pair<int, int>(run, orbit));
   }
 }

◆ mcheck()

bool MatacqProducer::mcheck ( const std::string & name )

private

Definition at line 1378 of file MatacqProducer.cc.

References mergeVDriftHistosByStation::name, and edm_modernize_messagelogger::stat.

                                                  {
   struct stat dummy;
   return 0 == stat(name.c_str(), &dummy);
   //   if(stat(name.c_str(), &dummy)==0){
   //     return true;
   //   } else{
   //     cout << "[Matacq " << now() << "] Failed to stat file '"
   //     << name.c_str() << "'. "
   //     << "Error " << errno << ": " << strerror(errno) << "\n";
   //     return false;
   //   }
 }

◆ mclose()

void MatacqProducer::mclose ( )

private

Definition at line 1404 of file MatacqProducer.cc.

References inFile_.

Referenced by mopen(), and ~MatacqProducer().

                             {
   if (inFile_ != nullptr)
     fclose(inFile_);
   inFile_ = nullptr;
 }

◆ meof()

bool MatacqProducer::meof ( )

private

Definition at line 1412 of file MatacqProducer.cc.

References inFile_.

                           {
   if (nullptr == inFile_)
     return true;
   return feof(inFile_) == 0;
 }

◆ misOpened()

bool MatacqProducer::misOpened ( )

private

Definition at line 1410 of file MatacqProducer.cc.

References inFile_.

Referenced by getMatacqFile().

1410 { return inFile_ != nullptr; }

MatacqProducer::inFile_

FILE_t inFile_

Definition: MatacqProducer.cc:237

◆ mopen()

bool MatacqProducer::mopen ( const std::string & name )

private

Definition at line 1391 of file MatacqProducer.cc.

References inFile_, inFileName_, mclose(), and mergeVDriftHistosByStation::name.

Referenced by getMatacqFile().

                                                 {
   if (inFile_ != nullptr)
     mclose();
   inFile_ = fopen(name.c_str(), "r");
   if (inFile_ != nullptr) {
     inFileName_ = name;
     return true;
   } else {
     inFileName_ = "";
     return false;
   }
 }

◆ mread()

bool MatacqProducer::mread	(	char *	buf,
		size_t	n,
		const char *	mess = `nullptr`,
		bool	peek = `false`
	)

private

Read a data block from input file. On failure file position is restored and if position restoring fails, file is rewind.

Parameters

buf	buffer to store read data
n	size of data block
mess	text to insert in the eventual error message.
peek	if true file position is restored after the data read

Returns: true on success, false on failure

Definition at line 1327 of file MatacqProducer.cc.

References visDQMUpload::buf, gather_cfg::cout, inFile_, lastOrb_, mrewind(), create_idmaps::n, now(), and verbosity_.

Referenced by getMatacqEvent(), getOrbitRange(), and MatacqProducer::PosEstimator::init().

                                                                            {
   if (nullptr == inFile_)
     return false;
   off_t pos = ftello(inFile_);
   bool rc = (pos != -1) && (1 == fread(buf, n, 1, inFile_));
   if (!rc) {
     if (verbosity_) {
       cout << "[Matacq " << now() << "] ";
       if (mess)
         cout << mess << ". ";
       cout << "Read failure from input matacq file.\n";
     }
     clearerr(inFile_);
   }
   if (peek || !rc) {  //need to restore file position
     if (0 != fseeko(inFile_, pos, SEEK_SET)) {
       if (verbosity_) {
         cout << "[Matacq " << now() << "] ";
         if (mess)
           cout << mess << ". ";
         cout << "Failed to restore file position of "
                 "before read error. Rewind file.\n";
       }
       //rewind(inFile_.get());
       mrewind();
       lastOrb_ = 0;
     }
   }
   return rc;
 }

◆ mrewind()

bool MatacqProducer::mrewind ( )

private

Definition at line 1371 of file MatacqProducer.cc.

References inFile_.

Referenced by getMatacqEvent(), MatacqProducer::PosEstimator::init(), mread(), and mseek().

                              {
   if (nullptr == inFile_)
     return false;
   clearerr(inFile_);
   return fseeko(inFile_, 0, SEEK_SET) != 0;
 }

◆ mseek()

bool MatacqProducer::mseek	(	filepos_t	offset,
		int	whence = `SEEK_SET`,
		const char *	mess = `nullptr`
	)

private

Move input file read pointer. On failure file is rewind.

Parameters

buf	buffer to store read data
n	size of data block
mess	text to insert in the eventual error message.

Returns: true on success, false on failure

Definition at line 1305 of file MatacqProducer.cc.

References gather_cfg::cout, inFile_, mrewind(), now(), HLT_IsoTrack_cff::offset, and verbosity_.

Referenced by getMatacqEvent(), getOrbitRange(), and MatacqProducer::PosEstimator::init().

                                                                      {
   if (nullptr == inFile_)
     return false;
   const int rc = fseeko(inFile_, offset, whence);
   if (rc != 0 && verbosity_) {
     cout << "[Matacq " << now() << "] ";
     if (mess)
       cout << mess << ". ";
     cout << "Random access error on input matacq file. "
             "Rewind file.\n";
     mrewind();
   }
   return rc == 0;
 }

◆ msize()

bool MatacqProducer::msize ( filepos_t & s )

private

Definition at line 1358 of file MatacqProducer.cc.

References visDQMUpload::buf, inFile_, and alignCSCRings::s.

Referenced by getMatacqEvent(), getOrbitRange(), and MatacqProducer::PosEstimator::init().

                                        {
   if (nullptr == inFile_)
     return false;
   struct stat buf;
   if (0 != fstat(fileno(inFile_), &buf)) {
     s = 0;
     return false;
   } else {
     s = buf.st_size;
     return true;
   }
 }

◆ mtell()

bool MatacqProducer::mtell ( filepos_t & pos )

private

Definition at line 1320 of file MatacqProducer.cc.

References inFile_.

Referenced by getMatacqEvent(), and getOrbitRange().

                                          {
   if (nullptr == inFile_)
     return false;
   pos = ftello(inFile_);
   return pos != -1;
 }

◆ newRun()

void MatacqProducer::newRun	(	int	prevRun,
		int	newRun
	)

private

Definition at line 1427 of file MatacqProducer.cc.

References eventSkipCounter_, logFile_, MatacqProducer::stats_t::nEvents, MatacqProducer::stats_t::nLaserEventsWithMatacq, MatacqProducer::stats_t::nNonLaserEventsWithMatacq, now(), runNumber_, and stats_.

Referenced by produce().

                                                    {
   runNumber_ = newRun;
   eventSkipCounter_ = 0;
   logFile_ << "[" << now() << "] Event count for run " << runNumber_ << ": "
            << "total: " << stats_.nEvents << ", "
            << "Laser event with Matacq data: " << stats_.nLaserEventsWithMatacq << ", "
            << "Non laser event (according to DCC header) with Matacq data: " << stats_.nNonLaserEventsWithMatacq << "\n"
            << flush;
 
   stats_.nEvents = 0;
   stats_.nLaserEventsWithMatacq = 0;
   stats_.nNonLaserEventsWithMatacq = 0;
 }

◆ produce()

void MatacqProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Produces the EDM products

Parameters

CMS	event
eventSetup	event conditions

Implements edm::one::EDProducerBase.

Definition at line 453 of file MatacqProducer.cc.

References addMatacqData(), gather_cfg::cout, disabled_, getRunNumber(), logTiming_, MatacqProducer::stats_t::nEvents, newRun(), now(), convertSQLiteXML::runNumber, runNumber_, startTime_, stats_, submitPVValidationJobs::t, timeLog_, timer_, and verbosity_.

                                                                              {
   if (verbosity_ >= 4)
     cout << "[Matacq " << now() << "] in MatacqProducer::produce" << endl;
   if (logTiming_) {
     timeval t;
     gettimeofday(&t, nullptr);
 
     timeLog_ << t.tv_sec << "." << setfill('0') << setw(3) << (t.tv_usec + 500) / 1000 << setfill(' ') << "\t"
              << (t.tv_usec - timer_.tv_usec) * 1. + (t.tv_sec - timer_.tv_sec) * 1.e6 << "\t";
     timer_ = t;
   }
 
   if (startTime_.tv_sec == 0)
     gettimeofday(&startTime_, nullptr);
   ++stats_.nEvents;
   if (disabled_)
     return;
   const uint32_t runNumber = getRunNumber(event);
   if (runNumber != runNumber_) {
     newRun(runNumber_, runNumber);
   }
   addMatacqData(event);
 
   if (logTiming_) {
     timeval t;
     gettimeofday(&t, nullptr);
     timeLog_ << (t.tv_usec - timer_.tv_usec) * 1. + (t.tv_sec - timer_.tv_sec) * 1.e6 << "\n";
     timer_ = t;
   }
 }

◆ runSubDir()

std::string MatacqProducer::runSubDir ( uint32_t runNumber )

staticprivate

Definition at line 1420 of file MatacqProducer.cc.

References convertSQLiteXML::runNumber, and units().

Referenced by getMatacqFile().

                                                       {
   int millions = runNumber / (1000 * 1000);
   int thousands = (runNumber - millions * 1000 * 1000) / 1000;
   int units = runNumber - millions * 1000 * 1000 - thousands * 1000;
   return fmt::sprintf("%03d/%03d/%03d", millions, thousands, units);
 }

Member Data Documentation

◆ bufferSize

const int MatacqProducer::bufferSize = 30000

staticprivate

Definition at line 239 of file MatacqProducer.cc.

Referenced by addMatacqData().

◆ data_

std::vector<unsigned char> MatacqProducer::data_

private

Definition at line 241 of file MatacqProducer.cc.

Referenced by addMatacqData(), and getMatacqEvent().

◆ digiInstanceName_

std::string MatacqProducer::digiInstanceName_

private

Instance name to use for the produced Matacq digi collection

Definition at line 188 of file MatacqProducer.cc.

Referenced by addMatacqData(), and MatacqProducer().

◆ disabled_

bool MatacqProducer::disabled_

private

Parameter to disable matacq data production. For timing purpose.

Definition at line 200 of file MatacqProducer.cc.

Referenced by produce().

◆ doOrbitOffset_

bool MatacqProducer::doOrbitOffset_

private

Switch for orbit ID correction.

See also: orbitOffsetFile_

Definition at line 264 of file MatacqProducer.cc.

Referenced by addMatacqData(), and MatacqProducer().

◆ eventSkipCounter_

int MatacqProducer::eventSkipCounter_

private

counter for event skipping

Definition at line 289 of file MatacqProducer.cc.

Referenced by addMatacqData(), getMatacqFile(), and newRun().

◆ fastRetrievalThresh_

int MatacqProducer::fastRetrievalThresh_

private

Definition at line 246 of file MatacqProducer.cc.

Referenced by getMatacqEvent().

◆ fileNames_

std::vector<std::string> MatacqProducer::fileNames_

private

Definition at line 184 of file MatacqProducer.cc.

Referenced by getMatacqFile().

◆ formatter_

MatacqDataFormatter MatacqProducer::formatter_

private

Definition at line 242 of file MatacqProducer.cc.

Referenced by addMatacqData().

◆ ignoreTriggerType_

bool MatacqProducer::ignoreTriggerType_

private

When true look for matacq data independently of trigger type.

Definition at line 231 of file MatacqProducer.cc.

Referenced by addMatacqData().

◆ inFile_

FILE_t MatacqProducer::inFile_

private

Stream of currently opened matacq file

Definition at line 237 of file MatacqProducer.cc.

Referenced by mclose(), meof(), misOpened(), mopen(), mread(), mrewind(), mseek(), msize(), and mtell().

◆ inFileName_

std::string MatacqProducer::inFileName_

private

Name of currently opened matacq file

Definition at line 268 of file MatacqProducer.cc.

Referenced by getMatacqEvent(), and mopen().

◆ inputRawCollection_

edm::InputTag MatacqProducer::inputRawCollection_

private

Name of the raw data collection the Matacq data must be merge to if merging is enabled.

Definition at line 217 of file MatacqProducer.cc.

◆ inputRawCollectionToken_

edm::EDGetTokenT<FEDRawDataCollection> MatacqProducer::inputRawCollectionToken_

private

EDM token to access the raw data collection the Matacq data must be merge to if merging is enabled.

Definition at line 222 of file MatacqProducer.cc.

Referenced by addMatacqData(), getCalibTriggerType(), getOrbitId(), and MatacqProducer().

◆ lastOrb_

int32_t MatacqProducer::lastOrb_

private

Definition at line 245 of file MatacqProducer.cc.

Referenced by getMatacqEvent(), getMatacqFile(), and mread().

◆ logFile_

std::ofstream MatacqProducer::logFile_

private

Log file

Definition at line 285 of file MatacqProducer.cc.

Referenced by MatacqProducer(), and newRun().

◆ logFileName_

std::string MatacqProducer::logFileName_

private

Log file name

Definition at line 281 of file MatacqProducer.cc.

Referenced by MatacqProducer().

◆ logTiming_

bool MatacqProducer::logTiming_

private

Switch for code timing.

Definition at line 308 of file MatacqProducer.cc.

Referenced by MatacqProducer(), produce(), and ~MatacqProducer().

◆ matacq_

MatacqRawEvent MatacqProducer::matacq_

private

Definition at line 233 of file MatacqProducer.cc.

Referenced by addMatacqData(), and getMatacqEvent().

◆ matacqFedId_

const int MatacqProducer::matacqFedId_ = 655

staticprivate

Definition at line 270 of file MatacqProducer.cc.

Referenced by addMatacqData().

◆ mergeRaw_

bool MatacqProducer::mergeRaw_

private

Switch for merging Matacq raw data with existing raw data collection.

Definition at line 227 of file MatacqProducer.cc.

Referenced by addMatacqData().

◆ onErrorDisablingEvtCnt_

int MatacqProducer::onErrorDisablingEvtCnt_

private

Number of events to skip in case of error

Definition at line 293 of file MatacqProducer.cc.

Referenced by getMatacqFile().

◆ openedFileRunNumber_

uint32_t MatacqProducer::openedFileRunNumber_

private

Definition at line 244 of file MatacqProducer.cc.

Referenced by getMatacqFile().

◆ orbitOffset_

std::map<uint32_t, uint32_t> MatacqProducer::orbitOffset_

private

Orbit offset table.

See also: orbitOffsetFile_

Definition at line 260 of file MatacqProducer.cc.

Referenced by addMatacqData(), and loadOrbitOffset().

◆ orbitOffsetFile_

std::string MatacqProducer::orbitOffsetFile_

private

File name of table with orbit offset between matacq event and DCC. Used to recover data suffering from orbit miss-synchonization

Definition at line 256 of file MatacqProducer.cc.

Referenced by loadOrbitOffset(), and MatacqProducer().

◆ orbitTolerance_

const int MatacqProducer::orbitTolerance_ = 80

staticprivate

Definition at line 243 of file MatacqProducer.cc.

Referenced by getMatacqEvent(), and getOrbitId().

◆ posEstim_

PosEstimator MatacqProducer::posEstim_

private

Definition at line 248 of file MatacqProducer.cc.

Referenced by getMatacqEvent(), getMatacqFile(), and MatacqProducer().

◆ produceDigis_

bool MatacqProducer::produceDigis_

private

Swictch for Matacq digi producion

Definition at line 208 of file MatacqProducer.cc.

Referenced by addMatacqData(), and MatacqProducer().

◆ produceRaw_

bool MatacqProducer::produceRaw_

private

Switch for Matacq FED raw data production

Definition at line 212 of file MatacqProducer.cc.

Referenced by addMatacqData(), and MatacqProducer().

◆ rawInstanceName_

std::string MatacqProducer::rawInstanceName_

private

Instance name to use for the produced Matacq raw data collection

Definition at line 192 of file MatacqProducer.cc.

Referenced by addMatacqData(), and MatacqProducer().

◆ runNumber_

uint32_t MatacqProducer::runNumber_

private

Number of the currently processed run

Definition at line 312 of file MatacqProducer.cc.

Referenced by newRun(), and produce().

◆ startTime_

timeval MatacqProducer::startTime_

private

Definition at line 250 of file MatacqProducer.cc.

Referenced by MatacqProducer(), produce(), and ~MatacqProducer().

◆ stats_

struct MatacqProducer::stats_t MatacqProducer::stats_

private

Referenced by addMatacqData(), newRun(), and produce().

◆ stats_init

const MatacqProducer::stats_t MatacqProducer::stats_init = {0, 0, 0}

staticprivate

Definition at line 278 of file MatacqProducer.cc.

◆ timeLog_

std::ofstream MatacqProducer::timeLog_

private

Output stream to log code timing

Definition at line 304 of file MatacqProducer.cc.

Referenced by MatacqProducer(), and produce().

◆ timeLogFile_

std::string MatacqProducer::timeLogFile_

private

Name of file to log timing

Definition at line 297 of file MatacqProducer.cc.

Referenced by MatacqProducer().

◆ timer_

timeval MatacqProducer::timer_

private

Buffer for timing

Definition at line 300 of file MatacqProducer.cc.

Referenced by MatacqProducer(), and produce().

◆ timing_

bool MatacqProducer::timing_

private

Parameter to switch module timing.

Definition at line 196 of file MatacqProducer.cc.

◆ verbosity_

int MatacqProducer::verbosity_

private

Verbosity level

Definition at line 204 of file MatacqProducer.cc.

Referenced by addMatacqData(), getMatacqEvent(), getMatacqFile(), MatacqProducer(), mread(), mseek(), MatacqProducer::PosEstimator::pos(), and produce().

Classes

Public Types

Public Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ FILE_t

◆ filepos_t

Member Enumeration Documentation

◆ calibTrigType_t

Constructor & Destructor Documentation

◆ MatacqProducer()

◆ ~MatacqProducer()

Member Function Documentation

◆ addMatacqData()

◆ getCalibTriggerType()

◆ getMatacqEvent()

◆ getMatacqFile()

◆ getOrbitId()

◆ getOrbitRange()

◆ getRunNumber()

◆ loadOrbitOffset()

◆ mcheck()

◆ mclose()

◆ meof()

◆ misOpened()

◆ mopen()

◆ mread()

◆ mrewind()

◆ mseek()

◆ msize()

◆ mtell()

◆ newRun()

◆ produce()

◆ runSubDir()

Member Data Documentation

◆ bufferSize

◆ data_

◆ digiInstanceName_

◆ disabled_

◆ doOrbitOffset_

◆ eventSkipCounter_

◆ fastRetrievalThresh_

◆ fileNames_

◆ formatter_

◆ ignoreTriggerType_

◆ inFile_

◆ inFileName_

◆ inputRawCollection_

◆ inputRawCollectionToken_

◆ lastOrb_

◆ logFile_

◆ logFileName_

◆ logTiming_

◆ matacq_

◆ matacqFedId_

◆ mergeRaw_

◆ onErrorDisablingEvtCnt_

◆ openedFileRunNumber_

◆ orbitOffset_

◆ orbitOffsetFile_

◆ orbitTolerance_

◆ posEstim_

◆ produceDigis_

◆ produceRaw_

◆ rawInstanceName_

◆ runNumber_

◆ startTime_

◆ stats_

◆ stats_init

◆ timeLog_

◆ timeLogFile_

◆ timer_

◆ timing_

◆ verbosity_