L1AnalysisL1UpgradeTfMuon.cc

Go to the documentation of this file.

#include "L1Trigger/L1TNtuples/interface/L1AnalysisL1UpgradeTfMuon.h"
#include <cmath>
L1Analysis::L1AnalysisL1UpgradeTfMuon::L1AnalysisL1UpgradeTfMuon() {}
 
L1Analysis::L1AnalysisL1UpgradeTfMuon::~L1AnalysisL1UpgradeTfMuon() {}
 
void L1Analysis::L1AnalysisL1UpgradeTfMuon::SetTfMuon(const l1t::RegionalMuonCandBxCollection& muon,
                                                      unsigned maxL1UpgradeTfMuon) {
  for (int ibx = muon.getFirstBX(); ibx <= muon.getLastBX(); ++ibx) {
    for (l1t::RegionalMuonCandBxCollection::const_iterator it = muon.begin(ibx);
         it != muon.end(ibx) && l1upgradetfmuon_.nTfMuons < maxL1UpgradeTfMuon;
         ++it) {
      if (it->hwPt() > 0) {
        l1upgradetfmuon_.tfMuonHwPt.push_back(it->hwPt());
        l1upgradetfmuon_.tfMuonHwEta.push_back(it->hwEta());
        l1upgradetfmuon_.tfMuonHwPhi.push_back(it->hwPhi());
        l1upgradetfmuon_.tfMuonGlobalPhi.push_back(
            l1t::MicroGMTConfiguration::calcGlobalPhi(it->hwPhi(), it->trackFinderType(), it->processor()));
        l1upgradetfmuon_.tfMuonHwSign.push_back(it->hwSign());
        l1upgradetfmuon_.tfMuonHwSignValid.push_back(it->hwSignValid());
        l1upgradetfmuon_.tfMuonHwQual.push_back(it->hwQual());
        l1upgradetfmuon_.tfMuonLink.push_back(it->link());
        l1upgradetfmuon_.tfMuonProcessor.push_back(it->processor());
        l1upgradetfmuon_.tfMuonTrackFinderType.push_back(it->trackFinderType());
        l1upgradetfmuon_.tfMuonHwHF.push_back(it->hwHF());
        l1upgradetfmuon_.tfMuonBx.push_back(ibx);
        std::map<std::string, int> decoded_track_address;
        if (it->trackFinderType() == l1t::tftype::bmtf) {
          int detSide = it->trackSubAddress(l1t::RegionalMuonCand::kWheelSide);
          int wheelNum = it->trackSubAddress(l1t::RegionalMuonCand::kWheelNum);
          int stat1 = it->trackSubAddress(l1t::RegionalMuonCand::kStat1);
          int stat2 = it->trackSubAddress(l1t::RegionalMuonCand::kStat2);
          int stat3 = it->trackSubAddress(l1t::RegionalMuonCand::kStat3);
          int stat4 = it->trackSubAddress(l1t::RegionalMuonCand::kStat4);
 
          int wheel = pow(-1, detSide) * wheelNum;
          l1upgradetfmuon_.tfMuonWh.push_back(wheel);
          l1upgradetfmuon_.tfMuonTrAdd.push_back(stat1);
          l1upgradetfmuon_.tfMuonTrAdd.push_back(stat2);
          l1upgradetfmuon_.tfMuonTrAdd.push_back(stat3);
          l1upgradetfmuon_.tfMuonTrAdd.push_back(stat4);
          decoded_track_address["wheel"] = wheel;
          decoded_track_address["station1"] = stat1;
          decoded_track_address["station2"] = stat2;
          decoded_track_address["station3"] = stat3;
          decoded_track_address["station4"] = stat4;
        } else if (it->trackFinderType() == l1t::omtf_neg || it->trackFinderType() == l1t::omtf_pos) {
          decoded_track_address["kLayers"] = it->trackSubAddress(l1t::RegionalMuonCand::kLayers);
          decoded_track_address["kWeight"] = it->trackSubAddress(l1t::RegionalMuonCand::kWeight);
        } else if (it->trackFinderType() == l1t::emtf_neg || it->trackFinderType() == l1t::emtf_pos) {
          decoded_track_address["kME1Seg"] = it->trackSubAddress(l1t::RegionalMuonCand::kME1Seg);
          decoded_track_address["kME1Ch"] = it->trackSubAddress(l1t::RegionalMuonCand::kME1Ch);
          decoded_track_address["kME2Seg"] = it->trackSubAddress(l1t::RegionalMuonCand::kME2Seg);
          decoded_track_address["kME2Ch"] = it->trackSubAddress(l1t::RegionalMuonCand::kME2Ch);
          decoded_track_address["kME3Seg"] = it->trackSubAddress(l1t::RegionalMuonCand::kME3Seg);
          decoded_track_address["kME3Ch"] = it->trackSubAddress(l1t::RegionalMuonCand::kME3Ch);
          decoded_track_address["kME4Seg"] = it->trackSubAddress(l1t::RegionalMuonCand::kME4Seg);
          decoded_track_address["kME4Ch"] = it->trackSubAddress(l1t::RegionalMuonCand::kME4Ch);
          decoded_track_address["kTrkNum"] = it->trackSubAddress(l1t::RegionalMuonCand::kTrkNum);
          decoded_track_address["kBX"] = it->trackSubAddress(l1t::RegionalMuonCand::kBX);
        }
        l1upgradetfmuon_.tfMuonDecodedTrAdd.push_back(decoded_track_address);
        l1upgradetfmuon_.tfMuonHwTrAdd.push_back(l1t::RegionalMuonRawDigiTranslator::generateRawTrkAddress(
            *it, true));  // TODO: We're assuming that we're dealing with Kalman muons here.
        l1upgradetfmuon_.nTfMuons++;
      }
    }
  }
}