L1TDTTPG Class Reference

#include <L1TDTTPG.h>

Inheritance diagram for L1TDTTPG:

Public Member Functions
	L1TDTTPG (const edm::ParameterSet &ps)
	~L1TDTTPG () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Member Functions
void	setMapPhLabel (MonitorElement *me)
void	setMapThLabel (MonitorElement *me)

Private Attributes
MonitorElement *	dttf_p_phi [3]
MonitorElement *	dttf_p_pt [3]
MonitorElement *	dttf_p_q [3]
MonitorElement *	dttf_p_qual [3]
MonitorElement *	dttpgphbestmap
MonitorElement *	dttpgphbestmapcorr
MonitorElement *	dttpgphbx [8]
MonitorElement *	dttpgphbxcomp
MonitorElement *	dttpgphmap
MonitorElement *	dttpgphmap2nd
MonitorElement *	dttpgphmapbx [3]
MonitorElement *	dttpgphmapcorr
MonitorElement *	dttpgphntrack
MonitorElement *	dttpgphquality [3]
MonitorElement *	dttpgphsector [3]
MonitorElement *	dttpgphstation [3]
MonitorElement *	dttpgphts2tag [3]
MonitorElement *	dttpgphwheel [3]
edm::InputTag	dttpgSource_
edm::EDGetTokenT< L1MuDTChambPhContainer >	dttpgSourcePhContainer_token_
edm::EDGetTokenT< L1MuDTChambThContainer >	dttpgSourceThContainer_token_
MonitorElement *	dttpgthbestmap
MonitorElement *	dttpgthbestmaph
MonitorElement *	dttpgthbx [3]
MonitorElement *	dttpgthmap
MonitorElement *	dttpgthmapbx [3]
MonitorElement *	dttpgthmaph
MonitorElement *	dttpgthntrack
MonitorElement *	dttpgthquality [3]
MonitorElement *	dttpgthsector [3]
MonitorElement *	dttpgthstation [3]
MonitorElement *	dttpgththeta [3]
MonitorElement *	dttpgthwheel [3]
std::ofstream	logFile_
bool	monitorDaemon_
int	nev_
std::string	outputFile_
std::string	trstring_
edm::EDGetTokenT< L1MuDTTrackContainer >	trToken_
bool	verbose_

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

Detailed Description

Definition at line 39 of file L1TDTTPG.h.

Constructor & Destructor Documentation

L1TDTTPG::L1TDTTPG

(

const edm::ParameterSet &

ps

)

Definition at line 93 of file L1TDTTPG.cc.

References gather_cfg::cout, dttpgSource_, edm::ParameterSet::getUntrackedParameter(), edm::InputTag::label(), outputFile_, edm::InputTag::process(), trstring_, trToken_, and verbose_.

    : dttpgSourcePhContainer_token_(consumes<L1MuDTChambPhContainer>(ps.getParameter<InputTag>("dttpgSource"))),
      dttpgSourceThContainer_token_(consumes<L1MuDTChambThContainer>(ps.getParameter<InputTag>("dttpgSource"))),
      dttpgSource_(ps.getParameter<InputTag>("dttpgSource")) {
  trstring_ = dttpgSource_.label() + ":" + "DATA" + ":" + dttpgSource_.process();
  trToken_ = consumes<L1MuDTTrackContainer>(trstring_);

  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);

  if (verbose_)
    cout << "L1TDTTPG: constructor...." << endl;

  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
  if (!outputFile_.empty()) {
    cout << "L1T Monitoring histograms will be saved to " << outputFile_.c_str() << endl;
  }

  bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
  if (disable) {
    outputFile_ = "";
  }
}

L1TDTTPG::~L1TDTTPG ( )

override

Definition at line 117 of file L1TDTTPG.cc.

{}

Member Function Documentation

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

overrideprotected

Definition at line 248 of file L1TDTTPG.cc.

References l1GtPatternGenerator_cfi::bx, L1MuDTChambPhDigi::code(), gather_cfg::cout, dttf_p_phi, dttf_p_pt, dttf_p_q, dttf_p_qual, dttpgphbestmap, dttpgphbestmapcorr, dttpgphbx, dttpgphbxcomp, dttpgphmap, dttpgphmap2nd, dttpgphmapbx, dttpgphmapcorr, dttpgphntrack, dttpgphquality, dttpgphsector, dttpgphstation, dttpgphts2tag, dttpgphwheel, dttpgSource_, dttpgSourcePhContainer_token_, dttpgSourceThContainer_token_, dttpgthbestmap, dttpgthbestmaph, dttpgthbx, dttpgthmap, dttpgthmapbx, dttpgthmaph, dttpgthntrack, dttpgthquality, dttpgthsector, dttpgthstation, dttpgththeta, dttpgthwheel, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), L1MuDTChambPhContainer::getContainer(), L1MuDTChambThContainer::getContainer(), L1MuDTTrackContainer::getContainer(), mps_fire::i, edm::HandleBase::isValid(), dqmiolumiharvest::j, dqmdumpme::k, edm::InputTag::label(), nev_, OrderedSet::t, trToken_, and verbose_.

                                                          {
  nev_++;
  if (verbose_)
    cout << "L1TDTTPG: analyze...." << endl;

  edm::Handle<L1MuDTChambPhContainer> myL1MuDTChambPhContainer;
  e.getByToken(dttpgSourcePhContainer_token_, myL1MuDTChambPhContainer);

  if (!myL1MuDTChambPhContainer.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuDTChambPhContainer with label " << dttpgSource_.label();
    return;
  }
  L1MuDTChambPhContainer::Phi_Container const *myPhContainer = myL1MuDTChambPhContainer->getContainer();

  edm::Handle<L1MuDTChambThContainer> myL1MuDTChambThContainer;
  e.getByToken(dttpgSourceThContainer_token_, myL1MuDTChambThContainer);

  if (!myL1MuDTChambThContainer.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuDTChambThContainer with label " << dttpgSource_.label();
    edm::LogInfo("DataNotFound") << "if this fails try to add DATA to the process name.";

    return;
  }
  L1MuDTChambThContainer::The_Container const *myThContainer = myL1MuDTChambThContainer->getContainer();

  int ndttpgphtrack = 0;
  int ndttpgthtrack = 0;
  int NumberOfSegmentsPhi[3] = {0, 0, 0};

  for (L1MuDTChambPhContainer::Phi_Container::const_iterator DTPhDigiItr = myPhContainer->begin();
       DTPhDigiItr != myPhContainer->end();
       ++DTPhDigiItr) {
    int bx = DTPhDigiItr->bxNum() - DTPhDigiItr->Ts2Tag();
    if (bx == -1)
      NumberOfSegmentsPhi[0]++;
    if (bx == 0)
      NumberOfSegmentsPhi[1]++;
    if (bx == 1)
      NumberOfSegmentsPhi[2]++;
  }
  /*Fill Histos for Segment counter for each bx separately */

  for (int k = 0; k < 3; k++) {
    dttpgphbx[k + 2]->Fill(NumberOfSegmentsPhi[k]);
  }
  int bxCounterDttfPhi = 0;  // = no. of BX's with non-zero data
  for (int k = 0; k < 3; k++) {
    if (NumberOfSegmentsPhi[k] > 0)
      bxCounterDttfPhi++;
  }

  /* the BX "code" */

  int bxCodePhi = 0;
  if (bxCounterDttfPhi == 0) {
    bxCodePhi = 0;
  } else if (bxCounterDttfPhi == 1) {
    for (int k = 0; k < 3; k++) {
      if (NumberOfSegmentsPhi[k] > 0)
        bxCodePhi = k + 2;
    }
  } else if (bxCounterDttfPhi == 2) {
    for (int k = 0; k < 3; k++) {
      if (NumberOfSegmentsPhi[k] == 0)
        bxCodePhi = 8 - k;
    }
  } else if (bxCounterDttfPhi == 3) {
    bxCodePhi = 10;
  }

  //The bx analyzer histo
  dttpgphbx[0]->Fill(bxCodePhi);

  const L1MuDTChambPhDigi *bestPhQualMap[5][12][4];
  memset(bestPhQualMap, 0, 240 * sizeof(L1MuDTChambPhDigi *));

  for (L1MuDTChambPhContainer::Phi_Container::const_iterator DTPhDigiItr = myPhContainer->begin();
       DTPhDigiItr != myPhContainer->end();
       ++DTPhDigiItr) {
    ndttpgphtrack++;

    int bxindex = DTPhDigiItr->bxNum() - DTPhDigiItr->Ts2Tag() + 1;

    dttpgphwheel[bxindex]->Fill(DTPhDigiItr->whNum());
    if (verbose_) {
      cout << "DTTPG phi wheel number " << DTPhDigiItr->whNum() << endl;
    }
    dttpgphstation[bxindex]->Fill(DTPhDigiItr->stNum());
    if (verbose_) {
      cout << "DTTPG phi station number " << DTPhDigiItr->stNum() << endl;
    }
    dttpgphsector[bxindex]->Fill(DTPhDigiItr->scNum());
    if (verbose_) {
      cout << "DTTPG phi sector number " << DTPhDigiItr->scNum() << endl;
    }
    dttpgphquality[bxindex]->Fill(DTPhDigiItr->code());
    if (verbose_) {
      cout << "DTTPG phi quality " << DTPhDigiItr->code() << endl;
    }
    dttpgphts2tag[bxindex]->Fill(DTPhDigiItr->Ts2Tag());
    if (verbose_) {
      cout << "DTTPG phi ts2tag " << DTPhDigiItr->Ts2Tag() << endl;
    }
    int ypos = DTPhDigiItr->scNum();
    int xpos = DTPhDigiItr->stNum() + 4 * (DTPhDigiItr->whNum() + 2);
    dttpgphmap->Fill(xpos, ypos);
    if (DTPhDigiItr->Ts2Tag())
      dttpgphmap2nd->Fill(xpos, ypos);
    dttpgphmapbx[bxindex]->Fill(xpos, ypos);
    if (DTPhDigiItr->code() > 3)
      dttpgphmapcorr->Fill(xpos, ypos);

    if (bestPhQualMap[DTPhDigiItr->whNum() + 2][DTPhDigiItr->scNum()][DTPhDigiItr->stNum() - 1] == nullptr ||
        bestPhQualMap[DTPhDigiItr->whNum() + 2][DTPhDigiItr->scNum()][DTPhDigiItr->stNum() - 1]->code() <
            DTPhDigiItr->code()) {
      bestPhQualMap[DTPhDigiItr->whNum() + 2][DTPhDigiItr->scNum()][DTPhDigiItr->stNum() - 1] = &(*DTPhDigiItr);
    }
  }

  for (int iwh = 0; iwh < 5; iwh++) {
    for (int isec = 0; isec < 12; isec++) {
      for (int ist = 0; ist < 4; ist++) {
        if (bestPhQualMap[iwh][isec][ist]) {
          int xpos = iwh * 4 + ist + 1;
          dttpgphbestmap->Fill(xpos, isec);
          if (bestPhQualMap[iwh][isec][ist]->code() > 3)
            dttpgphbestmapcorr->Fill(xpos, isec);
        }
      }
    }
  }

  int bestThQualMap[5][12][3];
  memset(bestThQualMap, 0, 180 * sizeof(int));
  //for( vector<L1MuDTChambThDigi>::const_iterator
  for (L1MuDTChambThContainer::The_Container::const_iterator DTThDigiItr = myThContainer->begin();
       DTThDigiItr != myThContainer->end();
       ++DTThDigiItr) {
    ndttpgthtrack++;

    int bxindex = DTThDigiItr->bxNum() + 1;

    dttpgthwheel[bxindex]->Fill(DTThDigiItr->whNum());
    if (verbose_) {
      cout << "DTTPG theta wheel number " << DTThDigiItr->whNum() << endl;
    }
    dttpgthstation[bxindex]->Fill(DTThDigiItr->stNum());
    if (verbose_) {
      cout << "DTTPG theta station number " << DTThDigiItr->stNum() << endl;
    }
    dttpgthsector[bxindex]->Fill(DTThDigiItr->scNum());
    if (verbose_) {
      cout << "DTTPG theta sector number " << DTThDigiItr->scNum() << endl;
    }
    dttpgthbx[bxindex]->Fill(DTThDigiItr->bxNum());
    if (verbose_) {
      cout << "DTTPG theta bx number " << DTThDigiItr->bxNum() << endl;
    }
    int thcode[7] = {0, 0, 0, 0, 0, 0, 0};
    for (int j = 0; j < 7; j++) {
      dttpgththeta[bxindex]->Fill(DTThDigiItr->position(j));
      if (verbose_) {
        cout << "DTTPG theta position " << DTThDigiItr->position(j) << endl;
      }
      thcode[j] = DTThDigiItr->code(j);
      dttpgthquality[bxindex]->Fill(thcode[j]);
      if (verbose_) {
        cout << "DTTPG theta quality " << DTThDigiItr->code(j) << endl;
      }
    }

    int ypos = DTThDigiItr->scNum();
    int xpos = DTThDigiItr->stNum() + 4 * (DTThDigiItr->whNum() + 2);
    int bestqual = 0;
    dttpgthmap->Fill(xpos, ypos);
    dttpgthmapbx[bxindex]->Fill(xpos, ypos);
    for (int pos = 0; pos < 7; pos++) {
      if (thcode[pos] > bestqual)
        bestqual = thcode[pos];
      if (thcode[pos] == 2)
        dttpgthmaph->Fill(xpos, ypos);
    }

    if (bestThQualMap[DTThDigiItr->whNum() + 2][DTThDigiItr->scNum()][DTThDigiItr->stNum() - 1] < bestqual) {
      bestThQualMap[DTThDigiItr->whNum() + 2][DTThDigiItr->scNum()][DTThDigiItr->stNum() - 1] = bestqual;
    }
  }

  for (int iwh = 0; iwh < 5; iwh++) {
    for (int isec = 0; isec < 12; isec++) {
      for (int ist = 0; ist < 3; ist++) {
        if (bestThQualMap[iwh][isec][ist]) {
          int xpos = iwh * 4 + ist + 1;
          dttpgthbestmap->Fill(xpos, isec);
          if (bestThQualMap[iwh][isec][ist] == 2)
            dttpgthbestmaph->Fill(xpos, isec);
        }
      }
    }
  }

  dttpgphntrack->Fill(ndttpgphtrack);
  if (verbose_) {
    cout << "DTTPG phi ntrack " << ndttpgphtrack << endl;
  }
  dttpgthntrack->Fill(ndttpgthtrack);
  if (verbose_) {
    cout << "DTTPG theta ntrack " << ndttpgthtrack << endl;
  }

  edm::Handle<L1MuDTTrackContainer> myL1MuDTTrackContainer;
  e.getByToken(trToken_, myL1MuDTTrackContainer);

  if (!myL1MuDTTrackContainer.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuDTTrackContainer with label " << dttpgSource_.label();
    return;
  }

  L1MuDTTrackContainer::TrackContainer const *t = myL1MuDTTrackContainer->getContainer();

  int NumberOfSegmentsOut[3] = {0, 0, 0};
  for (L1MuDTTrackContainer::TrackContainer::const_iterator i = t->begin(); i != t->end(); ++i) {
    if (i->bx() == -1)
      NumberOfSegmentsOut[0]++;
    if (i->bx() == 0)
      NumberOfSegmentsOut[1]++;
    if (i->bx() == 1)
      NumberOfSegmentsOut[2]++;
  }

  /*Fill Histos for Segment counter*/
  for (int k = 0; k < 3; k++) {
    dttpgphbx[k + 5]->Fill(NumberOfSegmentsOut[k]);
  }

  /*Bunch assigments*/

  int bxCounterDttfOut = 0;
  for (int k = 0; k < 3; k++) {
    if (NumberOfSegmentsOut[k] > 0)
      bxCounterDttfOut++;
  }

  int bxCodeOut = 0;
  if (bxCounterDttfOut == 0) {
    bxCodeOut = 0;
  } else if (bxCounterDttfOut == 1) {
    for (int k = 0; k < 3; k++) {
      if (NumberOfSegmentsOut[k] > 0)
        bxCodeOut = k + 2;
    }
  } else if (bxCounterDttfOut == 2) {
    for (int k = 0; k < 3; k++) {
      if (NumberOfSegmentsOut[k] == 0)
        bxCodeOut = 8 - k;
    }
  } else if (bxCounterDttfOut == 3) {
    bxCodeOut = 10;
  }

  //The bx analyzer histo
  dttpgphbx[1]->Fill(bxCodeOut);

  /*End Dttf Output Bunch analysis*/

  // the 2-DIM histo with phi.input vs. output
  dttpgphbxcomp->Fill(bxCodePhi, bxCodeOut);

  for (L1MuDTTrackContainer::TrackContainer::const_iterator i = t->begin(); i != t->end(); ++i) {
    if (verbose_) {
      std::cout << "bx = " << i->bx() << std::endl;
      std::cout << "quality (packed) = " << i->quality_packed() << std::endl;
      std::cout << "pt      (packed) = " << i->pt_packed() << std::endl;
      std::cout << "phi     (packed) = " << i->phi_packed() << std::endl;
      std::cout << "charge  (packed) = " << i->charge_packed() << std::endl;
    }

    int bxindex = i->bx() + 1;
    dttf_p_phi[bxindex]->Fill(i->phi_packed());
    dttf_p_qual[bxindex]->Fill(i->quality_packed());
    dttf_p_pt[bxindex]->Fill(i->pt_packed());
    dttf_p_q[bxindex]->Fill(i->charge_packed());
  }
}

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 123 of file L1TDTTPG.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dttf_p_phi, dttf_p_pt, dttf_p_q, dttf_p_qual, dttpgphbestmap, dttpgphbestmapcorr, dttpgphbx, dttpgphbxcomp, dttpgphmap, dttpgphmap2nd, dttpgphmapbx, dttpgphmapcorr, dttpgphntrack, dttpgphquality, dttpgphsector, dttpgphstation, dttpgphts2tag, dttpgphwheel, dttpgthbestmap, dttpgthbestmaph, dttpgthbx, dttpgthmap, dttpgthmapbx, dttpgthmaph, dttpgthntrack, dttpgthquality, dttpgthsector, dttpgthstation, dttpgththeta, dttpgthwheel, mps_fire::i, nev_, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), setMapPhLabel(), and setMapThLabel().

                                                                                             {
  nev_ = 0;
  ibooker.setCurrentFolder("L1T/L1TDTTPG");

  //hist1[0]
  dttpgphbx[0] = ibooker.book1D("BxEncoding_PHI", "Bunch encoding DTTF Phi", 11, 0, 11);
  //hist1[1]
  dttpgphbx[1] = ibooker.book1D("BxEncoding_OUT", "Bunch encoding DTTF Output", 11, 0, 11);

  //hist1[2]
  dttpgphbx[2] = ibooker.book1D("NumberOfSegmentsPHI_BunchNeg1", "Number of segments for bunch -1 Dttf Phi", 20, 0, 20);
  //hist1[3]
  dttpgphbx[3] = ibooker.book1D("NumberOfSegmentsPHI_Bunch0", "Number of segments for bunch 0 Dttf Phi", 20, 0, 20);
  //hist1[4]
  dttpgphbx[4] = ibooker.book1D("NumberOfSegmentsPHI_Bunch1", "Number of segments for bunch 1 Dttf Phi", 20, 0, 20);

  //hist1[5]
  dttpgphbx[5] =
      ibooker.book1D("NumberOfSegmentsOUT_BunchNeg1", "Number of segments for bunch -1 Dttf Output", 20, 0, 20);
  //hist1[6]
  dttpgphbx[6] = ibooker.book1D("NumberOfSegmentsOUT_Bunch0", "Number of segments for bunch 0 Dttf Output", 20, 0, 20);
  //hist1[7]
  dttpgphbx[7] = ibooker.book1D("NumberOfSegmentsOUT_Bunch1", "Number of segments for bunch 1 Dttf Output", 20, 0, 20);

  for (int i = 0; i < 2; i++) {
    dttpgphbx[i]->setBinLabel(1, "None");
    dttpgphbx[i]->setBinLabel(3, "Only bx=-1");
    dttpgphbx[i]->setBinLabel(4, "Only bx= 0");
    dttpgphbx[i]->setBinLabel(5, "Only bx=+1");
    dttpgphbx[i]->setBinLabel(7, "Bx=-1,0");
    dttpgphbx[i]->setBinLabel(8, "Bx=-1,1");
    dttpgphbx[i]->setBinLabel(9, "Bx= 0,1");
    dttpgphbx[i]->setBinLabel(11, "All bx");
  }

  dttpgphbxcomp = ibooker.book2D("BxEncoding_PHI_OUT", "Bunch encoding: DTTF Phi vs. Output", 11, 0, 11, 11, 0, 11);
  dttpgphbxcomp->setAxisTitle("DTTF (output)", 1);
  dttpgphbxcomp->setAxisTitle("PHI-TF", 2);
  for (int i = 1; i <= 2; i++) {
    dttpgphbxcomp->setBinLabel(1, "None", i);
    dttpgphbxcomp->setBinLabel(3, "Only bx=-1", i);
    dttpgphbxcomp->setBinLabel(4, "Only bx= 0", i);
    dttpgphbxcomp->setBinLabel(5, "Only bx=+1", i);
    dttpgphbxcomp->setBinLabel(7, "Bx=-1,0", i);
    dttpgphbxcomp->setBinLabel(8, "Bx=-1,1", i);
    dttpgphbxcomp->setBinLabel(9, "Bx= 0,1", i);
    dttpgphbxcomp->setBinLabel(11, "All bx", i);
  }

  dttpgphntrack = ibooker.book1D("DT_TPG_phi_ntrack", "DT TPG phi ntrack", 20, -0.5, 19.5);
  dttpgthntrack = ibooker.book1D("DT_TPG_theta_ntrack", "DT TPG theta ntrack", 20, -0.5, 19.5);

  for (int ibx = 0; ibx <= 2; ibx++) {
    ostringstream bxnum;
    bxnum << ibx - 1;
    string bxn;
    if (ibx < 2)
      bxn = bxnum.str();
    else
      bxn = "+" + bxnum.str();

    // Phi
    dttpgphwheel[ibx] =
        ibooker.book1D("DT_TPG_phi_wheel_number_" + bxn, "DT TPG phi wheel number " + bxn, 5, -2.5, 2.5);
    dttpgphsector[ibx] =
        ibooker.book1D("DT_TPG_phi_sector_number_" + bxn, "DT TPG phi sector number " + bxn, 12, -0.5, 11.5);
    dttpgphstation[ibx] =
        ibooker.book1D("DT_TPG_phi_station_number_" + bxn, "DT TPG phi station number " + bxn, 5, 0.5, 4.5);
    //  dttpgphphi[ibx] = dbe->book1D("DT_TPG_phi_"+bxn,
    //                    "DT TPG phi "+bxn, 100, -2100., 2100. ) ;
    //  dttpgphphiB[ibx] = dbe->book1D("DT_TPG_phiB_"+bxn,
    //                     "DT TPG phiB "+bxn, 100, -550., 550. ) ;
    dttpgphquality[ibx] = ibooker.book1D("DT_TPG_phi_quality_" + bxn, "DT TPG phi quality " + bxn, 8, -0.5, 7.5);
    dttpgphts2tag[ibx] = ibooker.book1D("DT_TPG_phi_Ts2Tag_" + bxn, "DT TPG phi Ts2Tag " + bxn, 2, -0.5, 1.5);
    //  dttpgphbxcnt[ibx] = dbe->book1D("DT_TPG_phi_BxCnt_"+bxn,
    //                  "DT TPG phi BxCnt "+bxn, 10, -0.5, 9.5 ) ;
    dttpgphmapbx[ibx] =
        ibooker.book2D("DT_TPG_phi_map_bx" + bxn, "Map of triggers per station (BX=" + bxn + ")", 20, 1, 21, 12, 0, 12);
    setMapPhLabel(dttpgphmapbx[ibx]);

    //Theta
    dttpgthbx[ibx] = ibooker.book1D("DT_TPG_theta_bx_" + bxn, "DT TPG theta bx " + bxn, 50, -24.5, 24.5);
    dttpgthwheel[ibx] =
        ibooker.book1D("DT_TPG_theta_wheel_number_" + bxn, "DT TPG theta wheel number " + bxn, 5, -2.5, 2.5);
    dttpgthsector[ibx] =
        ibooker.book1D("DT_TPG_theta_sector_number_" + bxn, "DT TPG theta sector number " + bxn, 12, -0.5, 11.5);
    dttpgthstation[ibx] =
        ibooker.book1D("DT_TPG_theta_station_number_" + bxn, "DT TPG theta station number " + bxn, 5, -0.5, 4.5);
    dttpgththeta[ibx] = ibooker.book1D("DT_TPG_theta_" + bxn, "DT TPG theta " + bxn, 20, -0.5, 19.5);
    dttpgthquality[ibx] = ibooker.book1D("DT_TPG_theta_quality_" + bxn, "DT TPG theta quality " + bxn, 8, -0.5, 7.5);
    dttpgthmapbx[ibx] = ibooker.book2D(
        "DT_TPG_theta_map_bx_" + bxn, "Map of triggers per station (BX=" + bxn + ")", 15, 1, 16, 12, 0, 12);
    setMapThLabel(dttpgthmapbx[ibx]);

    // Phi output
    dttf_p_phi[ibx] = ibooker.book1D("dttf_p_phi_" + bxn, "dttf phi output #phi " + bxn, 256, -0.5, 255.5);
    dttf_p_qual[ibx] = ibooker.book1D("dttf_p_qual_" + bxn, "dttf phi output qual " + bxn, 8, -0.5, 7.5);
    dttf_p_q[ibx] = ibooker.book1D("dttf_p_q_" + bxn, "dttf phi output q " + bxn, 2, -0.5, 1.5);
    dttf_p_pt[ibx] = ibooker.book1D("dttf_p_pt_" + bxn, "dttf phi output p_{t} " + bxn, 32, -0.5, 31.5);
  }

  dttpgphmap = ibooker.book2D("DT_TPG_phi_map", "Map of triggers per station", 20, 1, 21, 12, 0, 12);
  dttpgphmapcorr =
      ibooker.book2D("DT_TPG_phi_map_corr", "Map of correlated triggers per station", 20, 1, 21, 12, 0, 12);
  dttpgphmap2nd = ibooker.book2D("DT_TPG_phi_map_2nd", "Map of second tracks per station", 20, 1, 21, 12, 0, 12);
  dttpgphbestmap = ibooker.book2D("DT_TPG_phi_best_map", "Map of best triggers per station", 20, 1, 21, 12, 0, 12);
  dttpgphbestmapcorr =
      ibooker.book2D("DT_TPG_phi_best_map_corr", "Map of correlated best triggers per station", 20, 1, 21, 12, 0, 12);
  setMapPhLabel(dttpgphmap);
  setMapPhLabel(dttpgphmapcorr);
  setMapPhLabel(dttpgphmap2nd);
  setMapPhLabel(dttpgphbestmap);
  setMapPhLabel(dttpgphbestmapcorr);

  dttpgthmap = ibooker.book2D("DT_TPG_theta_map", "Map of triggers per station", 15, 1, 16, 12, 0, 12);
  dttpgthmaph = ibooker.book2D("DT_TPG_theta_map_h", "Map of H quality triggers per station", 15, 1, 16, 12, 0, 12);
  dttpgthbestmap = ibooker.book2D("DT_TPG_theta_best_map", "Map of besttriggers per station", 15, 1, 16, 12, 0, 12);
  dttpgthbestmaph =
      ibooker.book2D("DT_TPG_theta_best_map_h", "Map of H quality best triggers per station", 15, 1, 16, 12, 0, 12);
  setMapThLabel(dttpgthmap);
  setMapThLabel(dttpgthmaph);
  setMapThLabel(dttpgthbestmap);
  setMapThLabel(dttpgthbestmaph);
}

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 119 of file L1TDTTPG.cc.

                                                                  {
  //empty
}

void L1TDTTPG::setMapPhLabel ( MonitorElement *  me )

private

Definition at line 533 of file L1TDTTPG.cc.

References mps_fire::i, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                               {
  me->setAxisTitle("DTTF Sector", 2);
  for (int i = 0; i < 5; i++) {
    ostringstream wheel;
    wheel << i - 2;
    me->setBinLabel(1 + i * 4, "Wheel " + wheel.str(), 1);
  }
}

void L1TDTTPG::setMapThLabel ( MonitorElement *  me )

private

Definition at line 542 of file L1TDTTPG.cc.

References mps_fire::i, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                               {
  me->setAxisTitle("DTTF Sector", 2);
  for (int i = 0; i < 5; i++) {
    ostringstream wheel;
    wheel << i - 2;
    me->setBinLabel(1 + i * 3, "Wheel " + wheel.str(), 1);
  }
}

Member Data Documentation

MonitorElement* L1TDTTPG::dttf_p_phi[3]

private

Definition at line 92 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttf_p_pt[3]

private

Definition at line 93 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttf_p_q[3]

private

Definition at line 94 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttf_p_qual[3]

private

Definition at line 95 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbestmap

private

Definition at line 76 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbestmapcorr

private

Definition at line 77 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbx[8]

private

Definition at line 61 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphbxcomp

private

Definition at line 62 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmap

private

Definition at line 72 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmap2nd

private

Definition at line 74 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmapbx[3]

private

Definition at line 73 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphmapcorr

private

Definition at line 75 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphntrack

private

Definition at line 71 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphquality[3]

private

Definition at line 68 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphsector[3]

private

Definition at line 64 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphstation[3]

private

Definition at line 65 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphts2tag[3]

private

Definition at line 69 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgphwheel[3]

private

Definition at line 63 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag L1TDTTPG::dttpgSource_

private

Definition at line 104 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().

edm::EDGetTokenT<L1MuDTChambPhContainer> L1TDTTPG::dttpgSourcePhContainer_token_

private

Definition at line 102 of file L1TDTTPG.h.

Referenced by analyze().

edm::EDGetTokenT<L1MuDTChambThContainer> L1TDTTPG::dttpgSourceThContainer_token_

private

Definition at line 103 of file L1TDTTPG.h.

Referenced by analyze().

MonitorElement* L1TDTTPG::dttpgthbestmap

private

Definition at line 89 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthbestmaph

private

Definition at line 90 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthbx[3]

private

Definition at line 79 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthmap

private

Definition at line 86 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthmapbx[3]

private

Definition at line 87 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthmaph

private

Definition at line 88 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthntrack

private

Definition at line 85 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthquality[3]

private

Definition at line 84 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthsector[3]

private

Definition at line 81 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthstation[3]

private

Definition at line 82 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgththeta[3]

private

Definition at line 83 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TDTTPG::dttpgthwheel[3]

private

Definition at line 80 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

std::ofstream L1TDTTPG::logFile_

private

Definition at line 101 of file L1TDTTPG.h.

bool L1TDTTPG::monitorDaemon_

private

Definition at line 100 of file L1TDTTPG.h.

int L1TDTTPG::nev_

private

Definition at line 97 of file L1TDTTPG.h.

Referenced by analyze(), and bookHistograms().

std::string L1TDTTPG::outputFile_

private

Definition at line 98 of file L1TDTTPG.h.

Referenced by L1TDTTPG().

std::string L1TDTTPG::trstring_

private

Definition at line 106 of file L1TDTTPG.h.

Referenced by L1TDTTPG().

edm::EDGetTokenT<L1MuDTTrackContainer> L1TDTTPG::trToken_

private

Definition at line 107 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().

bool L1TDTTPG::verbose_

private

Definition at line 99 of file L1TDTTPG.h.

Referenced by analyze(), and L1TDTTPG().