#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 (L1MuDTChambPhContainer inCon, int bx, int wh, int sec, int st)

static int	nRPCHits (L1MuTMChambPhContainer 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::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::~DTRPCBxCorrection ( )

inline

Definition at line 35 of file DTRPCBxCorrection.h.

35 {};

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

                                                   {
   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
 }

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 run().

◆ deltaPhi()

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

static

Definition at line 195 of file DTRPCBxCorrection.cc.

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

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

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

◆ flipBit()

int DTRPCBxCorrection::flipBit ( int inv )

inlineprivate

Definition at line 50 of file DTRPCBxCorrection.h.

50 { return (inv ^ 1); };

Referenced by BxCorrection().

◆ getDTContainer()

L1MuDTChambPhContainer DTRPCBxCorrection::getDTContainer ( )

inline

Return Output PhContainer.

Definition at line 42 of file DTRPCBxCorrection.h.

42 { return m_dt_tsshifted; }

References m_dt_tsshifted.

◆ nRPCHits() [1/2]

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

static

Definition at line 219 of file DTRPCBxCorrection.cc.

                                                                                              {
   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;
 }

References l1GtPatternGenerator_cfi::bx, L1MuDTChambPhContainer::getContainer(), mps_fire::i, fileinputsource_cfi::sec, and findQualityFiles::size.

◆ nRPCHits() [2/2]

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

static

Definition at line 208 of file DTRPCBxCorrection.cc.

                                                                                              {
   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;
 }

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

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

◆ run()

void DTRPCBxCorrection::run ( const edm::EventSetup & c )

Definition at line 28 of file DTRPCBxCorrection.cc.

                                                   {
   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);
 }

References BxCorrection(), HltBtagPostValidation_cff::c, 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.

◆ sign()

int DTRPCBxCorrection::sign ( float inv )

private

Definition at line 200 of file DTRPCBxCorrection.cc.

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

Referenced by BxCorrection().

Member Data Documentation

◆ m_DphiWindow

int DTRPCBxCorrection::m_DphiWindow

private

Definition at line 62 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

◆ m_dt_tsshifted

L1MuDTChambPhContainer DTRPCBxCorrection::m_dt_tsshifted

private

Definition at line 57 of file DTRPCBxCorrection.h.

Referenced by getDTContainer(), and run().

◆ m_l1ttma_out

std::vector<L1MuDTChambPhDigi> DTRPCBxCorrection::m_l1ttma_out

private

Definition at line 59 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

◆ m_phiDTDigis

L1MuDTChambPhContainer DTRPCBxCorrection::m_phiDTDigis

private

Definition at line 55 of file DTRPCBxCorrection.h.

Referenced by BxCorrection().

◆ m_phiRPCDigis

L1MuDTChambPhContainer DTRPCBxCorrection::m_phiRPCDigis

private

Definition at line 56 of file DTRPCBxCorrection.h.

Referenced by BxCorrection().

◆ m_QualityLimit

int DTRPCBxCorrection::m_QualityLimit

private

Definition at line 61 of file DTRPCBxCorrection.h.

Referenced by BxCorrection(), and run().

◆ tmParamsHandle

edm::ESHandle<L1TTwinMuxParams> DTRPCBxCorrection::tmParamsHandle

Definition at line 39 of file DTRPCBxCorrection.h.

Referenced by run().

Public Member Functions

Static Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DTRPCBxCorrection()

◆ ~DTRPCBxCorrection()

Member Function Documentation

◆ BxCorrection()

◆ deltaPhi()

◆ flipBit()

◆ getDTContainer()

◆ nRPCHits() [1/2]

◆ nRPCHits() [2/2]

◆ run()

◆ sign()

Member Data Documentation

◆ m_DphiWindow

◆ m_dt_tsshifted

◆ m_l1ttma_out

◆ m_phiDTDigis

◆ m_phiRPCDigis

◆ m_QualityLimit

◆ tmParamsHandle