SAMuonCleaner Class Reference

#include <SAMuonCleaner.h>

Public Member Functions
std::vector< l1t::SAMuon >	cleanTFMuons (const std::vector< l1t::SAMuon > &muons)
	SAMuonCleaner ()=default
	~SAMuonCleaner ()=default

Private Member Functions
std::vector< l1t::SAMuon >	cleanTF (const std::vector< l1t::SAMuon > &tfMuons)
std::vector< l1t::SAMuon >	interTFClean (const std::vector< l1t::SAMuon > &bmtf, const std::vector< l1t::SAMuon > &omtf, const std::vector< l1t::SAMuon > &emtf)
void	overlapCleanTrack (l1t::SAMuon &source, const l1t::SAMuon &other, bool eq)
void	overlapCleanTrackInter (l1t::SAMuon &source, const l1t::SAMuon &other)
void	sort (std::vector< l1t::SAMuon > &in)
void	swap (std::vector< l1t::SAMuon > &, int i, int j)

Detailed Description

Definition at line 6 of file SAMuonCleaner.h.

Constructor & Destructor Documentation

SAMuonCleaner()

SAMuonCleaner::SAMuonCleaner ( )

default

~SAMuonCleaner()

SAMuonCleaner::~SAMuonCleaner ( )

default

Member Function Documentation

cleanTF()

std::vector< l1t::SAMuon > SAMuonCleaner::cleanTF ( const std::vector< l1t::SAMuon > &  tfMuons )

private

Definition at line 39 of file SAMuonCleaner.cc.

References mps_fire::i, dqmiolumiharvest::j, MillePedeFileConverter_cfg::out, overlapCleanTrack(), and source.

Referenced by cleanTFMuons().

                                                                               {
  std::vector<l1t::SAMuon> out;
  for (unsigned int i = 0; i < tfMuons.size(); ++i) {
    l1t::SAMuon source = tfMuons[i];
    for (unsigned int j = 0; j < tfMuons.size(); ++j) {
      if (i == j)
        continue;
      overlapCleanTrack(source, tfMuons[j], false);
    }
    if (source.stubs().size() > 1)
      out.push_back(source);
  }
  return out;
}

cleanTFMuons()

std::vector< l1t::SAMuon > SAMuonCleaner::cleanTFMuons

(

const std::vector< l1t::SAMuon > &

muons

)

Definition at line 119 of file SAMuonCleaner.cc.

References l1t::bmtf, cleanTF(), l1t::emtf_neg, l1t::emtf_pos, interTFClean(), amptDefaultParameters_cff::mu, DiMuonV_cfg::muons, l1t::omtf_neg, l1t::omtf_pos, MillePedeFileConverter_cfg::out, and sort().

Referenced by Phase2L1TGMTSAMuonGhostCleaner::prodMuons().

                                                                                  {
  std::vector<l1t::SAMuon> out;

  //split into collections
  std::vector<l1t::SAMuon> bmtf;
  std::vector<l1t::SAMuon> omtf_pos;
  std::vector<l1t::SAMuon> emtf_pos;
  std::vector<l1t::SAMuon> omtf_neg;
  std::vector<l1t::SAMuon> emtf_neg;

  for (const auto& mu : muons) {
    if (mu.tfType() == l1t::tftype::bmtf) {
      bmtf.push_back(mu);
    } else if (mu.tfType() == l1t::tftype::omtf_pos) {
      omtf_pos.push_back(mu);
    } else if (mu.tfType() == l1t::tftype::emtf_pos) {
      emtf_pos.push_back(mu);
    } else if (mu.tfType() == l1t::tftype::omtf_neg) {
      omtf_neg.push_back(mu);
    } else if (mu.tfType() == l1t::tftype::emtf_neg) {
      emtf_neg.push_back(mu);
    }
  }

  std::vector<l1t::SAMuon> omtf_cleaned = cleanTF(omtf_pos);
  std::vector<l1t::SAMuon> omtf_neg_cleaned = cleanTF(omtf_neg);
  omtf_cleaned.insert(omtf_cleaned.end(), omtf_neg_cleaned.begin(), omtf_neg_cleaned.end());
  sort(omtf_cleaned);

  std::vector<l1t::SAMuon> emtf_cleaned = cleanTF(emtf_pos);
  std::vector<l1t::SAMuon> emtf_neg_cleaned = cleanTF(emtf_neg);
  emtf_cleaned.insert(emtf_cleaned.end(), emtf_neg_cleaned.begin(), emtf_neg_cleaned.end());
  sort(emtf_cleaned);

  std::vector<l1t::SAMuon> cleaned = interTFClean(bmtf, omtf_cleaned, emtf_cleaned);
  sort(cleaned);
  return cleaned;
}

interTFClean()

std::vector< l1t::SAMuon > SAMuonCleaner::interTFClean ( const std::vector< l1t::SAMuon > &  bmtf, const std::vector< l1t::SAMuon > &  omtf, const std::vector< l1t::SAMuon > &  emtf  )

private

Definition at line 54 of file SAMuonCleaner.cc.

References l1t::bmtf, mps_fire::i, trackingPlots::other, MillePedeFileConverter_cfg::out, overlapCleanTrackInter(), and source.

Referenced by cleanTFMuons().

                                                                                       {
  std::vector<l1t::SAMuon> out = emtf;
  for (unsigned int i = 0; i < omtf.size(); ++i) {
    l1t::SAMuon source = omtf[i];
    for (const auto& other : emtf) {
      overlapCleanTrackInter(source, other);
    }
    if (source.stubs().size() > 1)
      out.push_back(source);
  }
  for (unsigned int i = 0; i < bmtf.size(); ++i) {
    l1t::SAMuon source = bmtf[i];
    for (const auto& other : omtf) {
      overlapCleanTrackInter(source, other);
    }
    if (source.stubs().size() > 1)
      out.push_back(source);
  }
  return out;
}

overlapCleanTrack()

void SAMuonCleaner::overlapCleanTrack ( l1t::SAMuon &  source, const l1t::SAMuon &  other, bool  eq  )

private

Definition at line 3 of file SAMuonCleaner.cc.

References runTheMatrix::keep, convertSQLiteXML::ok, trackingPlots::other, source, and SimL1Emulator_cff::stubs.

Referenced by cleanTF().

                                                                                        {
  int rank1 = source.hwQual();
  int rank2 = other.hwQual();
  bool keep = true;
  if ((eq && rank1 <= rank2) || ((!eq) && rank1 < rank2))
    keep = false;
  l1t::MuonStubRefVector stubs;
  for (const auto& s1 : source.stubs()) {
    bool ok = true;
    for (const auto& s2 : other.stubs()) {
      if ((*s1) == (*s2) && (!keep))
        ok = false;
    }
    if (ok) {
      stubs.push_back(s1);
    }
  }
  source.setStubs(stubs);
}

overlapCleanTrackInter()

void SAMuonCleaner::overlapCleanTrackInter ( l1t::SAMuon &  source, const l1t::SAMuon &  other  )

private

Definition at line 23 of file SAMuonCleaner.cc.

References runTheMatrix::keep, convertSQLiteXML::ok, trackingPlots::other, source, and SimL1Emulator_cff::stubs.

Referenced by interTFClean().

                                                                                    {
  bool keep = false;
  l1t::MuonStubRefVector stubs;
  for (const auto& s1 : source.stubs()) {
    bool ok = true;
    for (const auto& s2 : other.stubs()) {
      if ((*s1) == (*s2) && (!keep))
        ok = false;
    }
    if (ok) {
      stubs.push_back(s1);
    }
  }
  source.setStubs(stubs);
}

sort()

void SAMuonCleaner::sort ( std::vector< l1t::SAMuon > &  in )

private

Definition at line 96 of file SAMuonCleaner.cc.

References mps_fire::i, recoMuon::in, align_cfg::iteration, MillePedeFileConverter_cfg::out, reco::LeafCandidate::setP4(), swap(), HLT_2024v13_cff::track, and parallelization::uint.

Referenced by cleanTFMuons().

                                                 {
  l1t::SAMuon nullTrack;
  nullTrack.setP4(reco::Candidate::LorentzVector(0, 0, 0, 0.000001));
  while (in.size() < 32)
    in.push_back(nullTrack);

  for (uint iteration = 0; iteration < 16; ++iteration) {
    for (uint i = 0; i < 32; i = i + 2) {
      swap(in, i, i + 1);
    }
    for (uint i = 1; i < 31; i = i + 2) {
      swap(in, i, i + 1);
    }
  }

  std::vector<l1t::SAMuon> out;
  for (const auto& track : in) {
    if (!track.stubs().empty() && out.size() < 12)
      out.push_back(track);
  }
  in = out;
}

swap()

void SAMuonCleaner::swap ( std::vector< l1t::SAMuon > &  list, int  i, int  j  )

private

Definition at line 77 of file SAMuonCleaner.cc.

References mps_fire::i, dqmiolumiharvest::j, reco::LeafCandidate::pt(), HLT_2024v13_cff::pt1, HLT_2024v13_cff::pt2, and createJobs::tmp.

Referenced by sort().

                                                                 {
  const l1t::SAMuon& track1 = list[i];
  const l1t::SAMuon& track2 = list[j];

  int pt1 = track1.pt();
  int pt2 = track2.pt();
  bool swap = false;
  if (pt1 > pt2)
    swap = false;
  else
    swap = true;

  if (swap) {
    l1t::SAMuon tmp = list[i];
    list[i] = list[j];
    list[j] = tmp;
  }
}