DTLowQMatching.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: DTLowQMatching
//
//   DTLowQMatching
//
//
//   Author :
//   G. Flouris               U Ioannina    Nov. 2016
//   modifications: GKarathanasis  UAthens
//--------------------------------------------------
#include <iostream>
#include <iomanip>
#include <iterator>

#include "L1Trigger/L1TTwinMux/interface/DTLowQMatching.h"
#include "L1Trigger/L1TTwinMux/interface/DTRPCBxCorrection.h"
#include "L1Trigger/L1TTwinMux/interface/L1MuTMChambPhContainer.h"

using namespace std;

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


void DTLowQMatching::run( const edm::EventSetup& c) {

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

  m_DphiWindow   = tmParams.get_DphiWindowBxShift();

  Matching(0);
  Matching(1);
}

void DTLowQMatching::Matching(int track_seg){

  L1MuDTChambPhDigi * dtts=nullptr;
  L1MuDTChambPhDigi * rpcts1=nullptr;
  auto m_phiRPCDigis_tm=std::make_shared<L1MuTMChambPhContainer>();
  const std::vector<L1MuDTChambPhDigi> *phiChambVector;
  phiChambVector=m_phiRPCDigis.getContainer(); 
  m_phiRPCDigis_tm->setContainer(*phiChambVector );
  
  auto m_phiDTDigis_tm=std::make_shared<L1MuTMChambPhContainer>();
  const std::vector<L1MuDTChambPhDigi> *phiChambVectorDT;
  phiChambVectorDT=m_phiDTDigis->getContainer();
  m_phiDTDigis_tm->setContainer(*phiChambVectorDT );

  int bx=0, wheel=0, sector=0, station=1;
   //cout<<"LowQ Matching  "<<track_seg<<endl;
  for (wheel=-2;wheel<=2; wheel++ ){
   for (sector=0;sector<12; sector++ ){
    for (station=1; station<=4; station++){
     for(bx=-3; bx<=3; bx++){
         int matched = 0;
         for(int rpcbx=bx-1; rpcbx<=bx+1; rpcbx++){
            dtts=nullptr; rpcts1=nullptr; 
            dtts = m_phiDTDigis_tm->chPhiSegm(wheel,station,sector,bx ,track_seg);
            if(!dtts || dtts->code()>=2) continue;
            int nhits = 0;    
            nhits = DTRPCBxCorrection::nRPCHits(m_phiRPCDigis, rpcbx, wheel, sector, station);
            for(int hit=0; hit<nhits; hit++){
              rpcts1 = m_phiRPCDigis_tm->chPhiSegm(wheel, station, sector, rpcbx,hit);
              //If DT primitives with q<2 match with rpc hits do nothing else
              //'remove' the primitive by setting is quality to 7
              if(rpcts1 && DTRPCBxCorrection::deltaPhi(dtts->phi(),rpcts1->phi()) < m_DphiWindow) {
               matched++;
               continue;}
          }
        }//end of rpc bx
     if(matched == 0 && dtts && dtts->code()<2) {//dtts->setCode(7); 
       // int bx=dtts->bxNum(); cout<<bx<<endl;
       L1MuDTChambPhDigi dtts2(dtts->bxNum(),dtts->whNum(),dtts->scNum(),dtts->stNum(),dtts->phi(),dtts->phiB(),7,dtts->Ts2Tag(),dtts->BxCnt(),dtts->RpcBit());
       *dtts=dtts2;      
       } 
   }//end of dt bx
   }//end of station
   }//end of sc
  }//end of wheel
}