LCTQualityControl.h

Go to the documentation of this file.

#ifndef L1Trigger_CSCTriggerPrimitives_LCTQualityControl_h
#define L1Trigger_CSCTriggerPrimitives_LCTQualityControl_h

#include "DataFormats/CSCDigi/interface/CSCCLCTDigi.h"
#include "DataFormats/CSCDigi/interface/CSCALCTDigi.h"
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigi.h"
#include "L1Trigger/CSCTriggerPrimitives/interface/CSCBaseboard.h"

#include <vector>
#include <algorithm>

class LCTQualityControl : public CSCBaseboard {
public:
  // constructor
  LCTQualityControl(unsigned endcap,
                    unsigned station,
                    unsigned sector,
                    unsigned subsector,
                    unsigned chamber,
                    const edm::ParameterSet& conf);

  ~LCTQualityControl() override = default;

  // Check if the ALCT is valid
  void checkValidReadout(const CSCALCTDigi& alct) const;
  void checkValid(const CSCALCTDigi& alct, unsigned max_stubs = CSCConstants::MAX_ALCTS_PER_PROCESSOR) const;

  // Check if the CLCT is valid
  void checkValid(const CSCCLCTDigi& lct, unsigned max_stubs = CSCConstants::MAX_CLCTS_PER_PROCESSOR) const;

  // Check if the LCT is valid - TMB version
  void checkValid(const CSCCorrelatedLCTDigi& lct) const;

  // Check if the LCT is valid - MPC version
  void checkValid(const CSCCorrelatedLCTDigi& lct, const unsigned station, const unsigned ring) const;

  void checkRange(int parameter, int min_value, int max_value, const std::string& comment, unsigned& errors) const;

  template <class T>
  void reportErrors(const T& lct, const unsigned errors) const;

  // no more than 2 LCTs per BX in the readout
  void checkMultiplicityBX(const std::vector<CSCALCTDigi>& alcts) const;
  void checkMultiplicityBX(const std::vector<CSCCLCTDigi>& clcts) const;
  void checkMultiplicityBX(const std::vector<CSCCorrelatedLCTDigi>& lcts) const;
  template <class T>
  void checkMultiplicityBX(const std::vector<T>& lcts, unsigned nLCT) const;

  // for Phase-1 patterns
  int getSlopePhase1(unsigned pattern) const;

  /*
    CSC max strip & max wire
    Need to pass station and ring, because we want to reuse the LCTQualityControl in the
    MPC where the station and ring are considered arguments, not data members.
  */
  unsigned get_csc_max_wiregroup(unsigned station, unsigned ring) const;
  unsigned get_csc_max_halfstrip(unsigned station, unsigned ring) const;
  unsigned get_csc_max_quartstrip(unsigned station, unsigned ring) const;
  unsigned get_csc_max_eighthstrip(unsigned station, unsigned ring) const;

  // slope values
  std::pair<int, int> get_csc_clct_min_max_slope() const;

  // CLCT min, max CFEB numbers
  std::pair<unsigned, unsigned> get_csc_min_max_cfeb() const;

  // CSC min, max pattern
  std::pair<unsigned, unsigned> get_csc_min_max_pattern() const;
  std::pair<unsigned, unsigned> get_csc_min_max_pattern_run3() const;
  std::pair<unsigned, unsigned> get_csc_lct_min_max_pattern() const;

  /*
   CSC min, max quality
   Need to pass station and ring for the LCT implementation, because we want to
   reuse the LCTQualityControl in the MPC where the station and ring are considered
   arguments, not data members.
  */
  std::pair<unsigned, unsigned> get_csc_alct_min_max_quality() const;
  std::pair<unsigned, unsigned> get_csc_clct_min_max_quality() const;
  std::pair<unsigned, unsigned> get_csc_lct_min_max_quality(unsigned station, unsigned ring) const;

private:
  // min number of layers for a CLCT
  unsigned nplanes_clct_hit_pattern;
};

template <class T>
void LCTQualityControl::checkMultiplicityBX(const std::vector<T>& collection, unsigned nLCT) const {
  std::unordered_map<int, unsigned> freq;
  // check each BX
  for (const auto& p : collection) {
    if (p.isValid()) {
      freq[p.getBX()]++;

      // too many ALCTs, CLCTs or LCTs in this BX
      if (freq[p.getBX()] > nLCT) {
        edm::LogError("LCTQualityControl") << "Collection with more than " << nLCT << " in BX " << p.getBX();
        break;
      }
    }
  }
}

#endif