L1AnalysisGT.cc

Go to the documentation of this file.

#include <cstdint>
#include "L1Trigger/L1TNtuples/interface/L1AnalysisGT.h"
#include <cstdint>
 
L1Analysis::L1AnalysisGT::L1AnalysisGT() {}
 
L1Analysis::L1AnalysisGT::~L1AnalysisGT() {}
 
void L1Analysis::L1AnalysisGT::SetEvm(const L1GlobalTriggerEvmReadoutRecord* gtevmrr) {
  L1TcsWord tcsw = gtevmrr->tcsWord();
 
  //bx = tcsw.bxNr();
  //lumi = tcsw.luminositySegmentNr();
  //runn = tcsw.partRunNr();
  //eventn = tcsw.partTrigNr();
  //orbitn = tcsw.orbitNr();
  gt_.partrig_tcs = tcsw.partTrigNr();
  L1GtfeExtWord myGtfeExtWord = gtevmrr->gtfeWord();
  uint64_t gpsTime = myGtfeExtWord.gpsTime();
 
  gt_.gpsTimelo = gpsTime & 0xffffffff;
  gt_.gpsTimehi = (gpsTime >> 32) & 0xffffffff;
  gt_.bstMasterStatus = 0xffff & (myGtfeExtWord.bstMasterStatus());
  gt_.bstturnCountNumber = myGtfeExtWord.turnCountNumber();
  gt_.bstlhcFillNumber = myGtfeExtWord.lhcFillNumber();
  gt_.bstbeamMode = 0xffff & (myGtfeExtWord.beamMode());
  gt_.bstparticleTypeBeam1 = 0xffff & (myGtfeExtWord.particleTypeBeam1());
  gt_.bstparticleTypeBeam2 = 0xffff & (myGtfeExtWord.particleTypeBeam2());
  gt_.bstbeamMomentum = 0xffff & (myGtfeExtWord.beamMomentum());
  gt_.bsttotalIntensityBeam1 = myGtfeExtWord.totalIntensityBeam1();
  gt_.bsttotalIntensityBeam2 = myGtfeExtWord.totalIntensityBeam2();
}
 
void L1Analysis::L1AnalysisGT::Set(const L1GlobalTriggerReadoutRecord* gtrr) {
  for (int ibx = -1; ibx <= 1; ibx++) {
    const L1GtPsbWord psb = gtrr->gtPsbWord(0xbb0d, ibx);
    const L1GtPsbWord psb2 = gtrr->gtPsbWord(0xbb0e, ibx);
 
    // ------ ETT, ETM, HTT and HTM from PSB14:
 
    if (ibx == 0) {
      int psb_ett = psb2.aData(4);
      int ett_rank = psb_ett & 0xfff;
      gt_.RankETT = ett_rank;
      gt_.OvETT = (psb_ett >> 12) & 0x1;
 
      int psb_htt = psb2.bData(4);
      int htt_rank = psb_htt & 0xfff;
      gt_.RankHTT = htt_rank;
      gt_.OvHTT = (psb_htt >> 12) & 0x1;
 
      int psb_etmis = psb2.aData(5);
      int etmis_rank = psb_etmis & 0xfff;
      int psb_etmis_phi = psb2.bData(5) & 0x7F;
      gt_.RankETM = etmis_rank;
      gt_.PhiETM = psb_etmis_phi;
      gt_.OvETM = (psb_etmis >> 12) & 0x1;
 
      int psb_htmis = psb2.aData(3);
      int htmis_rank = (psb_htmis >> 5) & 0x7f;
      int htmis_phi = psb_htmis & 0x1F;
      gt_.RankHTM = htmis_rank;
      gt_.PhiHTM = htmis_phi;
      gt_.OvHTM = (psb_htmis >> 12) & 0x1;
    }
 
    // =---------------------------------------------
 
    std::vector<int> psbel;
    psbel.push_back(psb.aData(4));
    psbel.push_back(psb.aData(5));
    psbel.push_back(psb.bData(4));
    psbel.push_back(psb.bData(5));
    std::vector<int>::const_iterator ipsbel;
    for (ipsbel = psbel.begin(); ipsbel != psbel.end(); ipsbel++) {
      float rank = (*ipsbel) & 0x3f;
      if (rank > 0) {
        gt_.Bxel.push_back(ibx);
        gt_.Rankel.push_back(rank);
        gt_.Phiel.push_back(((*ipsbel) >> 10) & 0x1f);
        gt_.Etael.push_back((((*ipsbel >> 9) & 1) ? 10 - (((*ipsbel) >> 6) & 7) : (((*ipsbel) >> 6) & 7) + 11));
        gt_.Isoel.push_back(false);
      }
    }
    psbel.clear();
    psbel.push_back(psb.aData(6));
    psbel.push_back(psb.aData(7));
    psbel.push_back(psb.bData(6));
    psbel.push_back(psb.bData(7));
    for (ipsbel = psbel.begin(); ipsbel != psbel.end(); ipsbel++) {
      float rank = (*ipsbel) & 0x3f;
      if (rank > 0) {
        gt_.Bxel.push_back(ibx);
        gt_.Rankel.push_back(rank);
        gt_.Phiel.push_back(((*ipsbel) >> 10) & 0x1f);
        gt_.Etael.push_back((((*ipsbel >> 9) & 1) ? 10 - (((*ipsbel) >> 6) & 7) : (((*ipsbel) >> 6) & 7) + 11));
        gt_.Isoel.push_back(true);
      }
    }
 
    // central jets
    std::vector<int> psbjet;
    psbjet.push_back(psb.aData(2));
    psbjet.push_back(psb.aData(3));
    psbjet.push_back(psb.bData(2));
    psbjet.push_back(psb.bData(3));
    std::vector<int>::const_iterator ipsbjet;
    for (ipsbjet = psbjet.begin(); ipsbjet != psbjet.end(); ipsbjet++) {
      float rank = (*ipsbjet) & 0x3f;
      if (rank > 0) {
        gt_.Bxjet.push_back(ibx);
        gt_.Rankjet.push_back(rank);
        gt_.Phijet.push_back(((*ipsbjet) >> 10) & 0x1f);
        gt_.Etajet.push_back((((*ipsbjet >> 9) & 1) ? 10 - (((*ipsbjet) >> 6) & 7) : (((*ipsbjet) >> 6) & 7) + 11));
        gt_.Taujet.push_back(false);
        gt_.Fwdjet.push_back(false);
      }
    }
 
    // tau jets
    psbjet.clear();
    psbjet.push_back(psb2.aData(6));
    psbjet.push_back(psb2.aData(7));
    psbjet.push_back(psb2.bData(6));
    psbjet.push_back(psb2.bData(7));
    for (ipsbjet = psbjet.begin(); ipsbjet != psbjet.end(); ipsbjet++) {
      float rank = (*ipsbjet) & 0x3f;
      if (rank > 0) {
        gt_.Bxjet.push_back(ibx);
        gt_.Rankjet.push_back(rank);
        gt_.Phijet.push_back(((*ipsbjet) >> 10) & 0x1f);
        gt_.Etajet.push_back((((*ipsbjet >> 9) & 1) ? 10 - (((*ipsbjet) >> 6) & 7) : (((*ipsbjet) >> 6) & 7) + 11));
        gt_.Taujet.push_back(true);
        gt_.Fwdjet.push_back(false);
      }
    }
 
    // forward jets
    psbjet.clear();
    psbjet.push_back(psb.aData(0));
    psbjet.push_back(psb.aData(1));
    psbjet.push_back(psb.bData(0));
    psbjet.push_back(psb.bData(1));
    for (ipsbjet = psbjet.begin(); ipsbjet != psbjet.end(); ipsbjet++) {
      float rank = (*ipsbjet) & 0x3f;
      if (rank > 0) {
        gt_.Bxjet.push_back(ibx);
        gt_.Rankjet.push_back(rank);
        gt_.Phijet.push_back(((*ipsbjet) >> 10) & 0x1f);
        gt_.Etajet.push_back((((*ipsbjet >> 9) & 1) ? 3 - (((*ipsbjet) >> 6) & 7) : (((*ipsbjet) >> 6) & 7) + 18));
        gt_.Taujet.push_back(false);
        gt_.Fwdjet.push_back(true);
      }
    }
  }
  gt_.Nele = gt_.Bxel.size();
  gt_.Njet = gt_.Bxjet.size();
 
  L1GtFdlWord fdlWord = gtrr->gtFdlWord();
 
  gt_.tw1.resize(5, 0);
  gt_.tw2.resize(5, 0);
  gt_.tt.resize(5, 0);
 
  for (int iebx = 0; iebx < 5; iebx++) {
    DecisionWord gtDecisionWord = gtrr->decisionWord(iebx - 2);
 
    int dbitNumber = 0;
 
    DecisionWord::const_iterator GTdbitItr;
    for (GTdbitItr = gtDecisionWord.begin(); GTdbitItr != gtDecisionWord.end(); GTdbitItr++) {
      if (*GTdbitItr) {
        if (dbitNumber < 64) {
          gt_.tw1[iebx] |= (1LL << dbitNumber);
        } else {
          gt_.tw2[iebx] |= (1LL << (dbitNumber - 64));
        }
      }
      dbitNumber++;
    }
 
    dbitNumber = 0;
    TechnicalTriggerWord gtTTWord = gtrr->technicalTriggerWord(iebx - 2);
    TechnicalTriggerWord::const_iterator GTtbitItr;
    for (GTtbitItr = gtTTWord.begin(); GTtbitItr != gtTTWord.end(); GTtbitItr++) {
      if (*GTtbitItr) {
        gt_.tt[iebx] |= (1LL << dbitNumber);
      }
      dbitNumber++;
    }
  }
}