DTRPCBxCorrection Class Reference

#include <DTRPCBxCorrection.h>

Public Member Functions
	DTRPCBxCorrection (L1MuDTChambPhContainer, L1MuDTChambPhContainer)
L1MuDTChambPhContainer	getDTContainer ()
	Return Output PhContainer. More...
void	run (const edm::EventSetup &c)
	~DTRPCBxCorrection ()

Static Public Member Functions
static int	deltaPhi (int dt_phi, int rpc_strip)
static int	nRPCHits (L1MuTMChambPhContainer inCon, int bx, int wh, int sec, int st)
static int	nRPCHits (L1MuDTChambPhContainer inCon, int bx, int wh, int sec, int st)

Public Attributes
edm::ESHandle< L1TTwinMuxParams >	tmParamsHandle

Private Member Functions
void	BxCorrection (int track_seg)
int	flipBit (int inv)
int	sign (float)

Private Attributes
int	m_DphiWindow
L1MuDTChambPhContainer	m_dt_tsshifted
std::vector< L1MuDTChambPhDigi >	m_l1ttma_out
L1MuDTChambPhContainer	m_phiDTDigis
L1MuDTChambPhContainer	m_phiRPCDigis
int	m_QualityLimit

Detailed Description

Definition at line 32 of file DTRPCBxCorrection.h.

Constructor & Destructor Documentation

DTRPCBxCorrection::DTRPCBxCorrection	(	L1MuDTChambPhContainer	inphiDTDigis,
		L1MuDTChambPhContainer	inphiRPCDigis
	)

Definition at line 21 of file DTRPCBxCorrection.cc.

    : m_phiDTDigis(inphiDTDigis),
      m_phiRPCDigis(inphiRPCDigis){
          //  m_phiDTDigis=inphiDTDigis;
          //  m_phiRPCDigis=inphiRPCDigis;
      };

DTRPCBxCorrection::~DTRPCBxCorrection ( )

inline

Definition at line 35 of file DTRPCBxCorrection.h.

References HltBtagPostValidation_cff::c, and run().

{};

Member Function Documentation

void DTRPCBxCorrection::BxCorrection ( int  track_seg )

private

Concatanate all vectors in one

Shift primitives if vector>0

XXX: Source of discrepancies in order to send as second TS the two prims must come from different halves of the station this information does not exist in data 'simulate' this information by requiring different sign in phis and dphi>100

if there is no shift then put the original primitive

Definition at line 42 of file DTRPCBxCorrection.cc.

References l1GtPatternGenerator_cfi::bx, L1MuDTChambPhDigi::BxCnt(), L1MuDTChambPhDigi::bxNum(), L1MuTMChambPhContainer::chPhiSegm(), L1MuDTChambPhDigi::code(), dumpMFGeometry_cfg::delta, deltaPhi(), HLT_2018_cff::distance, flipBit(), L1MuDTChambPhContainer::getContainer(), m_DphiWindow, m_l1ttma_out, m_phiDTDigis, m_phiRPCDigis, m_QualityLimit, nhits, nRPCHits(), L1MuDTChambPhDigi::phi(), L1MuDTChambPhDigi::phiB(), L1MuDTChambPhDigi::scNum(), L1MuTMChambPhContainer::setContainer(), sign(), relativeConstraints::station, L1MuDTChambPhDigi::stNum(), L1MuDTChambPhDigi::Ts2Tag(), makeMuonMisalignmentScenario::wheel, and L1MuDTChambPhDigi::whNum().

Referenced by flipBit(), and run().

                                                  {
  L1MuTMChambPhContainer m_phiDTDigis_tm;
  //std::shared_ptr<L1MuTMChambPhContainer> m_phiDTDigis_tm (new L1MuTMChambPhContainer);
  const std::vector<L1MuDTChambPhDigi>* phiChambVectorDT;
  phiChambVectorDT = m_phiDTDigis.getContainer();
  m_phiDTDigis_tm.setContainer(*phiChambVectorDT);
  L1MuTMChambPhContainer m_phiRPCDigis_tm;
  //std::shared_ptr<L1MuTMChambPhContainer> m_phiRPCDigis_tm (new L1MuTMChambPhContainer);
  const std::vector<L1MuDTChambPhDigi>* phiChambVectorRPC;
  phiChambVectorRPC = m_phiRPCDigis.getContainer();
  m_phiRPCDigis_tm.setContainer(*phiChambVectorRPC);

  int ibx_dtm = 0, fbx_dtm = 0;
  int ibx_dtp = 0, fbx_dtp = 0;

  for (int wheel = -2; wheel <= 2; wheel++) {
    for (int sector = 0; sector < 12; sector++) {
      for (int station = 1; station <= 4; station++) {
        bool shifted[7] = {false, false, false, false, false, false, false};
        bool dups[7] = {false, false, false, false, false, false, false};
        bool secondTs[7] = {false, false, false, false, false, false, false};
        L1MuTMChambPhContainer shiftedPhiDTDigis;
        L1MuDTChambPhDigi* dtts_sh2nd = nullptr;
        for (int bx = 3; bx >= -3; bx--) {
          vector<int> delta_m, delta_p, delta_0;
          for (int rpcbx = bx - 1; rpcbx <= bx + 1; rpcbx++) {
            L1MuDTChambPhDigi* dtts = nullptr;
            L1MuDTChambPhDigi* rpcts1 = nullptr;
            dtts = m_phiDTDigis_tm.chPhiSegm(wheel, station, sector, bx, track_seg);

            if (!dtts)
              continue;
            int nhits = nRPCHits(m_phiRPCDigis_tm, rpcbx, wheel, sector, station);
            for (int hit = 0; hit < nhits; hit++) {
              rpcts1 = m_phiRPCDigis_tm.chPhiSegm(wheel, station, sector, rpcbx, hit);
              //Store in vectors the dphi of matched dt/rpc
              if (rpcts1 && dtts && dtts->code() < m_QualityLimit &&
                  deltaPhi(dtts->phi(), rpcts1->phi()) < m_DphiWindow) {
                if ((dtts->bxNum() - rpcbx) == -1) {
                  delta_m.push_back(deltaPhi(dtts->phi(), rpcts1->phi()));
                  ibx_dtm = dtts->bxNum();
                  fbx_dtm = rpcbx;
                }
                if ((dtts->bxNum() - rpcbx) == 0) {
                  delta_0.push_back(deltaPhi(dtts->phi(), rpcts1->phi()));
                }
                if ((dtts->bxNum() - rpcbx) == 1) {
                  delta_p.push_back(deltaPhi(dtts->phi(), rpcts1->phi()));
                  ibx_dtp = dtts->bxNum();
                  fbx_dtp = rpcbx;
                }
              }  //end if dtts and quality
            }
          }  //end of rpc bx and dtts, rpcts1, dttsnew go out of scope

          vector<int> delta;
          delta.insert(delta.end(), delta_0.begin(), delta_0.end());
          delta.insert(delta.end(), delta_p.begin(), delta_p.end());
          delta.insert(delta.end(), delta_m.begin(), delta_m.end());
          if (!delta.empty()) {
            L1MuDTChambPhDigi* dtts = nullptr;
            L1MuDTChambPhDigi* dttsnew = nullptr;
            L1MuDTChambPhDigi* dtts_sh = nullptr;
            std::vector<L1MuDTChambPhDigi> l1ttma_outsh;
            //Find the pair the min dphi(rpc,dt)
            unsigned int min_index = std::distance(delta.begin(), std::min_element(delta.begin(), delta.end())) + 0;
            int init_bx = 0, final_bx = 0;

            if (((delta_0.size() <= min_index) && (min_index < (delta_0.size() + delta_p.size())) &&
                 !delta_p.empty())) {
              init_bx = ibx_dtp;
              final_bx = fbx_dtp;
            } else if ((delta_0.size() + delta_p.size()) <= min_index && !delta_m.empty()) {
              init_bx = ibx_dtm;
              final_bx = fbx_dtm;
            } else
              continue;
            //Primitve to be shifted in place of dttsnew
            dtts = m_phiDTDigis_tm.chPhiSegm(wheel, station, sector, init_bx, track_seg);
            dttsnew = m_phiDTDigis_tm.chPhiSegm(wheel, station, sector, final_bx, track_seg);
            bool shift_1 = false;
            if (dtts && dtts->code() < m_QualityLimit && (!dttsnew || shifted[final_bx + 3] || dups[final_bx + 3])) {
              dtts_sh = new L1MuDTChambPhDigi(final_bx,
                                              dtts->whNum(),
                                              dtts->scNum(),
                                              dtts->stNum(),
                                              dtts->phi(),
                                              dtts->phiB(),
                                              dtts->code(),
                                              dtts->Ts2Tag(),
                                              dtts->BxCnt(),
                                              1);
              l1ttma_outsh.push_back(*dtts_sh);
              shifted[init_bx + 3] = true;
              shift_1 = true;
            }
            if (dtts && dtts->code() < m_QualityLimit && dttsnew)
              dups[init_bx + 3] = true;

            //dtts exists and qual lt m_QualityLimit and dttsnew exists and the previous (shift_1) prim was not shifted and there is empty space in second TS
            if (dtts && dtts->code() < m_QualityLimit && dttsnew && !shift_1 &&
                !m_phiDTDigis_tm.chPhiSegm(wheel, station, sector, final_bx, flipBit(track_seg))) {
              if (sign(dtts->phi()) != sign(dttsnew->phi()) && deltaPhi(-(dtts->phi()), dttsnew->phi()) > 100) {
                dtts_sh2nd = new L1MuDTChambPhDigi(final_bx,
                                                   dtts->whNum(),
                                                   dtts->scNum(),
                                                   dtts->stNum(),
                                                   dtts->phi(),
                                                   dtts->phiB(),
                                                   dtts->code(),
                                                   flipBit(track_seg),
                                                   dtts->BxCnt(),
                                                   1);
                secondTs[final_bx + 3] = true;
                dups[init_bx + 3] = false;
                shifted[init_bx + 3] = true;
              }
            }
            shiftedPhiDTDigis.setContainer(l1ttma_outsh);
          }
        }  //end of bx

        for (int bx = -3; bx <= 3; bx++) {
          L1MuDTChambPhDigi* dtts = nullptr;
          if (secondTs[bx + 3])
            if (dtts_sh2nd) {
              m_l1ttma_out.push_back(*dtts_sh2nd);
            }
          dtts = shiftedPhiDTDigis.chPhiSegm(wheel, station, sector, bx, track_seg);
          if (dtts) {
            m_l1ttma_out.push_back(*dtts);
            continue;
          }
          if (dups[bx + 3])
            continue;
          dtts = m_phiDTDigis_tm.chPhiSegm(wheel, station, sector, bx, track_seg);
          if (!shifted[bx + 3] && dtts) {
            m_l1ttma_out.push_back(*dtts);
          }
        }

      }  //end of station
    }    //end of sc
  }      //end of wheel
}

int DTRPCBxCorrection::deltaPhi ( int  dt_phi, int  rpc_strip  )

static

Definition at line 195 of file DTRPCBxCorrection.cc.

References funct::abs(), and HLT_2018_cff::delta_phi.

Referenced by BxCorrection(), getDTContainer(), and DTLowQMatching::Matching().

                                                          {
  int delta_phi = abs(dt_phi - rpc2dt_phi);
  return delta_phi;
}

int DTRPCBxCorrection::flipBit ( int  inv )

inlineprivate

Definition at line 50 of file DTRPCBxCorrection.h.

References BxCorrection().

Referenced by BxCorrection().

{ return (inv ^ 1); };

L1MuDTChambPhContainer DTRPCBxCorrection::getDTContainer ( )

inline

Return Output PhContainer.

Definition at line 42 of file DTRPCBxCorrection.h.

References l1GtPatternGenerator_cfi::bx, deltaPhi(), m_dt_tsshifted, nRPCHits(), fileinputsource_cfi::sec, and sign().

{ return m_dt_tsshifted; }

int DTRPCBxCorrection::nRPCHits ( L1MuTMChambPhContainer  inCon, int  bx, int  wh, int  sec, int  st  )

static

Definition at line 208 of file DTRPCBxCorrection.cc.

References L1MuTMChambPhContainer::getContainer(), mps_fire::i, and findQualityFiles::size.

Referenced by BxCorrection(), getDTContainer(), DTLowQMatching::Matching(), and IOPrinter::run().

                                                                                             {
  int size = 0;
  const std::vector<L1MuDTChambPhDigi>* vInCon = inCon.getContainer();
  for (auto& i : *vInCon) {
    if (bx == i.bxNum() && i.code() != 7 && i.whNum() == wh && i.scNum() == sec && i.stNum() == st)
      size++;
  }

  return size;
}

int DTRPCBxCorrection::nRPCHits ( L1MuDTChambPhContainer  inCon, int  bx, int  wh, int  sec, int  st  )

static

Definition at line 219 of file DTRPCBxCorrection.cc.

References L1MuDTChambPhContainer::getContainer(), mps_fire::i, and findQualityFiles::size.

                                                                                             {
  int size = 0;
  const std::vector<L1MuDTChambPhDigi>* vInCon = inCon.getContainer();
  for (auto& i : *vInCon) {
    if (bx == i.bxNum() && i.code() != 7 && i.whNum() == wh && i.scNum() == sec && i.stNum() == st)
      size++;
  }

  return size;
}

void DTRPCBxCorrection::run

(

const edm::EventSetup &

c

)

Definition at line 28 of file DTRPCBxCorrection.cc.

References BxCorrection(), edm::EventSetup::get(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), L1TTwinMuxParams::get_DphiWindowBxShift(), L1TTwinMuxParams::get_USERPCBXFORDTBELOWQUALITY(), m_DphiWindow, m_dt_tsshifted, m_l1ttma_out, m_QualityLimit, edm::ESHandle< T >::product(), L1MuDTChambPhContainer::setContainer(), and tmParamsHandle.

Referenced by ~DTRPCBxCorrection().

                                                  {
  const L1TTwinMuxParamsRcd& tmParamsRcd = c.get<L1TTwinMuxParamsRcd>();
  tmParamsRcd.get(tmParamsHandle);
  const L1TTwinMuxParams& tmParams = *tmParamsHandle.product();

  m_QualityLimit = tmParams.get_USERPCBXFORDTBELOWQUALITY();
  m_DphiWindow = tmParams.get_DphiWindowBxShift();

  BxCorrection(0);
  BxCorrection(1);

  m_dt_tsshifted.setContainer(m_l1ttma_out);
}

int DTRPCBxCorrection::sign ( float  inv )

private

Definition at line 200 of file DTRPCBxCorrection.cc.

Referenced by BxCorrection(), and getDTContainer().

                                     {
  if (inv < 0)
    return -1;
  if (inv > 0)
    return 1;
  return 0;
}

Member Data Documentation

int DTRPCBxCorrection::m_DphiWindow

private

Definition at line 62 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

L1MuDTChambPhContainer DTRPCBxCorrection::m_dt_tsshifted

private

Definition at line 57 of file DTRPCBxCorrection.h.

Referenced by getDTContainer(), and run().

std::vector<L1MuDTChambPhDigi> DTRPCBxCorrection::m_l1ttma_out

private

Definition at line 59 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

L1MuDTChambPhContainer DTRPCBxCorrection::m_phiDTDigis

private

Definition at line 55 of file DTRPCBxCorrection.h.

Referenced by BxCorrection().

L1MuDTChambPhContainer DTRPCBxCorrection::m_phiRPCDigis

private

Definition at line 56 of file DTRPCBxCorrection.h.

Referenced by BxCorrection().

int DTRPCBxCorrection::m_QualityLimit

private

Definition at line 61 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

edm::ESHandle<L1TTwinMuxParams> DTRPCBxCorrection::tmParamsHandle

Definition at line 39 of file DTRPCBxCorrection.h.

Referenced by run().