DTRPCBxCorrection Class Reference

#include <DTRPCBxCorrection.h>

Public Member Functions
	DTRPCBxCorrection (L1MuDTChambPhContainer, L1MuDTChambPhContainer)
L1MuDTChambPhContainer	getDTContainer ()
	Return Output PhContainer. More...
void	run (const L1TTwinMuxParams &)
	~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)

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 24 of file DTRPCBxCorrection.h.

Constructor & Destructor Documentation

DTRPCBxCorrection()

DTRPCBxCorrection::DTRPCBxCorrection	(	L1MuDTChambPhContainer	inphiDTDigis,
		L1MuDTChambPhContainer	inphiRPCDigis
	)

Definition at line 22 of file DTRPCBxCorrection.cc.

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

~DTRPCBxCorrection()

DTRPCBxCorrection::~DTRPCBxCorrection ( )

inline

Definition at line 27 of file DTRPCBxCorrection.h.

{};

Member Function Documentation

BxCorrection()

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 39 of file DTRPCBxCorrection.cc.

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

Referenced by 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
}

deltaPhi()

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

static

Definition at line 192 of file DTRPCBxCorrection.cc.

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

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

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

flipBit()

int DTRPCBxCorrection::flipBit ( int  inv )

inlineprivate

Definition at line 40 of file DTRPCBxCorrection.h.

Referenced by BxCorrection().

{ return (inv ^ 1); };

getDTContainer()

L1MuDTChambPhContainer DTRPCBxCorrection::getDTContainer ( )

inline

Return Output PhContainer.

Definition at line 32 of file DTRPCBxCorrection.h.

References m_dt_tsshifted.

{ return m_dt_tsshifted; }

nRPCHits() [1/2]

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

static

Definition at line 205 of file DTRPCBxCorrection.cc.

References simKBmtfDigis_cfi::bx, L1MuTMChambPhContainer::getContainer(), mps_fire::i, fileinputsource_cfi::sec, and findQualityFiles::size.

Referenced by BxCorrection(), 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;
}

nRPCHits() [2/2]

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

static

Definition at line 216 of file DTRPCBxCorrection.cc.

References simKBmtfDigis_cfi::bx, L1MuDTChambPhContainer::getContainer(), mps_fire::i, fileinputsource_cfi::sec, 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;
}

run()

void DTRPCBxCorrection::run

(

const L1TTwinMuxParams &

tmParams

)

Definition at line 29 of file DTRPCBxCorrection.cc.

References BxCorrection(), L1TTwinMuxParams::get_DphiWindowBxShift(), L1TTwinMuxParams::get_USERPCBXFORDTBELOWQUALITY(), m_DphiWindow, m_dt_tsshifted, m_l1ttma_out, m_QualityLimit, and L1MuDTChambPhContainer::setContainer().

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

  BxCorrection(0);
  BxCorrection(1);

  m_dt_tsshifted.setContainer(m_l1ttma_out);
}

sign()

int DTRPCBxCorrection::sign ( float  inv )

private

Definition at line 197 of file DTRPCBxCorrection.cc.

Referenced by BxCorrection().

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

Member Data Documentation

m_DphiWindow

int DTRPCBxCorrection::m_DphiWindow

private

Definition at line 52 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

m_dt_tsshifted

L1MuDTChambPhContainer DTRPCBxCorrection::m_dt_tsshifted

private

Definition at line 47 of file DTRPCBxCorrection.h.

Referenced by getDTContainer(), and run().

m_l1ttma_out

std::vector<L1MuDTChambPhDigi> DTRPCBxCorrection::m_l1ttma_out

private

Definition at line 49 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

m_phiDTDigis

L1MuDTChambPhContainer DTRPCBxCorrection::m_phiDTDigis

private

Definition at line 45 of file DTRPCBxCorrection.h.

Referenced by BxCorrection().

m_phiRPCDigis

L1MuDTChambPhContainer DTRPCBxCorrection::m_phiRPCDigis

private

Definition at line 46 of file DTRPCBxCorrection.h.

Referenced by BxCorrection().

m_QualityLimit

int DTRPCBxCorrection::m_QualityLimit

private

Definition at line 51 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().