Phase2L1TGMTFwdMuonTranslator Class Reference

Inheritance diagram for Phase2L1TGMTFwdMuonTranslator:

Public Member Functions
	Phase2L1TGMTFwdMuonTranslator (const edm::ParameterSet &)
	~Phase2L1TGMTFwdMuonTranslator () override=default
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	associateStubs (l1t::SAMuon &, const l1t::MuonStubRefVector &)
l1t::SAMuon	ConvertEMTFTrack (const l1t::phase2::EMTFTrack &track, const int bx_)
l1t::SAMuon	Convertl1tMuon (const l1t::RegionalMuonCand &mu, const int bx_, bool isDisplaced=false)
void	produce (edm::Event &, const edm::EventSetup &) override
l1t::MuonStubRefVector	selectLayerBX (const l1t::MuonStubRefVector &all, int bx, uint layer)

Private Attributes
edm::EDGetTokenT< l1t::phase2::EMTFTrackCollection >	emtfTrackToken_
edm::EDGetTokenT< l1t::RegionalMuonCandBxCollection >	omtfTrackToken_
edm::EDGetTokenT< l1t::MuonStubCollection >	stubToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 28 of file Phase2L1TGMTFwdMuonTranslator.cc.

Constructor & Destructor Documentation

Phase2L1TGMTFwdMuonTranslator()

Phase2L1TGMTFwdMuonTranslator::Phase2L1TGMTFwdMuonTranslator ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 60 of file Phase2L1TGMTFwdMuonTranslator.cc.

    : stubToken_(consumes<l1t::MuonStubCollection>(iConfig.getParameter<edm::InputTag>("stubs"))),
      omtfTrackToken_(consumes<l1t::RegionalMuonCandBxCollection>(iConfig.getParameter<edm::InputTag>("omtfTracks"))),
      emtfTrackToken_(consumes<l1t::phase2::EMTFTrackCollection>(iConfig.getParameter<edm::InputTag>("emtfTracks"))) {
  produces<std::vector<l1t::SAMuon> >("prompt").setBranchAlias("prompt");
  produces<std::vector<l1t::SAMuon> >("displaced").setBranchAlias("displaced");
}

~Phase2L1TGMTFwdMuonTranslator()

Phase2L1TGMTFwdMuonTranslator::~Phase2L1TGMTFwdMuonTranslator ( )

overridedefault

Member Function Documentation

associateStubs()

void Phase2L1TGMTFwdMuonTranslator::associateStubs ( l1t::SAMuon &  mu, const l1t::MuonStubRefVector &  stubs  )

private

Definition at line 219 of file Phase2L1TGMTFwdMuonTranslator.cc.

References spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, HLT_2024v14_cff::dPhi, mps_fire::i, nano_mu_digi_cff::layer, M_PI, amptDefaultParameters_cff::mu, selectLayerBX(), SimL1Emulator_cff::stubs, and parallelization::uint.

Referenced by produce().

                                                                                                   {
  for (unsigned int layer = 0; layer <= 4; ++layer) {
    l1t::MuonStubRefVector selectedStubs = selectLayerBX(stubs, 0, layer);
    int bestStubINT = -1;
    float dPhi = 1000.0;
    for (uint i = 0; i < selectedStubs.size(); ++i) {
      const l1t::MuonStubRef& stub = selectedStubs[i];
      float deltaPhi =
          (stub->quality() & 0x1) ? stub->offline_coord1() - mu.p4().phi() : stub->offline_coord2() - mu.p4().phi();
      if (deltaPhi > M_PI)
        deltaPhi = deltaPhi - 2 * M_PI;
      if (deltaPhi < -M_PI)
        deltaPhi = deltaPhi + 2 * M_PI;
      deltaPhi = fabs(deltaPhi);
      float deltaEta = (stub->etaQuality() == 0 || (stub->etaQuality() & 0x1))
                           ? fabs(stub->offline_eta1() - mu.p4().eta())
                           : fabs(stub->offline_eta2() - mu.p4().eta());
      if (deltaPhi < 0.3 && deltaEta < 0.3 && deltaPhi < dPhi) {
        dPhi = deltaPhi;
        bestStubINT = i;
      }
    }
    if (bestStubINT >= 0) {
      mu.addStub(selectedStubs[bestStubINT]);
    }
  }
}

ConvertEMTFTrack()

SAMuon Phase2L1TGMTFwdMuonTranslator::ConvertEMTFTrack ( const l1t::phase2::EMTFTrack &  track, const int  bx_  )

private

Definition at line 247 of file Phase2L1TGMTFwdMuonTranslator.cc.

References Phase2L1GMT::BITSETA, Phase2L1GMT::BITSPHI, Phase2L1GMT::BITSPT, Phase2L1GMT::BITSSAD0, emtf::phase2::tp::calcPhiGlobRadFromLoc(), emtf::phase2::tp::calcPhiLocRadFromInt(), emtf::phase2::tp::calcThetaRadFromInt(), ALCARECOTkAlJpsiMuMu_cff::charge, d0, l1t::emtf_neg, l1t::emtf_pos, PVValHelper::eta, dqm-mbProfile::log, Phase2L1GMT::LSBeta, Phase2L1GMT::LSBphi, Phase2L1GMT::LSBpt, phi, DiDispStaMuonMonitor_cfi::pt, l1t::SAMuon::setTF(), funct::tan(), HLT_2024v14_cff::track, leptonTimeLifeInfo_common_cff::track_phi, and HLTMuonOfflineAnalyzer_cfi::z0.

Referenced by produce().

                                                                                                     {
  // Convert EMTF Phi and Theta to Global Phi and Eta
  float track_phi =
      emtf::phase2::tp::calcPhiGlobRadFromLoc(track.sector(), emtf::phase2::tp::calcPhiLocRadFromInt(track.modelPhi()));
  float track_theta = emtf::phase2::tp::calcThetaRadFromInt(track.modelEta());
  float track_eta = -1 * std::log(std::tan(track_theta / 2));

  track_theta *= track.endcap();
  track_eta *= track.endcap();

  // Calculate Lorentz Vector
  // Muon mass taken from L1Trigger/L1TMuon/plugins/L1TMuonProducer.cc
  math::PtEtaPhiMLorentzVector p4(track.emtfPt() * LSBpt, track_eta, track_phi, 0.0);

  // Quantize Values
  ap_uint<BITSSAQUAL> qual = track.emtfQuality();  // Quality provided by EMTF to GMT
  int charge = track.emtfQ();                      // EMTF uses the same convention
  ap_uint<BITSGTPT> pt = track.emtfPt();           // Quantized by EMTF in the same units
  ap_int<BITSGTPHI> phi = round(track_phi / LSBphi);
  ap_int<BITSGTETA> eta = round(track_eta / LSBeta);
  ap_int<BITSSAZ0> z0 = track.emtfZ0();  // Quantized by EMTF in the same units
  ap_int<BITSSAD0> d0 = track.emtfD0();  // Quantized by EMTF in the same units

  //Here do not use the word format to GT but use the word format expected by GMT
  int bstart = 0;
  wordtype word(0);
  bstart = wordconcat<wordtype>(word, bstart, 1, 1);
  bstart = wordconcat<wordtype>(word, bstart, charge, 1);
  bstart = wordconcat<wordtype>(word, bstart, pt, BITSPT);
  bstart = wordconcat<wordtype>(word, bstart, phi, BITSPHI);
  bstart = wordconcat<wordtype>(word, bstart, eta, BITSETA);
  //  bstart = wordconcat<wordtype>(word, bstart, z0, BITSSAZ0); NOT YET SUPPORTED BY GMT
  bstart = wordconcat<wordtype>(word, bstart, d0, BITSSAD0);
  bstart = wordconcat<wordtype>(
      word, bstart, qual, 8);  //FOR NOW 8 bits to be efficienct with Ghost busting. THIS IS ***NOT*** THE FINAL QUALITY

  SAMuon samuon(p4, charge, pt.to_uint(), eta.to_int(), phi.to_int(), z0.to_int(), d0.to_int(), qual.to_uint());

  // +1=Positive Endcap and -1=Negative Endcap
  if (track.endcap() == 1)
    samuon.setTF(tftype::emtf_pos);
  else
    samuon.setTF(tftype::emtf_neg);

  samuon.setWord(word);

  return samuon;
}

Convertl1tMuon()

SAMuon Phase2L1TGMTFwdMuonTranslator::Convertl1tMuon ( const l1t::RegionalMuonCand &  mu, const int  bx_, bool  isDisplaced = false  )

private

Definition at line 160 of file Phase2L1TGMTFwdMuonTranslator.cc.

References Phase2L1GMT::BITSETA, Phase2L1GMT::BITSPHI, Phase2L1GMT::BITSPT, Phase2L1GMT::BITSSAD0, ALCARECOTkAlJpsiMuMu_cff::charge, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), d0, PVValHelper::eta, Phase2L1GMT::LSBeta, Phase2L1GMT::LSBphi, Phase2L1GMT::LSBpt, M_PI, amptDefaultParameters_cff::mu, phi, DiDispStaMuonMonitor_cfi::pt, l1t::SAMuon::setTF(), and HLTMuonOfflineAnalyzer_cfi::z0.

Referenced by produce().

                                                                                                                   {
  ap_uint<BITSSAQUAL> qual = mu.hwQual();
  int charge = mu.hwSign();

  ap_uint<BITSGTPT> pt = 0;
  if (!isDisplaced && mu.hwPt() > 0)
    pt = round(mu.hwPt() * 0.5 / LSBpt);  // Phase-1 LSB 0.5GeV
  if (isDisplaced && mu.hwPtUnconstrained() > 0)
    pt = round(mu.hwPtUnconstrained() * 1.0 / LSBpt);  // Phase-1 LSB 1.0GeV!!

  // BEWARE: THIS CONVERSION IS ONLY VALID FOR OMTF
  constexpr double p1phiLSB = 2 * M_PI / 576;
  // From the uGMTConfiguration of OMTF. OMTF send in local phi!!
  // all others correspond to 120 degree sectors = 192 in int-scale
  int globPhi = mu.processor() * 192 + mu.hwPhi();
  // first processor starts at CMS phi = 15 degrees (24 in int)... Handle wrap-around with %. Add 576 to make sure the number is positive
  globPhi = (globPhi + 600) % 576;
  ap_int<BITSGTPHI> phi = round(globPhi * p1phiLSB / LSBphi);     // Phase-1 LSB (2*pi/576)
  ap_int<BITSGTETA> eta = round(mu.hwEta() * 0.010875 / LSBeta);  // Phase-1 LSB 0.010875

  // FIXME: Below are not well defined in phase1 GMT
  // Using the version from Correlator for now
  ap_int<BITSSAZ0> z0 = 0;  // No tracks info in Phase 1
  // Use 2 bits with LSB = 30cm for BMTF and 25cm for EMTF currently, but subjet to change
  ap_int<BITSSAD0> d0 = mu.hwDXY();

  //Here do not use the word format to GT but use the word format expected by GMT
  int bstart = 0;
  wordtype word(0);
  bstart = wordconcat<wordtype>(word, bstart, 1, 1);
  bstart = wordconcat<wordtype>(word, bstart, charge, 1);
  bstart = wordconcat<wordtype>(word, bstart, pt, BITSPT);
  bstart = wordconcat<wordtype>(word, bstart, phi, BITSPHI);
  bstart = wordconcat<wordtype>(word, bstart, eta, BITSETA);
  //  bstart = wordconcat<wordtype>(word, bstart, z0, BITSSAZ0); NOT YET SUPPORTED BY GMT
  bstart = wordconcat<wordtype>(word, bstart, d0, BITSSAD0);
  bstart = wordconcat<wordtype>(
      word, bstart, qual, 8);  //FOR NOW 8 bits to be efficienct with Ghost busting. THIS IS ***NOT*** THE FINAL QUALITY

  // Calculate Lorentz Vector
  math::PtEtaPhiMLorentzVector p4(pt * LSBpt, eta * LSBeta, phi * LSBphi, 0.0);
  SAMuon samuon(p4, charge, pt.to_uint(), eta.to_int(), phi.to_int(), z0.to_int(), d0.to_int(), qual.to_uint());
  samuon.setTF(mu.trackFinderType());
  samuon.setWord(word);

  return samuon;
}  // -----  end of function Phase2L1TGMTFwdMuonTranslator::Convertl1tMuon  -----

fillDescriptions()

void Phase2L1TGMTFwdMuonTranslator::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 297 of file Phase2L1TGMTFwdMuonTranslator.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                               {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;

  // Input Collections
  desc.add<edm::InputTag>("stubs", edm::InputTag("gmtStubs"));
  desc.add<edm::InputTag>("emtfTracks", edm::InputTag("simEmtfDigisPhase2"));
  desc.add<edm::InputTag>("omtfTracks", edm::InputTag("simOmtfPhase2Digis"));

  // Register
  descriptions.add("Phase2L1TGMTFwdMuonTranslator", desc);
}

produce()

void Phase2L1TGMTFwdMuonTranslator::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 69 of file Phase2L1TGMTFwdMuonTranslator.cc.

References associateStubs(), BXVector< T >::at(), ConvertEMTFTrack(), Convertl1tMuon(), BPHMonitor_cff::displaced, emtfTrackToken_, mps_fire::i, iEvent, eostools::move(), amptDefaultParameters_cff::mu, omtfTrackToken_, BXVector< T >::size(), SimL1Emulator_cff::stubs, stubToken_, HLT_2024v14_cff::track, and parallelization::uint.

                                                                                         {
  using namespace edm;

  edm::Handle<l1t::MuonStubCollection> stubHandle;
  iEvent.getByToken(stubToken_, stubHandle);

  edm::Handle<l1t::phase2::EMTFTrackCollection> emtf_tracks;
  iEvent.getByToken(emtfTrackToken_, emtf_tracks);

  edm::Handle<l1t::RegionalMuonCandBxCollection> omtf_tracks;
  iEvent.getByToken(omtfTrackToken_, omtf_tracks);

  // Process Stubs
  l1t::MuonStubRefVector stubs;

  for (uint i = 0; i < stubHandle->size(); ++i) {
    l1t::MuonStubRef stub(stubHandle, i);

    if (stub->bxNum() == 0)
      stubs.push_back(stub);
  }

  // Collect Muons
  std::vector<SAMuon> prompt;

  //  TODO: Will receive hybrid stubs from OMTF/EMTF
  std::vector<SAMuon> displaced;

  // Convert OMTF Muons to SAMuons
  for (unsigned int i = 0; i < omtf_tracks->size(0); ++i) {
    const l1t::RegionalMuonCand& mu = omtf_tracks->at(0, i);
    // Since OMTF is using Phase-1 LSB, will convert to SAMuon locally
    // We should move to passing words in future
    l1t::SAMuon samuon;
    if (mu.hwPt() > 0) {
      samuon = Convertl1tMuon(mu, 0);
      //now associate the stubs
      associateStubs(samuon, stubs);
      prompt.push_back(samuon);
    }
    if (mu.hwPtUnconstrained() > 0) {
      samuon = Convertl1tMuon(mu, 0, true);
      //now associate the stubs
      associateStubs(samuon, stubs);
      displaced.push_back(samuon);
    }
  }

  // Convert EMTF++ Tracks to SAMuons
  for (unsigned int track_id = 0; track_id < emtf_tracks->size(); ++track_id) {
    const auto& track = emtf_tracks->at(track_id);

    // Short-Circuit: Only keep valid tracks that are in BX=0
    if ((track.valid() == 0) || (track.bx() != 0)) {
      continue;
    }

    // Short-Circuit: Only keep tracks with quality above 3 to avoid single hit tracks
    if (track.emtfQuality() <= 3) {
      continue;
    }

    // Short-Circuit: Only keep tracks with the max relevance score (127)
    if (track.emtfRels() != 127) {
      continue;
    }

    auto samuon = ConvertEMTFTrack(track, 0);

    //now associate the stubs
    associateStubs(samuon, stubs);

    // Add To Collections
    if (track.unconstrained()) {
      displaced.push_back(samuon);
    } else {
      prompt.push_back(samuon);
    }
  }

  // Output Prompt Muon Collection
  std::unique_ptr<std::vector<l1t::SAMuon> > prompt_ptr = std::make_unique<std::vector<l1t::SAMuon> >(prompt);
  std::unique_ptr<std::vector<l1t::SAMuon> > displaced_ptr = std::make_unique<std::vector<l1t::SAMuon> >(displaced);
  iEvent.put(std::move(prompt_ptr), "prompt");
  iEvent.put(std::move(displaced_ptr), "displaced");
}

selectLayerBX()

l1t::MuonStubRefVector Phase2L1TGMTFwdMuonTranslator::selectLayerBX ( const l1t::MuonStubRefVector &  all, int  bx, uint  layer  )

private

Definition at line 208 of file Phase2L1TGMTFwdMuonTranslator.cc.

References python.cmstools::all(), nano_mu_digi_cff::bx, nano_mu_digi_cff::layer, and MillePedeFileConverter_cfg::out.

Referenced by associateStubs().

                                                                                {
  l1t::MuonStubRefVector out;
  for (const auto& stub : all) {
    if (stub->bxNum() == bx && stub->tfLayer() == layer)
      out.push_back(stub);
  }
  return out;
}

Member Data Documentation

emtfTrackToken_

edm::EDGetTokenT<l1t::phase2::EMTFTrackCollection> Phase2L1TGMTFwdMuonTranslator::emtfTrackToken_

private

Definition at line 47 of file Phase2L1TGMTFwdMuonTranslator.cc.

Referenced by produce().

omtfTrackToken_

edm::EDGetTokenT<l1t::RegionalMuonCandBxCollection> Phase2L1TGMTFwdMuonTranslator::omtfTrackToken_

private

Definition at line 46 of file Phase2L1TGMTFwdMuonTranslator.cc.

Referenced by produce().

stubToken_

edm::EDGetTokenT<l1t::MuonStubCollection> Phase2L1TGMTFwdMuonTranslator::stubToken_

private

Definition at line 45 of file Phase2L1TGMTFwdMuonTranslator.cc.

Referenced by produce().