#include <SiStripSpyEventMatcher.h>

Classes
class	CountersWrapper

class	EventKey

struct	MatchingOutput

class	SpyDataCollections

Public Types
typedef edm::EventID	EventID

typedef std::set< EventID >	SpyEventList

Public Member Functions
void	getCollections (const edm::EventPrincipal &event, const uint32_t eventId, const uint8_t apvAddress, const SiStripFedCabling &cabling, MatchingOutput &matchingOutput)

void	getMatchedCollections (const uint32_t eventId, const uint8_t apvAddress, const SpyEventList *matchingEvents, const SiStripFedCabling &cabling, SpyDataCollections &collectionsToCreate)

void	initialize ()

const SpyEventList *	matchesForEvent (const uint32_t eventId, const uint8_t apvAddress) const

	SpyEventMatcher (const edm::ParameterSet &config)

virtual	~SpyEventMatcher ()

Private Types
typedef std::vector< uint32_t >	Counters

typedef std::shared_ptr < CountersWrapper >	CountersPtr

typedef edm::VectorInputSource	Source

Private Member Functions
void	addNextEventToMap (const edm::EventPrincipal &nextSpyEvent)

std::unique_ptr< Source >	constructSource (const edm::ParameterSet &sourceConfig)

void	operator() (const edm::EventPrincipal &event)

Static Private Member Functions
static void	findMatchingFeds (const uint32_t eventId, const uint8_t apvAddress, CountersPtr totalEventCounters, CountersPtr l1aCounters, CountersPtr apvAddresses, std::set< uint16_t > &matchingFeds)

static CountersPtr	getCounters (const edm::EventPrincipal &event, const edm::InputTag &tag, const bool mapKeyIsByFedID=true)

template<class T >
static const T *	getProduct (const edm::EventPrincipal &event, const edm::InputTag &tag)

static void	mergeMatchingData (const std::set< uint16_t > &matchingFeds, const FEDRawDataCollection &inputRawData, CountersPtr inputTotalEventCounters, CountersPtr inputL1ACounters, CountersPtr inputAPVAddresses, const edm::DetSetVector< SiStripRawDigi > inputScopeDigis, const edm::DetSetVector< SiStripRawDigi > inputPayloadDigis, const edm::DetSetVector< SiStripRawDigi > inputReorderedDigis, const edm::DetSetVector< SiStripRawDigi > inputVirginRawDigis, FEDRawDataCollection &outputRawData, std::vector< uint32_t > &outputTotalEventCounters, std::vector< uint32_t > &outputL1ACounters, std::vector< uint32_t > &outputAPVAddresses, std::vector< edm::DetSet< SiStripRawDigi > > outputScopeDigisVector, std::vector< edm::DetSet< SiStripRawDigi > > outputPayloadDigisVector, std::vector< edm::DetSet< SiStripRawDigi > > outputReorderedDigisVector, std::vector< edm::DetSet< SiStripRawDigi > > outputVirginRawDigisVector, const SiStripFedCabling &cabling)

Private Attributes
edm::InputTag	apvAddressesTag_

uint32_t	counterDiffMax_

std::map< EventKey, SpyEventList >	eventMatches_

std::unique_ptr < edm::EventPrincipal >	eventPrincipal_

edm::InputTag	l1aCountersTag_

edm::InputTag	payloadDigisTag_

std::unique_ptr < edm::ProcessConfiguration >	processConfiguration_

std::shared_ptr < edm::ProductRegistry >	productRegistry_

edm::InputTag	rawDataTag_

edm::InputTag	reorderedDigisTag_

edm::InputTag	scopeDigisTag_

std::unique_ptr < edm::VectorInputSource > const	source_

edm::InputTag	totalEventCountersTag_

edm::InputTag	virginRawDigisTag_

Static Private Attributes
static const char *const	mlLabel_ = "SpyEventMatcher"

Detailed Description

Definition at line 33 of file SiStripSpyEventMatcher.h.

Member Typedef Documentation

typedef std::vector<uint32_t> sistrip::SpyEventMatcher::Counters

private

Definition at line 103 of file SiStripSpyEventMatcher.h.

typedef std::shared_ptr<CountersWrapper> sistrip::SpyEventMatcher::CountersPtr

private

Definition at line 121 of file SiStripSpyEventMatcher.h.

typedef edm::EventID sistrip::SpyEventMatcher::EventID

Definition at line 69 of file SiStripSpyEventMatcher.h.

typedef edm::VectorInputSource sistrip::SpyEventMatcher::Source

private

Definition at line 124 of file SiStripSpyEventMatcher.h.

typedef std::set<EventID> sistrip::SpyEventMatcher::SpyEventList

Definition at line 70 of file SiStripSpyEventMatcher.h.

Constructor & Destructor Documentation

sistrip::SpyEventMatcher::SpyEventMatcher ( const edm::ParameterSet & config )

Definition at line 47 of file SiStripSpyEventMatcher.cc.

References edm::ParameterSet::emptyParameterSetID(), eventPrincipal_, processConfiguration_, productRegistry_, and source_.

       : rawDataTag_(config.getParameter<edm::InputTag>("RawSpyDataTag")),
         totalEventCountersTag_(config.getParameter<edm::InputTag>("SpyTotalEventCountersTag")),
         l1aCountersTag_(config.getParameter<edm::InputTag>("SpyL1ACountersTag")),
         apvAddressesTag_(config.getParameter<edm::InputTag>("SpyAPVAddressesTag")),
         scopeDigisTag_(config.getParameter<edm::InputTag>("SpyScopeDigisTag")),
         payloadDigisTag_(config.getParameter<edm::InputTag>("SpyPayloadDigisTag")),
         reorderedDigisTag_(config.getParameter<edm::InputTag>("SpyReorderedDigisTag")),
         virginRawDigisTag_(config.getParameter<edm::InputTag>("SpyVirginRawDigisTag")),
         counterDiffMax_(config.getParameter<uint32_t>("CounterDiffMaxAllowed")),
         productRegistry_(new edm::SignallingProductRegistry),
         source_(constructSource(config.getParameter<edm::ParameterSet>("SpySource"))),
         processConfiguration_(
             new edm::ProcessConfiguration(std::string("@MIXING"), edm::getReleaseVersion(), edm::getPassID())),
         eventPrincipal_() {
     // Use the empty parameter set for the parameter set ID of our "@MIXING" process.
     processConfiguration_->setParameterSetID(edm::ParameterSet::emptyParameterSetID());
     productRegistry_->setFrozen();
 
     eventPrincipal_ = std::make_unique<edm::EventPrincipal>(source_->productRegistry(),
                                                             std::make_shared<edm::BranchIDListHelper>(),
                                                             std::make_shared<edm::ThinnedAssociationsHelper>(),
                                                             *processConfiguration_,
                                                             nullptr);
   }

sistrip::SpyEventMatcher::~SpyEventMatcher ( )

virtual

Definition at line 45 of file SiStripSpyEventMatcher.cc.

45 {}

Member Function Documentation

void sistrip::SpyEventMatcher::addNextEventToMap ( const edm::EventPrincipal & nextSpyEvent )

private

Definition at line 104 of file SiStripSpyEventMatcher.cc.

References apvAddressesTag_, counterDiffMax_, edm::EventID::event(), eventMatches_, getCounters(), edm::EventPrincipal::id(), submitPVResolutionJobs::key, l1aCountersTag_, LogDebug, mlLabel_, contentValuesCheck::ss, and totalEventCountersTag_.

Referenced by initialize().

                                                                                {
     edm::EventID spyEventId = nextSpyEvent.id();
 
     CountersPtr totalEventCounters = getCounters(nextSpyEvent, totalEventCountersTag_);
     CountersPtr l1aCounters = getCounters(nextSpyEvent, l1aCountersTag_);
     CountersPtr apvAddresses = getCounters(nextSpyEvent, apvAddressesTag_, false);
     //loop over all FEDs. Maps should have same content and be in order so, avoid searches by iterating (and checking keys match)
     //add all possible event keys to the map
     std::vector<uint32_t>::const_iterator iTotalEventCount = totalEventCounters->begin();
     std::vector<uint32_t>::const_iterator iL1ACount = l1aCounters->begin();
     std::vector<uint32_t>::const_iterator iAPVAddress = apvAddresses->begin();
     //for debug
     std::map<EventKey, uint16_t> fedCounts;
     unsigned int fedid = 0;
     for (; ((iTotalEventCount != totalEventCounters->end()) && (iL1ACount != l1aCounters->end()) &&
             (iAPVAddress != apvAddresses->end()));
          (++iTotalEventCount, ++iL1ACount, ++iAPVAddress, ++fedid)) {
       if (*iAPVAddress == 0) {
         continue;
       }
 
       if (((*iTotalEventCount) > (*iL1ACount)) || ((*iL1ACount) - (*iTotalEventCount) > counterDiffMax_)) {
         LogWarning(mlLabel_) << "Spy event " << spyEventId.event() << " error in counter values for FED " << fedid
                              << ", totCount = " << *iTotalEventCount << ", L1Acount = " << *iL1ACount << std::endl;
 
         continue;
       }
 
       for (uint32_t eventId = (*iTotalEventCount) + 1; eventId <= (*iL1ACount) + 1; ++eventId) {
         EventKey key(eventId, *iAPVAddress);
         eventMatches_[key].insert(spyEventId);
         fedCounts[key]++;
       }
     }
 
     //for debug
     std::ostringstream ss;
     ss << "Spy event " << spyEventId.event() << " matches (eventID,apvAddress,nFEDs): ";
     for (std::map<EventKey, uint16_t>::const_iterator iEventFEDCount = fedCounts.begin();
          iEventFEDCount != fedCounts.end();
          ++iEventFEDCount) {
       ss << "(" << iEventFEDCount->first.eventId() << "," << uint16_t(iEventFEDCount->first.apvAddress()) << ","
          << iEventFEDCount->second << ") ";
     }
     LogDebug(mlLabel_) << ss.str();
   }

std::unique_ptr< SpyEventMatcher::Source > sistrip::SpyEventMatcher::constructSource ( const edm::ParameterSet & sourceConfig )

private

Definition at line 73 of file SiStripSpyEventMatcher.cc.

References getPayloadData::description, edm::VectorInputSourceFactory::get(), edm::VectorInputSourceFactory::makeVectorInputSource(), and productRegistry_.

                                                                                                            {
     const edm::VectorInputSourceFactory* sourceFactory = edm::VectorInputSourceFactory::get();
     edm::VectorInputSourceDescription description(productRegistry_, edm::PreallocationConfiguration());
     return sourceFactory->makeVectorInputSource(sourceConfig, description);
   }

void sistrip::SpyEventMatcher::findMatchingFeds	(	const uint32_t	eventId,
		const uint8_t	apvAddress,
		SpyEventMatcher::CountersPtr	totalEventCounters,
		SpyEventMatcher::CountersPtr	l1aCounters,
		SpyEventMatcher::CountersPtr	apvAddresses,
		std::set< uint16_t > &	matchingFeds
	)

staticprivate

Definition at line 277 of file SiStripSpyEventMatcher.cc.

Referenced by getCollections().

                                                                          {
     //loop over all FEDs. Maps should have same content and be in order so, avoid searches by iterating (and checking keys match)
     std::vector<uint32_t>::const_iterator iTotalEventCount = totalEventCounters->begin();
     std::vector<uint32_t>::const_iterator iL1ACount = l1aCounters->begin();
     std::vector<uint32_t>::const_iterator iAPVAddress = apvAddresses->begin();
     for (; ((iTotalEventCount != totalEventCounters->end()) && (iL1ACount != l1aCounters->end()) &&
             (iAPVAddress != apvAddresses->end()));
          (++iTotalEventCount, ++iL1ACount, ++iAPVAddress)) {
       if (*iAPVAddress == 0) {
         continue;
       }
       if ((eventId > *iTotalEventCount) && (eventId <= (*iL1ACount) + 1) && (*iAPVAddress == apvAddress)) {
         matchingFeds.insert(matchingFeds.end(), iTotalEventCount - totalEventCounters->begin());
       }
     }
   }

void sistrip::SpyEventMatcher::getCollections	(	const edm::EventPrincipal &	event,
		const uint32_t	eventId,
		const uint8_t	apvAddress,
		const SiStripFedCabling &	cabling,
		MatchingOutput &	matchingOutput
	)

Definition at line 171 of file SiStripSpyEventMatcher.cc.

Referenced by getMatchedCollections().

                                                            {
     //read the input collections from the event
     const FEDRawDataCollection* inputRawDataPtr = getProduct<FEDRawDataCollection>(event, rawDataTag_);
     if (!inputRawDataPtr) {
       throw cms::Exception(mlLabel_) << "Failed to get raw spy data with tag " << rawDataTag_ << " from spy event";
     }
     const FEDRawDataCollection& inputRawData = *inputRawDataPtr;
     CountersPtr inputTotalEventCounters = getCounters(event, totalEventCountersTag_);
     CountersPtr inputL1ACounters = getCounters(event, l1aCountersTag_);
     CountersPtr inputAPVAddresses = getCounters(event, apvAddressesTag_, false);
     const edm::DetSetVector<SiStripRawDigi>* inputScopeDigis =
         getProduct<edm::DetSetVector<SiStripRawDigi> >(event, scopeDigisTag_);
     const edm::DetSetVector<SiStripRawDigi>* inputPayloadDigis =
         getProduct<edm::DetSetVector<SiStripRawDigi> >(event, payloadDigisTag_);
     const edm::DetSetVector<SiStripRawDigi>* inputReorderedDigis =
         getProduct<edm::DetSetVector<SiStripRawDigi> >(event, reorderedDigisTag_);
     const edm::DetSetVector<SiStripRawDigi>* inputVirginRawDigis =
         getProduct<edm::DetSetVector<SiStripRawDigi> >(event, virginRawDigisTag_);
     //construct the output vectors if the digis were found and they do not exist
     if (inputScopeDigis && !mo.outputScopeDigisVector_.get())
       mo.outputScopeDigisVector_.reset(new std::vector<edm::DetSet<SiStripRawDigi> >);
     if (inputPayloadDigis && !mo.outputPayloadDigisVector_.get())
       mo.outputPayloadDigisVector_.reset(new std::vector<edm::DetSet<SiStripRawDigi> >);
     if (inputReorderedDigis && !mo.outputReorderedDigisVector_.get())
       mo.outputReorderedDigisVector_.reset(new std::vector<edm::DetSet<SiStripRawDigi> >);
     if (inputVirginRawDigis && !mo.outputVirginRawDigisVector_.get())
       mo.outputVirginRawDigisVector_.reset(new std::vector<edm::DetSet<SiStripRawDigi> >);
     //find matching FEDs
     std::set<uint16_t> matchingFeds;
     findMatchingFeds(eventId, apvAddress, inputTotalEventCounters, inputL1ACounters, inputAPVAddresses, matchingFeds);
     LogInfo(mlLabel_) << "Spy event " << event.id() << " has " << matchingFeds.size() << " matching FEDs";
     std::ostringstream ss;
     ss << "Matching FEDs for event " << event.id() << ": ";
     for (std::set<uint16_t>::const_iterator iFedId = matchingFeds.begin(); iFedId != matchingFeds.end(); ++iFedId) {
       ss << *iFedId << " ";
     }
     LogDebug(mlLabel_) << ss.str();
     //check there are no duplicates
     std::vector<uint16_t> duplicateFeds(std::min(mo.alreadyMergedFeds_.size(), matchingFeds.size()));
     std::vector<uint16_t>::iterator duplicatesBegin = duplicateFeds.begin();
     std::vector<uint16_t>::iterator duplicatesEnd = std::set_intersection(mo.alreadyMergedFeds_.begin(),
                                                                           mo.alreadyMergedFeds_.end(),
                                                                           matchingFeds.begin(),
                                                                           matchingFeds.end(),
                                                                           duplicatesBegin);
     if ((duplicatesEnd - duplicatesBegin) != 0) {
       std::ostringstream ss;
       ss << "Found a match for FEDs ";
       for (std::vector<uint16_t>::const_iterator iDup = duplicatesBegin; iDup != duplicatesEnd; ++iDup) {
         ss << *iDup << " ";
       }
       ss << ". Output SetSetVectors will be unusable!";
       LogError(mlLabel_) << ss.str();
     }
     //merge the matching data
     mergeMatchingData(matchingFeds,
                       inputRawData,
                       inputTotalEventCounters,
                       inputL1ACounters,
                       inputAPVAddresses,
                       inputScopeDigis,
                       inputPayloadDigis,
                       inputReorderedDigis,
                       inputVirginRawDigis,
                       mo.outputRawData_,
                       mo.outputTotalEventCounters_,
                       mo.outputL1ACounters_,
                       mo.outputAPVAddresses_,
                       mo.outputScopeDigisVector_.get(),
                       mo.outputPayloadDigisVector_.get(),
                       mo.outputReorderedDigisVector_.get(),
                       mo.outputVirginRawDigisVector_.get(),
                       cabling);
     mo.alreadyMergedFeds_.insert(matchingFeds.begin(), matchingFeds.end());
   }

SpyEventMatcher::CountersPtr sistrip::SpyEventMatcher::getCounters	(	const edm::EventPrincipal &	event,
		const edm::InputTag &	tag,
		const bool	mapKeyIsByFedID = `true`
	)

staticprivate

Definition at line 391 of file SiStripSpyEventMatcher.cc.

References edmPickEvents::event, Exception, sistrip::FED_ID_MAX, sistrip::SpyUtilities::fillFEDMajorities(), mlLabel_, and GlobalPosition_Frontier_DevDB_cff::tag.

Referenced by addNextEventToMap(), and getCollections().

                                                                                         {
     const std::vector<uint32_t>* vectorFromEvent = getProduct<std::vector<uint32_t> >(event, tag);
     if (vectorFromEvent) {
       //vector is from event so, will be deleted when the event is destroyed (and not before)
       return std::make_shared<CountersWrapper>(vectorFromEvent);
     } else {
       const std::map<uint32_t, uint32_t>* mapFromEvent = getProduct<std::map<uint32_t, uint32_t> >(event, tag);
       if (mapFromEvent) {
         std::vector<uint32_t>* newVector = new std::vector<uint32_t>(FED_ID_MAX + 1, 0);
         if (mapKeyIsByFedID) {
           for (std::map<uint32_t, uint32_t>::const_iterator iIdValue = mapFromEvent->begin();
                iIdValue != mapFromEvent->end();
                ++iIdValue) {
             newVector->at(iIdValue->first) = iIdValue->second;
           }
         } else {
           SpyUtilities::fillFEDMajorities(*mapFromEvent, *newVector);
         }
         //  std::cout << " -- Map " << tag << std::endl;
         //  for (uint32_t lIt= 0;
         //       lIt < newVector->size();
         //       lIt++) {
         //    std::cout << lIt << " " << newVector->at(lIt) << std::endl;
         //  }
         //vector was allocated here so, will need to be deleted when finished with
         CountersPtr newCountersPtr(new CountersWrapper(newVector, true));
         return newCountersPtr;
       } else {
         throw cms::Exception(mlLabel_) << "Unable to get product " << tag << " from spy event";
       }
     }
   }

void sistrip::SpyEventMatcher::getMatchedCollections	(	const uint32_t	eventId,
		const uint8_t	apvAddress,
		const SpyEventList *	matchingEvents,
		const SiStripFedCabling &	cabling,
		SpyDataCollections &	collectionsToCreate
	)

Definition at line 251 of file SiStripSpyEventMatcher.cc.

                                                                                        {
     if (!matchingEvents)
       return;
     size_t fileNameHash = 0U;
     FEDRawDataCollection outputRawData;
     MatchingOutput mo(outputRawData);
     source_->loopSpecified(
         *eventPrincipal_,
         fileNameHash,
         matchingEvents->begin(),
         matchingEvents->end(),
         [&](auto const& iE, auto const&) { this->getCollections(iE, eventId, apvAddress, cabling, mo); });
     collectionsToCreate = SpyDataCollections(mo.outputRawData_,
                                              mo.outputTotalEventCounters_,
                                              mo.outputL1ACounters_,
                                              mo.outputAPVAddresses_,
                                              mo.outputScopeDigisVector_.get(),
                                              mo.outputPayloadDigisVector_.get(),
                                              mo.outputReorderedDigisVector_.get(),
                                              mo.outputVirginRawDigisVector_.get());
   }

template<class T >

const T * sistrip::SpyEventMatcher::getProduct	(	const edm::EventPrincipal &	event,
		const edm::InputTag &	tag
	)

staticprivate

Definition at line 178 of file SiStripSpyEventMatcher.h.

References edmPickEvents::event, LogDebug, mlLabel_, and GlobalPosition_Frontier_DevDB_cff::tag.

                                                                                              {
     LogDebug(mlLabel_) << "Retrieving product " << tag;
     // Note: The third argument to getProductByTag can be a nullptr
     // as long as unscheduled execution of an EDProducer cannot occur
     // as a result of this function call (and with the current implementation
     // of SpyEventMatcher unscheduled execution never happens).
     const std::shared_ptr<const edm::Wrapper<T> > productWrapper = edm::getProductByTag<T>(event, tag, nullptr);
     if (productWrapper) {
       return productWrapper->product();
     } else {
       return nullptr;
     }
   }

void sistrip::SpyEventMatcher::initialize ( )

Definition at line 79 of file SiStripSpyEventMatcher.cc.

References addNextEventToMap(), eventMatches_, eventPrincipal_, LogDebug, SiStripPI::max, mlLabel_, source_, and contentValuesCheck::ss.

                                    {
     size_t fileNameHash = 0U;
     //add spy events to the map until there are none left
     source_->loopOverEvents(
         *eventPrincipal_,
         fileNameHash,
         std::numeric_limits<size_t>::max(),
         [this](auto const& iE, auto const&) {
           this->addNextEventToMap(iE);
           return true;
         },
         nullptr,
         nullptr,
         false);
     //debug
     std::ostringstream ss;
     ss << "Events with possible matches (eventID,apvAddress): ";
     for (std::map<EventKey, SpyEventList>::const_iterator iSpyEvent = eventMatches_.begin();
          iSpyEvent != eventMatches_.end();
          ++iSpyEvent) {
       ss << "(" << iSpyEvent->first.eventId() << "," << uint16_t(iSpyEvent->first.apvAddress()) << ") ";
     }
     LogDebug(mlLabel_) << ss.str();
   }

const SpyEventMatcher::SpyEventList * sistrip::SpyEventMatcher::matchesForEvent	(	const uint32_t	eventId,
		const uint8_t	apvAddress
	)		const

Definition at line 151 of file SiStripSpyEventMatcher.cc.

References eventMatches_, LogDebug, mlLabel_, and contentValuesCheck::ss.

                                                                                                         {
     EventKey eventKey(eventId, apvAddress);
     std::map<EventKey, SpyEventList>::const_iterator iMatch = eventMatches_.find(eventKey);
     if (iMatch == eventMatches_.end()) {
       LogDebug(mlLabel_) << "No match found for event " << eventId << " with APV address " << uint16_t(apvAddress);
       return nullptr;
     } else {
       std::ostringstream ss;
       ss << "Found matches to event " << eventId << " with address " << uint16_t(apvAddress) << " in spy events ";
       for (SpyEventList::const_iterator iMatchingSpyEvent = iMatch->second.begin();
            iMatchingSpyEvent != iMatch->second.end();
            ++iMatchingSpyEvent) {
         ss << iMatchingSpyEvent->event() << " ";
       }
       LogInfo(mlLabel_) << ss.str();
       return &(iMatch->second);
     }
   }

void sistrip::SpyEventMatcher::mergeMatchingData	(	const std::set< uint16_t > &	matchingFeds,
		const FEDRawDataCollection &	inputRawData,
		SpyEventMatcher::CountersPtr	inputTotalEventCounters,
		SpyEventMatcher::CountersPtr	inputL1ACounters,
		SpyEventMatcher::CountersPtr	inputAPVAddresses,
		const edm::DetSetVector< SiStripRawDigi > *	inputScopeDigis,
		const edm::DetSetVector< SiStripRawDigi > *	inputPayloadDigis,
		const edm::DetSetVector< SiStripRawDigi > *	inputReorderedDigis,
		const edm::DetSetVector< SiStripRawDigi > *	inputVirginRawDigis,
		FEDRawDataCollection &	outputRawData,
		std::vector< uint32_t > &	outputTotalEventCounters,
		std::vector< uint32_t > &	outputL1ACounters,
		std::vector< uint32_t > &	outputAPVAddresses,
		std::vector< edm::DetSet< SiStripRawDigi > > *	outputScopeDigisVector,
		std::vector< edm::DetSet< SiStripRawDigi > > *	outputPayloadDigisVector,
		std::vector< edm::DetSet< SiStripRawDigi > > *	outputReorderedDigisVector,
		std::vector< edm::DetSet< SiStripRawDigi > > *	outputVirginRawDigisVector,
		const SiStripFedCabling &	cabling
	)

staticprivate

Definition at line 299 of file SiStripSpyEventMatcher.cc.

References officialStyle::chan, FEDRawData::data(), edm::DetSetVector< T >::end(), sistrip::FEDCH_PER_FED, SiStripFedCabling::fedConnections(), FEDRawDataCollection::FEDData(), l1tstage2_dqm_sourceclient-live_cfg::fedId, sistrip::SpyUtilities::fedIndex(), edm::DetSetVector< T >::find(), sistrip::invalid_, LogDebug, mlLabel_, and FEDRawData::size().

Referenced by getCollections().

                                                                             {
     //reserve space in vectors
     if (inputScopeDigis) {
       outputScopeDigisVector->reserve(outputScopeDigisVector->size() +
                                       matchingFeds.size() *
                                           FEDCH_PER_FED);  //maximum number of channels on matching FEDs
     }
     if (inputPayloadDigis) {
       outputPayloadDigisVector->reserve(outputPayloadDigisVector->size() + matchingFeds.size() * FEDCH_PER_FED);
     }
     if (inputReorderedDigis) {
       outputReorderedDigisVector->reserve(outputReorderedDigisVector->size() + matchingFeds.size() * FEDCH_PER_FED);
     }
     if (inputVirginRawDigis) {
       outputVirginRawDigisVector->reserve(outputVirginRawDigisVector->size() +
                                           matchingFeds.size() * FEDCH_PER_FED /
                                               2);  //maximum number of dets on matching FEDs
     }
     //copy the data into output collections
     std::set<uint32_t> usedDetIds;
     for (std::set<uint16_t>::const_iterator iFedId = matchingFeds.begin(); iFedId != matchingFeds.end(); ++iFedId) {
       const uint32_t fedId = *iFedId;
       LogDebug(mlLabel_) << "Copying data for FED " << fedId;
       if (inputRawData.FEDData(fedId).size() && inputRawData.FEDData(fedId).data()) {
         outputRawData.FEDData(fedId) = inputRawData.FEDData(fedId);
       }
       outputTotalEventCounters[fedId] = (*inputTotalEventCounters)[fedId];
       outputL1ACounters[fedId] = (*inputL1ACounters)[fedId];
       outputAPVAddresses[fedId] = (*inputAPVAddresses)[fedId];
       for (uint8_t chan = 0; chan < FEDCH_PER_FED; ++chan) {
         uint32_t fedIndex = ((fedId & sistrip::invalid_) << 16) | (chan & sistrip::invalid_);
         ;
         if (inputScopeDigis) {
           edm::DetSetVector<SiStripRawDigi>::const_iterator iScopeDigis = inputScopeDigis->find(fedIndex);
           if (iScopeDigis != inputScopeDigis->end()) {
             outputScopeDigisVector->push_back(*iScopeDigis);
           }
         }
         if (inputPayloadDigis) {
           edm::DetSetVector<SiStripRawDigi>::const_iterator iPayloadDigis = inputPayloadDigis->find(fedIndex);
           if (iPayloadDigis != inputPayloadDigis->end()) {
             outputPayloadDigisVector->push_back(*iPayloadDigis);
           }
         }
         if (inputReorderedDigis) {
           edm::DetSetVector<SiStripRawDigi>::const_iterator iReorderedDigis = inputReorderedDigis->find(fedIndex);
           if (iReorderedDigis != inputReorderedDigis->end()) {
             outputReorderedDigisVector->push_back(*iReorderedDigis);
           }
         }
       }
       if (inputVirginRawDigis) {
         std::set<uint32_t> fedDetIds;
         auto conns = cabling.fedConnections(fedId);
         for (auto iConn = conns.begin(); iConn != conns.end(); ++iConn) {
           if (!iConn->isConnected())
             continue;
           const uint32_t detId = iConn->detId();
           if (usedDetIds.find(detId) != usedDetIds.end()) {
             LogError(mlLabel_) << "Duplicate DetID found " << detId << " skipping data for this Det from FED " << fedId;
             continue;
           }
           fedDetIds.insert(iConn->detId());
         }
         usedDetIds.insert(fedDetIds.begin(), fedDetIds.end());
         for (std::set<uint32_t>::const_iterator iDetId = fedDetIds.begin(); iDetId != fedDetIds.end(); ++iDetId) {
           edm::DetSetVector<SiStripRawDigi>::const_iterator iVirginRawDigis = inputVirginRawDigis->find(*iDetId);
           if (iVirginRawDigis != inputVirginRawDigis->end()) {
             outputVirginRawDigisVector->push_back(*iVirginRawDigis);
           }
         }
       }
     }
   }