L1MuDTSEU.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: L1MuDTSEU
//
//   Description: Single Extrapolation Unit
//
//
//
//   Author :
//   N. Neumeister            CERN EP
//
//--------------------------------------------------

//-----------------------
// This Class's Header --
//-----------------------

#include "L1Trigger/DTTrackFinder/src/L1MuDTSEU.h"

//---------------
// C++ Headers --
//---------------

#include <iostream>
#include <algorithm>

//-------------------------------
// Collaborating Class Headers --
//-------------------------------

#include "L1Trigger/DTTrackFinder/src/L1MuDTTFConfig.h"
#include "L1Trigger/DTTrackFinder/src/L1MuDTSectorProcessor.h"
#include "L1Trigger/DTTrackFinder/src/L1MuDTDataBuffer.h"
#include "L1Trigger/DTTrackFinder/src/L1MuDTTrackSegLoc.h"
#include "L1Trigger/DTTrackFinder/src/L1MuDTTrackSegPhi.h"
#include "L1Trigger/DTTrackFinder/src/L1MuDTExtrapolationUnit.h"
#include "L1Trigger/DTTrackFinder/src/L1MuDTEUX.h"
#include "L1Trigger/DTTrackFinder/src/L1MuDTERS.h"

using namespace std;

// --------------------------------
//       class L1MuDTSEU
//---------------------------------

//----------------
// Constructors --
//----------------

L1MuDTSEU::L1MuDTSEU(const L1MuDTSectorProcessor& sp, Extrapolation ext, unsigned int tsId) : 
     m_sp(sp), m_ext(ext), 
     m_startTS_Id(tsId), m_startTS(nullptr), m_EUXs(), m_ERS() {
 
  m_EUXs.reserve(12);

  for ( int target_ts = 0; target_ts < 12; target_ts++ ) {
    m_EUXs.push_back( new L1MuDTEUX(m_sp,*this,target_ts) );
  }

  m_ERS = new L1MuDTERS(*this);

}


//--------------
// Destructor --
//--------------

L1MuDTSEU::~L1MuDTSEU() {

  vector<L1MuDTEUX*>::iterator iter_eux;
  for ( iter_eux = m_EUXs.begin(); iter_eux != m_EUXs.end(); iter_eux++ ) {
    delete (*iter_eux);
    *iter_eux = nullptr;
  }

  m_startTS = nullptr;
  m_EUXs.clear();

  delete m_ERS;

}


//--------------
// Operations --
//--------------

//
// run SEU
//
void L1MuDTSEU::run(const edm::EventSetup& c) {

  if ( L1MuDTTFConfig::Debug(3) ) cout << "Run SEU " << m_ext << " " 
                                       << m_startTS_Id << endl;

  pair<int,int> ext_pair = L1MuDTExtrapolationUnit::which_ext(m_ext);
  int target = ext_pair.second;

  // check if it is a nextWheel or ownWheel SEU
  bool nextWheel = isNextWheelSEU();

  // relative addresses used
  //                                         nextWheel
  //           extrapolation               extrapolation
  //              address                    address
  //        +--------+--------+         +--------+--------+
  //    +   |        |        |     +   |        |        |
  //        |  4   5 |  6   7 |         |        |  6   7 |
  //    |   |        |        |     |   |        |        |
  //    |   +--------+--------+     |   +--------+--------+
  //        |........|        |         |........|        |
  //   phi  |..0...1.|  2   3 |    phi  |........|  2   3 |
  //        |........|        |         |........|        |
  //    |   +--------+--------+     |   +--------+--------+
  //    |   |        |        |     |   |        |        |
  //        |  8   9 | 10  11 |         |        | 10  11 |
  //    -   |        |        |     -   |        |        |
  //        +--------+--------+         +--------+--------+
  //
  //             -- eta --                   -- eta --  

  // loop over all 12 target addresses
  for ( int reladr = 0; reladr < 12; reladr++ ) {

    // for the nextWheel extrapolations only reladr: 2,3,6,7,10,11 
    if ( nextWheel && (reladr/2)%2 == 0 ) continue;

    const L1MuDTTrackSegPhi* target_ts = m_sp.data()->getTSphi(target, reladr);
    if ( target_ts && !target_ts->empty() ) {
      m_EUXs[reladr]->load(m_startTS, target_ts);
      m_EUXs[reladr]->run(c);
      if ( m_EUXs[reladr]->result() ) m_EXtable.set(reladr);
    }

  }

  if ( L1MuDTTFConfig::Debug(3) ) {
    int n_ext = numberOfExt();
    if ( n_ext > 0 ) cout << "number of successful EUX : " <<  n_ext << endl;
  }

  if ( m_ERS ) m_ERS->run();

  //  if ( m_ERS->address(0) != 15 ) m_QStable.set(m_ERS->address(0));
  //  if ( m_ERS->address(1) != 15 ) m_QStable.set(m_ERS->address(1));
  m_QStable = m_EXtable;

}


//
// reset SEU
//
void L1MuDTSEU::reset() {

  m_startTS = nullptr;
  vector<L1MuDTEUX*>::iterator iter_eux;
  for ( iter_eux = m_EUXs.begin(); iter_eux != m_EUXs.end(); iter_eux++ ) {
    (*iter_eux)->reset();
  }

  m_ERS->reset();
  
  m_EXtable.reset();
  m_QStable.reset();
  
}


//
// reset a single extrapolation
//
void L1MuDTSEU::reset(unsigned int relAdr) {

  m_EXtable.reset(relAdr);
  m_QStable.reset(relAdr);
  m_EUXs[relAdr]->reset();
//  m_ERS->reset();
  
}


//
// get number of successful extrapolations
//
int L1MuDTSEU::numberOfExt() const {

  int number = 0;
  vector<L1MuDTEUX*>::const_iterator iter_eux;
  for ( iter_eux = m_EUXs.begin(); iter_eux != m_EUXs.end(); iter_eux++ ) {
    if ( (*iter_eux)->result() ) number++;
  }

  return number;

}