/data/refman/pasoursint/CMSSW_5_3_3/src/DataFormats/L1CSCTrackFinder/src/TrackStub.cc

Go to the documentation of this file.

#include <DataFormats/L1CSCTrackFinder/interface/TrackStub.h>
#include <DataFormats/L1CSCTrackFinder/interface/CSCBitWidths.h>
#include <DataFormats/MuonDetId/interface/CSCTriggerNumbering.h>
#include <DataFormats/MuonDetId/interface/DTChamberId.h>
#include <DataFormats/MuonDetId/interface/CSCDetId.h>
#include <DataFormats/MuonDetId/interface/MuonSubdetId.h>

namespace csctf
{

  const double TrackStub::thePhiBinning = CSCTFConstants::SECTOR_RAD/(1<<CSCBitWidths::kGlobalPhiDataBitWidth);
  const double TrackStub::theEtaBinning = (CSCTFConstants::maxEta - CSCTFConstants::minEta)/(CSCTFConstants::etaBins);

  TrackStub::TrackStub(const CSCCorrelatedLCTDigi& aDigi,
                       const DetId& aDetId): CSCCorrelatedLCTDigi(aDigi),
                                             theDetId_(aDetId.rawId()),
                                             thePhi_(0),
                                             theEta_(0),
                                             link_(0)
  {}

  TrackStub::TrackStub(const CSCCorrelatedLCTDigi& aDigi,
                       const DetId& aDetId,
                       const unsigned& phi, const unsigned& eta): CSCCorrelatedLCTDigi(aDigi),
                                                                  theDetId_(aDetId.rawId()),
                                                                  thePhi_(phi),
                                                                  theEta_(eta),
                                                                  link_(0)
  {}



  TrackStub::TrackStub(const TrackStub& aTrackStub): CSCCorrelatedLCTDigi(aTrackStub),
                                                     theDetId_(aTrackStub.theDetId_),
                                                     thePhi_(aTrackStub.thePhi_),
                                                     theEta_(aTrackStub.theEta_),
                                                     link_(aTrackStub.link_)
  {}

  unsigned TrackStub::endcap() const
  {
    int e = 0;

    switch(DetId(theDetId_).subdetId())
      {
      case (MuonSubdetId::DT):
        e = (DTChamberId(theDetId_).wheel() > 0) ? 1 : 2;
        break;
      case (MuonSubdetId::CSC):
        e = CSCDetId(theDetId_).endcap();
        break;
      default:
        break;
      }

    return e;
  }

  unsigned TrackStub::station() const
  {
    int s = 0;

    switch(DetId(theDetId_).subdetId())
      {
      case (MuonSubdetId::DT):
        s = DTChamberId(theDetId_).station() + 4;
        break;
      case (MuonSubdetId::CSC):
        s = CSCDetId(theDetId_).station();
        break;
      default:
        break;
      }

    return s;
  }

  unsigned TrackStub::sector() const
  {
    int se = 0, temps = 0;

    switch(DetId(theDetId_).subdetId())
      {
      case (MuonSubdetId::DT): {
        temps = DTChamberId(theDetId_).sector();
        const unsigned int dt2csc[12] = {6,1,1,2,2,3,3,4,4,5,5,6};
        se = dt2csc[temps-1]; }
        break;
      case (MuonSubdetId::CSC):
        se = CSCTriggerNumbering::triggerSectorFromLabels(CSCDetId(theDetId_));
        break;
      default:
        break;
      }

    return se;
  }

  unsigned TrackStub::subsector() const
  {
    int ss = 0;

    switch(DetId(theDetId_).subdetId())
      {
      case (MuonSubdetId::DT): {
        ss = DTChamberId(theDetId_).sector();
        const unsigned int dt2csc_[12] = {2,1,2,1,2,1,2,1,2,1,2,1};
        ss = dt2csc_[ss-1]; }
        break;
      case (MuonSubdetId::CSC):
        ss = CSCTriggerNumbering::triggerSubSectorFromLabels(CSCDetId(theDetId_));
        break;
      default:
        break;
      }

    return ss;
  }

  unsigned TrackStub::cscid() const
  {
    if(DetId(theDetId_).subdetId() == MuonSubdetId::CSC)
      return CSCTriggerNumbering::triggerCscIdFromLabels(CSCDetId(theDetId_));

    return 0; // DT chambers obviously don't have a csc id :-D
  }

  unsigned TrackStub::cscidSeparateME1a() const
  {
    if(DetId(theDetId_).subdetId() != MuonSubdetId::CSC) return 0;
    CSCDetId id(theDetId_);
    unsigned normal_cscid = CSCTriggerNumbering::triggerCscIdFromLabels(id);
    if (id.station()==1 && id.ring()==4) return normal_cscid + 9; // 10,11,12 for ME1a
    return normal_cscid;
  }

  bool TrackStub::operator<(const TrackStub& rhs) const
  {
    return ( rhs.isValid() && ( (!(isValid())) || (getQuality() < rhs.getQuality()) ||
                                (getQuality() == rhs.getQuality() && cscid() < rhs.cscid()) ||
                              (getQuality() == rhs.getQuality() && cscid() == rhs.cscid() &&
                               (getTrknmb() == 2)) ) );
  }

  bool TrackStub::operator>(const TrackStub& rhs) const
  {
    return ( isValid() && ( (!(rhs.isValid())) || (getQuality() > rhs.getQuality()) ||
                            (getQuality() == rhs.getQuality() && cscid() > rhs.cscid()) ||
                            (getQuality() == rhs.getQuality() && cscid() == rhs.cscid() &&
                             (getTrknmb() == 1)) ) );
  }

}