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 45 of file L1TGMT.h.

Member Enumeration Documentation

ensubs

enum L1TGMT::ensubs

private

Enumerator
DTTF
RPCb
CSCTF
RPCf
GMT

Definition at line 62 of file L1TGMT.h.

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

Constructor & Destructor Documentation

L1TGMT()

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

override

Definition at line 36 of file L1TGMT.cc.

{}

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.

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

                                                        {
  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) {
      float dBx = Bx - bxnum_old_ + 3564.0 * (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(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]);
}

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.

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

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

phiconv_()

double L1TGMT::phiconv_ ( float  phi )

private

Definition at line 234 of file L1TGMT.cc.

References phi, and 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

bx_csc_rpc

MonitorElement* L1TGMT::bx_csc_rpc

private

Definition at line 93 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

bx_dt_csc

MonitorElement* L1TGMT::bx_dt_csc

private

Definition at line 94 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

bx_dt_rpc

MonitorElement* L1TGMT::bx_dt_rpc

private

Definition at line 92 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

bx_number

MonitorElement* L1TGMT::bx_number

private

Definition at line 75 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

bxnum_old_

int L1TGMT::bxnum_old_

private

Definition at line 108 of file L1TGMT.h.

Referenced by analyze().

dbx_chip

MonitorElement* L1TGMT::dbx_chip

private

Definition at line 76 of file L1TGMT.h.

Referenced by bookHistograms().

dist_eta_csc_rpc

MonitorElement* L1TGMT::dist_eta_csc_rpc

private

Definition at line 90 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

dist_eta_dt_csc

MonitorElement* L1TGMT::dist_eta_dt_csc

private

Definition at line 91 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

dist_eta_dt_rpc

MonitorElement* L1TGMT::dist_eta_dt_rpc

private

Definition at line 89 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

dist_phi_csc_rpc

MonitorElement* L1TGMT::dist_phi_csc_rpc

private

Definition at line 87 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

dist_phi_dt_csc

MonitorElement* L1TGMT::dist_phi_dt_csc

private

Definition at line 88 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

dist_phi_dt_rpc

MonitorElement* L1TGMT::dist_phi_dt_rpc

private

Definition at line 86 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

eta_dtcsc_and_rpc

MonitorElement* L1TGMT::eta_dtcsc_and_rpc

private

Definition at line 77 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

eta_dtcsc_only

MonitorElement* L1TGMT::eta_dtcsc_only

private

Definition at line 78 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

eta_rpc_only

MonitorElement* L1TGMT::eta_rpc_only

private

Definition at line 79 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

etaphi_dtcsc_and_rpc

MonitorElement* L1TGMT::etaphi_dtcsc_and_rpc

private

Definition at line 83 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

etaphi_dtcsc_only

MonitorElement* L1TGMT::etaphi_dtcsc_only

private

Definition at line 84 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

etaphi_rpc_only

MonitorElement* L1TGMT::etaphi_rpc_only

private

Definition at line 85 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

gmtSource_

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

private

Definition at line 104 of file L1TGMT.h.

Referenced by analyze().

l1muTrigscaleToken_

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

private

Definition at line 105 of file L1TGMT.h.

Referenced by bookHistograms(), and L1TGMT().

l1TrigptscaleToken_

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

private

Definition at line 106 of file L1TGMT.h.

Referenced by bookHistograms(), and L1TGMT().

logFile_

std::ofstream L1TGMT::logFile_

private

Definition at line 103 of file L1TGMT.h.

n_csctf_vs_dttf

MonitorElement* L1TGMT::n_csctf_vs_dttf

private

Definition at line 98 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

n_rpcb_vs_dttf

MonitorElement* L1TGMT::n_rpcb_vs_dttf

private

Definition at line 96 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

n_rpcf_vs_csctf

MonitorElement* L1TGMT::n_rpcf_vs_csctf

private

Definition at line 97 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

obnum_old_

int L1TGMT::obnum_old_

private

Definition at line 109 of file L1TGMT.h.

Referenced by analyze().

phi_dtcsc_and_rpc

MonitorElement* L1TGMT::phi_dtcsc_and_rpc

private

Definition at line 80 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

phi_dtcsc_only

MonitorElement* L1TGMT::phi_dtcsc_only

private

Definition at line 81 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

phi_rpc_only

MonitorElement* L1TGMT::phi_rpc_only

private

Definition at line 82 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

piconv_

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

staticprivate

Definition at line 112 of file L1TGMT.h.

Referenced by bookHistograms(), and phiconv_().

regional_triggers

MonitorElement* L1TGMT::regional_triggers

private

Definition at line 73 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_bits

MonitorElement* L1TGMT::subs_bits[5]

private

Definition at line 71 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_dbx

MonitorElement* L1TGMT::subs_dbx[4]

private

Definition at line 100 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_eta

MonitorElement* L1TGMT::subs_eta[5]

private

Definition at line 65 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_etaphi

MonitorElement* L1TGMT::subs_etaphi[5]

private

Definition at line 69 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_etaqty

MonitorElement* L1TGMT::subs_etaqty[5]

private

Definition at line 70 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_nbx

MonitorElement* L1TGMT::subs_nbx[5]

private

Definition at line 64 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_phi

MonitorElement* L1TGMT::subs_phi[5]

private

Definition at line 66 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_pt

MonitorElement* L1TGMT::subs_pt[5]

private

Definition at line 67 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

subs_qty

MonitorElement* L1TGMT::subs_qty[5]

private

Definition at line 68 of file L1TGMT.h.

Referenced by analyze(), and bookHistograms().

trsrc_old_

int L1TGMT::trsrc_old_

private

Definition at line 110 of file L1TGMT.h.

Referenced by analyze().

verbose_

const bool L1TGMT::verbose_

private

Definition at line 102 of file L1TGMT.h.

Referenced by analyze(), and L1TGMT().