MuonTriggerPrimitive.cc

Go to the documentation of this file.

#include "L1Trigger/L1TMuon/interface/deprecate/MuonTriggerPrimitive.h"

// the primitive types we can use
#include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigi.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambPhDigi.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambThDigi.h"
#include "DataFormats/RPCDigi/interface/RPCDigiL1Link.h"

// detector ID types
#include "DataFormats/MuonDetId/interface/DTChamberId.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"

using namespace L1TMuon;

namespace {
  const char subsystem_names[][4] = {"DT","CSC","RPC"};
}

//constructors from DT data
TriggerPrimitive::TriggerPrimitive(const DTChamberId& detid,
                   const L1MuDTChambPhDigi& digi_phi,
                   const int segment_number):
  _id(detid),
  _subsystem(TriggerPrimitive::kDT) {
  calculateDTGlobalSector(detid,_globalsector,_subsector);
  // fill in information from theta trigger
  _dt.theta_bti_group = -1;
  _dt.segment_number = segment_number;
  _dt.theta_code = -1;
  _dt.theta_quality = -1;
  // now phi trigger
  _dt.bx = digi_phi.bxNum();
  _dt.wheel = digi_phi.whNum();
  _dt.sector = digi_phi.scNum();
  _dt.station = digi_phi.stNum();
  _dt.radialAngle = digi_phi.phi();
  _dt.bendingAngle = digi_phi.phiB();
  _dt.qualityCode = digi_phi.code();
  _dt.Ts2TagCode = digi_phi.Ts2Tag();
  _dt.BxCntCode = digi_phi.BxCnt();
}

TriggerPrimitive::TriggerPrimitive(const DTChamberId& detid,
                   const L1MuDTChambThDigi& digi_th,
                   const int theta_bti_group):
  _id(detid),
  _subsystem(TriggerPrimitive::kDT) {
  calculateDTGlobalSector(detid,_globalsector,_subsector);
  // fill in information from theta trigger
  _dt.theta_bti_group = theta_bti_group;
  _dt.segment_number = digi_th.position(theta_bti_group);
  _dt.theta_code = digi_th.code(theta_bti_group);
  _dt.theta_quality = digi_th.quality(theta_bti_group);
  // now phi trigger
  _dt.bx = digi_th.bxNum();
  _dt.wheel = digi_th.whNum();
  _dt.sector = digi_th.scNum();
  _dt.station = digi_th.stNum();
  _dt.radialAngle = -1;
  _dt.bendingAngle = -1;
  _dt.qualityCode = -1;
  _dt.Ts2TagCode = -1;
  _dt.BxCntCode = -1;
}

TriggerPrimitive::TriggerPrimitive(const DTChamberId& detid,
                   const L1MuDTChambPhDigi& digi_phi,
                   const L1MuDTChambThDigi& digi_th,
                   const int theta_bti_group):
  _id(detid),
  _subsystem(TriggerPrimitive::kDT) {
  calculateDTGlobalSector(detid,_globalsector,_subsector);
  // fill in information from theta trigger
  _dt.theta_bti_group = theta_bti_group;
  _dt.segment_number = digi_th.position(theta_bti_group);
  _dt.theta_code = digi_th.code(theta_bti_group);
  _dt.theta_quality = digi_th.quality(theta_bti_group);
  // now phi trigger
  _dt.bx = digi_phi.bxNum();
  _dt.wheel = digi_phi.whNum();
  _dt.sector = digi_phi.scNum();
  _dt.station = digi_phi.stNum();
  _dt.radialAngle = digi_phi.phi();
  _dt.bendingAngle = digi_phi.phiB();
  _dt.qualityCode = digi_phi.code();
  _dt.Ts2TagCode = digi_phi.Ts2Tag();
  _dt.BxCntCode = digi_phi.BxCnt();
}

//constructor from CSC data
TriggerPrimitive::TriggerPrimitive(const CSCDetId& detid,
                   const CSCCorrelatedLCTDigi& digi):
  _id(detid),
  _subsystem(TriggerPrimitive::kCSC) {
  calculateCSCGlobalSector(detid,_globalsector,_subsector);
  _csc.trknmb  = digi.getTrknmb();
  _csc.valid   = digi.isValid();
  _csc.quality = digi.getQuality();
  _csc.keywire = digi.getKeyWG();
  _csc.strip   = digi.getStrip();
  _csc.pattern = digi.getPattern();
  _csc.bend    = digi.getBend();
  _csc.bx      = digi.getBX();
  _csc.mpclink = digi.getMPCLink();
  _csc.bx0     = digi.getBX0();
  _csc.syncErr = digi.getSyncErr();
  _csc.cscID   = digi.getCSCID();
}

// Create copy of TP1 with wire from TP2
TriggerPrimitive::TriggerPrimitive(const TriggerPrimitive& tp1, const TriggerPrimitive& tp2) {

  _subsystem    = TriggerPrimitive::kCSC;
  _id       = tp1.detId<CSCDetId>();
  _csc.trknmb   = tp1._csc.trknmb; 
  _csc.valid    = tp1._csc.valid; 
  _csc.quality  = tp1._csc.quality; 
  _csc.keywire  = tp2._csc.keywire; 
  _csc.strip    = tp1._csc.strip; 
  _csc.pattern  = tp1._csc.pattern; 
  _csc.bend     = tp1._csc.bend; 
  _csc.bx   = tp1._csc.bx; 
  _csc.mpclink  = tp1._csc.mpclink; 
  _csc.bx0  = tp1._csc.bx0; 
  _csc.syncErr  = tp1._csc.syncErr;
  _csc.cscID    = tp1._csc.cscID;

}


// constructor from RPC data
TriggerPrimitive::TriggerPrimitive(const RPCDetId& detid,
                   const unsigned strip,
                   const unsigned layer,
                   const uint16_t bx):
  _id(detid),
  _subsystem(TriggerPrimitive::kRPC) {
  calculateRPCGlobalSector(detid,_globalsector,_subsector);
  _rpc.strip = strip;
  _rpc.layer = layer;
  _rpc.bx = bx;
}
/*
TriggerPrimitive::TriggerPrimitive(const TriggerPrimitive& tp):
  _dt(tp._dt),
  _csc(tp._csc),
  _rpc(tp._rpc),
  _id(tp._id),
  _subsystem(tp._subsystem),
  _globalsector(tp._globalsector),
  _subsector(tp._subsector),
  _eta(tp._eta),
  _phi(tp._phi),
  _theta(tp._theta){
}*/

TriggerPrimitive& TriggerPrimitive::operator=(const TriggerPrimitive& tp) {
  this->_dt = tp._dt;
  this->_csc = tp._csc;
  this->_rpc = tp._rpc;
  this->_id = tp._id;
  this->_subsystem = tp._subsystem;
  this->_globalsector = tp._globalsector;
  this->_subsector = tp._subsector;
  this->_eta = tp._eta;
  this->_phi = tp._phi;
  return *this;
}


TriggerPrimitive::TriggerPrimitive(const TriggerPrimitive& tp):
  _dt(tp._dt),
  _csc(tp._csc),
  _rpc(tp._rpc),
  _id(tp._id),
  _subsystem(tp._subsystem),
  _globalsector(tp._globalsector),
  _subsector(tp._subsector),
  _eta(tp._eta),
  _phi(tp._phi),
  _rho(tp._rho),
  _theta(tp._theta){
}

bool TriggerPrimitive::operator==(const TriggerPrimitive& tp) const {
  return ( this->_dt.bx == tp._dt.bx &&
       this->_dt.wheel == tp._dt.wheel &&
       this->_dt.sector == tp._dt.sector &&
       this->_dt.station == tp._dt.station &&
       this->_dt.radialAngle == tp._dt.radialAngle &&
       this->_dt.bendingAngle == tp._dt.bendingAngle &&
       this->_dt.qualityCode == tp._dt.qualityCode &&
       this->_dt.Ts2TagCode == tp._dt.Ts2TagCode &&
       this->_dt.BxCntCode == tp._dt.BxCntCode &&
       this->_dt.theta_bti_group == tp._dt.theta_bti_group &&
       this->_dt.segment_number == tp._dt.segment_number &&
       this->_dt.theta_code == tp._dt.theta_code &&
       this->_dt.theta_quality == tp._dt.theta_quality &&
       this->_csc.trknmb == tp._csc.trknmb &&
       this->_csc.valid == tp._csc.valid &&
       this->_csc.quality == tp._csc.quality &&
       this->_csc.keywire == tp._csc.keywire &&
       this->_csc.strip == tp._csc.strip &&
       this->_csc.pattern == tp._csc.pattern &&
       this->_csc.bend == tp._csc.bend &&
       this->_csc.bx == tp._csc.bx &&
       this->_csc.mpclink == tp._csc.mpclink &&
       this->_csc.bx0 == tp._csc.bx0 &&
       this->_csc.syncErr == tp._csc.syncErr &&
       this->_csc.cscID == tp._csc.cscID &&
       this->_rpc.strip == tp._rpc.strip &&
       this->_rpc.layer == tp._rpc.layer &&
       this->_rpc.bx == tp._rpc.bx &&
       this->_id == tp._id &&
       this->_subsystem == tp._subsystem &&
       this->_globalsector == tp._globalsector &&
       this->_subsector == tp._subsector );
}

const int TriggerPrimitive::getBX() const {
  switch(_subsystem) {
  case kDT:
    return _dt.bx;
  case kCSC:
    return _csc.bx;
  case kRPC:
    return _rpc.bx;
  default:
    throw cms::Exception("Invalid Subsytem")
      << "The specified subsystem for this track stub is out of range"
      << std::endl;
  }
  return -1;
}

const int TriggerPrimitive::getStrip() const {
  switch(_subsystem) {
  case kDT:
    return -1;
  case kCSC:
    return _csc.strip;
  case kRPC:
    return _rpc.strip;
  default:
    throw cms::Exception("Invalid Subsytem")
      << "The specified subsystem for this track stub is out of range"
      << std::endl;
  }
  return -1;
}

const int TriggerPrimitive::getWire() const {
  switch(_subsystem) {
  case kDT:
    return -1;
  case kCSC:
    return _csc.keywire;
  case kRPC:
    return -1;
  default:
    throw cms::Exception("Invalid Subsytem")
      << "The specified subsystem for this track stub is out of range"
      << std::endl;
  }
  return -1;
}

const int TriggerPrimitive::getPattern() const {
  switch(_subsystem) {
  case kDT:
    return -1;
  case kCSC:
    return _csc.pattern;
  case kRPC:
    return -1;
  default:
    throw cms::Exception("Invalid Subsytem")
      << "The specified subsystem for this track stub is out of range"
      << std::endl;
  }
  return -1;
}
const int TriggerPrimitive::Id() const {
  switch(_subsystem) {
  case kDT:
    return -1;
  case kCSC:
    return _csc.cscID;
  case kRPC:
    return -1;
  default:
    throw cms::Exception("Invalid Subsytem")
      << "The specified subsystem for this track stub is out of range"
      << std::endl;
  }
  return -1;
}

void TriggerPrimitive::calculateDTGlobalSector(const DTChamberId& chid,
                           unsigned& global_sector,
                           unsigned& subsector ) {
}

void TriggerPrimitive::calculateCSCGlobalSector(const CSCDetId& chid,
                        unsigned& global_sector,
                        unsigned& subsector ) {
}

void TriggerPrimitive::calculateRPCGlobalSector(const RPCDetId& chid,
                        unsigned& global_sector,
                        unsigned& subsector ) {
}

void TriggerPrimitive::print(std::ostream& out) const {
  unsigned idx = (unsigned) _subsystem;
  out << subsystem_names[idx] << " Trigger Primitive" << std::endl;
  out << "eta: " << _eta << " phi: " << _phi
      << " bend: " << _theta << std::endl;
  switch(_subsystem) {
  case kDT:
    out << detId<DTChamberId>() << std::endl;
    out << "Local BX      : " << _dt.bx << std::endl;
    out << "Segment Nmb   : " << _dt.segment_number << std::endl;
    out << "Packed Phi    : " << _dt.radialAngle << std::endl;
    out << "Packed Bend   : " << _dt.bendingAngle << std::endl;
    out << "Quality Code  : " << _dt.qualityCode << std::endl;
    out << "Ts2Tag Code   : " << _dt.Ts2TagCode << std::endl;
    out << "BXCnt Code    : " << _dt.BxCntCode << std::endl;
    out << "Theta BTI Grp : " << _dt.theta_bti_group << std::endl;
    out << "Theta Code    : " << _dt.theta_code << std::endl;
    out << "Theta Quality : " << _dt.theta_quality << std::endl;
    break;
  case kCSC:
    out << detId<CSCDetId>() << std::endl;
    out << "Local BX      : " << _csc.bx << std::endl;
    out << "Segment Nmb   : " << _csc.trknmb << std::endl;
    out << "Segment Valid : " << _csc.valid << std::endl;
    out << "Quality Code  : " << _csc.quality << std::endl;
    out << "Key Wire Grp  : " << _csc.keywire << std::endl;
    out << "Half-Strip    : " << _csc.strip << std::endl;
    out << "CLCT Pattern  : " << _csc.pattern << std::endl;
    out << "Packed Bend   : " << _csc.bend << std::endl;
    out << "MPC Link      : " << _csc.mpclink << std::endl;
    out << "BX0           : " << _csc.bx0 << std::endl;
    out << "Sync Error    : " << _csc.syncErr << std::endl;
    out << "CSCID         : " << _csc.cscID << std::endl;
    break;
  case kRPC:
    out << detId<RPCDetId>() << std::endl;
    out << "Local BX      : " << _rpc.bx << std::endl;
    out << "Strip         : " << _rpc.strip << std::endl;
    out << "Layer         : " << _rpc.layer << std::endl;
    break;
  default:
    throw cms::Exception("Invalid Subsytem")
      << "The specified subsystem for this track stub is out of range"
      << std::endl;
  }
}