CSCCathodeLCTProcessor.h

Go to the documentation of this file.

#ifndef L1Trigger_CSCTriggerPrimitives_CSCCathodeLCTProcessor_h
#define L1Trigger_CSCTriggerPrimitives_CSCCathodeLCTProcessor_h

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/CSCDigi/interface/CSCComparatorDigiCollection.h"
#include "DataFormats/CSCDigi/interface/CSCCLCTDigi.h"
#include "DataFormats/CSCDigi/interface/CSCCLCTPreTriggerDigi.h"
#include "CondFormats/CSCObjects/interface/CSCDBL1TPParameters.h"
#include "L1Trigger/CSCCommonTrigger/interface/CSCConstants.h"

#include <vector>

class CSCGeometry;

class CSCCathodeLCTProcessor
{
 public:
  CSCCathodeLCTProcessor(unsigned endcap, unsigned station, unsigned sector,
             unsigned subsector, unsigned chamber,
             const edm::ParameterSet& conf,
             const edm::ParameterSet& comm,
             const edm::ParameterSet& ctmb);

  CSCCathodeLCTProcessor();

  void setConfigParameters(const CSCDBL1TPParameters* conf);

  void clear();

  void setCSCGeometry(const CSCGeometry *g) { csc_g = g; }

  std::vector<CSCCLCTDigi> run(const CSCComparatorDigiCollection* compdc);

  void run(const std::vector<int> halfstrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
       const std::vector<int> distrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);

  bool getDigis(const CSCComparatorDigiCollection* compdc);
  void getDigis(const CSCComparatorDigiCollection* compdc, const CSCDetId& id);

  CSCCLCTDigi bestCLCT[CSCConstants::MAX_CLCT_TBINS];

  CSCCLCTDigi secondCLCT[CSCConstants::MAX_CLCT_TBINS];

  std::vector<CSCCLCTDigi> readoutCLCTs();

  std::vector<CSCCLCTDigi> getCLCTs();

  std::vector<int> preTriggerBXs() const {return thePreTriggerBXs;}

  std::vector<CSCCLCTPreTriggerDigi> preTriggerDigis() const {return thePreTriggerDigis; }

  static void distripStagger(int stag_triad[CSCConstants::MAX_NUM_STRIPS_7CFEBS],
                 int stag_time[CSCConstants::MAX_NUM_STRIPS_7CFEBS],
                 int stag_digi[CSCConstants::MAX_NUM_STRIPS_7CFEBS],
                 int i_distrip, bool debug = false);

  void setRing(unsigned r) {theRing = r;}

  static const int pre_hit_pattern[2][CSCConstants::MAX_STRIPS_IN_PATTERN];
  static const int pattern[CSCConstants::NUM_CLCT_PATTERNS_PRE_TMB07][CSCConstants::MAX_STRIPS_IN_PATTERN+1];

  static const int pattern2007_offset[CSCConstants::MAX_HALFSTRIPS_IN_PATTERN];
  static const int pattern2007[CSCConstants::NUM_CLCT_PATTERNS][CSCConstants::MAX_HALFSTRIPS_IN_PATTERN+2];

  // we use these next ones to address the various bits inside the array that's
  // used to make the cathode LCTs.
  enum CLCT_INDICES {CLCT_PATTERN,
                     CLCT_BEND,
                     CLCT_STRIP,
                     CLCT_BX,
                     CLCT_STRIP_TYPE,
                     CLCT_QUALITY,
                     CLCT_CFEB,
                     CLCT_NUM_QUANTITIES= 7};

 private:
  int infoV;

  const unsigned theEndcap;
  const unsigned theStation;
  const unsigned theSector;
  const unsigned theSubsector;
  const unsigned theTrigChamber;

  const CSCGeometry* csc_g;

  // holders for easy access:
  unsigned int theRing;
  unsigned int theChamber;
  bool isME11;

  int numStrips;
  int stagger[CSCConstants::NUM_LAYERS];

  std::vector<CSCComparatorDigi> digiV[CSCConstants::NUM_LAYERS];
  std::vector<int> thePreTriggerBXs;
  std::vector<CSCCLCTPreTriggerDigi> thePreTriggerDigis;

  bool isMTCC;

  bool isTMB07;

  bool isSLHC;

  unsigned int alctClctOffset;

  unsigned int fifo_tbins,  fifo_pretrig; // only for test beam mode.
  unsigned int hit_persist, drift_delay;
  unsigned int nplanes_hit_pretrig, nplanes_hit_pattern;
  unsigned int pid_thresh_pretrig,  min_separation;
  unsigned int tmb_l1a_window_size;

  int start_bx_shift;

  bool smartME1aME1b, disableME1a, gangedME1a;

  int early_tbins;

  bool use_dead_time_zoning;
  unsigned int clct_state_machine_zone; // +- around a keystrip
  bool dynamic_state_machine_zone;  //use a pattern dependent zone

  unsigned int pretrig_trig_zone;

  bool use_corrected_bx;

  bool readout_earliest_2;

  static const unsigned int def_fifo_tbins,  def_fifo_pretrig;
  static const unsigned int def_hit_persist, def_drift_delay;
  static const unsigned int def_nplanes_hit_pretrig;
  static const unsigned int def_nplanes_hit_pattern;
  static const unsigned int def_pid_thresh_pretrig, def_min_separation;
  static const unsigned int def_tmb_l1a_window_size;

  void setDefaultConfigParameters();

  void checkConfigParameters();

  static const int cfeb_strips[2];

  //---------------- Methods common to all firmware versions ------------------
  void readComparatorDigis(std::vector<int>halfstrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
               std::vector<int> distrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);
  void readComparatorDigis(std::vector<int> halfstrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);
  void pulseExtension(
 const std::vector<int> time[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
 const int nStrips,
 unsigned int pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);

  //------------- Functions for idealized version for MC studies --------------
  std::vector<CSCCLCTDigi> findLCTs(
     const std::vector<int> strip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
     int stripType);
  bool preTrigger(
     const std::vector<int> strip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
     const int stripType, const int nStrips, int& first_bx);
  void getKeyStripData(
     const std::vector<int> strip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
     int keystrip_data[CSCConstants::NUM_HALF_STRIPS_7CFEBS][7],
     int nStrips, int first_bx, int& best_strip, int stripType);
  void getPattern(int pattern_num, int strip_value[CSCConstants::MAX_STRIPS_IN_PATTERN],
          int bx_time, int &quality, int &bend);
  bool hitIsGood(int hitTime, int BX);

  //-------------------- Functions for pre-2007 firmware ----------------------
  std::vector<CSCCLCTDigi> findLCTs(
  const std::vector<int> halfstrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
  const std::vector<int> distrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);
  bool preTrigger(
      const std::vector<int> strip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
   unsigned int pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
      const int stripType, const int nStrips,
      const int start_bx, int& first_bx);
  bool preTrigLookUp(const unsigned int pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
             const int stripType, const int nStrips,
             const unsigned int bx_time);
  void latchLCTs(const unsigned int pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
         int keyStrip[CSCConstants::MAX_CFEBS], unsigned int nhits[CSCConstants::MAX_CFEBS],
         const int stripType, const int nStrips, const int bx_time);
  void priorityEncode(const int h_keyStrip[CSCConstants::MAX_CFEBS],
              const unsigned int h_nhits[CSCConstants::MAX_CFEBS],
              const int d_keyStrip[CSCConstants::MAX_CFEBS],
              const unsigned int d_nhits[CSCConstants::MAX_CFEBS],
              int keystrip_data[2][7]);
  void getKeyStripData(const unsigned int h_pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
               const unsigned int d_pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
               int keystrip_data[2][7], const int first_bx);
  void getPattern(unsigned int pattern_num,
          const int strip_value[CSCConstants::MAX_STRIPS_IN_PATTERN],
          unsigned int& quality, unsigned int& bend);

  //--------------- Functions for 2007 version of the firmware ----------------
  std::vector<CSCCLCTDigi> findLCTs(
 const std::vector<int> halfstrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);
  bool preTrigger(
      const unsigned int pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
      const int start_bx, int& first_bx);
  bool ptnFinding(
      const unsigned int pulse[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
      const int nStrips, const unsigned int bx_time);
  void markBusyKeys(const int best_hstrip, const int best_patid,
            int quality[CSCConstants::NUM_HALF_STRIPS_7CFEBS]);

  unsigned int best_pid[CSCConstants::NUM_HALF_STRIPS_7CFEBS];
  unsigned int nhits[CSCConstants::NUM_HALF_STRIPS_7CFEBS];
  int first_bx_corrected[CSCConstants::NUM_HALF_STRIPS_7CFEBS];

  //--------------- Functions for SLHC studies ----------------

  std::vector<CSCCLCTDigi> findLCTsSLHC(
    const std::vector<int>  halfstrip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);

  bool ispretrig[CSCConstants::NUM_HALF_STRIPS_7CFEBS];

  //--------------------------- Auxiliary methods -----------------------------
  void dumpConfigParams() const;

  void dumpDigis(
      const std::vector<int> strip[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS],
      const int stripType, const int nStrips) const;

  //--------------------------- Methods for tests -----------------------------
  void testDistripStagger();
  void testLCTs();
  void printPatterns();
  void testPatterns();
  int findNumLayersHit(std::vector<int> stripsHit[CSCConstants::NUM_LAYERS][CSCConstants::NUM_HALF_STRIPS_7CFEBS]);
};

#endif