L1TGT Class Reference

#include <L1TGT.h>

Inheritance diagram for L1TGT:

Public Member Functions
	L1TGT (const edm::ParameterSet &ps)
	~L1TGT () override
Public Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
	DQMOneEDAnalyzer ()
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args...>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
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
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) final
	end section More...
void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
void	endLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override
Protected Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Private Types
enum	activeDAQ {   FDL = 0, PSB9, PSB13, PSB14,   PSB15, PSB19, PSB20, PSB21,   GMT }
enum	activeEVM { TCS, FDLEVM }
typedef std::vector< std::pair < int, int > >::const_iterator	CItVecPair

Private Member Functions
void	countPfsIndicesPerLs ()
bool	isActive (int word, int bit)
	book all histograms for the module More...

Private Attributes
MonitorElement *	algo_bits
MonitorElement *	algo_bits_corr
MonitorElement *	algo_bits_lumi
MonitorElement *	algo_tt_bits_corr
std::string	algoBitToName [128]
MonitorElement *	BST_beamMode
MonitorElement *	BST_beamMomentum
MonitorElement *	BST_intensityBeam1
MonitorElement *	BST_intensityBeam2
MonitorElement *	BST_lhcFillNumber
MonitorElement *	BST_MasterStatus
MonitorElement *	BST_turnCountNumber
MonitorElement *	dbx_module
MonitorElement *	event_lumi
MonitorElement *	event_number
MonitorElement *	event_type
MonitorElement *	evnum_trignum_lumi
MonitorElement *	gpsfreq
MonitorElement *	gpsfreqlum
MonitorElement *	gpsfreqwide
edm::EDGetTokenT < L1GlobalTriggerEvmReadoutRecord >	gtEvmSource_
	input tag for L1 GT EVM readout record More...
MonitorElement *	gtfe_bx
edm::EDGetTokenT < L1GlobalTriggerReadoutRecord >	gtSource_L1GT_
	input parameters More...
edm::EDGetTokenT < L1MuGMTReadoutCollection >	gtSource_L1MuGMT_
MonitorElement *	h_L1AlgoBX1
MonitorElement *	h_L1AlgoBX2
MonitorElement *	h_L1AlgoBX3
MonitorElement *	h_L1AlgoBX4
MonitorElement *	h_L1TechBX
edm::ESGetToken < L1GtTriggerMenu, L1GtTriggerMenuRcd >	l1gtTrigmenuToken_
std::map< std::string, bool >	l1TechTriggerDecision
std::map< std::string, bool >	l1TriggerDecision
std::string	m_histFolder
	histogram folder for L1 GT plots More...
MonitorElement *	m_monL1PfIndicesPerLs
MonitorElement *	m_monL1PrescaleFactorSet
MonitorElement *	m_monLsNrDiffTcsFdlEvm
MonitorElement *	m_monLsNrDiffTcsFdlEvmLs
MonitorElement *	m_monOrbitNrDiffTcsFdlEvm
MonitorElement *	m_monOrbitNrDiffTcsFdlEvmLs
int	m_nrEvJob
	number of events processed More...
int	m_nrEvRun
std::vector< std::pair< int, int > >	m_pairLsNumberPfIndex
bool	m_runInEndLumi
bool	m_runInEventLoop
	switches to choose the running of various methods More...
MonitorElement *	orbit_lumi
uint64_t	preGps_
uint64_t	preOrb_
MonitorElement *	setupversion_lumi
std::string	techBitToName [64]
std::map< std::string, bool > ::iterator	trig_iter
MonitorElement *	trigger_lumi
MonitorElement *	trigger_number
MonitorElement *	tt_bits
MonitorElement *	tt_bits_corr
MonitorElement *	tt_bits_lumi
bool	verbose_
	verbosity switch More...

Static Private Attributes
static const int	MaxLsNrDiffTcsFdlEvm = 24
static const int	MaxOrbitNrDiffTcsFdlEvm = 24

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Description: DQM for L1 Global Trigger.

Author

J. Berryhill, I. Mikulec

Vasile Mihai Ghete - HEPHY Vienna

Definition at line 49 of file L1TGT.h.

Member Typedef Documentation

typedef std::vector<std::pair<int, int> >::const_iterator L1TGT::CItVecPair

private

Definition at line 184 of file L1TGT.h.

Member Enumeration Documentation

enum L1TGT::activeDAQ

private

Enumerator
FDL
PSB9
PSB13
PSB14
PSB15
PSB19
PSB20
PSB21
GMT

Definition at line 81 of file L1TGT.h.

{ FDL = 0, PSB9, PSB13, PSB14, PSB15, PSB19, PSB20, PSB21, GMT };

enum L1TGT::activeEVM

private

Enumerator
TCS
FDLEVM

Definition at line 85 of file L1TGT.h.

{ TCS, FDLEVM };

Constructor & Destructor Documentation

L1TGT::L1TGT

(

const edm::ParameterSet &

ps

)

Definition at line 18 of file L1TGT.cc.

References edm::ParameterSet::getUntrackedParameter(), l1gtTrigmenuToken_, m_histFolder, and AlCaHLTBitMon_QueryRunRegistry::string.

    : gtSource_L1GT_(consumes<L1GlobalTriggerReadoutRecord>(ps.getParameter<edm::InputTag>("gtSource"))),
      gtSource_L1MuGMT_(consumes<L1MuGMTReadoutCollection>(ps.getParameter<edm::InputTag>("gtSource"))),
      gtEvmSource_(consumes<L1GlobalTriggerEvmReadoutRecord>(ps.getParameter<edm::InputTag>("gtEvmSource"))),
      m_runInEventLoop(ps.getUntrackedParameter<bool>("runInEventLoop", false)),
      m_runInEndLumi(ps.getUntrackedParameter<bool>("runInEndLumi", false)),
      verbose_(ps.getUntrackedParameter<bool>("verbose", false)),
      m_nrEvJob(0),
      m_nrEvRun(0),
      preGps_(0ULL),
      preOrb_(0ULL) {
  m_histFolder = ps.getUntrackedParameter<std::string>("HistFolder", "L1T/L1TGT");
  l1gtTrigmenuToken_ = esConsumes<edm::Transition::BeginRun>();
}

L1TGT::~L1TGT ( )

override

Definition at line 33 of file L1TGT.cc.

              {
  // empty
}

Member Function Documentation

void L1TGT::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  evSetup  )

overrideprotectedvirtual

get Global Trigger algo and technical trigger bit statistics

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 331 of file L1TGT.cc.

References accept(), L1GtfeWord::activeBoards(), algo_bits, algo_bits_corr, algo_bits_lumi, algo_tt_bits_corr, algoBitToName, EcalCondDBWriter_cfi::beam, L1GtfeExtWord::beamMode(), L1GtfeExtWord::beamMomentum(), BST_beamMode, BST_beamMomentum, BST_intensityBeam1, BST_intensityBeam2, BST_lhcFillNumber, BST_MasterStatus, BST_turnCountNumber, L1GtfeExtWord::bstMasterStatus(), L1TcsWord::bxNr(), L1GtPsbWord::bxNr(), L1GtFdlWord::bxNr(), L1GtfeWord::bxNr(), dbx_module, TauDecayModes::dec, event_lumi, event_number, event_type, L1TcsWord::eventNr(), evnum_trignum_lumi, FDL, dqm::impl::MonitorElement::Fill(), spr::find(), edm::Event::getByToken(), GMT, gpsfreq, gpsfreqlum, gpsfreqwide, L1GtfeExtWord::gpsTime(), gtEvmSource_, gtfe_bx, L1GtFdlWord::gtPrescaleFactorIndexAlgo(), gtSource_L1GT_, gtSource_L1MuGMT_, h_L1AlgoBX1, h_L1AlgoBX2, h_L1AlgoBX3, h_L1AlgoBX4, h_L1TechBX, mps_fire::i, isActive(), edm::HandleBase::isValid(), dqmiolumiharvest::j, l1TechTriggerDecision, l1TriggerDecision, L1GtfeExtWord::lhcFillNumber(), L1GtPsbWord::localBxNr(), L1GtFdlWord::localBxNr(), edm::EventBase::luminosityBlock(), L1TcsWord::luminositySegmentNr(), L1GtFdlWord::lumiSegmentNr(), m_monL1PrescaleFactorSet, m_monLsNrDiffTcsFdlEvm, m_monLsNrDiffTcsFdlEvmLs, m_monOrbitNrDiffTcsFdlEvm, m_monOrbitNrDiffTcsFdlEvmLs, m_nrEvJob, m_pairLsNumberPfIndex, funct::master(), MaxLsNrDiffTcsFdlEvm, MaxOrbitNrDiffTcsFdlEvm, orbit_lumi, L1TcsWord::orbitNr(), L1GtFdlWord::orbitNr(), edm::EventBase::orbitNumber(), L1TcsWord::partTrigNr(), preGps_, preOrb_, PSB9, L1GtfeWord::setupVersion(), setupversion_lumi, TCS, techBitToName, L1GtfeExtWord::totalIntensityBeam1(), L1GtfeExtWord::totalIntensityBeam2(), trig_iter, trigger_lumi, trigger_number, L1TcsWord::triggerType(), tt_bits, tt_bits_corr, tt_bits_lumi, L1GtfeExtWord::turnCountNumber(), and verbose_.

                                                                        {
  m_nrEvJob++;

  if (verbose_) {
    edm::LogInfo("L1TGT") << "L1TGT: analyze...." << std::endl;
  }

  // initialize Bx, orbit number, luminosity segment number to invalid value
  int tcsBx = -1;
  int gtfeEvmBx = -1;

  long long int orbitTcs = -1;
  int orbitEvmFdl = -1;

  int lsTcs = -1;
  int lsEvmFdl = -1;

  // get once only the LS block number, to be used in many histograms
  const int lsNumber = iEvent.luminosityBlock();

  // open EVM readout record if available
  edm::Handle<L1GlobalTriggerEvmReadoutRecord> gtEvmReadoutRecord;
  iEvent.getByToken(gtEvmSource_, gtEvmReadoutRecord);

  if (!gtEvmReadoutRecord.isValid()) {
    edm::LogInfo("L1TGT") << "can't find L1GlobalTriggerEvmReadoutRecord";
  } else {
    // get all info from the EVM record if available and fill the histograms

    const L1GtfeWord& gtfeEvmWord = gtEvmReadoutRecord->gtfeWord();
    const L1GtfeExtWord& gtfeEvmExtWord = gtEvmReadoutRecord->gtfeWord();

    gtfeEvmBx = gtfeEvmWord.bxNr();
    int gtfeEvmActiveBoards = gtfeEvmWord.activeBoards();

    if (isActive(gtfeEvmActiveBoards, TCS)) {
      // if TCS present in the record

      const L1TcsWord& tcsWord = gtEvmReadoutRecord->tcsWord();

      tcsBx = tcsWord.bxNr();
      orbitTcs = tcsWord.orbitNr();
      lsTcs = tcsWord.luminositySegmentNr();

      event_type->Fill(tcsWord.triggerType());
      orbit_lumi->Fill(lsNumber, orbitTcs);

      trigger_number->Fill(tcsWord.partTrigNr());
      event_number->Fill(tcsWord.eventNr());

      trigger_lumi->Fill(lsNumber, tcsWord.partTrigNr());
      event_lumi->Fill(lsNumber, tcsWord.eventNr());
      evnum_trignum_lumi->Fill(lsNumber,
                               static_cast<double>(tcsWord.eventNr()) / static_cast<double>(tcsWord.partTrigNr()));

      uint16_t master = gtfeEvmExtWord.bstMasterStatus();
      uint32_t turnCount = gtfeEvmExtWord.turnCountNumber();
      uint32_t lhcFill = gtfeEvmExtWord.lhcFillNumber();
      uint16_t beam = gtfeEvmExtWord.beamMode();
      uint16_t momentum = gtfeEvmExtWord.beamMomentum();
      uint32_t intensity1 = gtfeEvmExtWord.totalIntensityBeam1();
      uint32_t intensity2 = gtfeEvmExtWord.totalIntensityBeam2();

      BST_MasterStatus->Fill(lsNumber, static_cast<double>(master));
      BST_turnCountNumber->Fill(lsNumber, static_cast<double>(turnCount));
      BST_lhcFillNumber->Fill(static_cast<double>(lhcFill % 1000));
      BST_beamMode->Fill(lsNumber, static_cast<double>(beam));

      BST_beamMomentum->Fill(lsNumber, static_cast<double>(momentum));
      BST_intensityBeam1->Fill(lsNumber, static_cast<double>(intensity1));
      BST_intensityBeam2->Fill(lsNumber, static_cast<double>(intensity2));

      if (verbose_) {
        edm::LogInfo("L1TGT") << " check mode = " << beam << "    momentum " << momentum << " int2 " << intensity2
                              << std::endl;
      }

      uint64_t gpsr = gtfeEvmExtWord.gpsTime();
      uint64_t gpshi = (gpsr >> 32) & 0xffffffff;
      uint64_t gpslo = gpsr & 0xffffffff;
      uint64_t gps = gpshi * 1000000 + gpslo;
      //  edm::LogInfo("L1TGT") << "  gpsr = " << std::hex << gpsr << " hi=" << gpshi << " lo=" << gpslo << " gps=" << gps << std::endl;

      Long64_t delorb = orbitTcs - preOrb_;
      Long64_t delgps = gps - preGps_;
      Double_t freq = -1.;

      if (delgps > 0) {
        freq = ((Double_t)(delorb)) * 3564. / ((Double_t)(delgps));
      }

      if (delorb > 0) {
        gpsfreq->Fill(freq);
        gpsfreqwide->Fill(freq);
        gpsfreqlum->Fill(lsNumber, freq);
        if (verbose_) {
          if (freq > 200.) {
            edm::LogInfo("L1TGT") << " preOrb_ = " << preOrb_ << " orbitTcs=" << orbitTcs << " delorb=" << delorb
                                  << std::hex << " preGps_=" << preGps_ << " gps=" << gps << std::dec
                                  << " delgps=" << delgps << " freq=" << freq << std::endl;
          }
        }
      }

      preGps_ = gps;
      preOrb_ = orbitTcs;
    }

    // get info from FDL if active
    if (isActive(gtfeEvmActiveBoards, FDL)) {
      const L1GtFdlWord& fdlWord = gtEvmReadoutRecord->gtFdlWord();

      orbitEvmFdl = fdlWord.orbitNr();
      lsEvmFdl = fdlWord.lumiSegmentNr();
    }

    if ((orbitTcs >= 0) && (orbitEvmFdl >= 0)) {
      int diffOrbit = static_cast<float>(orbitTcs - orbitEvmFdl);
      edm::LogInfo("L1TGT") << "\n orbitTcs = " << orbitTcs << " orbitEvmFdl = " << orbitEvmFdl
                            << " diffOrbit = " << diffOrbit << " orbitEvent = " << iEvent.orbitNumber() << std::endl;

      if (diffOrbit >= MaxOrbitNrDiffTcsFdlEvm) {
        m_monOrbitNrDiffTcsFdlEvm->Fill(MaxOrbitNrDiffTcsFdlEvm);

      } else if (diffOrbit <= -MaxOrbitNrDiffTcsFdlEvm) {
        m_monOrbitNrDiffTcsFdlEvm->Fill(-MaxOrbitNrDiffTcsFdlEvm);

      } else {
        m_monOrbitNrDiffTcsFdlEvm->Fill(diffOrbit);
        m_monOrbitNrDiffTcsFdlEvmLs->Fill(lsNumber, diffOrbit);
      }

    } else {
      if (orbitTcs >= 0) {
        // EVM_FDL error
        m_monOrbitNrDiffTcsFdlEvm->Fill(MaxOrbitNrDiffTcsFdlEvm);
      } else if (orbitEvmFdl >= 0) {
        // TCS error
        m_monOrbitNrDiffTcsFdlEvm->Fill(-MaxOrbitNrDiffTcsFdlEvm);

      } else {
        // TCS and EVM_FDL error
        m_monOrbitNrDiffTcsFdlEvm->Fill(-MaxOrbitNrDiffTcsFdlEvm);
        m_monOrbitNrDiffTcsFdlEvm->Fill(MaxOrbitNrDiffTcsFdlEvm);
      }
    }

    if ((lsTcs >= 0) && (lsEvmFdl >= 0)) {
      int diffLs = static_cast<float>(lsTcs - lsEvmFdl);
      edm::LogInfo("L1TGT") << "\n lsTcs = " << lsTcs << " lsEvmFdl = " << lsEvmFdl << " diffLs = " << diffLs
                            << " lsEvent = " << lsNumber << std::endl;

      if (diffLs >= MaxLsNrDiffTcsFdlEvm) {
        m_monLsNrDiffTcsFdlEvm->Fill(MaxLsNrDiffTcsFdlEvm);

      } else if (diffLs <= -MaxLsNrDiffTcsFdlEvm) {
        m_monLsNrDiffTcsFdlEvm->Fill(-MaxLsNrDiffTcsFdlEvm);

      } else {
        m_monLsNrDiffTcsFdlEvm->Fill(diffLs);
        m_monLsNrDiffTcsFdlEvmLs->Fill(lsNumber, diffLs);
      }

    } else {
      if (lsTcs >= 0) {
        // EVM_FDL error
        m_monLsNrDiffTcsFdlEvm->Fill(MaxLsNrDiffTcsFdlEvm);
      } else if (lsEvmFdl >= 0) {
        // TCS error
        m_monLsNrDiffTcsFdlEvm->Fill(-MaxLsNrDiffTcsFdlEvm);

      } else {
        // TCS and EVM_FDL error
        m_monLsNrDiffTcsFdlEvm->Fill(-MaxLsNrDiffTcsFdlEvm);
        m_monLsNrDiffTcsFdlEvm->Fill(MaxLsNrDiffTcsFdlEvm);
      }
    }
  }

  // open GT DAQ readout record - exit if failed
  edm::Handle<L1GlobalTriggerReadoutRecord> gtReadoutRecord;
  iEvent.getByToken(gtSource_L1GT_, gtReadoutRecord);

  //edm::ESHandle<L1GtTriggerMenu> menuRcd;
  //evSetup.get<L1GtTriggerMenuRcd>().get(menuRcd);

  //const L1GtTriggerMenu* menu = menuRcd.product();

  if (!gtReadoutRecord.isValid()) {
    edm::LogInfo("L1TGT") << "can't find L1GlobalTriggerReadoutRecord";
    return;
  }

  if (gtReadoutRecord.isValid()) {
    unsigned int NmaxL1AlgoBit = gtReadoutRecord->decisionWord().size();
    unsigned int NmaxL1TechBit = gtReadoutRecord->technicalTriggerWord().size();

    const DecisionWord dWord = gtReadoutRecord->decisionWord();
    const TechnicalTriggerWord technicalTriggerWordBeforeMask = gtReadoutRecord->technicalTriggerWord();

    const std::vector<L1GtFdlWord>& m_gtFdlWord(gtReadoutRecord->gtFdlVector());
    int numberBxInEvent = m_gtFdlWord.size();
    int minBxInEvent = (numberBxInEvent + 1) / 2 - numberBxInEvent;

    for (unsigned int iBit = 0; iBit < NmaxL1AlgoBit; ++iBit) {
      bool accept = dWord[iBit];

      typedef std::map<std::string, bool>::value_type valType;
      trig_iter = l1TriggerDecision.find(algoBitToName[iBit]);
      if (trig_iter == l1TriggerDecision.end()) {
        l1TriggerDecision.insert(valType(algoBitToName[iBit], accept));
      } else {
        trig_iter->second = accept;
      }

      int ibx = 0;
      for (std::vector<L1GtFdlWord>::const_iterator itBx = m_gtFdlWord.begin(); itBx != m_gtFdlWord.end(); ++itBx) {
        const DecisionWord dWordBX = (*itBx).gtDecisionWord();
        bool accept = dWordBX[iBit];
        if (accept) {
          if (iBit < 32)
            h_L1AlgoBX1->Fill(iBit, minBxInEvent + ibx);
          else if (iBit < 64)
            h_L1AlgoBX2->Fill(iBit, minBxInEvent + ibx);
          else if (iBit < 96)
            h_L1AlgoBX3->Fill(iBit, minBxInEvent + ibx);
          else if (iBit < 128)
            h_L1AlgoBX4->Fill(iBit, minBxInEvent + ibx);
        }
        ibx++;
      }
    }

    for (unsigned int iBit = 0; iBit < NmaxL1TechBit; ++iBit) {
      bool accept = technicalTriggerWordBeforeMask[iBit];

      typedef std::map<std::string, bool>::value_type valType;
      trig_iter = l1TechTriggerDecision.find(techBitToName[iBit]);
      if (trig_iter == l1TechTriggerDecision.end())
        l1TechTriggerDecision.insert(valType(techBitToName[iBit], accept));
      else
        trig_iter->second = accept;

      int ibx = 0;
      for (std::vector<L1GtFdlWord>::const_iterator itBx = m_gtFdlWord.begin(); itBx != m_gtFdlWord.end(); ++itBx) {
        const DecisionWord dWordBX = (*itBx).gtTechnicalTriggerWord();
        bool accept = dWordBX[iBit];
        if (accept)
          h_L1TechBX->Fill(iBit, minBxInEvent + ibx);
        ibx++;
      }
    }
  }

  // initialize bx's to invalid value
  int gtfeBx = -1;
  int fdlBx[2] = {-1, -1};
  int psbBx[2][7] = {{-1, -1, -1, -1, -1, -1, -1}, {-1, -1, -1, -1, -1, -1, -1}};
  int gmtBx = -1;

  // get info from GTFE DAQ record
  const L1GtfeWord& gtfeWord = gtReadoutRecord->gtfeWord();
  gtfeBx = gtfeWord.bxNr();
  gtfe_bx->Fill(gtfeBx);
  setupversion_lumi->Fill(lsNumber, gtfeWord.setupVersion());
  int gtfeActiveBoards = gtfeWord.activeBoards();

  // look for GMT readout collection from the same source if GMT active
  if (isActive(gtfeActiveBoards, GMT)) {
    edm::Handle<L1MuGMTReadoutCollection> gmtReadoutCollection;
    iEvent.getByToken(gtSource_L1MuGMT_, gmtReadoutCollection);

    if (gmtReadoutCollection.isValid()) {
      gmtBx = gmtReadoutCollection->getRecord().getBxNr();
    }
  }

  // get info from FDL if active (including decision word)
  if (isActive(gtfeActiveBoards, FDL)) {
    const L1GtFdlWord& fdlWord = gtReadoutRecord->gtFdlWord();
    fdlBx[0] = fdlWord.bxNr();
    fdlBx[1] = fdlWord.localBxNr();

    const DecisionWord& gtDecisionWord = gtReadoutRecord->decisionWord();
    const TechnicalTriggerWord& gtTTWord = gtReadoutRecord->technicalTriggerWord();

    int dbitNumber = 0;
    DecisionWord::const_iterator GTdbitItr;
    algo_bits->Fill(-1.);  // fill underflow to normalize
    for (GTdbitItr = gtDecisionWord.begin(); GTdbitItr != gtDecisionWord.end(); GTdbitItr++) {
      if (*GTdbitItr) {
        algo_bits->Fill(dbitNumber);
        algo_bits_lumi->Fill(lsNumber, dbitNumber);
        int dbitNumber1 = 0;
        DecisionWord::const_iterator GTdbitItr1;
        for (GTdbitItr1 = gtDecisionWord.begin(); GTdbitItr1 != gtDecisionWord.end(); GTdbitItr1++) {
          if (*GTdbitItr1)
            algo_bits_corr->Fill(dbitNumber, dbitNumber1);
          dbitNumber1++;
        }
        int tbitNumber1 = 0;
        TechnicalTriggerWord::const_iterator GTtbitItr1;
        for (GTtbitItr1 = gtTTWord.begin(); GTtbitItr1 != gtTTWord.end(); GTtbitItr1++) {
          if (*GTtbitItr1)
            algo_tt_bits_corr->Fill(dbitNumber, tbitNumber1);
          tbitNumber1++;
        }
      }
      dbitNumber++;
    }

    int tbitNumber = 0;
    TechnicalTriggerWord::const_iterator GTtbitItr;
    tt_bits->Fill(-1.);  // fill underflow to normalize
    for (GTtbitItr = gtTTWord.begin(); GTtbitItr != gtTTWord.end(); GTtbitItr++) {
      if (*GTtbitItr) {
        tt_bits->Fill(tbitNumber);
        tt_bits_lumi->Fill(lsNumber, tbitNumber);
        int tbitNumber1 = 0;
        TechnicalTriggerWord::const_iterator GTtbitItr1;
        for (GTtbitItr1 = gtTTWord.begin(); GTtbitItr1 != gtTTWord.end(); GTtbitItr1++) {
          if (*GTtbitItr1)
            tt_bits_corr->Fill(tbitNumber, tbitNumber1);
          tbitNumber1++;
        }
      }
      tbitNumber++;
    }

    // fill the index of actual prescale factor set
    // the index for technical triggers and algorithm trigger is the same (constraint in L1 GT TS)
    // so we read only pfIndexAlgoTrig (uint16_t)

    const int pfIndexAlgoTrig = fdlWord.gtPrescaleFactorIndexAlgo();
    m_monL1PrescaleFactorSet->Fill(lsNumber, static_cast<float>(pfIndexAlgoTrig));

    //

    // check that the combination (lsNumber, pfIndex) is not already included
    // to avoid fake entries due to different event order

    std::pair<int, int> pairLsPfi = std::make_pair(lsNumber, pfIndexAlgoTrig);

    CItVecPair cIt = find(m_pairLsNumberPfIndex.begin(), m_pairLsNumberPfIndex.end(), pairLsPfi);

    if (cIt == m_pairLsNumberPfIndex.end()) {
      m_pairLsNumberPfIndex.push_back(pairLsPfi);
    }
  }

  // get info from active PSB's
  int ibit = PSB9;  // first psb
  // for now hardcode psb id's - TODO - get them from Vasile's board maps...
  int psbID[7] = {0xbb09, 0xbb0d, 0xbb0e, 0xbb0f, 0xbb13, 0xbb14, 0xbb15};
  for (int i = 0; i < 7; i++) {
    if (isActive(gtfeActiveBoards, ibit)) {
      L1GtPsbWord psbWord = gtReadoutRecord->gtPsbWord(psbID[i]);
      psbBx[0][i] = psbWord.bxNr();
      psbBx[1][i] = psbWord.localBxNr();
    }
    ibit++;
  }

  //fill the dbx histo
  if (gtfeEvmBx > -1)
    dbx_module->Fill(0., gtfeEvmBx - gtfeBx);
  if (tcsBx > -1)
    dbx_module->Fill(1., tcsBx - gtfeBx);
  for (int i = 0; i < 2; i++) {
    if (fdlBx[i] > -1)
      dbx_module->Fill(2. + i, fdlBx[i] - gtfeBx);
  }
  for (int j = 0; j < 7; j++) {
    for (int i = 0; i < 2; i++) {
      if (psbBx[i][j] > -1)
        dbx_module->Fill(4. + i + 2 * j, psbBx[i][j] - gtfeBx);
    }
  }
  if (gmtBx > -1)
    dbx_module->Fill(18., gmtBx - gtfeBx);
}

void L1TGT::beginLuminosityBlock ( const edm::LuminosityBlock &  , const edm::EventSetup &    )

inlinefinalprotected

end section

Definition at line 64 of file L1TGT.h.

{}

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

overrideprotectedvirtual

Implements DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 37 of file L1TGT.cc.

References algo_bits, algo_bits_corr, algo_bits_lumi, algo_tt_bits_corr, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), BST_beamMode, BST_beamMomentum, BST_intensityBeam1, BST_intensityBeam2, BST_lhcFillNumber, BST_MasterStatus, BST_turnCountNumber, dbx_module, event_lumi, event_number, event_type, evnum_trignum_lumi, edm::EventSetup::getData(), gpsfreq, gpsfreqlum, gpsfreqwide, L1GtTriggerMenu::gtAlgorithmMap(), gtfe_bx, L1GtTriggerMenu::gtTechnicalTriggerMap(), h_L1AlgoBX1, h_L1AlgoBX2, h_L1AlgoBX3, h_L1AlgoBX4, h_L1TechBX, l1gtTrigmenuToken_, m_histFolder, m_monL1PfIndicesPerLs, m_monL1PrescaleFactorSet, m_monLsNrDiffTcsFdlEvm, m_monLsNrDiffTcsFdlEvmLs, m_monOrbitNrDiffTcsFdlEvm, m_monOrbitNrDiffTcsFdlEvmLs, m_pairLsNumberPfIndex, MaxLsNrDiffTcsFdlEvm, MaxOrbitNrDiffTcsFdlEvm, relval_steps::menu, orbit_lumi, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), setupversion_lumi, trigger_lumi, trigger_number, tt_bits, tt_bits_corr, and tt_bits_lumi.

                                                                                                {
  ibooker.setCurrentFolder(m_histFolder);

  // book histograms
  const int TotalNrBinsLs = 1000;
  const double totalNrBinsLs = static_cast<double>(TotalNrBinsLs);

  ibooker.setCurrentFolder(m_histFolder);

  algo_bits = ibooker.book1D("algo_bits", "GT algorithm trigger bits", 128, -0.5, 127.5);
  algo_bits->setAxisTitle("Algorithm trigger bits", 1);

  algo_bits_corr =
      ibooker.book2D("algo_bits_corr", "GT algorithm trigger bit correlation", 128, -0.5, 127.5, 128, -0.5, 127.5);
  algo_bits_corr->setAxisTitle("Algorithm trigger bits", 1);
  algo_bits_corr->setAxisTitle("Algorithm trigger bits", 2);

  tt_bits = ibooker.book1D("tt_bits", "GT technical trigger bits", 64, -0.5, 63.5);
  tt_bits->setAxisTitle("Technical trigger bits", 1);

  tt_bits_corr = ibooker.book2D("tt_bits_corr", "GT technical trigger bit correlation", 64, -0.5, 63.5, 64, -0.5, 63.5);
  tt_bits_corr->setAxisTitle("Technical trigger bits", 1);
  tt_bits_corr->setAxisTitle("Technical trigger bits", 2);

  algo_tt_bits_corr = ibooker.book2D(
      "algo_tt_bits_corr", "GT algorithm - technical trigger bit correlation", 128, -0.5, 127.5, 64, -0.5, 63.5);
  algo_tt_bits_corr->setAxisTitle("Algorithm trigger bits", 1);
  algo_tt_bits_corr->setAxisTitle("Technical trigger bits", 2);

  algo_bits_lumi = ibooker.book2D(
      "algo_bits_lumi", "GT algorithm trigger bit rate per LS", TotalNrBinsLs, 0., totalNrBinsLs, 128, -0.5, 127.5);
  algo_bits_lumi->setAxisTitle("Luminosity segment", 1);
  algo_bits_lumi->setAxisTitle("Algorithm trigger bits", 2);

  tt_bits_lumi = ibooker.book2D(
      "tt_bits_lumi", "GT technical trigger bit rate per LS", TotalNrBinsLs, 0., totalNrBinsLs, 64, -0.5, 63.5);
  tt_bits_lumi->setAxisTitle("Luminosity segment", 1);
  tt_bits_lumi->setAxisTitle("Technical trigger bits", 2);

  event_type = ibooker.book1D("event_type", "GT event type", 10, -0.5, 9.5);
  event_type->setAxisTitle("Event type", 1);
  event_type->setBinLabel(2, "Physics", 1);
  event_type->setBinLabel(3, "Calibration", 1);
  event_type->setBinLabel(4, "Random", 1);
  event_type->setBinLabel(6, "Traced", 1);
  event_type->setBinLabel(7, "Test", 1);
  event_type->setBinLabel(8, "Error", 1);

  event_number = ibooker.book1D("event_number", "GT event number (from last resync)", 100, 0., 50000.);
  event_number->setAxisTitle("Event number", 1);

  event_lumi = ibooker.bookProfile(
      "event_lumi", "GT event number (from last resync) vs LS", TotalNrBinsLs, 0., totalNrBinsLs, 100, -0.1, 1.e15, "s");
  event_lumi->setAxisTitle("Luminosity segment", 1);
  event_lumi->setAxisTitle("Event number", 2);

  trigger_number = ibooker.book1D("trigger_number", "GT trigger number (from start run)", 100, 0., 50000.);
  trigger_number->setAxisTitle("Trigger number", 1);

  trigger_lumi = ibooker.bookProfile("trigger_lumi",
                                     "GT trigger number (from start run) vs LS",
                                     TotalNrBinsLs,
                                     0.,
                                     totalNrBinsLs,
                                     100,
                                     -0.1,
                                     1.e15,
                                     "s");
  trigger_lumi->setAxisTitle("Luminosity segment", 1);
  trigger_lumi->setAxisTitle("Trigger number", 2);

  evnum_trignum_lumi = ibooker.bookProfile(
      "evnum_trignum_lumi", "GT event/trigger number ratio vs LS", TotalNrBinsLs, 0., totalNrBinsLs, 100, -0.1, 2., "s");
  evnum_trignum_lumi->setAxisTitle("Luminosity segment", 1);
  evnum_trignum_lumi->setAxisTitle("Event/trigger number ratio", 2);

  orbit_lumi = ibooker.bookProfile(
      "orbit_lumi", "GT orbit number vs LS", TotalNrBinsLs, 0., totalNrBinsLs, 100, -0.1, 1.e15, "s");
  orbit_lumi->setAxisTitle("Luminosity segment", 1);
  orbit_lumi->setAxisTitle("Orbit number", 2);

  setupversion_lumi = ibooker.bookProfile(
      "setupversion_lumi", "GT setup version vs LS", TotalNrBinsLs, 0., totalNrBinsLs, 100, -0.1, 1.e10, "i");
  setupversion_lumi->setAxisTitle("Luminosity segment", 1);
  setupversion_lumi->setAxisTitle("Setup version", 2);

  gtfe_bx = ibooker.book1D("gtfe_bx", "GTFE Bx number", 3600, 0., 3600.);
  gtfe_bx->setAxisTitle("GTFE BX number", 1);

  dbx_module =
      ibooker.bookProfile("dbx_module", "delta Bx of GT modules wrt GTFE", 20, 0., 20., 100, -4000., 4000., "i");
  dbx_module->setAxisTitle("GT crate module", 1);
  dbx_module->setAxisTitle("Module Bx - GTFE Bx", 2);
  dbx_module->setBinLabel(1, "GTFEevm", 1);
  dbx_module->setBinLabel(2, "TCS", 1);
  dbx_module->setBinLabel(3, "FDL", 1);
  dbx_module->setBinLabel(4, "FDLloc", 1);
  dbx_module->setBinLabel(5, "PSB9", 1);
  dbx_module->setBinLabel(6, "PSB9loc", 1);
  dbx_module->setBinLabel(7, "PSB13", 1);
  dbx_module->setBinLabel(8, "PSB13loc", 1);
  dbx_module->setBinLabel(9, "PSB14", 1);
  dbx_module->setBinLabel(10, "PSB14loc", 1);
  dbx_module->setBinLabel(11, "PSB15", 1);
  dbx_module->setBinLabel(12, "PSB15loc", 1);
  dbx_module->setBinLabel(13, "PSB19", 1);
  dbx_module->setBinLabel(14, "PSB19loc", 1);
  dbx_module->setBinLabel(15, "PSB20", 1);
  dbx_module->setBinLabel(16, "PSB20loc", 1);
  dbx_module->setBinLabel(17, "PSB21", 1);
  dbx_module->setBinLabel(18, "PSB21loc", 1);
  dbx_module->setBinLabel(19, "GMT", 1);

  BST_MasterStatus =
      ibooker.book2D("BST_MasterStatus", "BST master status over LS", TotalNrBinsLs, 0., totalNrBinsLs, 6, -1., 5.);
  BST_MasterStatus->setAxisTitle("Luminosity segment", 1);
  BST_MasterStatus->setAxisTitle("BST master status", 2);
  BST_MasterStatus->setBinLabel(2, "Master Beam 1", 2);
  BST_MasterStatus->setBinLabel(3, "Master Beam 2", 2);

  BST_turnCountNumber =
      ibooker.book2D("BST_turnCountNumber", "BST turn count over LS", TotalNrBinsLs, 0., totalNrBinsLs, 250, 0., 4.3e9);
  BST_turnCountNumber->setAxisTitle("Luminosity segment", 1);
  BST_turnCountNumber->setAxisTitle("BST turn count number", 2);

  BST_lhcFillNumber = ibooker.book1D("BST_lhcFillNumber", "BST LHC fill number % 1000", 1000, 0., 1000.);
  BST_lhcFillNumber->setAxisTitle("BST LHC fill number modulo 1000");

  BST_beamMode = ibooker.book2D("BST_beamMode", "BST beam mode over LS", TotalNrBinsLs, 0., totalNrBinsLs, 25, 1., 26.);
  BST_beamMode->setAxisTitle("Luminosity segment", 1);
  BST_beamMode->setAxisTitle("Mode", 2);
  BST_beamMode->setBinLabel(1, "No mode", 2);
  BST_beamMode->setBinLabel(2, "Setup", 2);
  BST_beamMode->setBinLabel(3, "Inj pilot", 2);
  BST_beamMode->setBinLabel(4, "Inj intr", 2);
  BST_beamMode->setBinLabel(5, "Inj nomn", 2);
  BST_beamMode->setBinLabel(6, "Pre ramp", 2);
  BST_beamMode->setBinLabel(7, "Ramp", 2);
  BST_beamMode->setBinLabel(8, "Flat top", 2);
  BST_beamMode->setBinLabel(9, "Squeeze", 2);
  BST_beamMode->setBinLabel(10, "Adjust", 2);
  BST_beamMode->setBinLabel(11, "Stable", 2);
  BST_beamMode->setBinLabel(12, "Unstable", 2);
  BST_beamMode->setBinLabel(13, "Beam dump", 2);
  BST_beamMode->setBinLabel(14, "Ramp down", 2);
  BST_beamMode->setBinLabel(15, "Recovery", 2);
  BST_beamMode->setBinLabel(16, "Inj dump", 2);
  BST_beamMode->setBinLabel(17, "Circ dump", 2);
  BST_beamMode->setBinLabel(18, "Abort", 2);
  BST_beamMode->setBinLabel(19, "Cycling", 2);
  BST_beamMode->setBinLabel(20, "Warn beam dump", 2);
  BST_beamMode->setBinLabel(21, "No beam", 2);

  BST_beamMomentum =
      ibooker.book2D("BST_beamMomentum", "BST beam momentum", TotalNrBinsLs, 0., totalNrBinsLs, 100, 0., 7200.);
  BST_beamMomentum->setAxisTitle("Luminosity segment", 1);
  BST_beamMomentum->setAxisTitle("Beam momentum", 2);

  gpsfreq = ibooker.book1D("gpsfreq", "Clock frequency measured by GPS", 1000, 39.95, 40.2);
  gpsfreq->setAxisTitle("CMS clock frequency (MHz)");

  gpsfreqwide = ibooker.book1D("gpsfreqwide", "Clock frequency measured by GPS", 1000, -2., 200.);
  gpsfreqwide->setAxisTitle("CMS clock frequency (MHz)");

  gpsfreqlum = ibooker.book2D(
      "gpsfreqlum", "Clock frequency measured by GPS", TotalNrBinsLs, 0., totalNrBinsLs, 100, 39.95, 40.2);
  gpsfreqlum->setAxisTitle("Luminosity segment", 1);
  gpsfreqlum->setAxisTitle("CMS clock frequency (MHz)", 2);

  BST_intensityBeam1 =
      ibooker.book2D("BST_intensityBeam1", "Intensity beam 1", TotalNrBinsLs, 0., totalNrBinsLs, 1000, 0., 5000.);
  BST_intensityBeam1->setAxisTitle("Luminosity segment", 1);
  BST_intensityBeam1->setAxisTitle("Beam intensity", 2);

  BST_intensityBeam2 =
      ibooker.book2D("BST_intensityBeam2", "Intensity beam 2", TotalNrBinsLs, 0., totalNrBinsLs, 1000, 0., 5000.);
  BST_intensityBeam2->setAxisTitle("Luminosity segment", 1);
  BST_intensityBeam2->setAxisTitle("Beam intensity", 2);

  // prescale factor index monitoring

  m_monL1PrescaleFactorSet = ibooker.book2D(
      "L1PrescaleFactorSet", "Index of L1 prescale factor set", TotalNrBinsLs, 0., totalNrBinsLs, 25, 0., 25.);
  m_monL1PrescaleFactorSet->setAxisTitle("Luminosity segment", 1);
  m_monL1PrescaleFactorSet->setAxisTitle("L1 PF set index", 2);
  m_monL1PfIndicesPerLs =
      ibooker.book1D("L1PfIndicesPerLs", "Number of prescale factor indices used per LS", 10, 0., 10.);
  m_monL1PfIndicesPerLs->setAxisTitle("Number of PF indices used per LS", 1);
  m_monL1PfIndicesPerLs->setAxisTitle("Entries", 2);

  // TCS vs FDL common quantity monitoring

  ibooker.setCurrentFolder(m_histFolder + "/TCSvsEvmFDL");

  //    orbit number
  m_monOrbitNrDiffTcsFdlEvm = ibooker.book1D("OrbitNrDiffTcsFdlEvm",
                                             "Orbit number difference (TCS - EVM_FDL)",
                                             2 * MaxOrbitNrDiffTcsFdlEvm + 1,
                                             static_cast<float>(-(MaxOrbitNrDiffTcsFdlEvm + 1)),
                                             static_cast<float>(MaxOrbitNrDiffTcsFdlEvm + 1));
  m_monOrbitNrDiffTcsFdlEvm->setAxisTitle("Orbit number difference", 1);
  m_monOrbitNrDiffTcsFdlEvm->setAxisTitle("Entries/run", 2);
  m_monLsNrDiffTcsFdlEvm = ibooker.book1D("LsNrDiffTcsFdlEvm",
                                          "LS number difference (TCS - EVM_FDL)",
                                          2 * MaxLsNrDiffTcsFdlEvm + 1,
                                          static_cast<float>(-(MaxLsNrDiffTcsFdlEvm + 1)),
                                          static_cast<float>(MaxLsNrDiffTcsFdlEvm + 1));
  m_monLsNrDiffTcsFdlEvm->setAxisTitle("LS number difference", 1);
  m_monLsNrDiffTcsFdlEvm->setAxisTitle("Entries/run", 2);

  //    LS number
  m_monOrbitNrDiffTcsFdlEvmLs = ibooker.book2D("OrbitNrDiffTcsFdlEvmLs",
                                               "Orbit number difference (TCS - EVM_FDL)",
                                               TotalNrBinsLs,
                                               0.,
                                               totalNrBinsLs,
                                               2 * MaxOrbitNrDiffTcsFdlEvm + 1,
                                               static_cast<float>(-(MaxOrbitNrDiffTcsFdlEvm + 1)),
                                               static_cast<float>(MaxOrbitNrDiffTcsFdlEvm + 1));
  m_monOrbitNrDiffTcsFdlEvmLs->setAxisTitle("Luminosity segment", 1);
  m_monOrbitNrDiffTcsFdlEvmLs->setAxisTitle("Orbit number difference (TCS - EVM_FDL)", 2);

  m_monLsNrDiffTcsFdlEvmLs = ibooker.book2D("LsNrDiffTcsFdlEvmLs",
                                            "LS number difference (TCS - EVM_FDL)",
                                            TotalNrBinsLs,
                                            0.,
                                            totalNrBinsLs,
                                            2 * MaxLsNrDiffTcsFdlEvm + 1,
                                            static_cast<float>(-(MaxLsNrDiffTcsFdlEvm + 1)),
                                            static_cast<float>(MaxLsNrDiffTcsFdlEvm + 1));
  m_monLsNrDiffTcsFdlEvmLs->setAxisTitle("Luminosity segment", 1);
  m_monLsNrDiffTcsFdlEvmLs->setAxisTitle("LS number difference (TCS - EVM_FDL)", 2);

  ibooker.setCurrentFolder(m_histFolder);
  // clear bookkeeping for prescale factor change
  m_pairLsNumberPfIndex.clear();

  ibooker.setCurrentFolder(m_histFolder + "/PlotTrigsBx");

  //--------book AlgoBits/TechBits vs Bx Histogram-----------

  //edm::ESHandle<L1GtTriggerMenu> menuRcd;
  //evSetup.get<L1GtTriggerMenuRcd>().get(menuRcd);
  //menuRcd.product();
  const L1GtTriggerMenu* menu = &evSetup.getData(l1gtTrigmenuToken_);

  h_L1AlgoBX1 = ibooker.book2D("h_L1AlgoBX1", "L1 Algo Trigger BX (algo bit 0 to 31)", 32, -0.5, 31.5, 5, -2.5, 2.5);
  h_L1AlgoBX2 = ibooker.book2D("h_L1AlgoBX2", "L1 Algo Trigger BX (algo bit 32 to 63)", 32, 31.5, 63.5, 5, -2.5, 2.5);
  h_L1AlgoBX3 = ibooker.book2D("h_L1AlgoBX3", "L1 Algo Trigger BX (algo bit 64 to 95)", 32, 63.5, 95.5, 5, -2.5, 2.5);
  h_L1AlgoBX4 = ibooker.book2D("h_L1AlgoBX4", "L1 Algo Trigger BX (algo bit 96 to 127)", 32, 95.5, 127.5, 5, -2.5, 2.5);
  h_L1TechBX = ibooker.book2D("h_L1TechBX", "L1 Tech Trigger BX", 64, -0.5, 63.5, 5, -2.5, 2.5);

  for (CItAlgo algo = menu->gtAlgorithmMap().begin(); algo != menu->gtAlgorithmMap().end(); ++algo) {
    int itrig = (algo->second).algoBitNumber();
    //algoBitToName[itrig] = TString( (algo->second).algoName() );
    //const char* trigName =  (algo->second).algoName().c_str();
    if (itrig < 32) {
      //h_L1AlgoBX1->setBinLabel(itrig+1,trigName);
      h_L1AlgoBX1->setBinLabel(itrig + 1, std::to_string(itrig));
      h_L1AlgoBX1->setAxisTitle("Algorithm trigger bits", 1);
      h_L1AlgoBX1->setAxisTitle("BX (0=L1A)", 2);
    } else if (itrig < 64) {
      //h_L1AlgoBX2->setBinLabel(itrig+1-32,trigName);
      h_L1AlgoBX2->setBinLabel(itrig + 1 - 32, std::to_string(itrig));
      h_L1AlgoBX2->setAxisTitle("Algorithm trigger bits", 1);
      h_L1AlgoBX2->setAxisTitle("BX (0=L1A)", 2);
    } else if (itrig < 96) {
      //h_L1AlgoBX3->setBinLabel(itrig+1-64,trigName);
      h_L1AlgoBX3->setBinLabel(itrig + 1 - 64, std::to_string(itrig));
      h_L1AlgoBX3->setAxisTitle("Algorithm trigger bits", 1);
      h_L1AlgoBX3->setAxisTitle("BX (0=L1A)", 2);
    } else if (itrig < 128) {
      //h_L1AlgoBX4->setBinLabel(itrig+1-96,trigName);
      h_L1AlgoBX4->setBinLabel(itrig + 1 - 96, std::to_string(itrig));
      h_L1AlgoBX4->setAxisTitle("Algorithm trigger bits", 1);
      h_L1AlgoBX4->setAxisTitle("BX (0=L1A)", 2);
    }
  }

  // technical trigger bits
  for (CItAlgo techTrig = menu->gtTechnicalTriggerMap().begin(); techTrig != menu->gtTechnicalTriggerMap().end();
       ++techTrig) {
    int itrig = (techTrig->second).algoBitNumber();
    //techBitToName[itrig] = TString( (techTrig->second).algoName() );
    //const char* trigName =  (techTrig->second).algoName().c_str();
    h_L1TechBX->setBinLabel(itrig + 1, std::to_string(itrig));
    h_L1TechBX->setAxisTitle("Technical trigger bits", 1);
    h_L1TechBX->setAxisTitle("BX (0=L1A)", 2);
  }
}

void L1TGT::countPfsIndicesPerLs ( )

private

Definition at line 728 of file L1TGT.cc.

References dqm::impl::MonitorElement::Fill(), m_monL1PfIndicesPerLs, m_pairLsNumberPfIndex, dqm::impl::MonitorElement::Reset(), and verbose_.

Referenced by endLuminosityBlock().

                                 {
  if (verbose_) {
    edm::LogInfo("L1TGT") << "\n  Prescale factor indices used in a LS " << std::endl;

    for (CItVecPair cIt = m_pairLsNumberPfIndex.begin(); cIt != m_pairLsNumberPfIndex.end(); ++cIt) {
      edm::LogVerbatim("L1TGT") << "  lsNumber = " << (*cIt).first << " pfIndex = " << (*cIt).second << std::endl;
    }
    edm::LogVerbatim("L1TGT") << std::endl;
  }

  // reset the histogram...
  m_monL1PfIndicesPerLs->Reset();

  // sort the vector (for pairs: sort after first argument, then after the second argument)
  std::sort(m_pairLsNumberPfIndex.begin(), m_pairLsNumberPfIndex.end());

  int previousLsNumber = -1;
  int previousPfsIndex = -1;

  // count the number of pairs (lsNumber, pfIndex) per Ls
  // there are no duplicate entries, and pairs are sorted after both members
  // ... and fill again the histogram
  int pfsIndicesPerLs = 1;
  for (CItVecPair cIt = m_pairLsNumberPfIndex.begin(); cIt != m_pairLsNumberPfIndex.end(); ++cIt) {
    if ((*cIt).first == previousLsNumber) {
      if ((*cIt).second != previousPfsIndex) {
        previousPfsIndex = (*cIt).second;
        pfsIndicesPerLs++;
      }

    } else {
      // fill the histogram with the number of PF indices for the previous Ls
      if (previousLsNumber != -1) {
        m_monL1PfIndicesPerLs->Fill(pfsIndicesPerLs);
      }

      // new Ls
      previousLsNumber = (*cIt).first;
      previousPfsIndex = (*cIt).second;
      pfsIndicesPerLs = 1;
    }
  }
}

void L1TGT::dqmBeginRun ( const edm::Run &  iRrun, const edm::EventSetup &  evSetup  )

overrideprotectedvirtual

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 328 of file L1TGT.cc.

References m_nrEvRun.

{ m_nrEvRun = 0; }

void L1TGT::endLuminosityBlock ( const edm::LuminosityBlock &  iLumi, const edm::EventSetup &  evSetup  )

overrideprotected

Definition at line 715 of file L1TGT.cc.

References countPfsIndicesPerLs(), and m_runInEndLumi.

                                                                                            {
  if (m_runInEndLumi) {
    countPfsIndicesPerLs();
  }
}

bool L1TGT::isActive ( int  word, int  bit  )

private

book all histograms for the module

Definition at line 722 of file L1TGT.cc.

Referenced by analyze().

                                      {
  if (word & (1 << bit))
    return true;
  return false;
}

Member Data Documentation

MonitorElement* L1TGT::algo_bits

private

Definition at line 110 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::algo_bits_corr

private

Definition at line 111 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::algo_bits_lumi

private

Definition at line 115 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::algo_tt_bits_corr

private

Definition at line 114 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

std::string L1TGT::algoBitToName[128]

private

Definition at line 178 of file L1TGT.h.

Referenced by analyze().

MonitorElement* L1TGT::BST_beamMode

private

Definition at line 133 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::BST_beamMomentum

private

Definition at line 134 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::BST_intensityBeam1

private

Definition at line 135 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::BST_intensityBeam2

private

Definition at line 136 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::BST_lhcFillNumber

private

Definition at line 132 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::BST_MasterStatus

private

Definition at line 130 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::BST_turnCountNumber

private

Definition at line 131 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::dbx_module

private

Definition at line 128 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::event_lumi

private

Definition at line 120 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::event_number

private

Definition at line 119 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::event_type

private

Definition at line 117 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::evnum_trignum_lumi

private

Definition at line 123 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::gpsfreq

private

Definition at line 137 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::gpsfreqlum

private

Definition at line 139 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::gpsfreqwide

private

Definition at line 138 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<L1GlobalTriggerEvmReadoutRecord> L1TGT::gtEvmSource_

private

input tag for L1 GT EVM readout record

Definition at line 99 of file L1TGT.h.

Referenced by analyze().

MonitorElement* L1TGT::gtfe_bx

private

Definition at line 127 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> L1TGT::gtSource_L1GT_

private

input parameters

input tag for L1 GT DAQ readout record

Definition at line 95 of file L1TGT.h.

Referenced by analyze().

edm::EDGetTokenT<L1MuGMTReadoutCollection> L1TGT::gtSource_L1MuGMT_

private

Definition at line 96 of file L1TGT.h.

Referenced by analyze().

MonitorElement* L1TGT::h_L1AlgoBX1

private

Definition at line 159 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::h_L1AlgoBX2

private

Definition at line 160 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::h_L1AlgoBX3

private

Definition at line 161 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::h_L1AlgoBX4

private

Definition at line 162 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::h_L1TechBX

private

Definition at line 163 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

edm::ESGetToken<L1GtTriggerMenu, L1GtTriggerMenuRcd> L1TGT::l1gtTrigmenuToken_

private

Definition at line 100 of file L1TGT.h.

Referenced by bookHistograms(), and L1TGT().

std::map<std::string, bool> L1TGT::l1TechTriggerDecision

private

Definition at line 180 of file L1TGT.h.

Referenced by analyze().

std::map<std::string, bool> L1TGT::l1TriggerDecision

private

Definition at line 180 of file L1TGT.h.

Referenced by analyze().

std::string L1TGT::m_histFolder

private

histogram folder for L1 GT plots

Definition at line 173 of file L1TGT.h.

Referenced by bookHistograms(), and L1TGT().

MonitorElement* L1TGT::m_monL1PfIndicesPerLs

private

Definition at line 142 of file L1TGT.h.

Referenced by bookHistograms(), and countPfsIndicesPerLs().

MonitorElement* L1TGT::m_monL1PrescaleFactorSet

private

Definition at line 141 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::m_monLsNrDiffTcsFdlEvm

private

Definition at line 145 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::m_monLsNrDiffTcsFdlEvmLs

private

Definition at line 157 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::m_monOrbitNrDiffTcsFdlEvm

private

Definition at line 144 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::m_monOrbitNrDiffTcsFdlEvmLs

private

Definition at line 156 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

int L1TGT::m_nrEvJob

private

number of events processed

Definition at line 169 of file L1TGT.h.

Referenced by analyze().

int L1TGT::m_nrEvRun

private

Definition at line 170 of file L1TGT.h.

Referenced by dqmBeginRun().

std::vector<std::pair<int, int> > L1TGT::m_pairLsNumberPfIndex

private

Definition at line 183 of file L1TGT.h.

Referenced by analyze(), bookHistograms(), and countPfsIndicesPerLs().

bool L1TGT::m_runInEndLumi

private

Definition at line 104 of file L1TGT.h.

Referenced by endLuminosityBlock().

bool L1TGT::m_runInEventLoop

private

switches to choose the running of various methods

Definition at line 103 of file L1TGT.h.

const int L1TGT::MaxLsNrDiffTcsFdlEvm = 24

staticprivate

Definition at line 154 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

const int L1TGT::MaxOrbitNrDiffTcsFdlEvm = 24

staticprivate

Definition at line 153 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::orbit_lumi

private

Definition at line 124 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

uint64_t L1TGT::preGps_

private

Definition at line 175 of file L1TGT.h.

Referenced by analyze().

uint64_t L1TGT::preOrb_

private

Definition at line 176 of file L1TGT.h.

Referenced by analyze().

MonitorElement* L1TGT::setupversion_lumi

private

Definition at line 125 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

std::string L1TGT::techBitToName[64]

private

Definition at line 179 of file L1TGT.h.

Referenced by analyze().

std::map<std::string, bool>::iterator L1TGT::trig_iter

private

Definition at line 181 of file L1TGT.h.

Referenced by analyze().

MonitorElement* L1TGT::trigger_lumi

private

Definition at line 122 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::trigger_number

private

Definition at line 121 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::tt_bits

private

Definition at line 112 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::tt_bits_corr

private

Definition at line 113 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGT::tt_bits_lumi

private

Definition at line 116 of file L1TGT.h.

Referenced by analyze(), and bookHistograms().

bool L1TGT::verbose_

private

verbosity switch

Definition at line 107 of file L1TGT.h.

Referenced by analyze(), and countPfsIndicesPerLs().