L1TTwinMuxAlgorithm Class Reference

#include <L1TTwinMuxAlgorithm.h>

Public Member Functions
L1MuDTChambPhContainer	get_ph_tm_output ()
	Return Output PhContainer. More...
	L1TTwinMuxAlgorithm ()
void	run (edm::Handle< L1MuDTChambPhContainer > phiDigis, edm::Handle< L1MuDTChambThContainer > thetaDigis, edm::Handle< RPCDigiCollection > rpcDigis, const edm::EventSetup &c)
	~L1TTwinMuxAlgorithm ()

Private Member Functions
int	radialAngle (RPCDetId, const edm::EventSetup &, int)

Private Attributes
L1MuDTChambPhContainer	m_tm_phi_output
	Output PhContainer. More...
edm::ESHandle< L1TTwinMuxParams >	tmParamsHandle
	Event Setup Handler. More...

Detailed Description

Definition at line 33 of file L1TTwinMuxAlgorithm.h.

Constructor & Destructor Documentation

L1TTwinMuxAlgorithm::L1TTwinMuxAlgorithm ( )

inline

Definition at line 35 of file L1TTwinMuxAlgorithm.h.

{};

L1TTwinMuxAlgorithm::~L1TTwinMuxAlgorithm ( )

inline

Definition at line 36 of file L1TTwinMuxAlgorithm.h.

References EnergyCorrector::c, and run().

{};

Member Function Documentation

L1MuDTChambPhContainer L1TTwinMuxAlgorithm::get_ph_tm_output ( )

inline

Return Output PhContainer.

Definition at line 44 of file L1TTwinMuxAlgorithm.h.

References m_tm_phi_output, and radialAngle().

{  return m_tm_phi_output;}

int L1TTwinMuxAlgorithm::radialAngle ( RPCDetId  , const edm::EventSetup &  , int    )

private

Referenced by get_ph_tm_output().

void L1TTwinMuxAlgorithm::run	(	edm::Handle< L1MuDTChambPhContainer >	phiDigis,
		edm::Handle< L1MuDTChambThContainer >	thetaDigis,
		edm::Handle< RPCDigiCollection >	rpcDigis,
		const edm::EventSetup &	c
	)

ES Parameters

Align track segments that are coming in bx-1.

Clean RPC hits

Translate RPC digis to DT primitives.

Match low q DT primitives with RPC hits in dphiWindow

Correct(in bx) DT primitives by comparing them to RPC.

Add RPC primitives in case that there are no DT primitives.

Definition at line 33 of file L1TTwinMuxAlgorithm.cc.

References L1MuDTChambPhDigi::BxCnt(), L1MuDTChambPhDigi::bxNum(), L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), L1MuDTChambPhDigi::code(), gather_cfg::cout, edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), L1TTwinMuxParams::get_CorrectDTBxwRPC(), L1TTwinMuxParams::get_UseLowQDT(), L1TTwinMuxParams::get_UseOnlyDT(), L1TTwinMuxParams::get_UseOnlyRPC(), L1TTwinMuxParams::get_Verbose(), L1MuDTChambPhDigi::phi(), L1MuDTChambPhDigi::phiB(), IOPrinter::run(), L1MuDTChambPhDigi::scNum(), relativeConstraints::station, L1MuDTChambPhDigi::stNum(), L1MuDTChambPhDigi::Ts2Tag(), makeMuonMisalignmentScenario::wheel, and L1MuDTChambPhDigi::whNum().

Referenced by ~L1TTwinMuxAlgorithm().

                                                                                                            {


  const L1TTwinMuxParamsRcd& tmParamsRcd = c.get<L1TTwinMuxParamsRcd>();
  tmParamsRcd.get(tmParamsHandle);
  const L1TTwinMuxParams& tmParams = *tmParamsHandle.product();
  bool onlyRPC = tmParams.get_UseOnlyRPC();
  bool onlyDT = tmParams.get_UseOnlyDT();
  bool useLowQDT = tmParams.get_UseLowQDT();
  bool correctBX = tmParams.get_CorrectDTBxwRPC();
  bool verbose = tmParams.get_Verbose();



  auto alignedDTs=std::make_shared<AlignTrackSegments>(*inphiDigis);
  alignedDTs->run(c);
  L1MuDTChambPhContainer phiDigis = alignedDTs->getDTContainer();
  //if only DTs are required without bx correction
  //return the aligned track segments
  if(onlyDT && !correctBX && !useLowQDT) {
    m_tm_phi_output = phiDigis;
    return;
  }
  auto rpcHitCl=std::make_shared<RPCHitCleaner>(*rpcDigis);
  rpcHitCl->run(c);
  RPCDigiCollection rpcDigisCleaned = rpcHitCl->getRPCCollection();

  auto dt_from_rpc= std::make_shared<RPCtoDTTranslator>(rpcDigisCleaned);
  dt_from_rpc->run(c);
  L1MuDTChambPhContainer rpcPhiDigis = dt_from_rpc->getDTContainer();            //Primitves used for RPC->DT (only station 1 and 2)
  L1MuDTChambPhContainer rpcHitsPhiDigis = dt_from_rpc->getDTRPCHitsContainer(); //Primitves used for bx correction

  auto dtlowq= std::make_shared<DTLowQMatching>(&phiDigis, rpcHitsPhiDigis);
  dtlowq->run(c);

  if(onlyDT && !correctBX && useLowQDT) {
    m_tm_phi_output = phiDigis;
    if(verbose){
       IOPrinter ioPrinter;
       cout<<"======DT========"<<endl;
       ioPrinter.run(inphiDigis, m_tm_phi_output, rpcDigis, c);
       cout<<"======RPC========"<<endl;
       ioPrinter.run(&rpcHitsPhiDigis, m_tm_phi_output, &rpcDigisCleaned, c);
       cout<<"+++++++++++++++++++++++++++++++++++++++++++++++++"<<endl;
    }

    return;
  }


//  DTRPCBxCorrection *rpc_dt_bx = new DTRPCBxCorrection(phiDigis,rpcHitsPhiDigis);
  auto rpc_dt_bx=std::make_shared<DTRPCBxCorrection>(phiDigis,rpcHitsPhiDigis);
  rpc_dt_bx->run(c);

  L1MuDTChambPhContainer phiDigiscp = rpc_dt_bx->getDTContainer();

  std::vector<L1MuDTChambPhDigi> l1ttma_out;

  int bx=0, wheel=0, sector=0, station=1;

  for(bx=-3; bx<=3; bx++){
    for (wheel=-2;wheel<=2; wheel++ ){
      for (sector=0;sector<12; sector++ ){
        for (station=1; station<=4; station++){

           L1MuDTChambPhDigi const* dtts1 = phiDigiscp.chPhiSegm1(wheel,station,sector,bx);
           L1MuDTChambPhDigi const* dtts2 = phiDigiscp.chPhiSegm2(wheel,station,sector,bx);
           L1MuDTChambPhDigi const* rpcts1 = rpcPhiDigis.chPhiSegm1(wheel,station,sector,bx);

           if(!onlyRPC) {
              if(!dtts1 && !dtts2 && !rpcts1) continue;
              if(dtts1 && dtts1->code()!=7) {
                l1ttma_out.push_back(*dtts1);
              }
              if(dtts2 && dtts2->code()!=7) {
                l1ttma_out.push_back(*dtts2);
              }
              if(!onlyDT){
                if(!dtts1 && !dtts2 && rpcts1 && station<=2 ) {
                  l1ttma_out.emplace_back(rpcts1->bxNum() , rpcts1->whNum(), rpcts1->scNum(), rpcts1->stNum(),rpcts1->phi(), rpcts1->phiB(), rpcts1->code(), rpcts1->Ts2Tag(), rpcts1->BxCnt(),2);
                }
              }
            }

            else if(onlyRPC){
              if( rpcts1 && station<=2 ) {
                l1ttma_out.emplace_back(rpcts1->bxNum() , rpcts1->whNum(), rpcts1->scNum(), rpcts1->stNum(),rpcts1->phi(), rpcts1->phiB(), rpcts1->code(), rpcts1->Ts2Tag(), rpcts1->BxCnt(),2);
              }
            }

        }
      }
    }
  }

m_tm_phi_output.setContainer(l1ttma_out);

  if(verbose){
     IOPrinter ioPrinter;
     cout<<"======DT========"<<endl;
     ioPrinter.run(inphiDigis, m_tm_phi_output, rpcDigis, c);
     cout<<"======RPC========"<<endl;
     ioPrinter.run(&rpcHitsPhiDigis, m_tm_phi_output, &rpcDigisCleaned, c);
     cout<<"+++++++++++++++++++++++++++++++++++++++++++++++++"<<endl;
  }

}

Member Data Documentation

L1MuDTChambPhContainer L1TTwinMuxAlgorithm::m_tm_phi_output

private

Output PhContainer.

Definition at line 49 of file L1TTwinMuxAlgorithm.h.

Referenced by get_ph_tm_output().

edm::ESHandle< L1TTwinMuxParams > L1TTwinMuxAlgorithm::tmParamsHandle

private

Event Setup Handler.

Definition at line 52 of file L1TTwinMuxAlgorithm.h.