L1TCompare Class Reference

#include <L1TCompare.h>

Inheritance diagram for L1TCompare:

Classes
class	RctObject

Public Member Functions
	L1TCompare (const edm::ParameterSet &ps)
	~L1TCompare () 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 Types
typedef std::vector< L1TCompare::RctObject >	RctObjectCollection

Private Member Functions
bool	verbose () const

Static Private Member Functions
static bool	rctObjectComp (const RctObject &a, const RctObject &b)

Private Attributes
MonitorElement *	ecalTpgRctLeadingEmEta2_
MonitorElement *	ecalTpgRctLeadingEmEta_
MonitorElement *	ecalTpgRctLeadingEmPhi_
MonitorElement *	ecalTpgRctLeadingEmRank_
edm::InputTag	ecalTpgSource_
edm::EDGetTokenT< EcalTrigPrimDigiCollection >	ecalTpgSource_token_
edm::EDGetTokenT< L1GctJetCandCollection >	gctCenJetsToken_
edm::EDGetTokenT< L1GctEmCandCollection >	gctIsoEmCandsToken_
edm::EDGetTokenT< L1GctEmCandCollection >	gctNonIsoEmCandsToken_
edm::InputTag	gctSource_
std::ofstream	logFile_
bool	monitorDaemon_
int	nev_
std::string	outputFile_
MonitorElement *	rctGctLeadingIsoEmEta_
MonitorElement *	rctGctLeadingIsoEmPhi_
MonitorElement *	rctGctLeadingIsoEmRank_
MonitorElement *	rctGctLeadingNonIsoEmEta_
MonitorElement *	rctGctLeadingNonIsoEmPhi_
MonitorElement *	rctGctLeadingNonIsoEmRank_
edm::InputTag	rctSource_
edm::EDGetTokenT< L1CaloEmCollection >	rctSourceEm_token_
edm::EDGetTokenT< L1CaloRegionCollection >	rctSourceRctEmRgn_token_
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 62 of file L1TCompare.h.

Member Typedef Documentation

typedef std::vector<L1TCompare::RctObject> L1TCompare::RctObjectCollection

private

Definition at line 123 of file L1TCompare.h.

Constructor & Destructor Documentation

L1TCompare::L1TCompare

(

const edm::ParameterSet &

ps

)

Definition at line 90 of file L1TCompare.cc.

References gather_cfg::cout, gctCenJetsToken_, gctIsoEmCandsToken_, gctNonIsoEmCandsToken_, gctSource_, edm::ParameterSet::getUntrackedParameter(), edm::InputTag::label(), outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, verbose(), and verbose_.

    : rctSourceEm_token_(consumes<L1CaloEmCollection>(ps.getParameter<InputTag>("rctSource"))),
      rctSourceRctEmRgn_token_(consumes<L1CaloRegionCollection>(ps.getParameter<InputTag>("rctSource"))),
      rctSource_(ps.getParameter<InputTag>("rctSource")),
      gctSource_(ps.getParameter<InputTag>("gctSource")),
      ecalTpgSource_(ps.getParameter<edm::InputTag>("ecalTpgSource")),
      ecalTpgSource_token_(consumes<EcalTrigPrimDigiCollection>(ps.getParameter<edm::InputTag>("ecalTpgSource")))

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

  if (verbose())
    std::cout << "L1TCompare: constructor...." << std::endl;

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

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

  //set Token(-s)
  edm::InputTag gctCenJetsTag_(gctSource_.label(), "cenJets");
  edm::InputTag gctIsoEmCandsTag_(gctSource_.label(), "isoEm");
  edm::InputTag gctNonIsoEmCandsTag_(gctSource_.label(), "nonIsoEm");

  gctCenJetsToken_ = consumes<L1GctJetCandCollection>(gctCenJetsTag_);
  gctIsoEmCandsToken_ = consumes<L1GctEmCandCollection>(gctIsoEmCandsTag_);
  gctNonIsoEmCandsToken_ = consumes<L1GctEmCandCollection>(gctNonIsoEmCandsTag_);
}

L1TCompare::~L1TCompare ( )

override

Definition at line 125 of file L1TCompare.cc.

{}

Member Function Documentation

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

overrideprotected

Definition at line 191 of file L1TCompare.cc.

References edm::SortedCollection< T, SORT >::begin(), beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, ecalTpgRctLeadingEmEta2_, ecalTpgRctLeadingEmEta_, ecalTpgRctLeadingEmPhi_, ecalTpgRctLeadingEmRank_, ecalTpgSource_, ecalTpgSource_token_, edm::SortedCollection< T, SORT >::end(), dqm::impl::MonitorElement::Fill(), gctCenJetsToken_, gctIsoEmCandsToken_, gctNonIsoEmCandsToken_, edm::Event::getByToken(), h, triggerObjects_cff::id, edm::HandleBase::isValid(), edm::InputTag::label(), nev_, rctGctLeadingIsoEmEta_, rctGctLeadingIsoEmPhi_, rctGctLeadingIsoEmRank_, rctGctLeadingNonIsoEmEta_, rctGctLeadingNonIsoEmPhi_, rctGctLeadingNonIsoEmRank_, rctObjectComp(), rctSource_, rctSourceEm_token_, rctSourceRctEmRgn_token_, and verbose().

                                                            {
  ++nev_;
  if (verbose()) {
    std::cout << "L1TCompare: analyze...." << std::endl;
  }

  // L1E
  edm::Handle<L1EmParticleCollection> l1eIsoEm;
  edm::Handle<L1EmParticleCollection> l1eNonIsoEm;
  edm::Handle<L1JetParticleCollection> l1eCenJets;
  edm::Handle<L1JetParticleCollection> l1eForJets;
  edm::Handle<L1JetParticleCollection> l1eTauJets;
  //  edm::Handle < L1EtMissParticle > l1eEtMiss;
  edm::Handle<L1EtMissParticleCollection> l1eEtMiss;
  // RCT
  edm::Handle<L1CaloEmCollection> em;  // collection of L1CaloEmCands
  edm::Handle<L1CaloRegionCollection> rctEmRgn;

  // GCT
  edm::Handle<L1GctJetCandCollection> gctCenJets;
  edm::Handle<L1GctEmCandCollection> gctIsoEmCands;
  edm::Handle<L1GctEmCandCollection> gctNonIsoEmCands;

  e.getByToken(rctSourceEm_token_, em);

  if (!em.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1CaloEmCollection with label " << rctSource_.label();
    return;
  }

  e.getByToken(rctSourceRctEmRgn_token_, rctEmRgn);

  if (!rctEmRgn.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find "
                                 << "L1CaloRegionCollection with label " << rctSource_.label();
    return;
  }

  e.getByToken(gctCenJetsToken_, gctCenJets);
  e.getByToken(gctIsoEmCandsToken_, gctIsoEmCands);
  e.getByToken(gctNonIsoEmCandsToken_, gctNonIsoEmCands);

  if (!gctCenJets.isValid()) {
    std::cerr << "L1TGCT: could not find one of the classes?" << std::endl;
    return;
  }
  if (!gctIsoEmCands.isValid()) {
    std::cerr << "L1TGCT: could not find one of the classes?" << std::endl;
    return;
  }
  if (!gctNonIsoEmCands.isValid()) {
    std::cerr << "L1TGCT: could not find one of the classes?" << std::endl;
    return;
  }

  // GCT
  if (verbose()) {
    for (L1GctEmCandCollection::const_iterator iem = gctIsoEmCands->begin(); iem != gctIsoEmCands->end(); ++iem) {
      if (!iem->empty())
        std::cout << "GCT EM: " << iem->rank() << ", " << iem->etaIndex()
                  << "("
                  //<< int(iem->etaIndex()&0x3)*((iem->etaIndex()&0x4)?1:-1)
                  << "), " << iem->phiIndex() << std::endl;
    }
  }
  // rct phi: 0-17
  // rct eta: 0-21

  // Fill the RCT histograms

  // Regions
  RctObjectCollection rcj, rcj_iso, rcj_non_iso;
  for (L1CaloEmCollection::const_iterator iem = em->begin(); iem != em->end(); ++iem) {
    //   L1CaloRegionDetId id(false, iem->rctCrate(), iem->rctCard(),
    //           iem->rctRegion());
    L1CaloRegionDetId id(iem->rctCrate(), iem->rctCard(), iem->rctRegion());

    // RctObject h(id.gctEta(), id.gctPhi(), iem->rank());
    RctObject h(id.rctEta(), id.rctPhi(), iem->rank());
    if (!iem->isolated())
      rcj_non_iso.push_back(h);
    else
      rcj_iso.push_back(h);
    rcj.push_back(h);
  }
  // not so smart but ...
  std::sort(rcj.begin(), rcj.end(), rctObjectComp);
  std::sort(rcj_non_iso.begin(), rcj_non_iso.end(), rctObjectComp);
  std::sort(rcj_iso.begin(), rcj_iso.end(), rctObjectComp);
  if (verbose()) {
    for (RctObjectCollection::reverse_iterator ij = rcj_iso.rbegin();
         ij != rcj_iso.rend() && ij != rcj_iso.rbegin() + 8;
         ++ij) {
      std::cout << "RCT cj: " << ij->rank_ << ", " << ij->eta_ << ", " << ij->phi_ << std::endl;
    }
  }
  L1GctEmCandCollection::const_iterator lead_em = gctIsoEmCands->begin();
  if (!lead_em->empty()) {  // equivalent to rank == 0
    rctGctLeadingIsoEmEta_->Fill(lead_em->etaIndex(), rcj_iso.rbegin()->eta_);
    rctGctLeadingIsoEmPhi_->Fill(lead_em->phiIndex(), rcj_iso.rbegin()->phi_);
    rctGctLeadingIsoEmRank_->Fill(lead_em->rank(), rcj_iso.rbegin()->rank_);
  }

  // non-isolated
  if (verbose()) {
    for (L1GctEmCandCollection::const_iterator iem = gctNonIsoEmCands->begin(); iem != gctNonIsoEmCands->end(); ++iem) {
      if (!iem->empty())
        std::cout << "GCT EM non: " << iem->rank() << ", "
                  << iem->etaIndex()  //<< "("
                                      //<< int(iem->etaIndex()&0x3)*((iem->etaIndex()&0x4)?1:-1)
                  //<< ")"
                  << ", " << iem->phiIndex() << std::endl;
    }
  }
  if (verbose()) {
    for (RctObjectCollection::reverse_iterator ij = rcj_non_iso.rbegin();
         ij != rcj_non_iso.rend() && ij != rcj_non_iso.rbegin() + 8;
         ++ij) {
      std::cout << "RCT cj non: " << ij->rank_ << ", " << ij->eta_ << ", " << ij->phi_ << std::endl;
    }
  }
  lead_em = gctNonIsoEmCands->begin();
  if (!lead_em->empty()) {  // equivalent to rank != 0
    rctGctLeadingNonIsoEmEta_->Fill(lead_em->etaIndex(), rcj_non_iso.rbegin()->eta_);
    rctGctLeadingNonIsoEmPhi_->Fill(lead_em->phiIndex(), rcj_non_iso.rbegin()->phi_);
    rctGctLeadingNonIsoEmRank_->Fill(lead_em->rank(), rcj_non_iso.rbegin()->rank_);
  }

  // ECAL TPG's to RCT EM
  edm::Handle<EcalTrigPrimDigiCollection> eTP;
  e.getByToken(ecalTpgSource_token_, eTP);

  if (!eTP.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find EcalTrigPrimCollection with label " << ecalTpgSource_.label();
    return;
  }
  RctObjectCollection ecalobs;
  for (EcalTrigPrimDigiCollection::const_iterator ieTP = eTP->begin(); ieTP != eTP->end(); ieTP++) {
    ecalobs.push_back(RctObject(ieTP->id().ieta(), ieTP->id().iphi(), ieTP->compressedEt()));
  }
  std::sort(ecalobs.begin(), ecalobs.end(), rctObjectComp);
  if (verbose()) {
    for (RctObjectCollection::reverse_iterator ij = ecalobs.rbegin();
         ij != ecalobs.rend() && ij != ecalobs.rbegin() + 8;
         ++ij) {
      std::cout << "ECAL cj : " << ij->rank_ << ", " << ij->eta_ << ", " << ij->phi_ << std::endl;
    }
  }
  // abritrary cut
  if (rcj.rbegin()->rank_ > 4) {
    ecalTpgRctLeadingEmEta_->Fill(rcj.rbegin()->eta_, ecalobs.rbegin()->eta_);
    int e2 = (rcj.rbegin()->eta_ & 0x7UL) * ((rcj.rbegin()->eta_ & 0x8UL) ? 1 : -1);
    ecalTpgRctLeadingEmEta2_->Fill(e2, ecalobs.rbegin()->eta_);
    ecalTpgRctLeadingEmPhi_->Fill(rcj.rbegin()->phi_, ecalobs.rbegin()->phi_);
    ecalTpgRctLeadingEmRank_->Fill(rcj.rbegin()->rank_, ecalobs.rbegin()->rank_);
  }
  if (verbose()) {
    int seta = rcj.rbegin()->eta_;
    seta = (seta & 0x7UL) * (seta & 0x8 ? -1 : 1);
    std::cout << "ZZ: " << rcj.rbegin()->eta_ << " " << rcj.rbegin()->phi_ << " " << rcj.rbegin()->rank_ << " "
              << (++rcj.rbegin())->rank_ << " " << ecalobs.rbegin()->eta_ << " " << ecalobs.rbegin()->phi_ << " "
              << ecalobs.rbegin()->rank_ << " " << (++ecalobs.rbegin())->rank_ << " " << seta << " " << std::endl;
  }
}

void L1TCompare::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 131 of file L1TCompare.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), ecalTpgRctLeadingEmEta2_, ecalTpgRctLeadingEmEta_, ecalTpgRctLeadingEmPhi_, ecalTpgRctLeadingEmRank_, ETABINS, ETAMAX, ETAMIN, nev_, PHIBINS, PHIMAX, PHIMIN, R6BINS, R6MAX, R6MIN, rctGctLeadingIsoEmEta_, rctGctLeadingIsoEmPhi_, rctGctLeadingIsoEmRank_, rctGctLeadingNonIsoEmEta_, rctGctLeadingNonIsoEmPhi_, rctGctLeadingNonIsoEmRank_, dqm::impl::MonitorElement::setAxisTitle(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, TPETABINS, TPETAMAX, TPETAMIN, TPPHIBINS, TPPHIMAX, and TPPHIMIN.

                                                                                                         {
  nev_ = 0;

  ibooker.setCurrentFolder("L1T/Compare");

  // -------------------------------------------
  // RCT-GCT
  // -------------------------------------------
  // Isolated
  rctGctLeadingIsoEmRank_ =
      ibooker.book2D("rctGctLeadingIsoEmRank", "RCT-GCT: rank", R6BINS, R6MIN, R6MAX, R6BINS, R6MIN, R6MAX);
  rctGctLeadingIsoEmRank_->setAxisTitle(std::string("gct"), 1);
  rctGctLeadingIsoEmRank_->setAxisTitle(std::string("rct"), 2);
  rctGctLeadingIsoEmEta_ =
      ibooker.book2D("rctGctLeadingIsoEmEta", "RCT-GCT: #eta", ETABINS, ETAMIN, ETAMAX, ETABINS, ETAMIN, ETAMAX);
  rctGctLeadingIsoEmEta_->setAxisTitle(std::string("gct"), 1);
  rctGctLeadingIsoEmEta_->setAxisTitle(std::string("rct"), 2);

  rctGctLeadingIsoEmPhi_ =
      ibooker.book2D("rctGctLeadingIsoEmPhi", "RCT-GCT: #phi", PHIBINS, PHIMIN, PHIMAX, PHIBINS, PHIMIN, PHIMAX);
  rctGctLeadingIsoEmPhi_->setAxisTitle(std::string("gct"), 1);
  rctGctLeadingIsoEmPhi_->setAxisTitle(std::string("rct"), 2);
  // non-Isolated
  rctGctLeadingNonIsoEmRank_ =
      ibooker.book2D("rctGctLeadingNonIsoEmRank", "RCT-GCT: rank", R6BINS, R6MIN, R6MAX, R6BINS, R6MIN, R6MAX);
  rctGctLeadingNonIsoEmRank_->setAxisTitle(std::string("gct"), 1);
  rctGctLeadingNonIsoEmRank_->setAxisTitle(std::string("rct"), 2);

  rctGctLeadingNonIsoEmEta_ =
      ibooker.book2D("rctGctLeadingNonIsoEmEta", "RCT-GCT: #eta", ETABINS, ETAMIN, ETAMAX, ETABINS, ETAMIN, ETAMAX);
  rctGctLeadingNonIsoEmEta_->setAxisTitle(std::string("gct"), 1);
  rctGctLeadingNonIsoEmEta_->setAxisTitle(std::string("rct"), 2);

  rctGctLeadingNonIsoEmPhi_ =
      ibooker.book2D("rctGctLeadingNonIsoEmPhi", "RCT-GCT: #phi", PHIBINS, PHIMIN, PHIMAX, PHIBINS, PHIMIN, PHIMAX);
  rctGctLeadingNonIsoEmPhi_->setAxisTitle(std::string("gct"), 1);
  rctGctLeadingNonIsoEmPhi_->setAxisTitle(std::string("rct"), 2);
  // -------------------------------------------
  // ECAL TPG - RCT
  // -------------------------------------------
  ecalTpgRctLeadingEmRank_ =
      ibooker.book2D("ecalTpgRctLeadingEmRank", "ECAL TPG-RCT: rank", R6BINS, R6MIN, R6MAX, R6BINS, R6MIN, R6MAX);
  ecalTpgRctLeadingEmRank_->setAxisTitle(std::string("rct"), 1);
  ecalTpgRctLeadingEmRank_->setAxisTitle(std::string("ecal tp"), 2);

  ecalTpgRctLeadingEmEta_ =
      ibooker.book2D("ecalTpgRctLeadingEmEta", "ECAL TPG-RCT: #eta", 15, -0.5, 14.5, TPETABINS, TPETAMIN, TPETAMAX);
  ecalTpgRctLeadingEmEta_->setAxisTitle(std::string("rct"), 1);
  ecalTpgRctLeadingEmEta_->setAxisTitle(std::string("ecal tp"), 2);
  ecalTpgRctLeadingEmEta2_ =
      ibooker.book2D("ecalTpgRctLeadingEmEta2", "ECAL TPG-RCT: #eta (2)", 13, -6.5, 6.5, TPETABINS, TPETAMIN, TPETAMAX);
  ecalTpgRctLeadingEmEta2_->setAxisTitle(std::string("rct"), 1);
  ecalTpgRctLeadingEmEta2_->setAxisTitle(std::string("ecal tp"), 2);
  ecalTpgRctLeadingEmPhi_ = ibooker.book2D(
      "ecalTpgRctLeadingEmPhi", "ECAL TPG-RCT: #phi", PHIBINS, PHIMIN, PHIMAX, TPPHIBINS, TPPHIMIN, TPPHIMAX);
  ecalTpgRctLeadingEmPhi_->setAxisTitle(std::string("rct"), 1);
  ecalTpgRctLeadingEmPhi_->setAxisTitle(std::string("ecal tp"), 2);
  //}
}

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 127 of file L1TCompare.cc.

                                                                {
  //
}

static bool L1TCompare::rctObjectComp ( const RctObject &  a, const RctObject &  b  )

inlinestaticprivate

Definition at line 128 of file L1TCompare.h.

References L1TCompare::RctObject::eta_, L1TCompare::RctObject::phi_, and L1TCompare::RctObject::rank_.

Referenced by analyze().

                                                                    {
    // for equal rank I don't know what the appropriate sorting is.
    if (a.rank_ == b.rank_) {
      if (a.eta_ == b.eta_) {
        return a.phi_ < b.phi_;
      } else {
        return a.eta_ < b.eta_;
      }
    } else {
      return a.rank_ < b.rank_;
    }
  }

bool L1TCompare::verbose ( ) const

inlineprivate

Definition at line 100 of file L1TCompare.h.

References monitorDaemon_, and verbose_.

Referenced by analyze(), and L1TCompare().

{ return verbose_; };

Member Data Documentation

MonitorElement* L1TCompare::ecalTpgRctLeadingEmEta2_

private

Definition at line 93 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::ecalTpgRctLeadingEmEta_

private

Definition at line 92 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::ecalTpgRctLeadingEmPhi_

private

Definition at line 94 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::ecalTpgRctLeadingEmRank_

private

Definition at line 95 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag L1TCompare::ecalTpgSource_

private

Definition at line 108 of file L1TCompare.h.

Referenced by analyze().

edm::EDGetTokenT<EcalTrigPrimDigiCollection> L1TCompare::ecalTpgSource_token_

private

Definition at line 109 of file L1TCompare.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctJetCandCollection> L1TCompare::gctCenJetsToken_

private

Definition at line 112 of file L1TCompare.h.

Referenced by analyze(), and L1TCompare().

edm::EDGetTokenT<L1GctEmCandCollection> L1TCompare::gctIsoEmCandsToken_

private

Definition at line 113 of file L1TCompare.h.

Referenced by analyze(), and L1TCompare().

edm::EDGetTokenT<L1GctEmCandCollection> L1TCompare::gctNonIsoEmCandsToken_

private

Definition at line 114 of file L1TCompare.h.

Referenced by analyze(), and L1TCompare().

edm::InputTag L1TCompare::gctSource_

private

Definition at line 107 of file L1TCompare.h.

Referenced by L1TCompare().

std::ofstream L1TCompare::logFile_

private

Definition at line 102 of file L1TCompare.h.

bool L1TCompare::monitorDaemon_

private

Definition at line 100 of file L1TCompare.h.

Referenced by verbose().

int L1TCompare::nev_

private

Definition at line 97 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

std::string L1TCompare::outputFile_

private

Definition at line 98 of file L1TCompare.h.

Referenced by L1TCompare().

MonitorElement* L1TCompare::rctGctLeadingIsoEmEta_

private

Definition at line 83 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::rctGctLeadingIsoEmPhi_

private

Definition at line 84 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::rctGctLeadingIsoEmRank_

private

Definition at line 85 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::rctGctLeadingNonIsoEmEta_

private

Definition at line 87 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::rctGctLeadingNonIsoEmPhi_

private

Definition at line 88 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TCompare::rctGctLeadingNonIsoEmRank_

private

Definition at line 89 of file L1TCompare.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag L1TCompare::rctSource_

private

Definition at line 106 of file L1TCompare.h.

Referenced by analyze().

edm::EDGetTokenT<L1CaloEmCollection> L1TCompare::rctSourceEm_token_

private

Definition at line 104 of file L1TCompare.h.

Referenced by analyze().

edm::EDGetTokenT<L1CaloRegionCollection> L1TCompare::rctSourceRctEmRgn_token_

private

Definition at line 105 of file L1TCompare.h.

Referenced by analyze().

bool L1TCompare::verbose_

private

Definition at line 99 of file L1TCompare.h.

Referenced by L1TCompare(), and verbose().