L1AnalysisGMT.cc

Go to the documentation of this file.

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

L1Analysis::L1AnalysisGMT::L1AnalysisGMT() {}

L1Analysis::L1AnalysisGMT::~L1AnalysisGMT() {}

void L1Analysis::L1AnalysisGMT::Set(const L1MuGMTReadoutCollection* gmtrc,
                                    unsigned maxDTBX,
                                    unsigned maxCSC,
                                    unsigned maxRPC,
                                    unsigned maxGMT,
                                    bool physVal) {
  std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
  std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr;
  for (igmtrr = gmt_records.begin(); igmtrr != gmt_records.end(); igmtrr++) {
    std::vector<L1MuRegionalCand>::const_iterator iter1;
    std::vector<L1MuRegionalCand> rmc;

    if (igmtrr->getBxInEvent() == 0) {
      gmt_.EvBx = igmtrr->getBxNr();
    }

    //
    // DTBX Trigger
    //

    int iidt = 0;
    rmc = igmtrr->getDTBXCands();
    for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
      if ((unsigned)gmt_.Ndt < maxDTBX && !(*iter1).empty()) {
        gmt_.Bxdt.push_back((*iter1).bx());
        if (physVal) {
          gmt_.Etadt.push_back(float((*iter1).etaValue()));
          gmt_.Phidt.push_back(float((*iter1).phiValue()));
          gmt_.Ptdt.push_back(float((*iter1).ptValue()));
        } else {
          gmt_.Etadt.push_back(float((*iter1).eta_packed()));
          gmt_.Phidt.push_back(float((*iter1).phi_packed()));
          gmt_.Ptdt.push_back(float((*iter1).pt_packed()));
        }
        gmt_.Chadt.push_back((*iter1).chargeValue());
        if (!(*iter1).chargeValid())
          gmt_.Chadt.push_back(0);
        gmt_.FineEtadt.push_back((*iter1).isFineHalo());
        gmt_.Qualdt.push_back((*iter1).quality());
        gmt_.Dwdt.push_back((*iter1).getDataWord());
        gmt_.Chdt.push_back(iidt);

        gmt_.Ndt++;
      }
      iidt++;
    }

    //
    // CSC Trigger
    //

    rmc = igmtrr->getCSCCands();
    for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
      if ((unsigned)gmt_.Ncsc < maxCSC && !(*iter1).empty()) {
        gmt_.Bxcsc.push_back((*iter1).bx());
        if (physVal) {
          gmt_.Etacsc.push_back((*iter1).etaValue());
          gmt_.Phicsc.push_back((*iter1).phiValue());
          gmt_.Ptcsc.push_back((*iter1).ptValue());
        } else {
          gmt_.Etacsc.push_back((*iter1).eta_packed());
          gmt_.Phicsc.push_back((*iter1).phi_packed());
          gmt_.Ptcsc.push_back((*iter1).pt_packed());
        }
        gmt_.Chacsc.push_back((*iter1).chargeValue());
        if (!(*iter1).chargeValid())
          gmt_.Chacsc.push_back(0);
        gmt_.Qualcsc.push_back((*iter1).quality());
        gmt_.Dwcsc.push_back((*iter1).getDataWord());

        gmt_.Ncsc++;
      }
    }

    //
    // RPCb Trigger
    //

    rmc = igmtrr->getBrlRPCCands();
    for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
      if ((unsigned)gmt_.Nrpcb < maxRPC && !(*iter1).empty()) {
        gmt_.Bxrpcb.push_back((*iter1).bx());
        if (physVal) {
          gmt_.Etarpcb.push_back((*iter1).etaValue());
          gmt_.Phirpcb.push_back((*iter1).phiValue());
          gmt_.Ptrpcb.push_back((*iter1).ptValue());
        } else {
          gmt_.Etarpcb.push_back((*iter1).eta_packed());
          gmt_.Phirpcb.push_back((*iter1).phi_packed());
          gmt_.Ptrpcb.push_back((*iter1).pt_packed());
        }
        gmt_.Charpcb.push_back((*iter1).chargeValue());
        if (!(*iter1).chargeValid())
          gmt_.Charpcb.push_back(0);
        gmt_.Qualrpcb.push_back((*iter1).quality());
        gmt_.Dwrpcb.push_back((*iter1).getDataWord());

        gmt_.Nrpcb++;
      }
    }

    //
    // RPCf Trigger
    //

    rmc = igmtrr->getFwdRPCCands();
    for (iter1 = rmc.begin(); iter1 != rmc.end(); iter1++) {
      if ((unsigned)gmt_.Nrpcf < maxRPC && !(*iter1).empty()) {
        gmt_.Bxrpcf.push_back((*iter1).bx());
        if (physVal) {
          gmt_.Etarpcf.push_back((*iter1).etaValue());
          gmt_.Phirpcf.push_back((*iter1).phiValue());
          gmt_.Ptrpcf.push_back((*iter1).ptValue());
        } else {
          gmt_.Etarpcf.push_back((*iter1).eta_packed());
          gmt_.Phirpcf.push_back((*iter1).phi_packed());
          gmt_.Ptrpcf.push_back((*iter1).pt_packed());
        }
        gmt_.Charpcf.push_back((*iter1).chargeValue());
        if (!(*iter1).chargeValid())
          gmt_.Charpcf.push_back(0);
        gmt_.Qualrpcf.push_back((*iter1).quality());
        gmt_.Dwrpcf.push_back((*iter1).getDataWord());

        gmt_.Nrpcf++;
      }
    }

    //
    // GMT_. Trigger
    //

    std::vector<L1MuGMTExtendedCand>::const_iterator gmt_iter;
    std::vector<L1MuGMTExtendedCand> exc = igmtrr->getGMTCands();
    for (gmt_iter = exc.begin(); gmt_iter != exc.end(); gmt_iter++) {
      if ((unsigned)gmt_.N < maxGMT && !(*gmt_iter).empty()) {
        gmt_.CandBx.push_back((*gmt_iter).bx());
        if (physVal) {
          gmt_.Eta.push_back((*gmt_iter).etaValue());
          gmt_.Phi.push_back((*gmt_iter).phiValue());
          gmt_.Pt.push_back((*gmt_iter).ptValue());
        } else {
          gmt_.Eta.push_back((*gmt_iter).etaIndex());
          gmt_.Phi.push_back((*gmt_iter).phiIndex());
          gmt_.Pt.push_back((*gmt_iter).ptIndex());
        }
        gmt_.Cha.push_back((*gmt_iter).charge());
        if (!(*gmt_iter).charge_valid())
          gmt_.Cha.push_back(0);
        gmt_.Qual.push_back((*gmt_iter).quality());
        gmt_.Det.push_back((*gmt_iter).detector());
        gmt_.Rank.push_back((*gmt_iter).rank());
        gmt_.Isol.push_back((*gmt_iter).isol());
        gmt_.Mip.push_back((*gmt_iter).mip());
        gmt_.Dw.push_back((*gmt_iter).getDataWord());

        gmt_.IdxRPCb.push_back(-1);
        gmt_.IdxRPCf.push_back(-1);
        gmt_.IdxDTBX.push_back(-1);
        gmt_.IdxCSC.push_back(-1);

        if ((*gmt_iter).isMatchedCand() || (*gmt_iter).isRPC()) {
          if ((*gmt_iter).isFwd()) {
            gmt_.IdxRPCf.back() = (*gmt_iter).getRPCIndex();
          } else {
            gmt_.IdxRPCb.back() = (*gmt_iter).getRPCIndex();
          }
        }

        if ((*gmt_iter).isMatchedCand() || (!(*gmt_iter).isRPC())) {
          if ((*gmt_iter).isFwd()) {
            gmt_.IdxCSC.back() = (*gmt_iter).getDTCSCIndex();
          } else {
            gmt_.IdxDTBX.back() = (*gmt_iter).getDTCSCIndex();
          }
        }
        gmt_.N++;
      }
    }
  }
}