CSCMotherboard.h

Go to the documentation of this file.

#ifndef L1Trigger_CSCTriggerPrimitives_CSCMotherboard_h
#define L1Trigger_CSCTriggerPrimitives_CSCMotherboard_h

#include "L1Trigger/CSCTriggerPrimitives/interface/CSCAnodeLCTProcessor.h"
#include "L1Trigger/CSCTriggerPrimitives/interface/CSCCathodeLCTProcessor.h"
#include "L1Trigger/CSCTriggerPrimitives/interface/LCTContainer.h"
#include "L1Trigger/CSCTriggerPrimitives/interface/CSCALCTCrossCLCT.h"
#include "L1Trigger/CSCTriggerPrimitives/interface/CSCUpgradeAnodeLCTProcessor.h"
#include "L1Trigger/CSCTriggerPrimitives/interface/CSCUpgradeCathodeLCTProcessor.h"
#include "L1Trigger/CSCTriggerPrimitives/interface/LCTQualityAssignment.h"
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigi.h"
#include "DataFormats/CSCDigi/interface/CSCShowerDigi.h"
#include "CondFormats/CSCObjects/interface/CSCL1TPLookupTableCCLUT.h"
#include "CondFormats/CSCObjects/interface/CSCL1TPLookupTableME21ILT.h"
#include "CondFormats/CSCObjects/interface/CSCL1TPLookupTableME11ILT.h"

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

  ~CSCMotherboard() override = default;

  void run(const CSCWireDigiCollection* wiredc, const CSCComparatorDigiCollection* compdc);

  /*
    Returns vector of good correlated LCTs in the read-out time window.
    LCTs in the BX window [early_tbins,...,late_tbins] are considered
    good for physics. The central LCT BX is time bin 8.
    - tmb_l1a_window_size = 7 (Run-1, Run-2) -> [5, 6, 7, 8, 9, 10, 11]
    - tmb_l1a_window_size = 5 (Run-3)        ->    [6, 7, 8, 9, 10]
    - tmb_l1a_window_size = 3 (Run-4?)       ->       [7, 8, 9]

    Note, this function does not have an exact counterpart in the
    firmware. The reason is that the DAQ of LCTs is not correctly
    simulated in CMSSW - at least the simulation of the L1-accept.
    So, this function corresponds to both the trigger path and the
    DAQ path in the firmware. In general, the function will return
    LCTs that would not be used in the OMTF or EMTF emulator,
    because they are out-of-time relative for tracking purposes. For
    instance an LCT with BX5 would be read out by the DAQ, but would
    likely not be used by the EMTF.
  */
  virtual std::vector<CSCCorrelatedLCTDigi> readoutLCTs() const;

  // LCT selection: at most 2 in each BX
  void selectLCTs();

  std::vector<CSCShowerDigi> readoutShower() const;

  void clear();

  void setConfigParameters(const CSCDBL1TPParameters* conf);
  void setESLookupTables(const CSCL1TPLookupTableCCLUT* conf);
  void setESLookupTables(const CSCL1TPLookupTableME11ILT* conf);
  void setESLookupTables(const CSCL1TPLookupTableME21ILT* conf);

  std::unique_ptr<CSCAnodeLCTProcessor> alctProc;

  std::unique_ptr<CSCCathodeLCTProcessor> clctProc;

  // VK: change to protected, to allow inheritance
protected:
  // access to lookup tables via eventsetup
  const CSCL1TPLookupTableCCLUT* lookupTableCCLUT_;
  const CSCL1TPLookupTableME11ILT* lookupTableME11ILT_;
  const CSCL1TPLookupTableME21ILT* lookupTableME21ILT_;

  /* Containers for reconstructed ALCTs and CLCTs */
  std::vector<CSCALCTDigi> alctV;
  std::vector<CSCCLCTDigi> clctV;

  LCTContainer allLCTs_;

  /* Container with sorted and selected LCTs */
  std::vector<CSCCorrelatedLCTDigi> lctV;

  CSCShowerDigi showers_[CSCConstants::MAX_LCT_TBINS];

  // helper function to return ALCT/CLCT with correct central BX
  CSCALCTDigi getBXShiftedALCT(const CSCALCTDigi&) const;
  CSCCLCTDigi getBXShiftedCLCT(const CSCCLCTDigi&) const;

  unsigned int mpc_block_me1a;
  unsigned int alct_trig_enable, clct_trig_enable, match_trig_enable;
  unsigned int match_trig_window_size, tmb_l1a_window_size;

  bool drop_used_clcts;

  int early_tbins;

  bool readout_earliest_2;

  // when set to true, ignore CLCTs found in later BX's
  bool match_earliest_clct_only_;

  // encode special bits for high-multiplicity triggers
  std::vector<unsigned> showerSource_;
  unsigned thisShowerSource_;
  unsigned minbx_readout_;
  unsigned maxbx_readout_;

  bool ignoreAlctCrossClct_;

  /*
     Preferential index array in matching window, relative to the ALCT BX.
     Where the central match BX goes first,
     then the closest early, the closest late, etc.
  */
  std::vector<int> preferred_bx_match_;

  /* sort CLCT by bx if true, otherwise sort CLCT by quality+bending */
  bool sort_clct_bx_;

  static const unsigned int def_mpc_block_me1a;
  static const unsigned int def_alct_trig_enable, def_clct_trig_enable;
  static const unsigned int def_match_trig_enable, def_match_trig_window_size;
  static const unsigned int def_tmb_l1a_window_size;

  /* quality assignment */
  std::unique_ptr<LCTQualityAssignment> qualityAssignment_;

  /* quality control */
  std::unique_ptr<LCTQualityControl> qualityControl_;

  /*
    Helper class to check if an ALCT intersects with a CLCT. Normally
    this class should not be used. It is left in the code as a potential
    improvement for ME1/1 when unphysical LCTs are not desired. This
    function is not implemented in the firmware.
  */
  std::unique_ptr<CSCALCTCrossCLCT> cscOverlap_;

  void checkConfigParameters();

  /*sort CLCT by quality+bending and if CLCTs from different BX have
    same quality+bending, then rank CLCT by timing
   */
  void sortCLCTByQualBend(int alct_bx, std::vector<unsigned>& clctBxVector);
  /*
     For valid ALCTs in the trigger time window, look for CLCTs within the
     match-time window. Valid CLCTs are matched in-time. If a match was found
     for the best ALCT and best CLCT, also the second best ALCT and second
     best CLCT are sent to a correlation function "correlateLCTs" that will
     make the best-best pair and second-second pair (if applicable).
  */
  void matchALCTCLCT();

  /*
    This function matches maximum two ALCTs with maximum two CLCTs in
    a bunch crossing. The best ALCT is considered the one with the highest
    quality in a BX. Similarly for the best CLCT. If there is just one
    ALCT and just one CLCT, the correlated LCT is made from those two
    components. If there are exactly two ALCTs and two CLCTs, the best
    LCT and second best LCT are formed from the best ALCT-CLCT combination
    and the second best ALCT-CLCT combination. In case there is missing
    information (e.g. second best ALCT, but no second best CLCT), information
    is copied over.
   */
  void correlateLCTs(const CSCALCTDigi& bestALCT,
                     const CSCALCTDigi& secondALCT,
                     const CSCCLCTDigi& bestCLCT,
                     const CSCCLCTDigi& secondCLCT,
                     CSCCorrelatedLCTDigi& bLCT,
                     CSCCorrelatedLCTDigi& sLCT,
                     int type) const;

  /*
     This method calculates all the TMB words and then passes them to the
     constructor of correlated LCTs. The LCT data members are filled with
     information from the ALCT-CLCT combination.
  */
  void constructLCTs(
      const CSCALCTDigi& aLCT, const CSCCLCTDigi& cLCT, int type, int trknmb, CSCCorrelatedLCTDigi& lct) const;

  /*
    These functions copy valid ALCT/CLCT information to invalid the ALCT/CLCT
    if present, so that we always construct the maximum number of valid LCTs
  */
  void copyValidToInValidALCT(CSCALCTDigi&, CSCALCTDigi&) const;
  void copyValidToInValidCLCT(CSCCLCTDigi&, CSCCLCTDigi&) const;

  bool doesALCTCrossCLCT(const CSCALCTDigi&, const CSCCLCTDigi&) const;

  // CLCT pattern number: encodes the pattern number itself
  unsigned int encodePattern(const int clctPattern) const;

  void dumpConfigParams() const;

  /* match cathode shower and anode shower with and/or logic */
  void matchShowers(CSCShowerDigi* anode_showers, CSCShowerDigi* cathode_showers, bool andlogic);

  /* encode high multiplicity bits for Run-3 exotic triggers */
  void encodeHighMultiplicityBits();
};
#endif