L1AnalysisL1Upgrade.cc

Go to the documentation of this file.

#include "L1Trigger/L1TNtuples/interface/L1AnalysisL1Upgrade.h"

L1Analysis::L1AnalysisL1Upgrade::L1AnalysisL1Upgrade() {}

L1Analysis::L1AnalysisL1Upgrade::~L1AnalysisL1Upgrade() {}

void L1Analysis::L1AnalysisL1Upgrade::SetEm(const edm::Handle<l1t::EGammaBxCollection> em, unsigned maxL1Upgrade) {
  for (int ibx = em->getFirstBX(); ibx <= em->getLastBX(); ++ibx) {
    for (l1t::EGammaBxCollection::const_iterator it = em->begin(ibx);
         it != em->end(ibx) && l1upgrade_.nEGs < maxL1Upgrade;
         it++) {
      if (it->pt() > 0) {
        l1upgrade_.egEt.push_back(it->pt());
        l1upgrade_.egEta.push_back(it->eta());
        l1upgrade_.egPhi.push_back(it->phi());
        l1upgrade_.egIEt.push_back(it->hwPt());
        l1upgrade_.egIEta.push_back(it->hwEta());
        l1upgrade_.egIPhi.push_back(it->hwPhi());
        l1upgrade_.egIso.push_back(it->hwIso());
        l1upgrade_.egBx.push_back(ibx);
        l1upgrade_.egTowerIPhi.push_back(it->towerIPhi());
        l1upgrade_.egTowerIEta.push_back(it->towerIEta());
        l1upgrade_.egRawEt.push_back(it->rawEt());
        l1upgrade_.egIsoEt.push_back(it->isoEt());
        l1upgrade_.egFootprintEt.push_back(it->footprintEt());
        l1upgrade_.egNTT.push_back(it->nTT());
        l1upgrade_.egShape.push_back(it->shape());
        l1upgrade_.egTowerHoE.push_back(it->towerHoE());
        l1upgrade_.egHwQual.push_back(it->hwQual());
        l1upgrade_.nEGs++;
      }
    }
  }
}

void L1Analysis::L1AnalysisL1Upgrade::SetTau(const edm::Handle<l1t::TauBxCollection> tau, unsigned maxL1Upgrade) {
  for (int ibx = tau->getFirstBX(); ibx <= tau->getLastBX(); ++ibx) {
    for (l1t::TauBxCollection::const_iterator it = tau->begin(ibx);
         it != tau->end(ibx) && l1upgrade_.nTaus < maxL1Upgrade;
         it++) {
      if (it->pt() > 0) {
        l1upgrade_.tauEt.push_back(it->et());
        l1upgrade_.tauEta.push_back(it->eta());
        l1upgrade_.tauPhi.push_back(it->phi());
        l1upgrade_.tauIEt.push_back(it->hwPt());
        l1upgrade_.tauIEta.push_back(it->hwEta());
        l1upgrade_.tauIPhi.push_back(it->hwPhi());
        l1upgrade_.tauIso.push_back(it->hwIso());
        l1upgrade_.tauBx.push_back(ibx);
        l1upgrade_.tauTowerIPhi.push_back(it->towerIPhi());
        l1upgrade_.tauTowerIEta.push_back(it->towerIEta());
        l1upgrade_.tauRawEt.push_back(it->rawEt());
        l1upgrade_.tauIsoEt.push_back(it->isoEt());
        l1upgrade_.tauNTT.push_back(it->nTT());
        l1upgrade_.tauHasEM.push_back(it->hasEM());
        l1upgrade_.tauIsMerged.push_back(it->isMerged());
        l1upgrade_.tauHwQual.push_back(it->hwQual());
        l1upgrade_.nTaus++;
      }
    }
  }
}

void L1Analysis::L1AnalysisL1Upgrade::SetJet(const edm::Handle<l1t::JetBxCollection> jet, unsigned maxL1Upgrade) {
  for (int ibx = jet->getFirstBX(); ibx <= jet->getLastBX(); ++ibx) {
    for (l1t::JetBxCollection::const_iterator it = jet->begin(ibx);
         it != jet->end(ibx) && l1upgrade_.nJets < maxL1Upgrade;
         it++) {
      if (it->pt() > 0) {
        l1upgrade_.jetEt.push_back(it->et());
        l1upgrade_.jetEta.push_back(it->eta());
        l1upgrade_.jetPhi.push_back(it->phi());
        l1upgrade_.jetIEt.push_back(it->hwPt());
        l1upgrade_.jetIEta.push_back(it->hwEta());
        l1upgrade_.jetIPhi.push_back(it->hwPhi());
        l1upgrade_.jetBx.push_back(ibx);
        l1upgrade_.jetRawEt.push_back(it->rawEt());
        l1upgrade_.jetSeedEt.push_back(it->seedEt());
        l1upgrade_.jetTowerIEta.push_back(it->towerIEta());
        l1upgrade_.jetTowerIPhi.push_back(it->towerIPhi());
        l1upgrade_.jetPUEt.push_back(it->puEt());
        l1upgrade_.jetPUDonutEt0.push_back(it->puDonutEt(0));
        l1upgrade_.jetPUDonutEt1.push_back(it->puDonutEt(1));
        l1upgrade_.jetPUDonutEt2.push_back(it->puDonutEt(2));
        l1upgrade_.jetPUDonutEt3.push_back(it->puDonutEt(3));
        l1upgrade_.nJets++;
      }
    }
  }
}

void L1Analysis::L1AnalysisL1Upgrade::SetMuon(const edm::Handle<l1t::MuonBxCollection> muon, unsigned maxL1Upgrade) {
  for (int ibx = muon->getFirstBX(); ibx <= muon->getLastBX(); ++ibx) {
    for (l1t::MuonBxCollection::const_iterator it = muon->begin(ibx);
         it != muon->end(ibx) && l1upgrade_.nMuons < maxL1Upgrade;
         it++) {
      if (it->pt() > 0) {
        l1upgrade_.muonEt.push_back(it->et());
        l1upgrade_.muonEtUnconstrained.push_back(it->ptUnconstrained());
        l1upgrade_.muonEta.push_back(it->eta());
        l1upgrade_.muonPhi.push_back(it->phi());
        l1upgrade_.muonEtaAtVtx.push_back(it->etaAtVtx());
        l1upgrade_.muonPhiAtVtx.push_back(it->phiAtVtx());
        l1upgrade_.muonIEt.push_back(it->hwPt());
        l1upgrade_.muonIEtUnconstrained.push_back(it->hwPtUnconstrained());
        l1upgrade_.muonIEta.push_back(it->hwEta());
        l1upgrade_.muonIPhi.push_back(it->hwPhi());
        l1upgrade_.muonIEtaAtVtx.push_back(it->hwEtaAtVtx());
        l1upgrade_.muonIPhiAtVtx.push_back(it->hwPhiAtVtx());
        l1upgrade_.muonIDEta.push_back(it->hwDEtaExtra());
        l1upgrade_.muonIDPhi.push_back(it->hwDPhiExtra());
        l1upgrade_.muonChg.push_back(it->charge());
        l1upgrade_.muonIso.push_back(it->hwIso());
        l1upgrade_.muonQual.push_back(it->hwQual());
        l1upgrade_.muonDxy.push_back(it->hwDXY());
        l1upgrade_.muonTfMuonIdx.push_back(it->tfMuonIndex());
        l1upgrade_.muonBx.push_back(ibx);
        l1upgrade_.nMuons++;
      }
    }
  }
}

void L1Analysis::L1AnalysisL1Upgrade::SetSum(const edm::Handle<l1t::EtSumBxCollection> sums, unsigned maxL1Upgrade) {
  for (int ibx = sums->getFirstBX(); ibx <= sums->getLastBX(); ++ibx) {
    for (l1t::EtSumBxCollection::const_iterator it = sums->begin(ibx);
         it != sums->end(ibx) && l1upgrade_.nSums < maxL1Upgrade;
         it++) {
      int type = static_cast<int>(it->getType());
      l1upgrade_.sumType.push_back(type);
      l1upgrade_.sumEt.push_back(it->et());
      l1upgrade_.sumPhi.push_back(it->phi());
      l1upgrade_.sumIEt.push_back(it->hwPt());
      l1upgrade_.sumIPhi.push_back(it->hwPhi());
      l1upgrade_.sumBx.push_back(ibx);
      l1upgrade_.nSums++;
    }
  }
}