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::L1TTwinMuxAlgorithm ( )

inline

Definition at line 35 of file L1TTwinMuxAlgorithm.h.

{};

~L1TTwinMuxAlgorithm()

L1TTwinMuxAlgorithm::~L1TTwinMuxAlgorithm ( )

inline

Definition at line 36 of file L1TTwinMuxAlgorithm.h.

{};

Member Function Documentation

get_ph_tm_output()

L1MuDTChambPhContainer L1TTwinMuxAlgorithm::get_ph_tm_output ( )

inline

Return Output PhContainer.

Definition at line 44 of file L1TTwinMuxAlgorithm.h.

{ return m_tm_phi_output; }

References m_tm_phi_output.

radialAngle()

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

private

run()

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 32 of file L1TTwinMuxAlgorithm.cc.

                                                      {
  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();
 
  AlignTrackSegments alignedDTs{*inphiDigis};
  alignedDTs.run(c);
  L1MuDTChambPhContainer const& 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;
  }
  RPCHitCleaner rpcHitCl{*rpcDigis};
  rpcHitCl.run(c);
  RPCDigiCollection const& rpcDigisCleaned = rpcHitCl.getRPCCollection();
 
  RPCtoDTTranslator dt_from_rpc{rpcDigisCleaned};
  dt_from_rpc.run(c);
  L1MuDTChambPhContainer const& rpcPhiDigis =
      dt_from_rpc.getDTContainer();  //Primitves used for RPC->DT (only station 1 and 2)
  L1MuDTChambPhContainer const& rpcHitsPhiDigis =
      dt_from_rpc.getDTRPCHitsContainer();  //Primitves used for bx correction
 
  DTLowQMatching dtlowq{&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);
  DTRPCBxCorrection rpc_dt_bx{phiDigis, rpcHitsPhiDigis};
  rpc_dt_bx.run(c);
 
  L1MuDTChambPhContainer const& 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;
  }
}

References l1GtPatternGenerator_cfi::bx, L1MuDTChambPhDigi::BxCnt(), L1MuDTChambPhDigi::bxNum(), c, L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), L1MuDTChambPhDigi::code(), gather_cfg::cout, edm::eventsetup::EventSetupRecordImplementation< T >::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(), fakeTwinMuxParams_cff::useLowQDT, verbose, makeMuonMisalignmentScenario::wheel, and L1MuDTChambPhDigi::whNum().

Member Data Documentation

m_tm_phi_output

L1MuDTChambPhContainer L1TTwinMuxAlgorithm::m_tm_phi_output

private

Output PhContainer.

Definition at line 49 of file L1TTwinMuxAlgorithm.h.

Referenced by get_ph_tm_output().

tmParamsHandle

edm::ESHandle<L1TTwinMuxParams> L1TTwinMuxAlgorithm::tmParamsHandle

private

Event Setup Handler.

Definition at line 52 of file L1TTwinMuxAlgorithm.h.