CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | Static Private Attributes
L1TMP7ZeroSupp Class Reference

#include <L1TMP7ZeroSupp.h>

Inheritance diagram for L1TMP7ZeroSupp:
DQMEDAnalyzer edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 L1TMP7ZeroSupp (const edm::ParameterSet &ps)
 
 ~L1TMP7ZeroSupp () override
 
- Public Member Functions inherited from DQMEDAnalyzer
void accumulate (edm::Event const &ev, edm::EventSetup const &es) final
 
virtual void analyze (edm::Event const &, edm::EventSetup const &)
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual void bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer ()
 
 DQMEDAnalyzer (DQMEDAnalyzer const &)=delete
 
 DQMEDAnalyzer (DQMEDAnalyzer &&)=delete
 
void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override
 
void endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &run, edm::EventSetup const &setup) override
 
void endRunProduce (edm::Run &run, edm::EventSetup const &setup) override
 
 ~DQMEDAnalyzer () override=default
 
- Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
 EDProducer ()=default
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~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
 
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 () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- 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

void analyze (const edm::Event &e, const edm::EventSetup &c) override
 
void beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, const edm::Run &, const edm::EventSetup &) override
 
void dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Private Types

enum  binlabels {
  EVTS =0, EVTSGOOD, EVTSBAD, BLOCKS,
  ZSBLKSGOOD, ZSBLKSBAD, ZSBLKSBADFALSEPOS, ZSBLKSBADFALSENEG,
  BXBLOCKS, ZSBXBLKSGOOD, ZSBXBLKSBAD, ZSBXBLKSBADFALSEPOS,
  ZSBXBLKSBADFALSENEG, NBINLABELS
}
 
enum  ratioBinlabels {
  REVTS =0, RBLKS, RBLKSFALSEPOS, RBLKSFALSENEG,
  RBXBLKS, RBXBLKSFALSEPOS, RBXBLKSFALSENEG, RNBINLABELS
}
 

Private Member Functions

void bookCapIdHistograms (DQMStore::IBooker &ibooker, const unsigned int &id)
 

Private Attributes

int amc13HeaderSize_
 
int amc13TrailerSize_
 
int amcHeaderSize_
 
int amcTrailerSize_
 
MonitorElementcapIds_
 
int dataInvFlagMask_
 
std::vector< unsigned int > definedMaskCapIds_
 
std::map< unsigned int, MonitorElement * > errorSummaryDenMap_
 
std::map< unsigned int, MonitorElement * > errorSummaryNumMap_
 
edm::EDGetTokenT< FEDRawDataCollectionfedDataToken_
 
std::vector< int > fedIds_
 
std::vector< std::vector< int > > masks_
 
int maxFedReadoutSize_
 
std::string monitorDir_
 
int newZsFlagMask_
 
std::map< unsigned int, MonitorElement * > readoutSizeNoZSMap_
 
std::map< unsigned int, MonitorElement * > readoutSizeZSExpectedMap_
 
std::map< unsigned int, MonitorElement * > readoutSizeZSMap_
 
int slinkHeaderSize_
 
int slinkTrailerSize_
 
bool verbose_
 
std::map< unsigned int, MonitorElement * > zeroSuppValMap_
 
bool zsEnabled_
 
int zsFlagMask_
 

Static Private Attributes

static const unsigned int maxMasks_ = 16
 

Additional Inherited Members

- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase 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
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMTokenlumiToken_
 
edm::EDPutTokenT< DQMTokenrunToken_
 

Detailed Description

Definition at line 22 of file L1TMP7ZeroSupp.h.

Member Enumeration Documentation

Constructor & Destructor Documentation

L1TMP7ZeroSupp::L1TMP7ZeroSupp ( const edm::ParameterSet ps)

Definition at line 5 of file L1TMP7ZeroSupp.cc.

References gather_cfg::cout, TauDecayModes::dec, definedMaskCapIds_, edm::ParameterSet::exists(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, masks_, maxMasks_, AlCaHLTBitMon_QueryRunRegistry::string, verbose_, and ~L1TMP7ZeroSupp().

6  : fedDataToken_(consumes<FEDRawDataCollection>(ps.getParameter<edm::InputTag>("rawData"))),
7  zsEnabled_(ps.getUntrackedParameter<bool>("zsEnabled")),
8  fedIds_(ps.getParameter<std::vector<int>>("fedIds")),
9  slinkHeaderSize_(ps.getUntrackedParameter<int>("lenSlinkHeader")),
10  slinkTrailerSize_(ps.getUntrackedParameter<int>("lenSlinkTrailer")),
11  amc13HeaderSize_(ps.getUntrackedParameter<int>("lenAMC13Header")),
12  amc13TrailerSize_(ps.getUntrackedParameter<int>("lenAMC13Trailer")),
13  amcHeaderSize_(ps.getUntrackedParameter<int>("lenAMCHeader")),
14  amcTrailerSize_(ps.getUntrackedParameter<int>("lenAMCTrailer")),
15  newZsFlagMask_(ps.getUntrackedParameter<int>("newZsFlagMask")),
16  zsFlagMask_(ps.getUntrackedParameter<int>("zsFlagMask")),
17  dataInvFlagMask_(ps.getUntrackedParameter<int>("dataInvFlagMask")),
18  maxFedReadoutSize_(ps.getUntrackedParameter<int>("maxFEDReadoutSize")),
20  verbose_(ps.getUntrackedParameter<bool>("verbose"))
21 {
22  std::vector<int> zeroMask(6, 0);
23  masks_.reserve(maxMasks_);
24  for (unsigned int i = 0; i < maxMasks_; ++i) {
25  std::string maskCapIdStr{"maskCapId"+std::to_string(i)};
26  masks_.push_back(ps.getUntrackedParameter<std::vector<int>>(maskCapIdStr, zeroMask));
27  // which masks are defined?
28  if (ps.exists(maskCapIdStr)) {
29  definedMaskCapIds_.push_back(i);
30  }
31  }
32  if (verbose_) {
33  // check masks
34  std::cout << "masks" << std::endl;
35  for (unsigned int i = 0; i < maxMasks_; ++i) {
36  std::cout << "caption ID" << i << ":" << std::endl;
37  for (const auto& maskIt: masks_.at(i)) {
38  std::cout << std::hex << std::setw(8) << std::setfill('0') << maskIt << std::dec << std::endl;
39  }
40  }
41  std::cout << "----------" << std::endl;
42  }
43 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
edm::EDGetTokenT< FEDRawDataCollection > fedDataToken_
bool exists(std::string const &parameterName) const
checks if a parameter exists
std::string monitorDir_
std::vector< unsigned int > definedMaskCapIds_
std::vector< std::vector< int > > masks_
static const unsigned int maxMasks_
std::vector< int > fedIds_
L1TMP7ZeroSupp::~L1TMP7ZeroSupp ( )
overridedefault

Referenced by L1TMP7ZeroSupp().

Member Function Documentation

void L1TMP7ZeroSupp::analyze ( const edm::Event e,
const edm::EventSetup c 
)
overrideprotected

Definition at line 144 of file L1TMP7ZeroSupp.cc.

References amc13HeaderSize_, amc13TrailerSize_, amcHeaderSize_, amcTrailerSize_, groupFilesInBlocks::block, BLOCKS, BXBLOCKS, FEDHeader::bxID(), capIds_, FEDHeader::check(), FEDTrailer::check(), gather_cfg::cout, FEDTrailer::crc(), FEDRawData::data(), data, dataInvFlagMask_, TauDecayModes::dec, definedMaskCapIds_, end, errorSummaryDenMap_, errorSummaryNumMap_, EVTS, EVTSBAD, EVTSGOOD, FEDTrailer::evtStatus(), FEDRawDataCollection::FEDData(), fedDataToken_, l1t::stage2::layer2::fedId, fedIds_, MonitorElement::Fill(), FEDTrailer::fragmentLength(), edm::Event::getByToken(), RecoTauValidation_cfi::header, triggerObjects_cff::id, edm::HandleBase::isValid(), FEDHeader::lvl1ID(), masks_, maxMasks_, newZsFlagMask_, amc13::Packet::parse(), amc13::Packet::payload(), jets_cff::payload, RBLKS, RBLKSFALSENEG, RBLKSFALSEPOS, RBXBLKS, RBXBLKSFALSENEG, RBXBLKSFALSEPOS, readoutSizeNoZSMap_, readoutSizeZSExpectedMap_, readoutSizeZSMap_, REVTS, FEDRawData::size(), slinkHeaderSize_, slinkTrailerSize_, FEDHeader::sourceID(), command_line::start, FEDHeader::triggerType(), FEDTrailer::ttsBits(), verbose_, FEDHeader::version(), zeroSuppValMap_, ZSBLKSBAD, ZSBLKSBADFALSENEG, ZSBLKSBADFALSEPOS, ZSBLKSGOOD, ZSBXBLKSBAD, ZSBXBLKSBADFALSENEG, ZSBXBLKSBADFALSEPOS, ZSBXBLKSGOOD, zsEnabled_, and zsFlagMask_.

144  {
145 
146  if (verbose_) edm::LogInfo("L1TDQM") << "L1TMP7ZeroSupp: analyze..." << std::endl;
147 
149  e.getByToken(fedDataToken_, feds);
150 
151  if (!feds.isValid()) {
152  edm::LogError("L1TDQM") << "Cannot analyse: no FEDRawDataCollection found";
153  return;
154  }
155 
158  for (const auto &id: definedMaskCapIds_) {
159  zeroSuppValMap_[id]->Fill(EVTS);
160  errorSummaryDenMap_[id]->Fill(REVTS);
161  }
162 
163  std::map<unsigned int, bool> evtGood;
164  evtGood[maxMasks_] = true;
165  for (const auto &id: definedMaskCapIds_) {
166  evtGood[id] = true;
167  }
168  unsigned valid_count = 0;
169  for (const auto& fedId: fedIds_) {
170  const FEDRawData& l1tRcd = feds->FEDData(fedId);
171 
172  unsigned int fedDataSize = l1tRcd.size();
173  std::map<unsigned int, unsigned int> readoutSizeNoZSMap;
174  std::map<unsigned int, unsigned int> readoutSizeZSMap;
175  std::map<unsigned int, unsigned int> readoutSizeZSExpectedMap;
176  readoutSizeNoZSMap[maxMasks_] = 0;
177  readoutSizeZSMap[maxMasks_] = 0;
178  readoutSizeZSExpectedMap[maxMasks_] = 0;
179  for (const auto &id: definedMaskCapIds_) {
180  readoutSizeNoZSMap[id] = 0;
181  readoutSizeZSMap[id] = 0;
182  readoutSizeZSExpectedMap[id] = 0;
183  }
184 
185  edm::LogInfo("L1TDQM") << "Found FEDRawDataCollection with ID " << fedId << " and size " << l1tRcd.size();
186 
188  if (l1tRcd.size() > 0) {
189  edm::LogError("L1TDQM") << "Cannot analyse: invalid L1T raw data (size = " << l1tRcd.size() << ") for ID " << fedId << ".";
190  }
191  continue;
192  } else {
193  valid_count++;
194  }
195 
196  const unsigned char *data = l1tRcd.data();
197  FEDHeader header(data);
198 
199  if (header.check()) {
200  edm::LogInfo("L1TDQM") << "Found SLink header:" << " Trigger type " << header.triggerType() << " L1 event ID " << header.lvl1ID() << " BX Number " << header.bxID() << " FED source " << header.sourceID() << " FED version " << header.version();
201  } else {
202  edm::LogWarning("L1TDQM") << "Did not find a SLink header!";
203  }
204 
205  FEDTrailer trailer(data + (l1tRcd.size() - slinkTrailerSize_));
206 
207  if (trailer.check()) {
208  edm::LogInfo("L1TDQM") << "Found SLink trailer:" << " Length " << trailer.fragmentLength() << " CRC " << trailer.crc() << " Status " << trailer.evtStatus() << " Throttling bits " << trailer.ttsBits();
209  } else {
210  edm::LogWarning("L1TDQM") << "Did not find a SLink trailer!";
211  }
212 
213  amc13::Packet packet;
214  if (!packet.parse(
215  (const uint64_t*) data,
216  (const uint64_t*) (data + slinkHeaderSize_),
217  (l1tRcd.size() - slinkHeaderSize_ - slinkTrailerSize_) / 8,
218  header.lvl1ID(),
219  header.bxID())) {
220  edm::LogError("L1TDQM") << "Could not extract AMC13 Packet.";
221  return;
222  }
223 
224  for (auto& amc: packet.payload()) {
225  if (amc.size() == 0)
226  continue;
227 
228  auto payload64 = amc.data();
229  auto start = (const uint32_t*) payload64.get();
230  // Want to have payload size in 32 bit words, but AMC measures
231  // it in 64 bit words -> factor 2.
232  const uint32_t * end = start + (amc.size() * 2);
233 
234  auto payload = std::make_unique<l1t::MP7Payload>(start, end, false);
235 
236  // getBlock() returns a non-null unique_ptr on success
237  std::unique_ptr<l1t::Block> block;
238  while ((block = payload->getBlock()) != nullptr) {
239  if (verbose_) {
240  std::cout << ">>> check zero suppression for block <<<" << std::endl
241  << "hdr: " << std::hex << std::setw(8) << std::setfill('0') << block->header().raw() << std::dec
242  << " (ID " << block->header().getID() << ", size " << block->header().getSize()
243  << ", CapID 0x" << std::hex << std::setw(2) << std::setfill('0') << block->header().getCapID()
244  << ", flags 0x" << std::hex << std::setw(2) << std::setfill('0') << block->header().getFlags()
245  << ")" << std::dec << std::endl;
246  for (const auto& word: block->payload()) {
247  std::cout << "data: " << std::hex << std::setw(8) << std::setfill('0') << word << std::dec << std::endl;
248  }
249  }
250 
251  unsigned int blockCapId = block->header().getCapID();
252  unsigned int blockSize = block->header().getSize() * 4; // times 4 to get the size in byte
253  unsigned int blockHeaderSize = sizeof(block->header().raw());
254  unsigned int blockHeaderFlags = block->header().getFlags();
255  bool newZsFlagSet = ((blockHeaderFlags & newZsFlagMask_) != 0); // use the per-BX ZS
256  bool blockZsFlagSet = newZsFlagSet ? true : ((blockHeaderFlags & zsFlagMask_) != 0); // ZS validation flag for whole block
257  bool dataInvertFlagSet = newZsFlagSet && ((blockHeaderFlags & dataInvFlagMask_) != 0); // invert the data before applying the mask
258 
259  capIds_->Fill(blockCapId);
260 
261  // fill the denominator histograms
262  bool capIdDefined = false;
267  if (zeroSuppValMap_.find(blockCapId) != zeroSuppValMap_.end()) {
268  capIdDefined = true;
269  zeroSuppValMap_[blockCapId]->Fill(BLOCKS);
270  errorSummaryDenMap_[blockCapId]->Fill(RBLKS);
271  errorSummaryDenMap_[blockCapId]->Fill(RBLKSFALSEPOS);
272  errorSummaryDenMap_[blockCapId]->Fill(RBLKSFALSENEG);
273  }
274 
275  auto totalBlockSize = blockHeaderSize;
276  if (!newZsFlagSet) {
277  totalBlockSize += blockSize;
278  }
279  auto totalBlockSizeExpected = totalBlockSize;
280  auto totalBlockSizeNoZS = blockHeaderSize + blockSize;
281 
282  auto bxBlocks = block->getBxBlocks(6, newZsFlagSet); // 6 32 bit MP7 payload words per BX
283 
284  // check all BX blocks
285  bool allToSuppress = true;
286  for (const auto& bxBlock: bxBlocks) {
287  bool toSuppress = false;
288  bool bxZsFlagSet = ((bxBlock.header().getFlags() & zsFlagMask_) != 0); // ZS validation flag
289 
290  // check if this bxblock should be suppressed
291  unsigned int wordcounter = 0;
292  unsigned int wordsum = 0;
293  for (const auto& word: bxBlock.payload()) {
294  if (dataInvertFlagSet) {
295  wordsum += masks_[blockCapId].at(wordcounter%6) & (~word);
296  } else {
297  wordsum += masks_[blockCapId].at(wordcounter%6) & word;
298  }
299  if (verbose_) {
300  std::cout << "word: " << std::hex << std::setw(8) << std::setfill('0') << word << std::dec
301  << ", maskword" << wordcounter%6 << ": " << std::hex << std::setw(8) << std::setfill('0')
302  << masks_[blockCapId].at(wordcounter%6) << std::dec << ", wordsum: " << wordsum << std::endl;
303  }
304  if (wordsum > 0) {
305  if (verbose_) {
306  std::cout << "wordsum not 0: this BX block should be kept" << std::endl;
307  }
308  break;
309  }
310  ++wordcounter;
311  }
312  // the sum of payload words must be 0 for correct ZS
313  if (wordsum == 0 && zsEnabled_) {
314  toSuppress = true;
315  if (verbose_) {
316  std::cout << "wordsum == 0: this BX block should be zero suppressed" << std::endl;
317  }
318  }
319  // update the overall block status
320  allToSuppress = allToSuppress && toSuppress;
321 
322  // only fill the BX related things for the per-BX ZS
323  if (newZsFlagSet) {
324  // the ZS flag of the block is the AND of all BX block ZS flags
325  blockZsFlagSet = blockZsFlagSet && bxZsFlagSet;
326 
327  // fill the BX related bins of the denominator histogram
332  if (capIdDefined) {
333  zeroSuppValMap_[blockCapId]->Fill(BXBLOCKS);
334  errorSummaryDenMap_[blockCapId]->Fill(RBXBLKS);
335  errorSummaryDenMap_[blockCapId]->Fill(RBXBLKSFALSEPOS);
336  errorSummaryDenMap_[blockCapId]->Fill(RBXBLKSFALSENEG);
337  }
338 
339  unsigned int totalBxBlockSize = bxBlock.getSize() * 4 + sizeof(bxBlock.header().raw()); // times 4 to get the size in byte
340  // check if zero suppression flag agrees for the BX block
341  if (toSuppress && bxZsFlagSet) {
342  if (verbose_) std::cout << "GOOD BX block with ZS flag true" << std::endl;
344  if (capIdDefined) {
345  zeroSuppValMap_[blockCapId]->Fill(ZSBXBLKSGOOD);
346  }
347  } else if (!toSuppress && !bxZsFlagSet) {
348  if (verbose_) std::cout << "GOOD BX block with ZS flag false" << std::endl;
349  totalBlockSize += totalBxBlockSize;
350  totalBlockSizeExpected += totalBxBlockSize;
352  if (capIdDefined) {
353  zeroSuppValMap_[blockCapId]->Fill(ZSBXBLKSGOOD);
354  }
355  } else if (!toSuppress && bxZsFlagSet) {
356  if (verbose_) std::cout << "BAD BX block with ZS flag true" << std::endl;
357  totalBlockSizeExpected += totalBxBlockSize;
362  evtGood[maxMasks_] = false;
363  if (capIdDefined) {
364  zeroSuppValMap_[blockCapId]->Fill(ZSBXBLKSBAD);
365  zeroSuppValMap_[blockCapId]->Fill(ZSBXBLKSBADFALSEPOS);
366  errorSummaryNumMap_[blockCapId]->Fill(RBXBLKS);
367  errorSummaryNumMap_[blockCapId]->Fill(RBXBLKSFALSEPOS);
368  evtGood[blockCapId] = false;
369  }
370  } else {
371  if (verbose_) std::cout << "BAD BX block with ZS flag false" << std::endl;
372  totalBlockSize += totalBxBlockSize;
377  evtGood[maxMasks_] = false;
378  if (capIdDefined) {
379  zeroSuppValMap_[blockCapId]->Fill(ZSBXBLKSBAD);
380  zeroSuppValMap_[blockCapId]->Fill(ZSBXBLKSBADFALSENEG);
381  errorSummaryNumMap_[blockCapId]->Fill(RBXBLKS);
382  errorSummaryNumMap_[blockCapId]->Fill(RBXBLKSFALSENEG);
383  evtGood[blockCapId] = false;
384  }
385  }
386  }
387  }
388 
389  readoutSizeNoZSMap[maxMasks_] += totalBlockSizeNoZS;
390  if (capIdDefined) {
391  readoutSizeNoZSMap[blockCapId] += totalBlockSizeNoZS;
392  }
393 
394  // check if zero suppression flag agrees for the whole block
395  if (allToSuppress && blockZsFlagSet) {
396  if (verbose_) std::cout << "GOOD block with ZS flag true" << std::endl;
398  if (capIdDefined) {
399  zeroSuppValMap_[blockCapId]->Fill(ZSBLKSGOOD);
400  }
401  } else if (!allToSuppress && !blockZsFlagSet) {
402  if (verbose_) std::cout << "GOOD block with ZS flag false" << std::endl;
404  readoutSizeZSMap[maxMasks_] += totalBlockSize;
405  readoutSizeZSExpectedMap[maxMasks_] += totalBlockSizeExpected;
406  if (capIdDefined) {
407  zeroSuppValMap_[blockCapId]->Fill(ZSBLKSGOOD);
408  readoutSizeZSMap[blockCapId] += totalBlockSize;
409  readoutSizeZSExpectedMap[blockCapId] += totalBlockSizeExpected;
410  }
411  } else if (!allToSuppress && blockZsFlagSet) {
412  if (verbose_) std::cout << "BAD block with ZS flag true" << std::endl;
417  readoutSizeZSExpectedMap[maxMasks_] += totalBlockSizeExpected;
418  evtGood[maxMasks_] = false;
419  if (capIdDefined) {
420  zeroSuppValMap_[blockCapId]->Fill(ZSBLKSBAD);
421  zeroSuppValMap_[blockCapId]->Fill(ZSBLKSBADFALSEPOS);
422  errorSummaryNumMap_[blockCapId]->Fill(RBLKS);
423  errorSummaryNumMap_[blockCapId]->Fill(RBLKSFALSEPOS);
424  readoutSizeZSExpectedMap[blockCapId] += totalBlockSizeExpected;
425  evtGood[blockCapId] = false;
426  }
427  } else {
428  if (verbose_) std::cout << "BAD block with ZS flag false" << std::endl;
433  readoutSizeZSMap[maxMasks_] += totalBlockSize;
434  evtGood[maxMasks_] = false;
435  if (capIdDefined) {
436  zeroSuppValMap_[blockCapId]->Fill(ZSBLKSBAD);
437  zeroSuppValMap_[blockCapId]->Fill(ZSBLKSBADFALSENEG);
438  errorSummaryNumMap_[blockCapId]->Fill(RBLKS);
439  errorSummaryNumMap_[blockCapId]->Fill(RBLKSFALSENEG);
440  readoutSizeZSMap[blockCapId] += totalBlockSize;
441  evtGood[blockCapId] = false;
442  }
443  }
444  }
445  }
446  if (verbose_) {
447  std::cout << "FED data size: " << fedDataSize << " bytes" << std::endl;
448  std::cout << "Payload size no ZS: " << readoutSizeNoZSMap[maxMasks_] << " bytes" << std::endl;
449  std::cout << "Payload size ZS: " << readoutSizeZSMap[maxMasks_] << " bytes" << std::endl;
450  std::cout << "Payload size expected ZS: " << readoutSizeZSExpectedMap[maxMasks_] << " bytes" << std::endl;
451  std::cout << "Filled readout size ZS with headers: " << readoutSizeZSMap[maxMasks_] + fedDataSize - readoutSizeNoZSMap[maxMasks_] << " bytes" << std::endl;
452  std::cout << "Filled expected readout size ZS with headers: " << readoutSizeZSExpectedMap[maxMasks_] + fedDataSize - readoutSizeNoZSMap[maxMasks_] << " bytes" << std::endl;
453  }
454  readoutSizeNoZSMap_[maxMasks_]->Fill(fedDataSize);
455  readoutSizeZSMap_[maxMasks_]->Fill(readoutSizeZSMap[maxMasks_] + fedDataSize - readoutSizeNoZSMap[maxMasks_]);
456  readoutSizeZSExpectedMap_[maxMasks_]->Fill(readoutSizeZSExpectedMap[maxMasks_] + fedDataSize - readoutSizeNoZSMap[maxMasks_]);
457  for (const auto &id: definedMaskCapIds_) {
458  readoutSizeNoZSMap_[id]->Fill(readoutSizeNoZSMap[id]);
459  readoutSizeZSMap_[id]->Fill(readoutSizeZSMap[id]);
460  readoutSizeZSExpectedMap_[id]->Fill(readoutSizeZSExpectedMap[id]);
461  }
462  }
463 
464  if (evtGood[maxMasks_]) {
466  } else {
469  }
470  for (const auto &id: definedMaskCapIds_) {
471  if (evtGood[id]) {
472  zeroSuppValMap_[id]->Fill(EVTSGOOD);
473  } else {
474  zeroSuppValMap_[id]->Fill(EVTSBAD);
475  errorSummaryNumMap_[id]->Fill(REVTS);
476  }
477  }
478 }
Definition: start.py:1
std::vector< amc::Packet > payload() const
Definition: AMC13Spec.h:80
std::map< unsigned int, MonitorElement * > readoutSizeZSExpectedMap_
std::map< unsigned int, MonitorElement * > zeroSuppValMap_
edm::EDGetTokenT< FEDRawDataCollection > fedDataToken_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:519
std::map< unsigned int, MonitorElement * > readoutSizeNoZSMap_
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
void Fill(long long x)
std::map< unsigned int, MonitorElement * > errorSummaryNumMap_
std::vector< unsigned int > definedMaskCapIds_
MonitorElement * capIds_
const FEDRawData & FEDData(int fedid) const
retrieve data for fed
std::map< unsigned int, MonitorElement * > errorSummaryDenMap_
std::vector< std::vector< int > > masks_
#define end
Definition: vmac.h:39
static const unsigned int maxMasks_
bool isValid() const
Definition: HandleBase.h:74
unsigned long long uint64_t
Definition: Time.h:15
bool parse(const uint64_t *start, const uint64_t *data, unsigned int size, unsigned int lv1, unsigned int bx, bool legacy_mc=false, bool mtf7_mode=false)
Definition: AMC13Spec.cc:70
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:28
std::vector< int > fedIds_
Definition: AMCSpec.h:8
std::map< unsigned int, MonitorElement * > readoutSizeZSMap_
void L1TMP7ZeroSupp::beginLuminosityBlock ( const edm::LuminosityBlock ,
const edm::EventSetup  
)
overrideprotected

Definition at line 72 of file L1TMP7ZeroSupp.cc.

72 {}
void L1TMP7ZeroSupp::bookCapIdHistograms ( DQMStore::IBooker ibooker,
const unsigned int &  id 
)
private

Definition at line 88 of file L1TMP7ZeroSupp.cc.

References BLOCKS, DQMStore::IBooker::book1D(), BXBLOCKS, errorSummaryDenMap_, errorSummaryNumMap_, EVTS, EVTSBAD, EVTSGOOD, triggerObjects_cff::id, maxFedReadoutSize_, maxMasks_, NBINLABELS, RBLKS, RBLKSFALSENEG, RBLKSFALSEPOS, RBXBLKS, RBXBLKSFALSENEG, RBXBLKSFALSEPOS, readoutSizeNoZSMap_, readoutSizeZSExpectedMap_, readoutSizeZSMap_, REVTS, RNBINLABELS, AlCaHLTBitMon_QueryRunRegistry::string, zeroSuppValMap_, ZSBLKSBAD, ZSBLKSBADFALSENEG, ZSBLKSBADFALSEPOS, ZSBLKSGOOD, ZSBXBLKSBAD, ZSBXBLKSBADFALSENEG, ZSBXBLKSBADFALSEPOS, and ZSBXBLKSGOOD.

Referenced by bookHistograms().

88  {
89  std::string summaryTitleText = "Zero suppression validation summary";
90  std::string sizeTitleText;
91  if (id == maxMasks_) {
92  sizeTitleText = "FED readout ";
93  } else {
94  summaryTitleText = summaryTitleText+", caption id "+std::to_string(id);
95  sizeTitleText = "cumulated caption id "+std::to_string(id)+" block ";
96  }
97 
98  zeroSuppValMap_[id] = ibooker.book1D("zeroSuppVal", summaryTitleText, NBINLABELS, 0, NBINLABELS);
99  zeroSuppValMap_[id]->setAxisTitle("ZS status", 1);
100  zeroSuppValMap_[id]->setBinLabel(EVTS+1, "events", 1);
101  zeroSuppValMap_[id]->setBinLabel(EVTSGOOD+1, "good events", 1);
102  zeroSuppValMap_[id]->setBinLabel(EVTSBAD+1, "bad events", 1);
103  zeroSuppValMap_[id]->setBinLabel(BLOCKS+1, "blocks", 1);
104  zeroSuppValMap_[id]->setBinLabel(ZSBLKSGOOD+1, "good blocks", 1);
105  zeroSuppValMap_[id]->setBinLabel(ZSBLKSBAD+1, "bad blocks", 1);
106  zeroSuppValMap_[id]->setBinLabel(ZSBLKSBADFALSEPOS+1, "false pos.", 1);
107  zeroSuppValMap_[id]->setBinLabel(ZSBLKSBADFALSENEG+1, "false neg.", 1);
108  zeroSuppValMap_[id]->setBinLabel(BXBLOCKS+1, "BX blocks", 1);
109  zeroSuppValMap_[id]->setBinLabel(ZSBXBLKSGOOD+1, "good BX blocks", 1);
110  zeroSuppValMap_[id]->setBinLabel(ZSBXBLKSBAD+1, "bad BX blocks", 1);
111  zeroSuppValMap_[id]->setBinLabel(ZSBXBLKSBADFALSEPOS+1, "BX false pos.", 1);
112  zeroSuppValMap_[id]->setBinLabel(ZSBXBLKSBADFALSENEG+1, "BX false neg.", 1);
113 
114  errorSummaryNumMap_[id] = ibooker.book1D("errorSummaryNum", summaryTitleText, RNBINLABELS, 0, RNBINLABELS);
115  errorSummaryNumMap_[id]->setBinLabel(REVTS+1, "bad events", 1);
116  errorSummaryNumMap_[id]->setBinLabel(RBLKS+1, "bad blocks", 1);
117  errorSummaryNumMap_[id]->setBinLabel(RBLKSFALSEPOS+1, "false pos.", 1);
118  errorSummaryNumMap_[id]->setBinLabel(RBLKSFALSENEG+1, "false neg.", 1);
119  errorSummaryNumMap_[id]->setBinLabel(RBXBLKS+1, "bad BX blocks", 1);
120  errorSummaryNumMap_[id]->setBinLabel(RBXBLKSFALSEPOS+1, "BX false pos.", 1);
121  errorSummaryNumMap_[id]->setBinLabel(RBXBLKSFALSENEG+1, "BX false neg.", 1);
122 
123  errorSummaryDenMap_[id] = ibooker.book1D("errorSummaryDen", "denominators", RNBINLABELS, 0, RNBINLABELS);
124  errorSummaryDenMap_[id]->setBinLabel(REVTS+1, "# events", 1);
125  errorSummaryDenMap_[id]->setBinLabel(RBLKS+1, "# blocks", 1);
126  errorSummaryDenMap_[id]->setBinLabel(RBLKSFALSEPOS+1, "# blocks", 1);
127  errorSummaryDenMap_[id]->setBinLabel(RBLKSFALSENEG+1, "# blocks", 1);
128  errorSummaryDenMap_[id]->setBinLabel(RBXBLKS+1, "# BX blocks", 1);
129  errorSummaryDenMap_[id]->setBinLabel(RBXBLKSFALSEPOS+1, "# BX blocks", 1);
130  errorSummaryDenMap_[id]->setBinLabel(RBXBLKSFALSENEG+1, "# BX blocks", 1);
131  // Setting canExtend to false is needed to get the correct behaviour when running multithreaded.
132  // Otherwise, when merging the histgrams of the threads, TH1::Merge sums bins that have the same label in one bin.
133  // This needs to come after the calls to setBinLabel.
134  errorSummaryDenMap_[id]->getTH1F()->GetXaxis()->SetCanExtend(false);
135 
136  readoutSizeNoZSMap_[id] = ibooker.book1D("readoutSize", sizeTitleText + "size", 100, 0, maxFedReadoutSize_);
137  readoutSizeNoZSMap_[id]->setAxisTitle("size (byte)", 1);
138  readoutSizeZSMap_[id] = ibooker.book1D("readoutSizeZS", sizeTitleText + "size with zero suppression", 100, 0, maxFedReadoutSize_);
139  readoutSizeZSMap_[id]->setAxisTitle("size (byte)", 1);
140  readoutSizeZSExpectedMap_[id] = ibooker.book1D("readoutSizeZSExpected", "Expected " + sizeTitleText + "size with zero suppression", 100, 0, maxFedReadoutSize_);
141  readoutSizeZSExpectedMap_[id]->setAxisTitle("size (byte)", 1);
142 }
std::map< unsigned int, MonitorElement * > readoutSizeZSExpectedMap_
std::map< unsigned int, MonitorElement * > zeroSuppValMap_
std::map< unsigned int, MonitorElement * > readoutSizeNoZSMap_
std::map< unsigned int, MonitorElement * > errorSummaryNumMap_
std::map< unsigned int, MonitorElement * > errorSummaryDenMap_
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:118
static const unsigned int maxMasks_
std::map< unsigned int, MonitorElement * > readoutSizeZSMap_
void L1TMP7ZeroSupp::bookHistograms ( DQMStore::IBooker ibooker,
const edm::Run ,
const edm::EventSetup  
)
overrideprotected

Definition at line 74 of file L1TMP7ZeroSupp.cc.

References DQMStore::IBooker::book1D(), bookCapIdHistograms(), capIds_, definedMaskCapIds_, maxMasks_, monitorDir_, MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

74  {
75  // overall summary
78  capIds_ = ibooker.book1D("capIds", "Caption ids found in data", maxMasks_, 0, maxMasks_);
79  capIds_->setAxisTitle("caption id", 1);
80 
81  // per caption id subdirectories
82  for (const auto &id: definedMaskCapIds_) {
83  ibooker.setCurrentFolder(monitorDir_+"/CapId"+std::to_string(id));
84  bookCapIdHistograms(ibooker, id);
85  }
86 }
std::string monitorDir_
void bookCapIdHistograms(DQMStore::IBooker &ibooker, const unsigned int &id)
std::vector< unsigned int > definedMaskCapIds_
MonitorElement * capIds_
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:118
static const unsigned int maxMasks_
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:279
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
void L1TMP7ZeroSupp::dqmBeginRun ( const edm::Run r,
const edm::EventSetup c 
)
overrideprotected

Definition at line 70 of file L1TMP7ZeroSupp.cc.

70 {}
void L1TMP7ZeroSupp::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 47 of file L1TMP7ZeroSupp.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addOptionalUntracked(), edm::ParameterSetDescription::addUntracked(), mps_fire::i, maxMasks_, and AlCaHLTBitMon_QueryRunRegistry::string.

47  {
49  desc.add<edm::InputTag>("rawData");
50  desc.add<std::vector<int>>("fedIds")->setComment("FED ids to analyze.");
51  desc.addUntracked<bool>("zsEnabled", true)->setComment("MP7 zero suppression is enabled.");
52  desc.addUntracked<int>("lenSlinkHeader", 8)->setComment("Number of Slink header bytes.");
53  desc.addUntracked<int>("lenSlinkTrailer", 8)->setComment("Number of Slink trailer bytes.");
54  desc.addUntracked<int>("lenAMC13Header", 8)->setComment("Number of AMC13 header bytes.");
55  desc.addUntracked<int>("lenAMC13Trailer", 8)->setComment("Number of AMC13 trailer bytes.");
56  desc.addUntracked<int>("lenAMCHeader", 8)->setComment("Number of AMC header bytes.");
57  desc.addUntracked<int>("lenAMCTrailer", 0)->setComment("Number of AMC trailer bytes.");
58  desc.addUntracked<int>("zsFlagMask", 0x1)->setComment("Zero suppression flag mask.");
59  desc.addUntracked<int>("newZsFlagMask", 0x2)->setComment("Per-BX zero suppression flag mask.");
60  desc.addUntracked<int>("dataInvFlagMask", 0x1)->setComment("Data inversion flag mask.");
61  desc.addUntracked<int>("maxFEDReadoutSize", 10000)->setComment("Maximal FED readout size histogram x-axis value.");
62  for (unsigned int i = 0; i < maxMasks_; ++i) {
63  desc.addOptionalUntracked<std::vector<int>>("maskCapId"+std::to_string(i))->setComment("ZS mask for caption id "+std::to_string(i)+".");
64  }
65  desc.addUntracked<std::string>("monitorDir", "")->setComment("Target directory in the DQM file. Will be created if not existing.");
66  desc.addUntracked<bool>("verbose", false);
67  descriptions.add("l1tMP7ZeroSupp", desc);
68 }
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
static const unsigned int maxMasks_
ParameterDescriptionBase * add(U const &iLabel, T const &value)
void add(std::string const &label, ParameterSetDescription const &psetDescription)
ParameterDescriptionBase * addOptionalUntracked(U const &iLabel, T const &value)

Member Data Documentation

int L1TMP7ZeroSupp::amc13HeaderSize_
private

Definition at line 53 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

int L1TMP7ZeroSupp::amc13TrailerSize_
private

Definition at line 54 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

int L1TMP7ZeroSupp::amcHeaderSize_
private

Definition at line 55 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

int L1TMP7ZeroSupp::amcTrailerSize_
private

Definition at line 56 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

MonitorElement* L1TMP7ZeroSupp::capIds_
private

Definition at line 76 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and bookHistograms().

int L1TMP7ZeroSupp::dataInvFlagMask_
private

Definition at line 59 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

std::vector<unsigned int> L1TMP7ZeroSupp::definedMaskCapIds_
private

Definition at line 68 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), bookHistograms(), and L1TMP7ZeroSupp().

std::map<unsigned int, MonitorElement*> L1TMP7ZeroSupp::errorSummaryDenMap_
private

Definition at line 72 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and bookCapIdHistograms().

std::map<unsigned int, MonitorElement*> L1TMP7ZeroSupp::errorSummaryNumMap_
private

Definition at line 71 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and bookCapIdHistograms().

edm::EDGetTokenT<FEDRawDataCollection> L1TMP7ZeroSupp::fedDataToken_
private

Definition at line 45 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

std::vector<int> L1TMP7ZeroSupp::fedIds_
private

Definition at line 47 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

std::vector<std::vector<int> > L1TMP7ZeroSupp::masks_
private

Definition at line 48 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and L1TMP7ZeroSupp().

int L1TMP7ZeroSupp::maxFedReadoutSize_
private

Definition at line 61 of file L1TMP7ZeroSupp.h.

Referenced by bookCapIdHistograms().

const unsigned int L1TMP7ZeroSupp::maxMasks_ = 16
staticprivate
std::string L1TMP7ZeroSupp::monitorDir_
private

Definition at line 63 of file L1TMP7ZeroSupp.h.

Referenced by bookHistograms().

int L1TMP7ZeroSupp::newZsFlagMask_
private

Definition at line 57 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

std::map<unsigned int, MonitorElement*> L1TMP7ZeroSupp::readoutSizeNoZSMap_
private

Definition at line 73 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and bookCapIdHistograms().

std::map<unsigned int, MonitorElement*> L1TMP7ZeroSupp::readoutSizeZSExpectedMap_
private

Definition at line 75 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and bookCapIdHistograms().

std::map<unsigned int, MonitorElement*> L1TMP7ZeroSupp::readoutSizeZSMap_
private

Definition at line 74 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and bookCapIdHistograms().

int L1TMP7ZeroSupp::slinkHeaderSize_
private

Definition at line 51 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

int L1TMP7ZeroSupp::slinkTrailerSize_
private

Definition at line 52 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

bool L1TMP7ZeroSupp::verbose_
private

Definition at line 64 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and L1TMP7ZeroSupp().

std::map<unsigned int, MonitorElement*> L1TMP7ZeroSupp::zeroSuppValMap_
private

Definition at line 70 of file L1TMP7ZeroSupp.h.

Referenced by analyze(), and bookCapIdHistograms().

bool L1TMP7ZeroSupp::zsEnabled_
private

Definition at line 46 of file L1TMP7ZeroSupp.h.

Referenced by analyze().

int L1TMP7ZeroSupp::zsFlagMask_
private

Definition at line 58 of file L1TMP7ZeroSupp.h.

Referenced by analyze().