#include <L1TGMT.h>

Inheritance diagram for L1TGMT:

Public Member Functions
	L1TGMT (const edm::ParameterSet &ps)

	~L1TGMT () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Types
enum	ensubs { DTTF = 0, RPCb, CSCTF, RPCf, GMT }

Private Member Functions
void	book_ (const edm::EventSetup &c)

double	phiconv_ (float phi)

Private Attributes
MonitorElement *	bx_csc_rpc

MonitorElement *	bx_dt_csc

MonitorElement *	bx_dt_rpc

MonitorElement *	bx_number

int	bxnum_old_

MonitorElement *	dbx_chip

MonitorElement *	dist_eta_csc_rpc

MonitorElement *	dist_eta_dt_csc

MonitorElement *	dist_eta_dt_rpc

MonitorElement *	dist_phi_csc_rpc

MonitorElement *	dist_phi_dt_csc

MonitorElement *	dist_phi_dt_rpc

MonitorElement *	eta_dtcsc_and_rpc

MonitorElement *	eta_dtcsc_only

MonitorElement *	eta_rpc_only

MonitorElement *	etaphi_dtcsc_and_rpc

MonitorElement *	etaphi_dtcsc_only

MonitorElement *	etaphi_rpc_only

const edm::EDGetTokenT< L1MuGMTReadoutCollection >	gmtSource_

edm::ESGetToken< L1MuTriggerScales, L1MuTriggerScalesRcd >	l1muTrigscaleToken_

edm::ESGetToken< L1MuTriggerPtScale, L1MuTriggerPtScaleRcd >	l1TrigptscaleToken_

std::ofstream	logFile_

MonitorElement *	n_csctf_vs_dttf

MonitorElement *	n_rpcb_vs_dttf

MonitorElement *	n_rpcf_vs_csctf

int	obnum_old_

MonitorElement *	phi_dtcsc_and_rpc

MonitorElement *	phi_dtcsc_only

MonitorElement *	phi_rpc_only

MonitorElement *	regional_triggers

MonitorElement *	subs_bits [5]

MonitorElement *	subs_dbx [4]

MonitorElement *	subs_eta [5]

MonitorElement *	subs_etaphi [5]

MonitorElement *	subs_etaqty [5]

MonitorElement *	subs_nbx [5]

MonitorElement *	subs_phi [5]

MonitorElement *	subs_pt [5]

MonitorElement *	subs_qty [5]

int	trsrc_old_

const bool	verbose_

Static Private Attributes
static const double	piconv_ = 180. / acos(-1.)

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 46 of file L1TGMT.h.

Member Enumeration Documentation

◆ ensubs

enum L1TGMT::ensubs

private

Enumerator
DTTF
RPCb
CSCTF
RPCf
GMT

Definition at line 63 of file L1TGMT.h.

63 { DTTF = 0, RPCb, CSCTF, RPCf, GMT };

Constructor & Destructor Documentation

◆ L1TGMT()

L1TGMT::L1TGMT ( const edm::ParameterSet & ps )

Definition at line 23 of file L1TGMT.cc.

     : verbose_(ps.getUntrackedParameter<bool>("verbose", false))  // verbosity switch
       ,
       gmtSource_(consumes<L1MuGMTReadoutCollection>(ps.getParameter<InputTag>("gmtSource"))),
       bxnum_old_(0),
       obnum_old_(0),
       trsrc_old_(0) {
   if (verbose_)
     cout << "L1TGMT: constructor...." << endl;
   l1muTrigscaleToken_ = esConsumes<edm::Transition::BeginRun>();
   l1TrigptscaleToken_ = esConsumes<edm::Transition::BeginRun>();
 }

References gather_cfg::cout, l1muTrigscaleToken_, l1TrigptscaleToken_, and verbose_.

◆ ~L1TGMT()

L1TGMT::~L1TGMT ( )

override

Definition at line 36 of file L1TGMT.cc.

36 {}

Member Function Documentation

◆ analyze()

void L1TGMT::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 38 of file L1TGMT.cc.

                                                         {
   if (verbose_)
     cout << "L1TGMT: analyze...." << endl;
  
   edm::Handle<L1MuGMTReadoutCollection> pCollection;
   e.getByToken(gmtSource_, pCollection);
  
   if (!pCollection.isValid()) {
     edm::LogInfo("DataNotFound") << "can't find L1MuGMTReadoutCollection";
     return;
   }
  
   // remember the bx of 1st candidate of each system (9=none)
   int bx1st[4] = {9, 9, 9, 9};
  
   // get GMT readout collection
   L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
   // get record vector
   vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
   // loop over records of individual bx's
   vector<L1MuGMTReadoutRecord>::const_iterator RRItr;
  
   for (RRItr = gmt_records.begin(); RRItr != gmt_records.end(); RRItr++) {
     vector<L1MuRegionalCand> INPCands[4] = {
         RRItr->getDTBXCands(), RRItr->getBrlRPCCands(), RRItr->getCSCCands(), RRItr->getFwdRPCCands()};
     vector<L1MuGMTExtendedCand> GMTCands = RRItr->getGMTCands();
  
     vector<L1MuRegionalCand>::const_iterator INPItr;
     vector<L1MuGMTExtendedCand>::const_iterator GMTItr;
     vector<L1MuGMTExtendedCand>::const_iterator GMTItr2;
  
     int BxInEvent = RRItr->getBxInEvent();
  
     // count non-empty candidates in this bx
     int nSUBS[5] = {0, 0, 0, 0, 0};
     for (int i = 0; i < 4; i++) {
       for (INPItr = INPCands[i].begin(); INPItr != INPCands[i].end(); ++INPItr) {
         if (!INPItr->empty()) {
           nSUBS[i]++;
           if (bx1st[i] == 9)
             bx1st[i] = BxInEvent;
         }
       }
       subs_nbx[i]->Fill(float(nSUBS[i]), float(BxInEvent));
     }
  
     for (GMTItr = GMTCands.begin(); GMTItr != GMTCands.end(); ++GMTItr) {
       if (!GMTItr->empty())
         nSUBS[GMT]++;
     }
     subs_nbx[GMT]->Fill(float(nSUBS[GMT]), float(BxInEvent));
  
     // from here care only about the L1A bunch crossing
     if (BxInEvent != 0)
       continue;
  
     // get the absolute bx number of the L1A
     int Bx = RRItr->getBxNr();
  
     bx_number->Fill(double(Bx));
  
     for (int i = 0; i < 4; i++) {
       for (INPItr = INPCands[i].begin(); INPItr != INPCands[i].end(); ++INPItr) {
         if (INPItr->empty())
           continue;
         subs_eta[i]->Fill(INPItr->etaValue());
         subs_phi[i]->Fill(phiconv_(INPItr->phiValue()));
         subs_pt[i]->Fill(INPItr->ptValue());
         subs_qty[i]->Fill(INPItr->quality());
         subs_etaphi[i]->Fill(INPItr->etaValue(), phiconv_(INPItr->phiValue()));
         subs_etaqty[i]->Fill(INPItr->etaValue(), INPItr->quality());
         int word = INPItr->getDataWord();
         for (int j = 0; j < 32; j++) {
           if (word & (1 << j))
             subs_bits[i]->Fill(float(j));
         }
       }
     }
  
     for (GMTItr = GMTCands.begin(); GMTItr != GMTCands.end(); ++GMTItr) {
       if (GMTItr->empty())
         continue;
       subs_eta[GMT]->Fill(GMTItr->etaValue());
       subs_phi[GMT]->Fill(phiconv_(GMTItr->phiValue()));
       subs_pt[GMT]->Fill(GMTItr->ptValue());
       subs_qty[GMT]->Fill(GMTItr->quality());
       subs_etaphi[GMT]->Fill(GMTItr->etaValue(), phiconv_(GMTItr->phiValue()));
       subs_etaqty[GMT]->Fill(GMTItr->etaValue(), GMTItr->quality());
       int word = GMTItr->getDataWord();
       for (int j = 0; j < 32; j++) {
         if (word & (1 << j))
           subs_bits[GMT]->Fill(float(j));
       }
  
       if (GMTItr->isMatchedCand()) {
         if (GMTItr->quality() > 3) {
           eta_dtcsc_and_rpc->Fill(GMTItr->etaValue());
           phi_dtcsc_and_rpc->Fill(phiconv_(GMTItr->phiValue()));
           etaphi_dtcsc_and_rpc->Fill(GMTItr->etaValue(), phiconv_(GMTItr->phiValue()));
         }
       } else if (GMTItr->isRPC()) {
         if (GMTItr->quality() > 3) {
           eta_rpc_only->Fill(GMTItr->etaValue());
           phi_rpc_only->Fill(phiconv_(GMTItr->phiValue()));
           etaphi_rpc_only->Fill(GMTItr->etaValue(), phiconv_(GMTItr->phiValue()));
         }
       } else {
         if (GMTItr->quality() > 3) {
           eta_dtcsc_only->Fill(GMTItr->etaValue());
           phi_dtcsc_only->Fill(phiconv_(GMTItr->phiValue()));
           etaphi_dtcsc_only->Fill(GMTItr->etaValue(), phiconv_(GMTItr->phiValue()));
         }
  
         if (GMTItr != GMTCands.end()) {
           for (GMTItr2 = GMTCands.begin(); GMTItr2 != GMTCands.end(); ++GMTItr2) {
             if (GMTItr2 == GMTItr)
               continue;
             if (GMTItr2->empty())
               continue;
             if (GMTItr2->isRPC()) {
               if (GMTItr->isFwd()) {
                 dist_eta_csc_rpc->Fill(GMTItr->etaValue() - GMTItr2->etaValue());
                 dist_phi_csc_rpc->Fill(phiconv_(GMTItr->phiValue()) - phiconv_(GMTItr2->phiValue()));
               } else {
                 dist_eta_dt_rpc->Fill(GMTItr->etaValue() - GMTItr2->etaValue());
                 dist_phi_dt_rpc->Fill(phiconv_(GMTItr->phiValue()) - phiconv_(GMTItr2->phiValue()));
               }
             } else {
               if (!(GMTItr->isFwd()) && GMTItr2->isFwd()) {
                 dist_eta_dt_csc->Fill(GMTItr->etaValue() - GMTItr2->etaValue());
                 dist_phi_dt_csc->Fill(phiconv_(GMTItr->phiValue()) - phiconv_(GMTItr2->phiValue()));
               } else if (GMTItr->isFwd() && !(GMTItr2->isFwd())) {
                 dist_eta_dt_csc->Fill(GMTItr2->etaValue() - GMTItr->etaValue());
                 dist_phi_dt_csc->Fill(phiconv_(GMTItr->phiValue()) - phiconv_(GMTItr2->phiValue()));
               }
             }
           }
         }
       }
     }
  
     n_rpcb_vs_dttf->Fill(float(nSUBS[DTTF]), float(nSUBS[RPCb]));
     n_rpcf_vs_csctf->Fill(float(nSUBS[CSCTF]), float(nSUBS[RPCf]));
     n_csctf_vs_dttf->Fill(float(nSUBS[DTTF]), float(nSUBS[CSCTF]));
  
     regional_triggers->Fill(-1.);  // fill underflow for normalization
     if (nSUBS[GMT])
       regional_triggers->Fill(0.);  // fill all muon bin
     int ioff = 1;
     for (int i = 0; i < 4; i++) {
       if (nSUBS[i])
         regional_triggers->Fill(float(5 * i + nSUBS[i] + ioff));
     }
     if (nSUBS[DTTF] && (nSUBS[RPCb] || nSUBS[RPCf]))
       regional_triggers->Fill(22.);
     if (nSUBS[DTTF] && nSUBS[CSCTF])
       regional_triggers->Fill(23.);
     if (nSUBS[CSCTF] && (nSUBS[RPCb] || nSUBS[RPCf]))
       regional_triggers->Fill(24.);
     if (nSUBS[DTTF] && nSUBS[CSCTF] && (nSUBS[RPCb] || nSUBS[RPCf]))
       regional_triggers->Fill(25.);
  
     // fill only if previous event corresponds to previous trigger
     //    if( (Ev - evnum_old_) == 1 && bxnum_old_ > -1 ) {
     // assume getting all events in a sequence (usefull only from reco data)
     if (bxnum_old_ > -1) {
       int dBx = Bx - bxnum_old_ + 3564 * (e.orbitNumber() - obnum_old_);
       for (int id = 0; id < 4; id++) {
         if (trsrc_old_ & (1 << id)) {
           for (int i = 0; i < 4; i++) {
             if (nSUBS[i])
               subs_dbx[i]->Fill(float(dBx), float(id));
           }
         }
       }
     }
  
     // save quantities for the next event
     bxnum_old_ = Bx;
     obnum_old_ = e.orbitNumber();
     trsrc_old_ = 0;
     for (int i = 0; i < 4; i++) {
       if (nSUBS[i])
         trsrc_old_ |= (1 << i);
     }
   }
  
   if (bx1st[DTTF] < 9 && bx1st[RPCb] < 9)
     bx_dt_rpc->Fill(bx1st[DTTF], bx1st[RPCb]);
   if (bx1st[CSCTF] < 9 && bx1st[RPCf] < 9)
     bx_csc_rpc->Fill(bx1st[CSCTF], bx1st[RPCf]);
   if (bx1st[DTTF] < 9 && bx1st[CSCTF] < 9)
     bx_dt_csc->Fill(bx1st[DTTF], bx1st[CSCTF]);
 }

References bx_csc_rpc, bx_dt_csc, bx_dt_rpc, bx_number, bxnum_old_, gather_cfg::cout, CSCTF, dist_eta_csc_rpc, dist_eta_dt_csc, dist_eta_dt_rpc, dist_phi_csc_rpc, dist_phi_dt_csc, dist_phi_dt_rpc, DTTF, MillePedeFileConverter_cfg::e, eta_dtcsc_and_rpc, eta_dtcsc_only, eta_rpc_only, etaphi_dtcsc_and_rpc, etaphi_dtcsc_only, etaphi_rpc_only, dqm::impl::MonitorElement::Fill(), L1MuGMTReadoutCollection::getRecords(), GMT, gmtSource_, mps_fire::i, edm::HandleBase::isValid(), dqmiolumiharvest::j, n_csctf_vs_dttf, n_rpcb_vs_dttf, n_rpcf_vs_csctf, obnum_old_, phi_dtcsc_and_rpc, phi_dtcsc_only, phi_rpc_only, phiconv_(), edm::Handle< T >::product(), regional_triggers, RPCb, RPCf, subs_bits, subs_dbx, subs_eta, subs_etaphi, subs_etaqty, subs_nbx, subs_phi, subs_pt, subs_qty, trsrc_old_, and verbose_.

◆ book_()

void L1TGMT::book_ ( const edm::EventSetup & c )

private

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 241 of file L1TGMT.cc.

                                                                                            {
   std::string subs[5] = {"DTTF", "RPCb", "CSCTF", "RPCf", "GMT"};
  
   const L1MuTriggerScales* scales = &c.getData(l1muTrigscaleToken_);
   const L1MuTriggerPtScale* scalept = &c.getData(l1TrigptscaleToken_);
  
   ibooker.setCurrentFolder("L1T/L1TGMT");
  
   int nqty = 8;
   double qtymin = -0.5;
   double qtymax = 7.5;
  
   float phiscale[145];
   int nphiscale;
   {
     int nbins = scales->getPhiScale()->getNBins();
     if (nbins > 144)
       nbins = 144;
     for (int j = 0; j <= nbins; j++) {
       phiscale[j] = piconv_ * scales->getPhiScale()->getValue(j);
     }
     nphiscale = nbins;
   }
  
   float qscale[9];
   {
     for (int j = 0; j < 9; j++) {
       qscale[j] = -0.5 + j;
     }
   }
  
   // pt scale first bin reserved for empty muon
   float ptscale[32];
   int nptscale;
   {
     int nbins = scalept->getPtScale()->getNBins() - 1;
     if (nbins > 31)
       nbins = 31;
     for (int j = 1; j <= nbins; j++) {
       ptscale[j - 1] = scalept->getPtScale()->getValue(j);
     }
     ptscale[nbins] = ptscale[nbins - 1] + 10.;  // make reasonable size last bin
     nptscale = nbins;
   }
  
   float etascale[5][66];
   int netascale[5];
   // DTTF eta scale
   {
     int nbins = scales->getRegionalEtaScale(DTTF)->getNBins();
     if (nbins > 65)
       nbins = 65;
     for (int j = 0; j <= nbins; j++) {
       etascale[DTTF][j] = scales->getRegionalEtaScale(DTTF)->getValue(j);
     }
     netascale[DTTF] = nbins;
   }
   // RPCb etascale
   {
     int nbins = scales->getRegionalEtaScale(RPCb)->getNBins();
     if (nbins > 65)
       nbins = 65;
     for (int j = 0; j <= nbins; j++) {
       etascale[RPCb][j] = scales->getRegionalEtaScale(RPCb)->getValue(j);
     }
     netascale[RPCb] = nbins;
   }
   // CSCTF etascale
   // special case - need to mirror 2*32 bins
   {
     int nbins = scales->getRegionalEtaScale(CSCTF)->getNBins();
     if (nbins > 32)
       nbins = 32;
  
     int i = 0;
     for (int j = nbins; j >= 0; j--, i++) {
       etascale[CSCTF][i] = (-1) * scales->getRegionalEtaScale(CSCTF)->getValue(j);
     }
     for (int j = 0; j <= nbins; j++, i++) {
       etascale[CSCTF][i] = scales->getRegionalEtaScale(CSCTF)->getValue(j);
     }
     netascale[CSCTF] = i - 1;
   }
   // RPCf etascale
   {
     int nbins = scales->getRegionalEtaScale(RPCf)->getNBins();
     if (nbins > 65)
       nbins = 65;
     for (int j = 0; j <= nbins; j++) {
       etascale[RPCf][j] = scales->getRegionalEtaScale(RPCf)->getValue(j);
     }
     netascale[RPCf] = nbins;
   }
   // GMT etascale
   {
     int nbins = scales->getGMTEtaScale()->getNBins();
     if (nbins > 32)
       nbins = 32;
  
     int i = 0;
     for (int j = nbins; j > 0; j--, i++) {
       etascale[GMT][i] = (-1) * scales->getGMTEtaScale()->getValue(j);
     }
     for (int j = 0; j <= nbins; j++, i++) {
       etascale[GMT][i] = scales->getGMTEtaScale()->getValue(j);
     }
     netascale[GMT] = i - 1;
   }
  
   std::string hname("");
   std::string htitle("");
  
   for (int i = 0; i < 5; i++) {
     hname = subs[i] + "_nbx";
     htitle = subs[i] + " multiplicity in bx";
     subs_nbx[i] = ibooker.book2D(hname.data(), htitle.data(), 4, 1., 5., 5, -2.5, 2.5);
     subs_nbx[i]->setAxisTitle(subs[i] + " candidates", 1);
     subs_nbx[i]->setAxisTitle("bx wrt L1A", 2);
  
     hname = subs[i] + "_eta";
     htitle = subs[i] + " eta value";
     subs_eta[i] = ibooker.book1D(hname.data(), htitle.data(), netascale[i], etascale[i]);
     subs_eta[i]->setAxisTitle("eta", 1);
  
     hname = subs[i] + "_phi";
     htitle = subs[i] + " phi value";
     subs_phi[i] = ibooker.book1D(hname.data(), htitle.data(), nphiscale, phiscale);
     subs_phi[i]->setAxisTitle("phi (deg)", 1);
  
     hname = subs[i] + "_pt";
     htitle = subs[i] + " pt value";
     subs_pt[i] = ibooker.book1D(hname.data(), htitle.data(), nptscale, ptscale);
     subs_pt[i]->setAxisTitle("L1 pT (GeV)", 1);
  
     hname = subs[i] + "_qty";
     htitle = subs[i] + " qty value";
     subs_qty[i] = ibooker.book1D(hname.data(), htitle.data(), nqty, qtymin, qtymax);
     subs_qty[i]->setAxisTitle(subs[i] + " quality", 1);
  
     hname = subs[i] + "_etaphi";
     htitle = subs[i] + " phi vs eta";
     subs_etaphi[i] = ibooker.book2D(hname.data(), htitle.data(), netascale[i], etascale[i], nphiscale, phiscale);
     subs_etaphi[i]->setAxisTitle("eta", 1);
     subs_etaphi[i]->setAxisTitle("phi (deg)", 2);
  
     hname = subs[i] + "_etaqty";
     htitle = subs[i] + " qty vs eta";
     subs_etaqty[i] = ibooker.book2D(hname.data(), htitle.data(), netascale[i], etascale[i], nqty, qscale);
     subs_etaqty[i]->setAxisTitle("eta", 1);
     subs_etaqty[i]->setAxisTitle(subs[i] + " quality", 2);
  
     hname = subs[i] + "_bits";
     htitle = subs[i] + " bit population";
     subs_bits[i] = ibooker.book1D(hname.data(), htitle.data(), 32, -0.5, 31.5);
     subs_bits[i]->setAxisTitle("bit number", 1);
   }
  
   regional_triggers = ibooker.book1D("Regional_trigger", "Muon trigger contribution", 27, 0., 27.);
   regional_triggers->setAxisTitle("regional trigger", 1);
   int ib = 1;
   regional_triggers->setBinLabel(ib++, "All muons", 1);
   ib++;
   regional_triggers->setBinLabel(ib++, "DT 1mu", 1);
   regional_triggers->setBinLabel(ib++, "DT 2mu", 1);
   regional_triggers->setBinLabel(ib++, "DT 3mu", 1);
   regional_triggers->setBinLabel(ib++, "DT 4mu", 1);
   ib++;
   regional_triggers->setBinLabel(ib++, "RPCb 1mu", 1);
   regional_triggers->setBinLabel(ib++, "RPCb 2mu", 1);
   regional_triggers->setBinLabel(ib++, "RPCb 3mu", 1);
   regional_triggers->setBinLabel(ib++, "RPCb 4mu", 1);
   ib++;
   regional_triggers->setBinLabel(ib++, "CSC 1mu", 1);
   regional_triggers->setBinLabel(ib++, "CSC 2mu", 1);
   regional_triggers->setBinLabel(ib++, "CSC 3mu", 1);
   regional_triggers->setBinLabel(ib++, "CSC 4mu", 1);
   ib++;
   regional_triggers->setBinLabel(ib++, "RPCf 1mu", 1);
   regional_triggers->setBinLabel(ib++, "RPCf 2mu", 1);
   regional_triggers->setBinLabel(ib++, "RPCf 3mu", 1);
   regional_triggers->setBinLabel(ib++, "RPCf 4mu", 1);
   ib++;
   regional_triggers->setBinLabel(ib++, "DT & RPC", 1);
   regional_triggers->setBinLabel(ib++, "DT & CSC", 1);
   regional_triggers->setBinLabel(ib++, "CSC & RPC", 1);
   regional_triggers->setBinLabel(ib++, "DT & CSC & RPC", 1);
  
   bx_number = ibooker.book1D("Bx_Number", "Bx number ROP chip", 3564, 0., 3564.);
   bx_number->setAxisTitle("bx number", 1);
  
   dbx_chip = ibooker.bookProfile("dbx_Chip", "bx count difference wrt ROP chip", 5, 0., 5., 100, -4000., 4000., "i");
   dbx_chip->setAxisTitle("chip name", 1);
   dbx_chip->setAxisTitle("delta bx", 2);
   dbx_chip->setBinLabel(1, "IND", 1);
   dbx_chip->setBinLabel(2, "INB", 1);
   dbx_chip->setBinLabel(3, "INC", 1);
   dbx_chip->setBinLabel(4, "INF", 1);
   dbx_chip->setBinLabel(5, "SRT", 1);
  
   eta_dtcsc_and_rpc =
       ibooker.book1D("eta_DTCSC_and_RPC", "eta of confirmed GMT candidates", netascale[GMT], etascale[GMT]);
   eta_dtcsc_and_rpc->setAxisTitle("eta", 1);
  
   eta_dtcsc_only =
       ibooker.book1D("eta_DTCSC_only", "eta of unconfirmed DT/CSC candidates", netascale[GMT], etascale[GMT]);
   eta_dtcsc_only->setAxisTitle("eta", 1);
  
   eta_rpc_only = ibooker.book1D("eta_RPC_only", "eta of unconfirmed RPC candidates", netascale[GMT], etascale[GMT]);
   eta_rpc_only->setAxisTitle("eta", 1);
  
   phi_dtcsc_and_rpc = ibooker.book1D("phi_DTCSC_and_RPC", "phi of confirmed GMT candidates", nphiscale, phiscale);
   phi_dtcsc_and_rpc->setAxisTitle("phi (deg)", 1);
  
   phi_dtcsc_only = ibooker.book1D("phi_DTCSC_only", "phi of unconfirmed DT/CSC candidates", nphiscale, phiscale);
   phi_dtcsc_only->setAxisTitle("phi (deg)", 1);
  
   phi_rpc_only = ibooker.book1D("phi_RPC_only", "phi of unconfirmed RPC candidates", nphiscale, phiscale);
   phi_rpc_only->setAxisTitle("phi (deg)", 1);
  
   etaphi_dtcsc_and_rpc = ibooker.book2D("etaphi_DTCSC_and_RPC",
                                         "eta vs phi map of confirmed GMT candidates",
                                         netascale[GMT],
                                         etascale[GMT],
                                         nphiscale,
                                         phiscale);
   etaphi_dtcsc_and_rpc->setAxisTitle("eta", 1);
   etaphi_dtcsc_and_rpc->setAxisTitle("phi (deg)", 2);
  
   etaphi_dtcsc_only = ibooker.book2D("etaphi_DTCSC_only",
                                      "eta vs phi map of unconfirmed DT/CSC candidates",
                                      netascale[GMT],
                                      etascale[GMT],
                                      nphiscale,
                                      phiscale);
   etaphi_dtcsc_only->setAxisTitle("eta", 1);
   etaphi_dtcsc_only->setAxisTitle("phi (deg)", 2);
  
   etaphi_rpc_only = ibooker.book2D("etaphi_RPC_only",
                                    "eta vs phi map of unconfirmed RPC candidates",
                                    netascale[GMT],
                                    etascale[GMT],
                                    nphiscale,
                                    phiscale);
   etaphi_rpc_only->setAxisTitle("eta", 1);
   etaphi_rpc_only->setAxisTitle("phi (deg)", 2);
  
   dist_phi_dt_rpc = ibooker.book1D("dist_phi_DT_RPC", "Dphi between DT and RPC candidates", 100, -125., 125.);
   dist_phi_dt_rpc->setAxisTitle("delta phi (deg)", 1);
  
   dist_phi_csc_rpc = ibooker.book1D("dist_phi_CSC_RPC", "Dphi between CSC and RPC candidates", 100, -125., 125.);
   dist_phi_csc_rpc->setAxisTitle("delta phi (deg)", 1);
  
   dist_phi_dt_csc = ibooker.book1D("dist_phi_DT_CSC", "Dphi between DT and CSC candidates", 100, -125., 125.);
   dist_phi_dt_csc->setAxisTitle("delta phi (deg)", 1);
  
   dist_eta_dt_rpc = ibooker.book1D("dist_eta_DT_RPC", "Deta between DT and RPC candidates", 40, -1., 1.);
   dist_eta_dt_rpc->setAxisTitle("delta eta", 1);
  
   dist_eta_csc_rpc = ibooker.book1D("dist_eta_CSC_RPC", "Deta between CSC and RPC candidates", 40, -1., 1.);
   dist_eta_csc_rpc->setAxisTitle("delta eta", 1);
  
   dist_eta_dt_csc = ibooker.book1D("dist_eta_DT_CSC", "Deta between DT and CSC candidates", 40, -1., 1.);
   dist_eta_dt_csc->setAxisTitle("delta eta", 1);
  
   n_rpcb_vs_dttf = ibooker.book2D("n_RPCb_vs_DTTF", "n cands RPCb vs DTTF", 5, -0.5, 4.5, 5, -0.5, 4.5);
   n_rpcb_vs_dttf->setAxisTitle("DTTF candidates", 1);
   n_rpcb_vs_dttf->setAxisTitle("barrel RPC candidates", 2);
  
   n_rpcf_vs_csctf = ibooker.book2D("n_RPCf_vs_CSCTF", "n cands RPCf vs CSCTF", 5, -0.5, 4.5, 5, -0.5, 4.5);
   n_rpcf_vs_csctf->setAxisTitle("CSCTF candidates", 1);
   n_rpcf_vs_csctf->setAxisTitle("endcap RPC candidates", 2);
  
   n_csctf_vs_dttf = ibooker.book2D("n_CSCTF_vs_DTTF", "n cands CSCTF vs DTTF", 5, -0.5, 4.5, 5, -0.5, 4.5);
   n_csctf_vs_dttf->setAxisTitle("DTTF candidates", 1);
   n_csctf_vs_dttf->setAxisTitle("CSCTF candidates", 2);
  
   bx_dt_rpc = ibooker.book2D("bx_DT_vs_RPC", "1st bx DT vs. RPC", 5, -2.5, 2.5, 5, -2.5, 2.5);
   bx_dt_rpc->setAxisTitle("bx of 1st DTTF candidate", 1);
   bx_dt_rpc->setAxisTitle("bx of 1st RPCb candidate", 2);
  
   bx_csc_rpc = ibooker.book2D("bx_CSC_vs_RPC", "1st bx CSC vs. RPC", 5, -2.5, 2.5, 5, -2.5, 2.5);
   bx_csc_rpc->setAxisTitle("bx of 1st CSCTF candidate", 1);
   bx_csc_rpc->setAxisTitle("bx of 1st RPCf candidate", 2);
  
   bx_dt_csc = ibooker.book2D("bx_DT_vs_CSC", "1st bx DT vs. CSC", 5, -2.5, 2.5, 5, -2.5, 2.5);
   bx_dt_csc->setAxisTitle("bx of 1st DTTF candidate", 1);
   bx_dt_csc->setAxisTitle("bx of 1st CSCTF candidate", 2);
  
   for (int i = 0; i < 4; i++) {
     hname = subs[i] + "_dbx";
     htitle = "dBx " + subs[i] + " to previous event";
     subs_dbx[i] = ibooker.book2D(hname.data(), htitle.data(), 1000, 0., 1000., 4, 0., 4.);
     for (int j = 0; j < 4; j++) {
       subs_dbx[i]->setBinLabel((j + 1), subs[j], 2);
     }
   }
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), bx_csc_rpc, bx_dt_csc, bx_dt_rpc, bx_number, c, CSCTF, dbx_chip, dist_eta_csc_rpc, dist_eta_dt_csc, dist_eta_dt_rpc, dist_phi_csc_rpc, dist_phi_dt_csc, dist_phi_dt_rpc, DTTF, eta_dtcsc_and_rpc, eta_dtcsc_only, eta_rpc_only, etaphi_dtcsc_and_rpc, etaphi_dtcsc_only, etaphi_rpc_only, L1MuTriggerScales::getGMTEtaScale(), L1MuScale::getNBins(), L1MuTriggerScales::getPhiScale(), L1MuTriggerPtScale::getPtScale(), L1MuTriggerScales::getRegionalEtaScale(), L1MuScale::getValue(), GMT, mps_fire::i, cuy::ib, dqmiolumiharvest::j, l1muTrigscaleToken_, l1TrigptscaleToken_, n_csctf_vs_dttf, n_rpcb_vs_dttf, n_rpcf_vs_csctf, LaserClient_cfi::nbins, phi_dtcsc_and_rpc, phi_dtcsc_only, phi_rpc_only, piconv_, emtf::ptscale, regional_triggers, RPCb, RPCf, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, subs_bits, subs_dbx, subs_eta, subs_etaphi, subs_etaqty, subs_nbx, subs_phi, subs_pt, and subs_qty.

◆ phiconv_()

double L1TGMT::phiconv_ ( float phi )

private

Definition at line 234 of file L1TGMT.cc.

                                  {
   double phiout = double(phi);
   phiout *= piconv_;
   phiout += 0.001;  // add a small value to get off the bin edge
   return phiout;
 }

References phi, and piconv_.

Referenced by analyze().

Member Data Documentation

◆ bx_csc_rpc

MonitorElement* L1TGMT::bx_csc_rpc

private

Definition at line 94 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ bx_dt_csc

MonitorElement* L1TGMT::bx_dt_csc

private

Definition at line 95 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ bx_dt_rpc

MonitorElement* L1TGMT::bx_dt_rpc

private

Definition at line 93 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ bx_number

MonitorElement* L1TGMT::bx_number

private

Definition at line 76 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ bxnum_old_

int L1TGMT::bxnum_old_

private

Definition at line 109 of file L1TGMT.h.

Referenced by analyze().

◆ dbx_chip

MonitorElement* L1TGMT::dbx_chip

private

Definition at line 77 of file L1TGMT.h.

Referenced by bookHistograms().

◆ dist_eta_csc_rpc

MonitorElement* L1TGMT::dist_eta_csc_rpc

private

Definition at line 91 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ dist_eta_dt_csc

MonitorElement* L1TGMT::dist_eta_dt_csc

private

Definition at line 92 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ dist_eta_dt_rpc

MonitorElement* L1TGMT::dist_eta_dt_rpc

private

Definition at line 90 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ dist_phi_csc_rpc

MonitorElement* L1TGMT::dist_phi_csc_rpc

private

Definition at line 88 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ dist_phi_dt_csc

MonitorElement* L1TGMT::dist_phi_dt_csc

private

Definition at line 89 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ dist_phi_dt_rpc

MonitorElement* L1TGMT::dist_phi_dt_rpc

private

Definition at line 87 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ eta_dtcsc_and_rpc

MonitorElement* L1TGMT::eta_dtcsc_and_rpc

private

Definition at line 78 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ eta_dtcsc_only

MonitorElement* L1TGMT::eta_dtcsc_only

private

Definition at line 79 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ eta_rpc_only

MonitorElement* L1TGMT::eta_rpc_only

private

Definition at line 80 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ etaphi_dtcsc_and_rpc

MonitorElement* L1TGMT::etaphi_dtcsc_and_rpc

private

Definition at line 84 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ etaphi_dtcsc_only

MonitorElement* L1TGMT::etaphi_dtcsc_only

private

Definition at line 85 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ etaphi_rpc_only

MonitorElement* L1TGMT::etaphi_rpc_only

private

Definition at line 86 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ gmtSource_

const edm::EDGetTokenT<L1MuGMTReadoutCollection> L1TGMT::gmtSource_

private

Definition at line 105 of file L1TGMT.h.

Referenced by analyze().

◆ l1muTrigscaleToken_

edm::ESGetToken<L1MuTriggerScales, L1MuTriggerScalesRcd> L1TGMT::l1muTrigscaleToken_

private

Definition at line 106 of file L1TGMT.h.

Referenced by bookHistograms(), and L1TGMT().

◆ l1TrigptscaleToken_

edm::ESGetToken<L1MuTriggerPtScale, L1MuTriggerPtScaleRcd> L1TGMT::l1TrigptscaleToken_

private

Definition at line 107 of file L1TGMT.h.

Referenced by bookHistograms(), and L1TGMT().

◆ logFile_

std::ofstream L1TGMT::logFile_

private

Definition at line 104 of file L1TGMT.h.

◆ n_csctf_vs_dttf

MonitorElement* L1TGMT::n_csctf_vs_dttf

private

Definition at line 99 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ n_rpcb_vs_dttf

MonitorElement* L1TGMT::n_rpcb_vs_dttf

private

Definition at line 97 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ n_rpcf_vs_csctf

MonitorElement* L1TGMT::n_rpcf_vs_csctf

private

Definition at line 98 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ obnum_old_

int L1TGMT::obnum_old_

private

Definition at line 110 of file L1TGMT.h.

Referenced by analyze().

◆ phi_dtcsc_and_rpc

MonitorElement* L1TGMT::phi_dtcsc_and_rpc

private

Definition at line 81 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ phi_dtcsc_only

MonitorElement* L1TGMT::phi_dtcsc_only

private

Definition at line 82 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ phi_rpc_only

MonitorElement* L1TGMT::phi_rpc_only

private

Definition at line 83 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ piconv_

const double L1TGMT::piconv_ = 180. / acos(-1.)

staticprivate

Definition at line 113 of file L1TGMT.h.

Referenced by bookHistograms(), and phiconv_().

◆ regional_triggers

MonitorElement* L1TGMT::regional_triggers

private

Definition at line 74 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_bits

MonitorElement* L1TGMT::subs_bits[5]

private

Definition at line 72 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_dbx

MonitorElement* L1TGMT::subs_dbx[4]

private

Definition at line 101 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_eta

MonitorElement* L1TGMT::subs_eta[5]

private

Definition at line 66 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_etaphi

MonitorElement* L1TGMT::subs_etaphi[5]

private

Definition at line 70 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_etaqty

MonitorElement* L1TGMT::subs_etaqty[5]

private

Definition at line 71 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_nbx

MonitorElement* L1TGMT::subs_nbx[5]

private

Definition at line 65 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_phi

MonitorElement* L1TGMT::subs_phi[5]

private

Definition at line 67 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_pt

MonitorElement* L1TGMT::subs_pt[5]

private

Definition at line 68 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ subs_qty

MonitorElement* L1TGMT::subs_qty[5]

private

Definition at line 69 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

◆ trsrc_old_

int L1TGMT::trsrc_old_

private

Definition at line 111 of file L1TGMT.h.

Referenced by analyze().

◆ verbose_

const bool L1TGMT::verbose_

private

Definition at line 103 of file L1TGMT.h.

Referenced by analyze(), and L1TGMT().

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ ensubs

Constructor & Destructor Documentation

◆ L1TGMT()

◆ ~L1TGMT()

Member Function Documentation

◆ analyze()

◆ book_()

◆ bookHistograms()

◆ phiconv_()

Member Data Documentation

◆ bx_csc_rpc

◆ bx_dt_csc

◆ bx_dt_rpc

◆ bx_number

◆ bxnum_old_

◆ dbx_chip

◆ dist_eta_csc_rpc

◆ dist_eta_dt_csc

◆ dist_eta_dt_rpc

◆ dist_phi_csc_rpc

◆ dist_phi_dt_csc

◆ dist_phi_dt_rpc

◆ eta_dtcsc_and_rpc

◆ eta_dtcsc_only

◆ eta_rpc_only

◆ etaphi_dtcsc_and_rpc

◆ etaphi_dtcsc_only

◆ etaphi_rpc_only

◆ gmtSource_

◆ l1muTrigscaleToken_

◆ l1TrigptscaleToken_

◆ logFile_

◆ n_csctf_vs_dttf

◆ n_rpcb_vs_dttf

◆ n_rpcf_vs_csctf

◆ obnum_old_

◆ phi_dtcsc_and_rpc

◆ phi_dtcsc_only

◆ phi_rpc_only

◆ piconv_

◆ regional_triggers

◆ subs_bits

◆ subs_dbx

◆ subs_eta

◆ subs_etaphi

◆ subs_etaqty

◆ subs_nbx

◆ subs_phi

◆ subs_pt

◆ subs_qty

◆ trsrc_old_

◆ verbose_