L1TdeGCT Class Reference

#include <L1TdeGCT.h>

Inheritance diagram for L1TdeGCT:

Public Member Functions
	L1TdeGCT (const edm::ParameterSet &)
	~L1TdeGCT () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
int	verbose ()

Private Attributes
std::string	cLabel [nGctColl_]
int	colCount [nGctColl_]
int	colCount_stage1Layer2 [nStage1Layer2Coll_]
edm::EDGetTokenT < L1DataEmulRecord >	DEsource_
MonitorElement *	deword [nGctColl_]
MonitorElement *	deword_stage1layer2 [nStage1Layer2Coll_]
MonitorElement *	dword [nGctColl_]
MonitorElement *	dword_stage1layer2 [nStage1Layer2Coll_]
std::string	errLabel [nerr] = {"Agree", "Loc. Agree", "L.Disagree", "Data only", "Emul only"}
MonitorElement *	errortype [nGctColl_]
MonitorElement *	errortype_stage1layer2 [nStage1Layer2Coll_]
MonitorElement *	eta [nGctColl_]
MonitorElement *	eta_stage1layer2 [nStage1Layer2Coll_]
MonitorElement *	etaData [nGctColl_]
MonitorElement *	etaData_stage1layer2 [nStage1Layer2Coll_]
const double	etaMaxim = 21.5
const double	etaMinim = -0.5
const int	etaNBins = 22
MonitorElement *	etaphi [nGctColl_]
MonitorElement *	etaphi_stage1layer2 [nStage1Layer2Coll_]
MonitorElement *	eword [nGctColl_]
MonitorElement *	eword_stage1layer2 [nStage1Layer2Coll_]
bool	hasRecord_
std::string	histFile_
std::string	histFolder_
bool	m_stage1_layer2_
MonitorElement *	masked [nGctColl_]
MonitorElement *	masked_stage1layer2 [nStage1Layer2Coll_]
bool	monitorDaemon_
const int	nbit = 32
int	nWithCol [nGctColl_]
int	nWithCol_stage1Layer2 [nStage1Layer2Coll_]
MonitorElement *	phi [nGctColl_]
MonitorElement *	phi_stage1layer2 [nStage1Layer2Coll_]
MonitorElement *	phiData [nGctColl_]
MonitorElement *	phiData_stage1layer2 [nStage1Layer2Coll_]
const double	phiMaxim = 17.5
const double	phiMinim = -0.5
const int	phiNBins = 18
MonitorElement *	rnk [nGctColl_]
MonitorElement *	rnk_stage1layer2 [nStage1Layer2Coll_]
MonitorElement *	rnkData [nGctColl_]
MonitorElement *	rnkData_stage1layer2 [nStage1Layer2Coll_]
std::string	sLabel [nStage1Layer2Coll_]
MonitorElement *	sysncand [2]
MonitorElement *	sysrates
int	verbose_

Static Private Attributes
static const int	nerr = 5
static const int	nGctColl_ = dedefs::GCThfbit - dedefs::GCTisolaem + 1
static const int	nStage1Layer2Coll_ = dedefs::GCTisotaujets - dedefs::GCTisolaem + 1

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 27 of file L1TdeGCT.h.

Constructor & Destructor Documentation

L1TdeGCT::L1TdeGCT ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 10 of file L1TdeGCT.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and verbose.

                                                 {
  verbose_ = iConfig.getUntrackedParameter<int>("VerboseFlag", 0);

  if (verbose())
    std::cout << "L1TdeGCT::L1TdeGCT()...\n" << std::flush;

  DEsource_ = consumes<L1DataEmulRecord>(iConfig.getParameter<edm::InputTag>("DataEmulCompareSource"));
  histFolder_ = iConfig.getUntrackedParameter<std::string>("HistFolder");

  histFile_ = iConfig.getUntrackedParameter<std::string>("HistFile", "");
  if (iConfig.getUntrackedParameter<bool>("disableROOToutput", true))
    histFile_ = "";

  if (!histFile_.empty()) {
    edm::LogInfo("OutputRootFile") << "L1TEmulator GCT specific histograms will be saved to " << histFile_.c_str()
                                   << std::endl;
  }

  hasRecord_ = true;

  if (verbose())
    std::cout << "L1TdeGCT::L1TdeGCT()...done.\n" << std::flush;
  m_stage1_layer2_ = iConfig.getParameter<bool>("stage1_layer2_");
}

L1TdeGCT::~L1TdeGCT ( )

override

Definition at line 35 of file L1TdeGCT.cc.

{}

Member Function Documentation

void L1TdeGCT::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprotectedvirtual

get the comparison results

get the de candidates

— Fill histograms(me) —

(note see L1Trigger/HardwareValidation/interface/DEtrait.h)

bitset loop

event based rate

object based rate

get the de candidates

— Fill histograms(me) —

(note see L1Trigger/HardwareValidation/interface/DEtrait.h)

bitset loop

event based rate

object based rate

Reimplemented from DQMEDAnalyzer.

Definition at line 336 of file L1TdeGCT.cc.

References edm::HandleBase::clear(), gather_cfg::cout, TauDecayModes::dec, PVValHelper::eta, fastmath::etaphi(), dedefs::GCT, dedefs::GCTisolaem, edm::Event::getByToken(), mps_fire::i, edm::HandleBase::isValid(), dqmiolumiharvest::j, LogDebug, RPCpg::rate(), L1DataEmulDigi::reset(), dedefs::SystLabel, and verbose.

                                                                          {
  if (!hasRecord_)
    return;

  if (verbose())
    std::cout << "L1TdeGCT::analyze()  start\n" << std::flush;

  edm::Handle<L1DataEmulRecord> deRecord;
  iEvent.getByToken(DEsource_, deRecord);

  if (!deRecord.isValid()) {
    edm::LogInfo("DataNotFound") << "Cannot find L1DataEmulRecord"
                                 << " Please verify that comparator was successfully executed."
                                 << " Emulator DQM for GCT will be skipped!" << std::endl;
    hasRecord_ = false;
    return;
  }

  if (m_stage1_layer2_ == false) {
    bool isComp = deRecord->get_isComp(GCT);
    if (!isComp) {
      if (verbose())
        std::cout << "[L1TdeGCT] Gct information not generated in de-record."
                  << " Skiping event!\n"
                  << std::flush;
      return;
    }

    int DEncand[2] = {0};
    for (int j = 0; j < 2; j++)
      DEncand[j] = deRecord->getNCand(GCT, j);

    if (verbose())
      std::cout << "[L1TdeGCT] ncands"
                << " data: " << DEncand[0] << " emul: " << DEncand[1] << std::endl;

    L1DEDigiCollection deColl;
    deColl = deRecord->getColl();

    // extract the GCT comparison digis
    L1DEDigiCollection gctColl;

    gctColl.reserve(20);
    gctColl.clear();

    for (L1DEDigiCollection::const_iterator it = deColl.begin(); it != deColl.end(); it++)
      if (!it->empty())
        if (it->sid() == GCT)
          gctColl.push_back(*it);

    if (verbose()) {
      std::cout << "[L1TdeGCT] record has " << gctColl.size() << " gct de digis\n" << std::flush;
      for (L1DEDigiCollection::const_iterator it = gctColl.begin(); it != gctColl.end(); it++)
        std::cout << "\t" << *it << std::endl;
    }

    const int nullVal = L1DataEmulDigi().reset();


    // d|e candidate loop
    for (L1DEDigiCollection::const_iterator it = gctColl.begin(); it != gctColl.end(); it++) {
      // sid should be GCT
      int sid = it->sid();
      // cid: GCTisolaem, GCTnoisoem, GCTcenjets, GCTforjets, GCTtaujets
      int cid = it->cid();

      if (verbose())
        std::cout << "[L1TdeGCT] processing digi "
                  << " sys:" << sid << " type:" << cid << " \n\t" << *it << "\n"
                  << std::flush;

      //assert(cid==GCT);
      if (sid != GCT || it->empty()) {
        LogDebug("L1TdeGCT") << "consistency check failure, non-gct digis!";
        continue;
      }

      int type = it->type();
      double phiv = it->x1();
      double etav = it->x2();
      float rankarr[2];
      it->rank(rankarr);
      float rnkv = rankarr[0];

      double wei;

      unsigned int mask = (~0x0);

      // shift coll type for starting at zero
      int ccid = cid - dedefs::GCTisolaem;
      if (ccid < 0 || ccid >= nGctColl_) {
        LogDebug("L1TdeGCT") << "consistency check failure, col type outbounds:" << ccid << "\n";
        ccid = 0;
      }

      //type: 0:agree 1:loc.agree, 2:loc.disagree, 3:data.only, 4:emul.only
      if (type < 4)
        sysncand[0]->Fill(ccid);
      if (type < 5 && type != 3)
        sysncand[1]->Fill(ccid);

      errortype[ccid]->Fill(type);

      wei = (type == 0) ? 0. : 1.;
      if (etav != nullVal && phiv != nullVal)
        etaphi[ccid]->Fill(etav, phiv, wei);
      if (etav != nullVal)
        eta[ccid]->Fill(etav, wei);
      if (phiv != nullVal)
        phi[ccid]->Fill(phiv, wei);
      rnk[ccid]->Fill(rnkv, wei);

      //exclude e-only cands (only data)
      wei = (type == 4) ? 0. : 1.;
      if (etav != nullVal)
        etaData[ccid]->Fill(etav, wei);
      if (phiv != nullVal)
        phiData[ccid]->Fill(phiv, wei);
      rnkData[ccid]->Fill(rnkv, wei);

      // GCT trigger bits
      unsigned int word[2];
      it->data(word);
      std::bitset<32> dbits(word[0]);
      std::bitset<32> ebits(word[1]);
      unsigned int dexor = ((word[0]) ^ (word[1]));
      //disagreeing bits
      std::bitset<32> debits(dexor);
      //disagreeing bits after masking
      std::bitset<32> dembits(((dexor) & (mask)));

      if (verbose())
        std::cout << "l1degct"
                  << " sid:" << sid << " cid:" << cid << "\n"
                  << " data:0x" << std::hex << word[0] << std::dec << " bitset:" << dbits << "\n"
                  << " emul:0x" << std::hex << word[1] << std::dec << " bitset:" << ebits << "\n"
                  << "  xor:0x" << std::hex << dexor << std::dec << " bitset:" << debits
                  << " bitset:" << ((dbits) ^ (ebits)) << "\n"
                  << std::flush;

      for (int ibit = 0; ibit < 32; ibit++) {
        //comparison gives no info if there's only 1 candidate
        wei = (type == 3 || type == 4) ? 0. : 1.;
        if (dbits[ibit])
          dword[ccid]->Fill(ibit, wei);
        if (ebits[ibit])
          eword[ccid]->Fill(ibit, wei);
        if (debits[ibit])
          deword[ccid]->Fill(ibit, wei);
        //if(dembits[ibit])masked[sid]->Fill(ibit,wei);
      }
    }

    //error rates per GCT trigger object type
    int hasCol[nGctColl_] = {0};
    int nagree[nGctColl_] = {0};
    for (L1DEDigiCollection::const_iterator it = gctColl.begin(); it != gctColl.end(); it++) {
      int ccid = it->cid() - dedefs::GCTisolaem;
      ccid = (ccid < 0 || ccid >= nGctColl_) ? 0 : ccid;
      hasCol[ccid]++;
      if (!it->type())
        nagree[ccid]++;
    }
    for (int i = 0; i < nGctColl_; i++) {
      if (!hasCol[i])
        continue;
      //nWithCol[i]++;
      //if(nagree[i]<hasCol[i]) colCount[i]++;
      nWithCol[i] += hasCol[i];  //#of objects
      colCount[i] += nagree[i];  //#of agreements
    }
    for (int i = 0; i < nGctColl_; i++) {
      int ibin = i + 1;
      double rate = nWithCol[i] ? 1. - 1. * colCount[i] / nWithCol[i] : 0.;
      sysrates->setBinContent(ibin, rate);
      if (verbose()) {
        std::cout << "[L1TDEMON] analyze rate computation\t\n"
                  << " colid:" << i << "(so far)"
                  << " nWithCol: " << nWithCol[i] << " colCount: " << colCount[i] << "(this event)"
                  << "hasCol: " << hasCol[i] << " nagree: " << nagree[i] << " rate:" << sysrates->getBinContent(ibin)
                  << "\n"
                  << std::flush;
        if (rate > 1. || rate < 0.)
          std::cout << "problem, error rate for " << SystLabel[i] << " is " << sysrates->getBinContent(ibin) << "\n"
                    << std::flush;
      }
    }
  }

  if (m_stage1_layer2_ == true) {
    bool isComp = deRecord->get_isComp(GCT);
    if (!isComp) {
      if (verbose())
        std::cout << "[L1TdeGCT] Gct information not generated in de-record."
                  << " Skiping event!\n"
                  << std::flush;
      return;
    }

    int DEncand[2] = {0};
    for (int j = 0; j < 2; j++)
      DEncand[j] = deRecord->getNCand(GCT, j);

    if (verbose())
      std::cout << "[L1TdeGCT] ncands"
                << " data: " << DEncand[0] << " emul: " << DEncand[1] << std::endl;

    L1DEDigiCollection deColl;
    deColl = deRecord->getColl();

    // extract the GCT comparison digis
    L1DEDigiCollection stage1layer2Coll;

    stage1layer2Coll.reserve(21);
    stage1layer2Coll.clear();

    for (L1DEDigiCollection::const_iterator it = deColl.begin(); it != deColl.end(); it++)
      if (!it->empty())
        if (it->sid() == GCT)
          stage1layer2Coll.push_back(*it);

    if (verbose()) {
      std::cout << "[L1TdeSTAGE1LAYER2] record has " << stage1layer2Coll.size() << " stage1layer2 de digis\n"
                << std::endl;
      for (L1DEDigiCollection::const_iterator it = stage1layer2Coll.begin(); it != stage1layer2Coll.end(); it++)
        std::cout << "\t" << *it << std::endl;
    }

    const int nullVal = L1DataEmulDigi().reset();


    // d|e candidate loop
    for (L1DEDigiCollection::const_iterator it = stage1layer2Coll.begin(); it != stage1layer2Coll.end(); it++) {
      // sid should be GCT
      int sid = it->sid();

      int cid = it->cid();

      if (verbose())
        std::cout << "[L1TdeStage1Layer2] processing digi "
                  << " sys:" << sid << " type:" << cid << " \n\t" << *it << "\n"
                  << std::endl;

      if (sid != GCT || it->empty()) {
        LogDebug("L1TdeGCT") << "consistency check failure, non-stage1layer2 digis!";
        continue;
      }

      int type = it->type();
      double phiv = it->x1();
      double etav = it->x2();
      float rankarr[2];
      it->rank(rankarr);
      float rnkv = rankarr[0];

      double wei;

      unsigned int mask = (~0x0);

      // shift coll type for starting at zero
      int ccid = cid - dedefs::GCTisolaem;
      if (ccid < 0 || ccid >= nStage1Layer2Coll_) {
        LogDebug("L1TdeGCT") << "consistency check failure, col type outbounds:" << ccid << "\n";
        ccid = 0;
      }

      //type: 0:agree 1:loc.agree, 2:loc.disagree, 3:data.only, 4:emul.only
      if (type < 4)
        sysncand[0]->Fill(ccid);
      if (type < 5 && type != 3)
        sysncand[1]->Fill(ccid);
      errortype_stage1layer2[ccid]->Fill(type);
      wei = (type == 0) ? 0. : 1.;
      if (etav != nullVal && phiv != nullVal)
        etaphi_stage1layer2[ccid]->Fill(etav, phiv, wei);
      if (etav != nullVal)
        eta_stage1layer2[ccid]->Fill(etav, wei);
      if (phiv != nullVal)
        phi_stage1layer2[ccid]->Fill(phiv, wei);
      rnk_stage1layer2[ccid]->Fill(rnkv, wei);

      //exclude e-only cands (only data)
      wei = (type == 4) ? 0. : 1.;
      if (etav != nullVal)
        etaData_stage1layer2[ccid]->Fill(etav, wei);
      if (phiv != nullVal)
        phiData_stage1layer2[ccid]->Fill(phiv, wei);
      rnkData_stage1layer2[ccid]->Fill(rnkv, wei);

      // GCT trigger bits
      unsigned int word_stage1layer2[2];
      it->data(word_stage1layer2);
      std::bitset<32> dbits(word_stage1layer2[0]);
      std::bitset<32> ebits(word_stage1layer2[1]);
      unsigned int dexor = ((word_stage1layer2[0]) ^ (word_stage1layer2[1]));
      //disagreeing bits
      std::bitset<32> debits(dexor);
      //disagreeing bits after masking
      std::bitset<32> dembits(((dexor) & (mask)));

      if (verbose())
        std::cout << "l1degct"
                  << " sid:" << sid << " cid:" << cid << "\n"
                  << " data:0x" << std::hex << word_stage1layer2[0] << std::dec << " bitset:" << dbits << "\n"
                  << " emul:0x" << std::hex << word_stage1layer2[1] << std::dec << " bitset:" << ebits << "\n"
                  << "  xor:0x" << std::hex << dexor << std::dec << " bitset:" << debits
                  << " bitset:" << ((dbits) ^ (ebits)) << "\n"
                  << std::flush;

      for (int ibit = 0; ibit < 32; ibit++) {
        //comparison gives no info if there's only 1 candidate
        wei = (type == 3 || type == 4) ? 0. : 1.;
        if (dbits[ibit])
          dword_stage1layer2[ccid]->Fill(ibit, wei);
        if (ebits[ibit])
          eword_stage1layer2[ccid]->Fill(ibit, wei);
        if (debits[ibit])
          deword_stage1layer2[ccid]->Fill(ibit, wei);
        //if(dembits[ibit])masked[sid]->Fill(ibit,wei);
      }
    }
    //error rates per GCT trigger object type
    int hasCol[nStage1Layer2Coll_] = {0};
    int nagree[nStage1Layer2Coll_] = {0};
    for (L1DEDigiCollection::const_iterator it = stage1layer2Coll.begin(); it != stage1layer2Coll.end(); it++) {
      int ccid = it->cid() - dedefs::GCTisolaem;
      ccid = (ccid < 0 || ccid >= nStage1Layer2Coll_) ? 0 : ccid;
      hasCol[ccid]++;
      if (!it->type())
        nagree[ccid]++;
    }
    for (int i = 0; i < nStage1Layer2Coll_; i++) {
      if (!hasCol[i])
        continue;
      //nWithCol[i]++;
      //if(nagree[i]<hasCol[i]) colCount[i]++;
      nWithCol_stage1Layer2[i] += hasCol[i];  //#of objects
      colCount_stage1Layer2[i] += nagree[i];  //#of agreements
    }
    for (int i = 0; i < nStage1Layer2Coll_; i++) {
      int ibin = i + 1;
      double rate = nWithCol_stage1Layer2[i] ? 1. - 1. * colCount_stage1Layer2[i] / nWithCol_stage1Layer2[i] : 0.;
      sysrates->setBinContent(ibin, rate);
      if (verbose()) {
        std::cout << "[L1TDEMON] analyze rate computation\t\n"
                  << " colid:" << i << "(so far)"
                  << " nWithCol: " << nWithCol_stage1Layer2[i] << " colCount: " << colCount_stage1Layer2[i]
                  << "(this event)"
                  << "hasCol: " << hasCol[i] << " nagree: " << nagree[i] << " rate:" << sysrates->getBinContent(ibin)
                  << "\n"
                  << std::flush;
        if (rate > 1. || rate < 0.)
          std::cout << "problem, error rate for " << SystLabel[i] << " is " << sysrates->getBinContent(ibin) << "\n"
                    << std::flush;
      }
    }
  }
}

void L1TdeGCT::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 37 of file L1TdeGCT.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), gather_cfg::cout, PVValHelper::eta, fastmath::etaphi(), mps_fire::i, dqmiolumiharvest::j, dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, and verbose.

                                                                                           {
  int rnkNBins = 63;
  double rnkMinim = 0.5;
  double rnkMaxim = 63.5;

  if (verbose())
    std::cout << "L1TdeGCT::beginRun()  start\n" << std::flush;

  ibooker.setCurrentFolder(histFolder_);
  //for Legacy GCT
  if (m_stage1_layer2_ == false) {
    sysrates = ibooker.book1D("sysrates", "RATE OF COMPARISON FAILURES", nGctColl_, 0, nGctColl_);

    for (int j = 0; j < 2; j++) {
      std::string lbl("sysncand");
      lbl += (j == 0 ? "Data" : "Emul");
      std::string title("GCT OBJECT MULTIPLICITY ");
      title += (j == 0 ? "(DATA)" : "(EMULATOR)");
      sysncand[j] = ibooker.book1D(lbl.data(), title.data(), nGctColl_, 0, nGctColl_);
    }

    for (int j = 0; j < nGctColl_; j++) {
      ibooker.setCurrentFolder(std::string(histFolder_ + "/" + cLabel[j]));

      std::string lbl("");
      std::string title("");
      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "ErrorFlag";
      title += cLabel[j];
      title += " ErrorFlag";
      errortype[j] = ibooker.book1D(lbl.data(), title.data(), nerr, 0, nerr);

      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "Eta";
      title += cLabel[j];
      title += " ETA OF COMPARISON FAILURES";
      eta[j] = ibooker.book1D(lbl.data(), title.data(), etaNBins, etaMinim, etaMaxim);
      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "Phi";
      title += cLabel[j];
      title += " PHI OF COMPARISON FAILURES";
      phi[j] = ibooker.book1D(lbl.data(), title.data(), phiNBins, phiMinim, phiMaxim);

      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "Etaphi";
      title += cLabel[j];
      title += " ETA PHI OF COMPARISON FAILURES";
      etaphi[j] = ibooker.book2D(lbl.data(), title.data(), etaNBins, etaMinim, etaMaxim, phiNBins, phiMinim, phiMaxim);
      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "Rank";
      title += cLabel[j];
      title += " RANK OF COMPARISON FAILURES";
      rnk[j] = ibooker.book1D(lbl.data(), title.data(), rnkNBins, rnkMinim, rnkMaxim);
      //
      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "Eta";
      lbl += "Data";
      title += cLabel[j];
      title += " ETA (DATA)";
      etaData[j] = ibooker.book1D(lbl.data(), title.data(), etaNBins, etaMinim, etaMaxim);
      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "Phi";
      lbl += "Data";
      title += cLabel[j];
      title += " PHI (DATA)";
      phiData[j] = ibooker.book1D(lbl.data(), title.data(), phiNBins, phiMinim, phiMaxim);

      lbl.clear();
      title.clear();
      lbl += cLabel[j];
      lbl += "Rank";
      lbl += "Data";
      title += cLabel[j];
      title += " RANK (DATA)";
      rnkData[j] = ibooker.book1D(lbl.data(), title.data(), rnkNBins, rnkMinim, rnkMaxim);
      lbl.clear();
      lbl += cLabel[j];
      lbl += "Dword";
      dword[j] = ibooker.book1D(lbl.data(), lbl.data(), nbit, 0, nbit);
      lbl.clear();
      lbl += cLabel[j];
      lbl += "Eword";
      eword[j] = ibooker.book1D(lbl.data(), lbl.data(), nbit, 0, nbit);
      lbl.clear();
      lbl += cLabel[j];
      lbl += "DEword";
      deword[j] = ibooker.book1D(lbl.data(), lbl.data(), nbit, 0, nbit);
      //lbl.clear();
      //lbl+=cLabel[j];lbl+="Masked";
      //masked[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
    }

    for (int i = 0; i < nGctColl_; i++) {
      sysrates->setBinLabel(i + 1, cLabel[i]);
      sysncand[0]->setBinLabel(i + 1, cLabel[i]);
      sysncand[1]->setBinLabel(i + 1, cLabel[i]);
    }

    for (int i = 0; i < nGctColl_; i++) {
      for (int j = 0; j < nerr; j++) {
        errortype[i]->setBinLabel(j + 1, errLabel[j]);
      }
    }

    for (int i = 0; i < nGctColl_; i++) {
      etaphi[i]->setAxisTitle("GCT #eta", 1);
      etaphi[i]->setAxisTitle("GCT #phi", 2);
      eta[i]->setAxisTitle("GCT #eta");
      phi[i]->setAxisTitle("GCT #phi");
      rnk[i]->setAxisTitle("Rank");
      etaData[i]->setAxisTitle("GCT #eta");
      phiData[i]->setAxisTitle("GCT #phi");
      rnkData[i]->setAxisTitle("Rank");
      dword[i]->setAxisTitle("trigger data word bit");
      eword[i]->setAxisTitle("trigger data word bit");
      deword[i]->setAxisTitle("trigger data word bit");
    }

    for (int i = 0; i < nGctColl_; i++) {
      colCount[i] = 0;
      nWithCol[i] = 0;
    }
  }

  // for stage 1 layer 2

  if (m_stage1_layer2_ == true) {
    sysrates = ibooker.book1D("sysrates", "RATE OF COMPARISON FAILURES", nStage1Layer2Coll_, 0, nStage1Layer2Coll_);

    for (int j = 0; j < 2; j++) {
      std::string lbl("sysncand");
      lbl += (j == 0 ? "Data" : "Emul");
      std::string title("Stage1Layer2 OBJECT MULTIPLICITY ");
      title += (j == 0 ? "(DATA)" : "(EMULATOR)");
      sysncand[j] = ibooker.book1D(lbl.data(), title.data(), nStage1Layer2Coll_, 0, nStage1Layer2Coll_);
    }

    for (int j = 0; j < nStage1Layer2Coll_; j++) {
      ibooker.setCurrentFolder(std::string(histFolder_ + "/" + sLabel[j]));

      if (sLabel[j] == "MHT") {
        rnkNBins = 127;
        rnkMinim = 0.5;
        rnkMaxim = 127.5;
      }

      if ((sLabel[j] == "HT") || (sLabel[j] == "ET") || (sLabel[j] == "MET")) {
        rnkNBins = 4096;
        rnkMinim = -0.5;
        rnkMaxim = 4095.5;
      }

      if (sLabel[j] == "Stage1HFSums") {
        rnkNBins = 9;
        rnkMinim = -0.5;
        rnkMaxim = 8.5;
      }

      if (sLabel[j] == "IsoTauJet") {
        rnkNBins = 64;
        rnkMinim = -0.5;
        rnkMaxim = 63.5;
      }

      std::string lbl("");
      std::string title("");
      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "ErrorFlag";
      title += sLabel[j];
      title += " ErrorFlag";
      errortype_stage1layer2[j] = ibooker.book1D(lbl.data(), title.data(), nerr, 0, nerr);

      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "Eta";
      title += sLabel[j];
      title += " ETA OF COMPARISON FAILURES";
      eta_stage1layer2[j] = ibooker.book1D(lbl.data(), title.data(), etaNBins, etaMinim, etaMaxim);
      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "Phi";
      title += sLabel[j];
      title += " PHI OF COMPARISON FAILURES";
      phi_stage1layer2[j] = ibooker.book1D(lbl.data(), title.data(), phiNBins, phiMinim, phiMaxim);

      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "Etaphi";
      title += sLabel[j];
      title += " ETA PHI OF COMPARISON FAILURES";
      etaphi_stage1layer2[j] =
          ibooker.book2D(lbl.data(), title.data(), etaNBins, etaMinim, etaMaxim, phiNBins, phiMinim, phiMaxim);
      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "Rank";
      title += sLabel[j];
      title += " RANK OF COMPARISON FAILURES";
      rnk_stage1layer2[j] = ibooker.book1D(lbl.data(), title.data(), rnkNBins, rnkMinim, rnkMaxim);
      //
      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "Eta";
      lbl += "Data";
      title += sLabel[j];
      title += " ETA (DATA)";
      etaData_stage1layer2[j] = ibooker.book1D(lbl.data(), title.data(), etaNBins, etaMinim, etaMaxim);
      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "Phi";
      lbl += "Data";
      title += sLabel[j];
      title += " PHI (DATA)";
      phiData_stage1layer2[j] = ibooker.book1D(lbl.data(), title.data(), phiNBins, phiMinim, phiMaxim);

      lbl.clear();
      title.clear();
      lbl += sLabel[j];
      lbl += "Rank";
      lbl += "Data";
      title += sLabel[j];
      title += " RANK (DATA)";
      rnkData_stage1layer2[j] = ibooker.book1D(lbl.data(), title.data(), rnkNBins, rnkMinim, rnkMaxim);
      lbl.clear();
      lbl += sLabel[j];
      lbl += "Dword";
      dword_stage1layer2[j] = ibooker.book1D(lbl.data(), lbl.data(), nbit, 0, nbit);
      lbl.clear();
      lbl += sLabel[j];
      lbl += "Eword";
      eword_stage1layer2[j] = ibooker.book1D(lbl.data(), lbl.data(), nbit, 0, nbit);
      lbl.clear();
      lbl += sLabel[j];
      lbl += "DEword";
      deword_stage1layer2[j] = ibooker.book1D(lbl.data(), lbl.data(), nbit, 0, nbit);
      //lbl.clear();
      //lbl+=cLabel[j];lbl+="Masked";
      //masked_stage1layer2[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
    }

    for (int i = 0; i < nStage1Layer2Coll_; i++) {
      sysrates->setBinLabel(i + 1, sLabel[i]);
      sysncand[0]->setBinLabel(i + 1, sLabel[i]);
      sysncand[1]->setBinLabel(i + 1, sLabel[i]);
    }

    for (int i = 0; i < nStage1Layer2Coll_; i++) {
      for (int j = 0; j < nerr; j++) {
        errortype_stage1layer2[i]->setBinLabel(j + 1, errLabel[j]);
      }
    }

    for (int i = 0; i < nStage1Layer2Coll_; i++) {
      etaphi_stage1layer2[i]->setAxisTitle("Stage1Layer2 #eta", 1);
      etaphi_stage1layer2[i]->setAxisTitle("Stage1Layer2 #phi", 2);
      eta_stage1layer2[i]->setAxisTitle("Stage1Layer2 #eta");
      phi_stage1layer2[i]->setAxisTitle("Stage1Layer2 #phi");
      rnk_stage1layer2[i]->setAxisTitle("Rank");
      etaData_stage1layer2[i]->setAxisTitle("Stage1Layer2 #eta");
      phiData_stage1layer2[i]->setAxisTitle("Stage1Layer2 #phi");
      rnkData_stage1layer2[i]->setAxisTitle("Rank");
      dword_stage1layer2[i]->setAxisTitle("trigger data word bit");
      eword_stage1layer2[i]->setAxisTitle("trigger data word bit");
      deword_stage1layer2[i]->setAxisTitle("trigger data word bit");
    }

    for (int i = 0; i < nStage1Layer2Coll_; i++) {
      colCount_stage1Layer2[i] = 0;
      nWithCol_stage1Layer2[i] = 0;
    }
  }

  if (verbose())
    std::cout << "L1TdeGCT::beginJob()  end.\n" << std::flush;
}

int L1TdeGCT::verbose ( )

inlineprivate

Definition at line 44 of file L1TdeGCT.h.

References verbose_.

{ return verbose_; }

Member Data Documentation

std::string L1TdeGCT::cLabel[nGctColl_]

private

Initial value:

= {
      "IsoEM", "NoisoEM", "CenJet", "ForJet", "TauJet", "HT", "MET", "ET", "MHT", "HFSums", "HFCnts"}

Definition at line 77 of file L1TdeGCT.h.

int L1TdeGCT::colCount[nGctColl_]

private

Definition at line 62 of file L1TdeGCT.h.

int L1TdeGCT::colCount_stage1Layer2[nStage1Layer2Coll_]

private

Definition at line 65 of file L1TdeGCT.h.

edm::EDGetTokenT<L1DataEmulRecord> L1TdeGCT::DEsource_

private

Definition at line 39 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::deword[nGctColl_]

private

Definition at line 111 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::deword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 115 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::dword[nGctColl_]

private

Definition at line 109 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::dword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 113 of file L1TdeGCT.h.

std::string L1TdeGCT::errLabel[nerr] = {"Agree", "Loc. Agree", "L.Disagree", "Data only", "Emul only"}

private

Definition at line 81 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::errortype[nGctColl_]

private

Definition at line 87 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::errortype_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 97 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::eta[nGctColl_]

private

Definition at line 90 of file L1TdeGCT.h.

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

MonitorElement* L1TdeGCT::eta_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 100 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::etaData[nGctColl_]

private

Definition at line 93 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::etaData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 103 of file L1TdeGCT.h.

const double L1TdeGCT::etaMaxim = 21.5

private

Definition at line 74 of file L1TdeGCT.h.

const double L1TdeGCT::etaMinim = -0.5

private

Definition at line 73 of file L1TdeGCT.h.

const int L1TdeGCT::etaNBins = 22

private

Definition at line 72 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::etaphi[nGctColl_]

private

Definition at line 89 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::etaphi_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 99 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::eword[nGctColl_]

private

Definition at line 110 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::eword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 114 of file L1TdeGCT.h.

bool L1TdeGCT::hasRecord_

private

Definition at line 40 of file L1TdeGCT.h.

std::string L1TdeGCT::histFile_

private

Definition at line 47 of file L1TdeGCT.h.

std::string L1TdeGCT::histFolder_

private

Definition at line 50 of file L1TdeGCT.h.

bool L1TdeGCT::m_stage1_layer2_

private

Definition at line 52 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::masked[nGctColl_]

private

Definition at line 112 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::masked_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 116 of file L1TdeGCT.h.

bool L1TdeGCT::monitorDaemon_

private

Definition at line 55 of file L1TdeGCT.h.

const int L1TdeGCT::nbit = 32

private

Definition at line 76 of file L1TdeGCT.h.

const int L1TdeGCT::nerr = 5

staticprivate

Definition at line 75 of file L1TdeGCT.h.

const int L1TdeGCT::nGctColl_ = dedefs::GCThfbit - dedefs::GCTisolaem + 1

staticprivate

Definition at line 58 of file L1TdeGCT.h.

const int L1TdeGCT::nStage1Layer2Coll_ = dedefs::GCTisotaujets - dedefs::GCTisolaem + 1

staticprivate

Definition at line 59 of file L1TdeGCT.h.

int L1TdeGCT::nWithCol[nGctColl_]

private

Definition at line 63 of file L1TdeGCT.h.

int L1TdeGCT::nWithCol_stage1Layer2[nStage1Layer2Coll_]

private

Definition at line 66 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::phi[nGctColl_]

private

Definition at line 91 of file L1TdeGCT.h.

Referenced by Particle.Particle::__str__(), and ntupleDataFormat.Track::phiPull().

MonitorElement* L1TdeGCT::phi_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 101 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::phiData[nGctColl_]

private

Definition at line 94 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::phiData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 104 of file L1TdeGCT.h.

const double L1TdeGCT::phiMaxim = 17.5

private

Definition at line 71 of file L1TdeGCT.h.

const double L1TdeGCT::phiMinim = -0.5

private

Definition at line 70 of file L1TdeGCT.h.

const int L1TdeGCT::phiNBins = 18

private

Definition at line 69 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::rnk[nGctColl_]

private

Definition at line 92 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::rnk_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 102 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::rnkData[nGctColl_]

private

Definition at line 95 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::rnkData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 105 of file L1TdeGCT.h.

std::string L1TdeGCT::sLabel[nStage1Layer2Coll_]

private

Initial value:

= {
      "IsoEM", "NoisoEM", "CenJet", "ForJet", "TauJet", "HT", "MET", "ET", "MHT", "Stage1HFSums", "HFCnts", "IsoTauJet"}

Definition at line 79 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::sysncand[2]

private

Definition at line 85 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::sysrates

private

Definition at line 84 of file L1TdeGCT.h.

int L1TdeGCT::verbose_

private

Definition at line 43 of file L1TdeGCT.h.

Referenced by verbose().