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 26 of file L1TdeGCT.h.

Constructor & Destructor Documentation

L1TdeGCT()

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::~L1TdeGCT ( )

override

Definition at line 35 of file L1TdeGCT.cc.

{}

Member Function Documentation

analyze()

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 gather_cfg::cout, TauDecayModes::dec, PVValHelper::eta, fastmath::etaphi(), dedefs::GCT, dedefs::GCTisolaem, L1DataEmulRecord::get_isComp(), L1DataEmulRecord::getColl(), L1DataEmulRecord::getNCand(), mps_fire::i, iEvent, 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;
      }
    }
  }
}

bookHistograms()

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, metDiagnosticParameterSet_cfi::etaNBins, fastmath::etaphi(), mps_fire::i, dqmiolumiharvest::j, metDiagnosticParameterSet_cfi::phiNBins, 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;
}

verbose()

int L1TdeGCT::verbose ( )

inlineprivate

Definition at line 43 of file L1TdeGCT.h.

References verbose_.

{ return verbose_; }

Member Data Documentation

cLabel

std::string L1TdeGCT::cLabel[nGctColl_]

private

Initial value:

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

Definition at line 76 of file L1TdeGCT.h.

colCount

int L1TdeGCT::colCount[nGctColl_]

private

Definition at line 61 of file L1TdeGCT.h.

colCount_stage1Layer2

int L1TdeGCT::colCount_stage1Layer2[nStage1Layer2Coll_]

private

Definition at line 64 of file L1TdeGCT.h.

DEsource_

edm::EDGetTokenT<L1DataEmulRecord> L1TdeGCT::DEsource_

private

Definition at line 38 of file L1TdeGCT.h.

deword

MonitorElement* L1TdeGCT::deword[nGctColl_]

private

Definition at line 110 of file L1TdeGCT.h.

deword_stage1layer2

MonitorElement* L1TdeGCT::deword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 114 of file L1TdeGCT.h.

dword

MonitorElement* L1TdeGCT::dword[nGctColl_]

private

Definition at line 108 of file L1TdeGCT.h.

dword_stage1layer2

MonitorElement* L1TdeGCT::dword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 112 of file L1TdeGCT.h.

errLabel

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

private

Definition at line 80 of file L1TdeGCT.h.

errortype

MonitorElement* L1TdeGCT::errortype[nGctColl_]

private

Definition at line 86 of file L1TdeGCT.h.

errortype_stage1layer2

MonitorElement* L1TdeGCT::errortype_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 96 of file L1TdeGCT.h.

eta

MonitorElement* L1TdeGCT::eta[nGctColl_]

private

Definition at line 89 of file L1TdeGCT.h.

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

eta_stage1layer2

MonitorElement* L1TdeGCT::eta_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 99 of file L1TdeGCT.h.

etaData

MonitorElement* L1TdeGCT::etaData[nGctColl_]

private

Definition at line 92 of file L1TdeGCT.h.

etaData_stage1layer2

MonitorElement* L1TdeGCT::etaData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 102 of file L1TdeGCT.h.

etaMaxim

const double L1TdeGCT::etaMaxim = 21.5

private

Definition at line 73 of file L1TdeGCT.h.

etaMinim

const double L1TdeGCT::etaMinim = -0.5

private

Definition at line 72 of file L1TdeGCT.h.

etaNBins

const int L1TdeGCT::etaNBins = 22

private

Definition at line 71 of file L1TdeGCT.h.

etaphi

MonitorElement* L1TdeGCT::etaphi[nGctColl_]

private

Definition at line 88 of file L1TdeGCT.h.

etaphi_stage1layer2

MonitorElement* L1TdeGCT::etaphi_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 98 of file L1TdeGCT.h.

eword

MonitorElement* L1TdeGCT::eword[nGctColl_]

private

Definition at line 109 of file L1TdeGCT.h.

eword_stage1layer2

MonitorElement* L1TdeGCT::eword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 113 of file L1TdeGCT.h.

hasRecord_

bool L1TdeGCT::hasRecord_

private

Definition at line 39 of file L1TdeGCT.h.

histFile_

std::string L1TdeGCT::histFile_

private

Definition at line 46 of file L1TdeGCT.h.

histFolder_

std::string L1TdeGCT::histFolder_

private

Definition at line 49 of file L1TdeGCT.h.

m_stage1_layer2_

bool L1TdeGCT::m_stage1_layer2_

private

Definition at line 51 of file L1TdeGCT.h.

masked

MonitorElement* L1TdeGCT::masked[nGctColl_]

private

Definition at line 111 of file L1TdeGCT.h.

masked_stage1layer2

MonitorElement* L1TdeGCT::masked_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 115 of file L1TdeGCT.h.

monitorDaemon_

bool L1TdeGCT::monitorDaemon_

private

Definition at line 54 of file L1TdeGCT.h.

nbit

const int L1TdeGCT::nbit = 32

private

Definition at line 75 of file L1TdeGCT.h.

nerr

const int L1TdeGCT::nerr = 5

staticprivate

Definition at line 74 of file L1TdeGCT.h.

nGctColl_

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

staticprivate

Definition at line 57 of file L1TdeGCT.h.

nStage1Layer2Coll_

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

staticprivate

Definition at line 58 of file L1TdeGCT.h.

nWithCol

int L1TdeGCT::nWithCol[nGctColl_]

private

Definition at line 62 of file L1TdeGCT.h.

nWithCol_stage1Layer2

int L1TdeGCT::nWithCol_stage1Layer2[nStage1Layer2Coll_]

private

Definition at line 65 of file L1TdeGCT.h.

phi

MonitorElement* L1TdeGCT::phi[nGctColl_]

private

Definition at line 90 of file L1TdeGCT.h.

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

phi_stage1layer2

MonitorElement* L1TdeGCT::phi_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 100 of file L1TdeGCT.h.

phiData

MonitorElement* L1TdeGCT::phiData[nGctColl_]

private

Definition at line 93 of file L1TdeGCT.h.

phiData_stage1layer2

MonitorElement* L1TdeGCT::phiData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 103 of file L1TdeGCT.h.

phiMaxim

const double L1TdeGCT::phiMaxim = 17.5

private

Definition at line 70 of file L1TdeGCT.h.

phiMinim

const double L1TdeGCT::phiMinim = -0.5

private

Definition at line 69 of file L1TdeGCT.h.

phiNBins

const int L1TdeGCT::phiNBins = 18

private

Definition at line 68 of file L1TdeGCT.h.

rnk

MonitorElement* L1TdeGCT::rnk[nGctColl_]

private

Definition at line 91 of file L1TdeGCT.h.

rnk_stage1layer2

MonitorElement* L1TdeGCT::rnk_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 101 of file L1TdeGCT.h.

rnkData

MonitorElement* L1TdeGCT::rnkData[nGctColl_]

private

Definition at line 94 of file L1TdeGCT.h.

rnkData_stage1layer2

MonitorElement* L1TdeGCT::rnkData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 104 of file L1TdeGCT.h.

sLabel

std::string L1TdeGCT::sLabel[nStage1Layer2Coll_]

private

Initial value:

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

Definition at line 78 of file L1TdeGCT.h.

sysncand

MonitorElement* L1TdeGCT::sysncand[2]

private

Definition at line 84 of file L1TdeGCT.h.

sysrates

MonitorElement* L1TdeGCT::sysrates

private

Definition at line 83 of file L1TdeGCT.h.

verbose_

int L1TdeGCT::verbose_

private

Definition at line 42 of file L1TdeGCT.h.

Referenced by verbose().