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 L1TTwinMuxParams &, const RPCGeometry &)
	~L1TTwinMuxAlgorithm ()

Private Attributes
L1MuDTChambPhContainer	m_tm_phi_output
	Output PhContainer. More...

Detailed Description

Definition at line 27 of file L1TTwinMuxAlgorithm.h.

Constructor & Destructor Documentation

L1TTwinMuxAlgorithm()

L1TTwinMuxAlgorithm::L1TTwinMuxAlgorithm ( )

inline

Definition at line 29 of file L1TTwinMuxAlgorithm.h.

{};

~L1TTwinMuxAlgorithm()

L1TTwinMuxAlgorithm::~L1TTwinMuxAlgorithm ( )

inline

Definition at line 30 of file L1TTwinMuxAlgorithm.h.

{};

Member Function Documentation

get_ph_tm_output()

L1MuDTChambPhContainer L1TTwinMuxAlgorithm::get_ph_tm_output ( )

inline

Return Output PhContainer.

Definition at line 39 of file L1TTwinMuxAlgorithm.h.

{ return m_tm_phi_output; }

References m_tm_phi_output.

Referenced by L1TTwinMuxProducer::produce().

run()

void L1TTwinMuxAlgorithm::run	(	edm::Handle< L1MuDTChambPhContainer >	phiDigis,
		edm::Handle< L1MuDTChambThContainer >	thetaDigis,
		edm::Handle< RPCDigiCollection >	rpcDigis,
		const L1TTwinMuxParams &	tmParams,
		const RPCGeometry &	rpcGeometry
	)

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.

                                                              {
  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();
  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();
  RPCDigiCollection const& rpcDigisCleaned = rpcHitCl.getRPCCollection();
 
  RPCtoDTTranslator dt_from_rpc{rpcDigisCleaned};
  dt_from_rpc.run(rpcGeometry);
  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(tmParams);
 
  if (onlyDT && !correctBX && useLowQDT) {
    m_tm_phi_output = phiDigis;
    if (verbose) {
      IOPrinter ioPrinter;
      cout << "======DT========" << endl;
      ioPrinter.run(inphiDigis, m_tm_phi_output, rpcDigis, rpcGeometry);
      cout << "======RPC========" << endl;
      ioPrinter.run(&rpcHitsPhiDigis, m_tm_phi_output, &rpcDigisCleaned, rpcGeometry);
      cout << "+++++++++++++++++++++++++++++++++++++++++++++++++" << endl;
    }
 
    return;
  }
 
  //  DTRPCBxCorrection *rpc_dt_bx = new DTRPCBxCorrection(phiDigis,rpcHitsPhiDigis);
  DTRPCBxCorrection rpc_dt_bx{phiDigis, rpcHitsPhiDigis};
  rpc_dt_bx.run(tmParams);
 
  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, rpcGeometry);
    cout << "======RPC========" << endl;
    ioPrinter.run(&rpcHitsPhiDigis, m_tm_phi_output, &rpcDigisCleaned, rpcGeometry);
    cout << "+++++++++++++++++++++++++++++++++++++++++++++++++" << endl;
  }
}

References l1GtPatternGenerator_cfi::bx, L1MuDTChambPhDigi::BxCnt(), L1MuDTChambPhDigi::bxNum(), L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), L1MuDTChambPhDigi::code(), gather_cfg::cout, 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().

Referenced by L1TTwinMuxProducer::produce().

Member Data Documentation

m_tm_phi_output

L1MuDTChambPhContainer L1TTwinMuxAlgorithm::m_tm_phi_output

private

Output PhContainer.

Definition at line 43 of file L1TTwinMuxAlgorithm.h.

Referenced by get_ph_tm_output().