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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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 >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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::L1TdeGCT ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 10 of file L1TdeGCT.cc.

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

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

~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.

                                                                          {
  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 = 1.;
 
      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 (it->type() < 4)
        sysncand[0]->Fill(ccid);
      if (it->type() < 5 && it->type() != 3)
        sysncand[1]->Fill(ccid);
 
      errortype[ccid]->Fill(type);
 
      wei = 1.;
      if (!type)
        wei = 0.;
      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 = 1.;
      if (type == 4)
        wei = 0.;
      if (etav != nullVal)
        etaData[ccid]->Fill(etav, wei);
      if (phiv != nullVal)
        phiData[ccid]->Fill(phiv, wei);
      rnkData[ccid]->Fill(rnkv, wei);
      wei = 1;
 
      // 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++) {
        wei = 1.;
        //comparison gives no info if there's only 1 candidate
        if (type == 3 || type == 4)
          wei = 0.;
        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);
      }
      wei = 1;
    }
 
    //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 = 1.;
 
      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 (it->type() < 4)
        sysncand[0]->Fill(ccid);
      if (it->type() < 5 && it->type() != 3)
        sysncand[1]->Fill(ccid);
      errortype_stage1layer2[ccid]->Fill(type);
      wei = 1.;
      if (!type)
        wei = 0.;
      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 = 1.;
      if (type == 4)
        wei = 0.;
      if (etav != nullVal)
        etaData_stage1layer2[ccid]->Fill(etav, wei);
      if (phiv != nullVal)
        phiData_stage1layer2[ccid]->Fill(phiv, wei);
      rnkData_stage1layer2[ccid]->Fill(rnkv, wei);
      wei = 1;
 
      // 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++) {
        wei = 1.;
        //comparison gives no info if there's only 1 candidate
        if (type == 3 || type == 4)
          wei = 0.;
        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);
      }
      wei = 1;
    }
    //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;
      }
    }
  }
}

References gather_cfg::cout, TauDecayModes::dec, PVValHelper::eta, fastmath::etaphi(), HcalObjRepresent::Fill(), 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.

bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 37 of file L1TdeGCT.cc.

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

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, overlapproblemtsosanalyzer_cfi::title, and verbose.

verbose()

int L1TdeGCT::verbose ( )

inlineprivate

Definition at line 44 of file L1TdeGCT.h.

{ return verbose_; }

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

colCount

int L1TdeGCT::colCount[nGctColl_]

private

Definition at line 62 of file L1TdeGCT.h.

colCount_stage1Layer2

int L1TdeGCT::colCount_stage1Layer2[nStage1Layer2Coll_]

private

Definition at line 65 of file L1TdeGCT.h.

DEsource_

edm::EDGetTokenT<L1DataEmulRecord> L1TdeGCT::DEsource_

private

Definition at line 39 of file L1TdeGCT.h.

deword

MonitorElement* L1TdeGCT::deword[nGctColl_]

private

Definition at line 111 of file L1TdeGCT.h.

deword_stage1layer2

MonitorElement* L1TdeGCT::deword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 115 of file L1TdeGCT.h.

dword

MonitorElement* L1TdeGCT::dword[nGctColl_]

private

Definition at line 109 of file L1TdeGCT.h.

dword_stage1layer2

MonitorElement* L1TdeGCT::dword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 113 of file L1TdeGCT.h.

errLabel

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

private

Definition at line 81 of file L1TdeGCT.h.

errortype

MonitorElement* L1TdeGCT::errortype[nGctColl_]

private

Definition at line 87 of file L1TdeGCT.h.

errortype_stage1layer2

MonitorElement* L1TdeGCT::errortype_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 97 of file L1TdeGCT.h.

eta

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().

eta_stage1layer2

MonitorElement* L1TdeGCT::eta_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 100 of file L1TdeGCT.h.

etaData

MonitorElement* L1TdeGCT::etaData[nGctColl_]

private

Definition at line 93 of file L1TdeGCT.h.

etaData_stage1layer2

MonitorElement* L1TdeGCT::etaData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 103 of file L1TdeGCT.h.

etaMaxim

const double L1TdeGCT::etaMaxim = 21.5

private

Definition at line 74 of file L1TdeGCT.h.

etaMinim

const double L1TdeGCT::etaMinim = -0.5

private

Definition at line 73 of file L1TdeGCT.h.

etaNBins

const int L1TdeGCT::etaNBins = 22

private

Definition at line 72 of file L1TdeGCT.h.

etaphi

MonitorElement* L1TdeGCT::etaphi[nGctColl_]

private

Definition at line 89 of file L1TdeGCT.h.

etaphi_stage1layer2

MonitorElement* L1TdeGCT::etaphi_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 99 of file L1TdeGCT.h.

eword

MonitorElement* L1TdeGCT::eword[nGctColl_]

private

Definition at line 110 of file L1TdeGCT.h.

eword_stage1layer2

MonitorElement* L1TdeGCT::eword_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 114 of file L1TdeGCT.h.

hasRecord_

bool L1TdeGCT::hasRecord_

private

Definition at line 40 of file L1TdeGCT.h.

histFile_

std::string L1TdeGCT::histFile_

private

Definition at line 47 of file L1TdeGCT.h.

histFolder_

std::string L1TdeGCT::histFolder_

private

Definition at line 50 of file L1TdeGCT.h.

m_stage1_layer2_

bool L1TdeGCT::m_stage1_layer2_

private

Definition at line 52 of file L1TdeGCT.h.

masked

MonitorElement* L1TdeGCT::masked[nGctColl_]

private

Definition at line 112 of file L1TdeGCT.h.

masked_stage1layer2

MonitorElement* L1TdeGCT::masked_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 116 of file L1TdeGCT.h.

monitorDaemon_

bool L1TdeGCT::monitorDaemon_

private

Definition at line 55 of file L1TdeGCT.h.

nbit

const int L1TdeGCT::nbit = 32

private

Definition at line 76 of file L1TdeGCT.h.

nerr

const int L1TdeGCT::nerr = 5

staticprivate

Definition at line 75 of file L1TdeGCT.h.

nGctColl_

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

staticprivate

Definition at line 58 of file L1TdeGCT.h.

nStage1Layer2Coll_

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

staticprivate

Definition at line 59 of file L1TdeGCT.h.

nWithCol

int L1TdeGCT::nWithCol[nGctColl_]

private

Definition at line 63 of file L1TdeGCT.h.

nWithCol_stage1Layer2

int L1TdeGCT::nWithCol_stage1Layer2[nStage1Layer2Coll_]

private

Definition at line 66 of file L1TdeGCT.h.

phi

MonitorElement* L1TdeGCT::phi[nGctColl_]

private

Definition at line 91 of file L1TdeGCT.h.

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

phi_stage1layer2

MonitorElement* L1TdeGCT::phi_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 101 of file L1TdeGCT.h.

phiData

MonitorElement* L1TdeGCT::phiData[nGctColl_]

private

Definition at line 94 of file L1TdeGCT.h.

phiData_stage1layer2

MonitorElement* L1TdeGCT::phiData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 104 of file L1TdeGCT.h.

phiMaxim

const double L1TdeGCT::phiMaxim = 17.5

private

Definition at line 71 of file L1TdeGCT.h.

phiMinim

const double L1TdeGCT::phiMinim = -0.5

private

Definition at line 70 of file L1TdeGCT.h.

phiNBins

const int L1TdeGCT::phiNBins = 18

private

Definition at line 69 of file L1TdeGCT.h.

rnk

MonitorElement* L1TdeGCT::rnk[nGctColl_]

private

Definition at line 92 of file L1TdeGCT.h.

rnk_stage1layer2

MonitorElement* L1TdeGCT::rnk_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 102 of file L1TdeGCT.h.

rnkData

MonitorElement* L1TdeGCT::rnkData[nGctColl_]

private

Definition at line 95 of file L1TdeGCT.h.

rnkData_stage1layer2

MonitorElement* L1TdeGCT::rnkData_stage1layer2[nStage1Layer2Coll_]

private

Definition at line 105 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 79 of file L1TdeGCT.h.

sysncand

MonitorElement* L1TdeGCT::sysncand[2]

private

Definition at line 85 of file L1TdeGCT.h.

sysrates

MonitorElement* L1TdeGCT::sysrates

private

Definition at line 84 of file L1TdeGCT.h.

verbose_

int L1TdeGCT::verbose_

private

Definition at line 43 of file L1TdeGCT.h.

Referenced by verbose().