L1Analysis::L1AnalysisGT Class Reference

#include <L1AnalysisGT.h>

Public Member Functions
L1AnalysisGTDataFormat *	getData ()
	L1AnalysisGT ()
void	Reset ()
void	Set (const L1GlobalTriggerReadoutRecord *gtrr)
void	SetEvm (const L1GlobalTriggerEvmReadoutRecord *gtevmrr)
	~L1AnalysisGT ()

Private Attributes
L1AnalysisGTDataFormat	gt_

Detailed Description

Definition at line 17 of file L1AnalysisGT.h.

Constructor & Destructor Documentation

L1AnalysisGT()

L1Analysis::L1AnalysisGT::L1AnalysisGT

(

)

Definition at line 5 of file L1AnalysisGT.cc.

{}

~L1AnalysisGT()

L1Analysis::L1AnalysisGT::~L1AnalysisGT

(

)

Definition at line 7 of file L1AnalysisGT.cc.

{}

Member Function Documentation

getData()

L1AnalysisGTDataFormat* L1Analysis::L1AnalysisGT::getData ( void  )

inline

Definition at line 25 of file L1AnalysisGT.h.

{ return &gt_; }

References gt_.

Reset()

void L1Analysis::L1AnalysisGT::Reset ( )

inline

Definition at line 24 of file L1AnalysisGT.h.

{ gt_.Reset(); }

References gt_, and L1Analysis::L1AnalysisGTDataFormat::Reset().

Set()

void L1Analysis::L1AnalysisGT::Set

(

const L1GlobalTriggerReadoutRecord *

gtrr

)

get Global Trigger algo and technical triger bit statistics

Definition at line 34 of file L1AnalysisGT.cc.

                                                                         {
  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++;
    }
  }
}

References L1GtPsbWord::aData(), L1GtPsbWord::bData(), L1GlobalTriggerReadoutRecord::decisionWord(), L1GlobalTriggerReadoutRecord::gtFdlWord(), L1GlobalTriggerReadoutRecord::gtPsbWord(), L1DTConfigFromDB_cfi::LL, L1GlobalTriggerReadoutRecord::technicalTriggerWord(), and testProducerWithPsetDescEmpty_cfi::x1.

SetEvm()

void L1Analysis::L1AnalysisGT::SetEvm

(

const L1GlobalTriggerEvmReadoutRecord *

gtevmrr

)

Definition at line 9 of file L1AnalysisGT.cc.

                                                                                  {
  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();
}

References L1GtfeExtWord::beamMode(), L1GtfeExtWord::beamMomentum(), L1GtfeExtWord::bstMasterStatus(), L1GtfeExtWord::gpsTime(), L1GlobalTriggerEvmReadoutRecord::gtfeWord(), L1GtfeExtWord::lhcFillNumber(), L1GtfeExtWord::particleTypeBeam1(), L1GtfeExtWord::particleTypeBeam2(), L1TcsWord::partTrigNr(), L1GlobalTriggerEvmReadoutRecord::tcsWord(), L1GtfeExtWord::totalIntensityBeam1(), L1GtfeExtWord::totalIntensityBeam2(), and L1GtfeExtWord::turnCountNumber().

Member Data Documentation

gt_

L1AnalysisGTDataFormat L1Analysis::L1AnalysisGT::gt_

private

Definition at line 28 of file L1AnalysisGT.h.

Referenced by getData(), and Reset().