L1TGMT Class Reference

#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

enum L1TGMT::ensubs

private

Enumerator
DTTF
RPCb
CSCTF
RPCf
GMT

Definition at line 63 of file L1TGMT.h.

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

Constructor & Destructor Documentation

L1TGMT::L1TGMT

(

const edm::ParameterSet &

ps

)

Definition at line 23 of file L1TGMT.cc.

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

    : 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>();
}

L1TGMT::~L1TGMT ( )

override

Definition at line 36 of file L1TGMT.cc.

{}

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 38 of file L1TGMT.cc.

References SplitLinear::begin, 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, eta_dtcsc_and_rpc, eta_dtcsc_only, eta_rpc_only, etaphi_dtcsc_and_rpc, etaphi_dtcsc_only, etaphi_rpc_only, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), 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_, edm::EventBase::orbitNumber(), 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_.

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

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

private

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 241 of file L1TGMT.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), bx_csc_rpc, bx_dt_csc, bx_dt_rpc, bx_number, 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, edm::EventSetup::getData(), 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, hlt_dqm_clientPB-live_cfg::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.

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

double L1TGMT::phiconv_ ( float  phi )

private

Definition at line 234 of file L1TGMT.cc.

References piconv_.

Referenced by analyze().

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

Member Data Documentation

MonitorElement* L1TGMT::bx_csc_rpc

private

Definition at line 94 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::bx_dt_csc

private

Definition at line 95 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::bx_dt_rpc

private

Definition at line 93 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::bx_number

private

Definition at line 76 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

int L1TGMT::bxnum_old_

private

Definition at line 109 of file L1TGMT.h.

Referenced by analyze().

MonitorElement* L1TGMT::dbx_chip

private

Definition at line 77 of file L1TGMT.h.

Referenced by bookHistograms().

MonitorElement* L1TGMT::dist_eta_csc_rpc

private

Definition at line 91 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::dist_eta_dt_csc

private

Definition at line 92 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::dist_eta_dt_rpc

private

Definition at line 90 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::dist_phi_csc_rpc

private

Definition at line 88 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::dist_phi_dt_csc

private

Definition at line 89 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::dist_phi_dt_rpc

private

Definition at line 87 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::eta_dtcsc_and_rpc

private

Definition at line 78 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::eta_dtcsc_only

private

Definition at line 79 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::eta_rpc_only

private

Definition at line 80 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::etaphi_dtcsc_and_rpc

private

Definition at line 84 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::etaphi_dtcsc_only

private

Definition at line 85 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::etaphi_rpc_only

private

Definition at line 86 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 105 of file L1TGMT.h.

Referenced by analyze().

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

private

Definition at line 106 of file L1TGMT.h.

Referenced by bookHistograms(), and L1TGMT().

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

private

Definition at line 107 of file L1TGMT.h.

Referenced by bookHistograms(), and L1TGMT().

std::ofstream L1TGMT::logFile_

private

Definition at line 104 of file L1TGMT.h.

MonitorElement* L1TGMT::n_csctf_vs_dttf

private

Definition at line 99 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::n_rpcb_vs_dttf

private

Definition at line 97 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::n_rpcf_vs_csctf

private

Definition at line 98 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

int L1TGMT::obnum_old_

private

Definition at line 110 of file L1TGMT.h.

Referenced by analyze().

MonitorElement* L1TGMT::phi_dtcsc_and_rpc

private

Definition at line 81 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::phi_dtcsc_only

private

Definition at line 82 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::phi_rpc_only

private

Definition at line 83 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

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

staticprivate

Definition at line 113 of file L1TGMT.h.

Referenced by bookHistograms(), and phiconv_().

MonitorElement* L1TGMT::regional_triggers

private

Definition at line 74 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_bits[5]

private

Definition at line 72 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_dbx[4]

private

Definition at line 101 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_eta[5]

private

Definition at line 66 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_etaphi[5]

private

Definition at line 70 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_etaqty[5]

private

Definition at line 71 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_nbx[5]

private

Definition at line 65 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_phi[5]

private

Definition at line 67 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_pt[5]

private

Definition at line 68 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGMT::subs_qty[5]

private

Definition at line 69 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

int L1TGMT::trsrc_old_

private

Definition at line 111 of file L1TGMT.h.

Referenced by analyze().

const bool L1TGMT::verbose_

private

Definition at line 103 of file L1TGMT.h.

Referenced by analyze(), and L1TGMT().