DTBtiChip.h

Go to the documentation of this file.

//-------------------------------------------------
//
//
//--------------------------------------------------
#ifndef DT_BTI_CHIP_H
#define DT_BTI_CHIP_H

//------------------------------------
// Collaborating Class Declarations --
//------------------------------------
class DTBtiHit;
class DTBtiTrig;
class DTBtiTrigData;
class DTDigi;

//----------------------
// Base Class Headers --
//----------------------
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/MuonDetId/interface/DTBtiId.h"
#include "L1Trigger/DTUtilities/interface/DTTrigGeom.h"
#include "L1TriggerConfig/DTTPGConfig/interface/DTConfig.h"
#include "L1TriggerConfig/DTTPGConfig/interface/BitArray.h"
#include "L1TriggerConfig/DTTPGConfig/interface/DTConfigBti.h"
#include "L1Trigger/DTBti/interface/DTBtiCard.h"

//---------------
// C++ Headers --
//---------------
#include <vector>
#include <memory>
#include <array>

//              ---------------------
//              -- Class Interface --
//              ---------------------

class DTBtiChip {

  public:

  //DTBtiChip(DTTrigGeom* geom, int supl, int n);

  DTBtiChip(DTBtiCard* card, DTTrigGeom* geom, int supl, int n, DTConfigBti* _config );

  DTBtiChip(DTBtiChip&& bti) = delete;
  DTBtiChip(DTBtiChip const& bti) = delete;

  ~DTBtiChip();

  DTBtiChip& operator=(DTBtiChip&& bti) = delete;
  DTBtiChip& operator=(DTBtiChip const& bti) = delete;

  void add_digi(int cell, const DTDigi* digi);

  void add_digi_clock(int cell, int clock_digi);

  std::vector<const DTDigi*> const& get_CellDigis(int cell) const { return _digis[cell];}

  void run();

  void clear();

  void addTrig(int step, std::unique_ptr<DTBtiTrig> btitrig);

  // Public const methods

  inline int number() const { return _id.bti(); }

  inline int superlayer() const { return _id.superlayer(); }

  inline LocalPoint localPosition() const { return _geom->localPosition(_id); }

  inline GlobalPoint CMSPosition() const { return _geom->CMSPosition(_id); }

  int nCellHit() const;

  int nTrig(int step) const;

  //  // Return the trigger vector
  std::vector<std::unique_ptr<DTBtiTrig>> const& trigList(int step) const;

  DTBtiTrig const* trigger(int step, unsigned n) const;

  DTBtiTrigData triggerData(int step, unsigned n) const;

  //inline DTConfig* config() const { return _geom->config(); }

  inline DTConfigBti* config() const { return _config; }


  inline DTTrigGeom* geom() const { return _geom; }

  inline DTBtiId id() const { return _id; }

  inline int wheel() const { return _geom->wheel(); }

  inline int station() const { return _geom->station(); }

  inline int sector() const { return _geom->sector(); }

  void init_clock();    // initialization from clocks



 private:

  void init(); // initialization
  void tick(); // next step (80 MHz)
  // return current step (40MHz)
  inline int currentStep() const { return (int)(((float)(_curStep)+0.5)/2); }
  inline int currentIntStep() const { return _curStep; } // int. step (80MHz)
  void computeSums();                                 // compute sums and diffs
  void sum(const int s, const int a, const int b);    //   "     a sum and dif
  void computeEqs();                                  // compute X and K equat.
  void findTrig();                                    // find triggers
  int keepTrig(const int eq, const int acp, const int code); // find  a trigger
  int keepTrigPatt(int flag, const int eq, const int pattType, int hlflag);//SV
  bool matchEq(float eqA, float eqB, int AC);
  void acceptMask(BitArray<80> * BitArrPtr,int k,int accep);   
  void doLTS();                                   // adjacent LTRIG suppression
  int store(const int eq, const int code, const int K, const int X, 
             float KeqAB=0., float KeqBC=0., float KeqCD=0., 
             float KeqAC=0., float KeqBD=0., float KeqAD=0.); 
  void eraseTrigger(int step, unsigned n); // Erase the requested trigger
  void setSnap();
  void reSumSet(); //remainder table compute
  int reSum(int a, int b) { return reSumAr[a][b+2];}
  int reSum23(int a, int b) { return reSumAr23[a][b+2];}

 private:

  // parent card
  DTBtiCard* _card;

  DTTrigGeom*  _geom;
  DTConfigBti* _config;

  DTBtiId _id;

  // input data from DTDigis
  std::array<std::vector<const DTDigi*>,9> _digis;
  // input data from clock digis
  std::array<std::vector<int >,9> _digis_clock;

  // output data (ordered by step number)
  std::array<std::vector<std::unique_ptr<DTBtiTrig>>,DTConfig::NSTEPL - DTConfig::NSTEPF+1>  _trigs;

  // internal use variables
  int _curStep;                      // current step
  std::array<std::vector<DTBtiHit*>,9> _hits;    // current hits in cells
  std::array<int,9> _thisStepUsedTimes;         // current used times in cells (JTRIG)
  std::array<DTBtiHit*,9> _thisStepUsedHit; // link to currently used hits
  int _nStepUsedHits;                // number of currently used hits
  std::array<float,25> _sums, _difs;        // time sums and differences 
  float _Keq[32][6];                 // The K equations
  float _Xeq[32][2];                 // The X equations
  int _MinKAcc;                      // min K value accepted by DTBtiChip 
  int _MaxKAcc;                      // max K value accepted by DTBtiChip
  int _MinKleftTraco;            // K limits for left traco
  int _MaxKleftTraco;
  int _MinKcenterTraco;              // K limits for center traco
  int _MaxKcenterTraco;
  int _MinKrightTraco;               // K limits for right traco
  int _MaxKrightTraco;

  float _XeqAC_patt0, _XeqBD_patt0;  // special pattern 0 X equations
  float _XeqAB_patt0, _XeqCD_patt0;

  float _KTR[32][2];                 //
  float _JTR[32][3];                 //
  int init_done;                     // initialization flag
  std::array<int,9> _busyStart_clock;            // SV - busy wire flag

  //snap register variables
  int ST43, RE43, ST23, RE23, ST, ST2, ST3, ST4, ST5, ST7;
  //remainder table
  int reSumAr[3][5];
  int reSumAr23[3][5];

};
#endif