L1TGCT Class Reference

#include <L1TGCT.h>

Inheritance diagram for L1TGCT:

Public Member Functions
	L1TGCT (const edm::ParameterSet &ps)
virtual	~L1TGCT ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
virtual void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)
virtual void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
virtual void	dqmBeginRun (const edm::Run &, const edm::EventSetup &)
Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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 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 Attributes
int	filterTriggerType_
edm::InputTag	gctCenJetsSource_
edm::EDGetTokenT < L1GctJetCandCollection >	gctCenJetsSourceToken_
edm::InputTag	gctEnergySumsSource_
edm::EDGetTokenT < L1GctHFRingEtSumsCollection >	gctEnergySumsSourceToken_
edm::InputTag	gctForJetsSource_
edm::EDGetTokenT < L1GctJetCandCollection >	gctForJetsSourceToken_
edm::InputTag	gctIsoEmSource_
edm::EDGetTokenT < L1GctEmCandCollection >	gctIsoEmSourceToken_
edm::InputTag	gctIsoTauJetsSource_
edm::EDGetTokenT < L1GctJetCandCollection >	gctIsoTauJetsSourceToken_
edm::InputTag	gctNonIsoEmSource_
edm::EDGetTokenT < L1GctEmCandCollection >	gctNonIsoEmSourceToken_
edm::InputTag	gctTauJetsSource_
edm::EDGetTokenT < L1GctJetCandCollection >	gctTauJetsSourceToken_
edm::EDGetTokenT < L1GctEtHadCollection >	l1EtHadToken_
edm::EDGetTokenT < L1GctEtMissCollection >	l1EtMissToken_
edm::EDGetTokenT < L1GctEtTotalCollection >	l1EtTotalToken_
MonitorElement *	l1GctAllEmOccRankBx_
MonitorElement *	l1GctAllJetsEtEtaPhi_
MonitorElement *	l1GctAllJetsOccEtaPhi_
MonitorElement *	l1GctAllJetsOccRankBx_
MonitorElement *	l1GctCenJetsEtEtaPhi_
MonitorElement *	l1GctCenJetsOccEtaPhi_
MonitorElement *	l1GctCenJetsRank_
MonitorElement *	l1GctEtHad_
MonitorElement *	l1GctEtHadOccBx_
MonitorElement *	l1GctEtHadOf_
MonitorElement *	l1GctEtMiss_
MonitorElement *	l1GctEtMissHtMissCorr_
MonitorElement *	l1GctEtMissHtMissCorrPhi_
MonitorElement *	l1GctEtMissOccBx_
MonitorElement *	l1GctEtMissOf_
MonitorElement *	l1GctEtMissPhi_
MonitorElement *	l1GctEtTotal_
MonitorElement *	l1GctEtTotalEtHadCorr_
MonitorElement *	l1GctEtTotalOccBx_
MonitorElement *	l1GctEtTotalOf_
MonitorElement *	l1GctForJetsEtEtaPhi_
MonitorElement *	l1GctForJetsOccEtaPhi_
MonitorElement *	l1GctForJetsRank_
MonitorElement *	l1GctHFRing1ETSumNegEta_
MonitorElement *	l1GctHFRing1ETSumPosEta_
MonitorElement *	l1GctHFRing1PosEtaNegEta_
MonitorElement *	l1GctHFRing1TowerCountNegEta_
MonitorElement *	l1GctHFRing1TowerCountPosEta_
MonitorElement *	l1GctHFRing1TowerCountPosEtaNegEta_
MonitorElement *	l1GctHFRing2ETSumNegEta_
MonitorElement *	l1GctHFRing2ETSumPosEta_
MonitorElement *	l1GctHFRing2PosEtaNegEta_
MonitorElement *	l1GctHFRing2TowerCountNegEta_
MonitorElement *	l1GctHFRing2TowerCountPosEta_
MonitorElement *	l1GctHFRing2TowerCountPosEtaNegEta_
MonitorElement *	l1GctHFRingETSumOccBx_
MonitorElement *	l1GctHFRingRatioNegEta_
MonitorElement *	l1GctHFRingRatioPosEta_
MonitorElement *	l1GctHFRingTowerCountOccBx_
MonitorElement *	l1GctHtMiss_
MonitorElement *	l1GctHtMissOccBx_
MonitorElement *	l1GctHtMissOf_
MonitorElement *	l1GctHtMissPhi_
MonitorElement *	l1GctIsoEmOccEtaPhi_
MonitorElement *	l1GctIsoEmRank_
MonitorElement *	l1GctIsoEmRankEtaPhi_
MonitorElement *	l1GctIsoTauJetsEtEtaPhi_
MonitorElement *	l1GctIsoTauJetsOccEtaPhi_
MonitorElement *	l1GctIsoTauJetsRank_
MonitorElement *	l1GctNonIsoEmOccEtaPhi_
MonitorElement *	l1GctNonIsoEmRank_
MonitorElement *	l1GctNonIsoEmRankEtaPhi_
MonitorElement *	l1GctTauJetsEtEtaPhi_
MonitorElement *	l1GctTauJetsOccEtaPhi_
MonitorElement *	l1GctTauJetsRank_
edm::EDGetTokenT < L1GctHFBitCountsCollection >	l1HFCountsToken_
edm::EDGetTokenT < L1GctHtMissCollection >	l1HtMissToken_
std::ofstream	logFile_
MonitorElement *	lumisecId_
bool	m_stage1_layer2_
	filter TriggerType More...
bool	monitorDaemon_
int	nev_
std::string	outputFile_
MonitorElement *	runId_
MonitorElement *	triggerType_
bool	verbose_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 124 of file L1TGCT.h.

Constructor & Destructor Documentation

L1TGCT::L1TGCT

(

const edm::ParameterSet &

ps

)

Definition at line 60 of file L1TGCT.cc.

References gctCenJetsSourceToken_, gctEnergySumsSourceToken_, gctForJetsSourceToken_, gctIsoEmSourceToken_, gctIsoTauJetsSourceToken_, gctNonIsoEmSourceToken_, gctTauJetsSourceToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), l1EtHadToken_, l1EtMissToken_, l1EtTotalToken_, l1HFCountsToken_, l1HtMissToken_, m_stage1_layer2_, outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

                                         :
  gctCenJetsSource_(ps.getParameter<edm::InputTag>("gctCentralJetsSource")),
  gctForJetsSource_(ps.getParameter<edm::InputTag>("gctForwardJetsSource")),
  gctTauJetsSource_(ps.getParameter<edm::InputTag>("gctTauJetsSource")),
  gctIsoTauJetsSource_(ps.getParameter<edm::InputTag>("gctIsoTauJetsSource")),
  gctEnergySumsSource_(ps.getParameter<edm::InputTag>("gctEnergySumsSource")),
  gctIsoEmSource_(ps.getParameter<edm::InputTag>("gctIsoEmSource")),
  gctNonIsoEmSource_(ps.getParameter<edm::InputTag>("gctNonIsoEmSource")),
  m_stage1_layer2_(ps.getParameter<bool>("stage1_layer2_")),
  filterTriggerType_ (ps.getParameter< int >("filterTriggerType"))
{

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

  if (verbose_)
    edm::LogInfo("L1TGCT") << "L1TGCT: constructor...." << std::endl;

  outputFile_ = ps.getUntrackedParameter < std::string > ("outputFile", "");
  if (outputFile_.size() != 0) {
    edm::LogInfo("L1TGCT") << "L1T Monitoring histograms will be saved to "
                           << outputFile_ << std::endl;
  }

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

  //set Token(-s)
  gctIsoEmSourceToken_ = consumes<L1GctEmCandCollection>(ps.getParameter<edm::InputTag>("gctIsoEmSource"));
  gctNonIsoEmSourceToken_ = consumes<L1GctEmCandCollection>(ps.getParameter<edm::InputTag>("gctNonIsoEmSource"));
  gctCenJetsSourceToken_ = consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctCentralJetsSource"));
  gctForJetsSourceToken_ = consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctForwardJetsSource"));
  gctTauJetsSourceToken_ = consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctTauJetsSource"));
  if(m_stage1_layer2_ == true){
      gctIsoTauJetsSourceToken_=consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctIsoTauJetsSource"));
  }
  gctEnergySumsSourceToken_ = consumes<L1GctHFRingEtSumsCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
  l1HFCountsToken_ = consumes<L1GctHFBitCountsCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
  l1EtMissToken_ = consumes<L1GctEtMissCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
  l1HtMissToken_ = consumes<L1GctHtMissCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
  l1EtHadToken_ = consumes<L1GctEtHadCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
  l1EtTotalToken_ = consumes<L1GctEtTotalCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
}

L1TGCT::~L1TGCT ( )

virtual

Definition at line 106 of file L1TGCT.cc.

{
}

Member Function Documentation

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

protectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 236 of file L1TGCT.cc.

References edm::EventBase::experimentType(), MonitorElement::Fill(), filterTriggerType_, gctCenJetsSource_, gctCenJetsSourceToken_, gctEnergySumsSource_, gctEnergySumsSourceToken_, gctForJetsSource_, gctForJetsSourceToken_, gctIsoEmSource_, gctIsoEmSourceToken_, gctIsoTauJetsSource_, gctIsoTauJetsSourceToken_, gctNonIsoEmSource_, gctNonIsoEmSourceToken_, gctTauJetsSource_, gctTauJetsSourceToken_, edm::Event::getByToken(), i, edm::EventBase::isRealData(), edm::HandleBase::isValid(), l1EtHadToken_, l1EtMissToken_, l1EtTotalToken_, l1GctAllEmOccRankBx_, l1GctAllJetsEtEtaPhi_, l1GctAllJetsOccEtaPhi_, l1GctAllJetsOccRankBx_, l1GctCenJetsEtEtaPhi_, l1GctCenJetsOccEtaPhi_, l1GctCenJetsRank_, l1GctEtHad_, l1GctEtHadOccBx_, l1GctEtHadOf_, l1GctEtMiss_, l1GctEtMissHtMissCorr_, l1GctEtMissHtMissCorrPhi_, l1GctEtMissOccBx_, l1GctEtMissOf_, l1GctEtMissPhi_, l1GctEtTotal_, l1GctEtTotalEtHadCorr_, l1GctEtTotalOccBx_, l1GctEtTotalOf_, l1GctForJetsEtEtaPhi_, l1GctForJetsOccEtaPhi_, l1GctForJetsRank_, l1GctHFRing1ETSumNegEta_, l1GctHFRing1ETSumPosEta_, l1GctHFRing1PosEtaNegEta_, l1GctHFRing1TowerCountNegEta_, l1GctHFRing1TowerCountPosEta_, l1GctHFRing1TowerCountPosEtaNegEta_, l1GctHFRing2ETSumNegEta_, l1GctHFRing2ETSumPosEta_, l1GctHFRing2PosEtaNegEta_, l1GctHFRing2TowerCountNegEta_, l1GctHFRing2TowerCountPosEta_, l1GctHFRing2TowerCountPosEtaNegEta_, l1GctHFRingETSumOccBx_, l1GctHFRingRatioNegEta_, l1GctHFRingRatioPosEta_, l1GctHFRingTowerCountOccBx_, l1GctHtMiss_, l1GctHtMissOccBx_, l1GctHtMissOf_, l1GctHtMissPhi_, l1GctIsoEmOccEtaPhi_, l1GctIsoEmRank_, l1GctIsoEmRankEtaPhi_, l1GctIsoTauJetsEtEtaPhi_, l1GctIsoTauJetsOccEtaPhi_, l1GctIsoTauJetsRank_, l1GctNonIsoEmOccEtaPhi_, l1GctNonIsoEmRank_, l1GctNonIsoEmRankEtaPhi_, l1GctTauJetsEtEtaPhi_, l1GctTauJetsOccEtaPhi_, l1GctTauJetsRank_, l1HFCountsToken_, l1HtMissToken_, m_stage1_layer2_, objects.METAnalyzer::met, Vispa.Plugins.EdmBrowser.EdmDataAccessor::ne(), nev_, HLT_25ns14e33_v1_cff::triggerType, triggerType_, and verbose_.

{
  nev_++;
  if (verbose_) {
    edm::LogInfo("L1TGCT") << "L1TGCT: analyze...." << std::endl;
  }

  
  // filter according trigger type
  //  enum ExperimentType {
  //        Undefined          =  0,
  //        PhysicsTrigger     =  1,
  //        CalibrationTrigger =  2,
  //        RandomTrigger      =  3,
  //        Reserved           =  4,
  //        TracedEvent        =  5,
  //        TestTrigger        =  6,
  //        ErrorTrigger       = 15

  // fill a histogram with the trigger type, for normalization fill also last bin
  // ErrorTrigger + 1
  double triggerType = static_cast<double> (e.experimentType()) + 0.001;
  double triggerTypeLast = static_cast<double> (edm::EventAuxiliary::ExperimentType::ErrorTrigger)
                          + 0.001;
  triggerType_->Fill(triggerType);
  triggerType_->Fill(triggerTypeLast + 1);

  // filter only if trigger type is greater than 0, negative values disable filtering
  if (filterTriggerType_ >= 0) {

      // now filter, for real data only
      if (e.isRealData()) {
          if (!(e.experimentType() == filterTriggerType_)) {

              edm::LogInfo("L1TGCT") << "\n Event of TriggerType "
                      << e.experimentType() << " rejected" << std::endl;
              return;

          }
      }

  }

  // Get all the collections
  edm::Handle < L1GctEmCandCollection > l1IsoEm;
  edm::Handle < L1GctEmCandCollection > l1NonIsoEm;
  edm::Handle < L1GctJetCandCollection > l1CenJets;
  edm::Handle < L1GctJetCandCollection > l1ForJets;
  edm::Handle < L1GctJetCandCollection > l1TauJets;
  if(m_stage1_layer2_ == true) {
    edm::Handle < L1GctJetCandCollection > l1IsoTauJets;
    e.getByToken(gctIsoTauJetsSourceToken_, l1IsoTauJets);
  }
  edm::Handle < L1GctHFRingEtSumsCollection > l1HFSums; 
  edm::Handle < L1GctHFBitCountsCollection > l1HFCounts;
  edm::Handle < L1GctEtMissCollection >  l1EtMiss;
  edm::Handle < L1GctHtMissCollection >  l1HtMiss;
  edm::Handle < L1GctEtHadCollection >   l1EtHad;
  edm::Handle < L1GctEtTotalCollection > l1EtTotal;

  e.getByToken(gctIsoEmSourceToken_, l1IsoEm);
  e.getByToken(gctNonIsoEmSourceToken_, l1NonIsoEm);
  e.getByToken(gctCenJetsSourceToken_, l1CenJets);
  e.getByToken(gctForJetsSourceToken_, l1ForJets);
  e.getByToken(gctTauJetsSourceToken_, l1TauJets);
  e.getByToken(gctEnergySumsSourceToken_, l1HFSums);
  e.getByToken(l1HFCountsToken_, l1HFCounts);
  e.getByToken(l1EtMissToken_, l1EtMiss);
  e.getByToken(l1HtMissToken_, l1HtMiss);
  e.getByToken(l1EtHadToken_, l1EtHad);
  e.getByToken(l1EtTotalToken_, l1EtTotal);

  // Fill histograms

  // Central jets
  if (l1CenJets.isValid()) {
    for (L1GctJetCandCollection::const_iterator cj = l1CenJets->begin();cj != l1CenJets->end(); cj++) {
      // only plot central BX
      if (cj->bx()==0) {
        l1GctCenJetsRank_->Fill(cj->rank());
        // only plot eta and phi maps for non-zero candidates
        if (cj->rank()) {
          l1GctAllJetsEtEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi(),cj->rank());
          l1GctAllJetsOccEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi());
          l1GctCenJetsEtEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi(),cj->rank());
          l1GctCenJetsOccEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi());
        }
      }
      if (cj->rank()) l1GctAllJetsOccRankBx_->Fill(cj->bx(),cj->rank()); // for all BX
    }
  } else {    
    edm::LogWarning("DataNotFound") << " Could not find l1CenJets label was " << gctCenJetsSource_ ;
  }

  // Forward jets
  if (l1ForJets.isValid()) {
    for (L1GctJetCandCollection::const_iterator fj = l1ForJets->begin(); fj != l1ForJets->end(); fj++) {
      // only plot central BX
      if (fj->bx()==0) {
        l1GctForJetsRank_->Fill(fj->rank());
        // only plot eta and phi maps for non-zero candidates
        if (fj->rank()) {
          l1GctAllJetsEtEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi(),fj->rank());
          l1GctAllJetsOccEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi());
          l1GctForJetsEtEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi(),fj->rank());
          l1GctForJetsOccEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi());    
        }
      }
      if (fj->rank()) l1GctAllJetsOccRankBx_->Fill(fj->bx(),fj->rank()); // for all BX
    }
  } else {    
    edm::LogWarning("DataNotFound") << " Could not find l1ForJets label was " << gctForJetsSource_ ;
  }

  // Tau jets
  if (l1TauJets.isValid()) {
    for (L1GctJetCandCollection::const_iterator tj = l1TauJets->begin(); tj != l1TauJets->end(); tj++) {
      // only plot central BX
      if (tj->bx()==0) {
        l1GctTauJetsRank_->Fill(tj->rank());
        // only plot eta and phi maps for non-zero candidates
        if (tj->rank()) {
          l1GctTauJetsEtEtaPhi_->Fill(tj->regionId().ieta(),tj->regionId().iphi(),tj->rank());
          l1GctTauJetsOccEtaPhi_->Fill(tj->regionId().ieta(),tj->regionId().iphi());
        }
      }
      if (tj->rank()) l1GctAllJetsOccRankBx_->Fill(tj->bx(),tj->rank()); // for all BX
    }
  } else {    
    edm::LogWarning("DataNotFound") << " Could not find l1TauJets label was " << gctTauJetsSource_ ;
  }

   // IsoTau jets
  if (m_stage1_layer2_ == true) {
   edm::Handle < L1GctJetCandCollection > l1IsoTauJets;
   e.getByToken(gctIsoTauJetsSourceToken_, l1IsoTauJets);
   if (l1IsoTauJets.isValid()) {
    for (L1GctJetCandCollection::const_iterator itj = l1IsoTauJets->begin(); itj != l1IsoTauJets->end(); itj++) {
      // only plot central BX
      if (itj->bx()==0) {
        l1GctIsoTauJetsRank_->Fill(itj->rank());
        // only plot eta and phi maps for non-zero candidates
        if (itj->rank()) {
          l1GctIsoTauJetsEtEtaPhi_->Fill(itj->regionId().ieta(),itj->regionId().iphi(),itj->rank());
          l1GctIsoTauJetsOccEtaPhi_->Fill(itj->regionId().ieta(),itj->regionId().iphi());
        }
      }
      if (itj->rank()) l1GctAllJetsOccRankBx_->Fill(itj->bx(),itj->rank()); // for all BX
    }
   } else {    
    edm::LogWarning("DataNotFound") << " Could not find l1IsoTauJets label was " << gctIsoTauJetsSource_ ;
   }
  }

  
  // Missing ET
  if (l1EtMiss.isValid()) { 
    for (L1GctEtMissCollection::const_iterator met = l1EtMiss->begin(); met != l1EtMiss->end(); met++) {
      // only plot central BX
      if (met->bx()==0) {
        if (met->overFlow() == 0 && met->et() > 0) {
          //Avoid problems with met=0 candidates affecting MET_PHI plots
          l1GctEtMiss_->Fill(met->et());
          l1GctEtMissPhi_->Fill(met->phi());
        }
        l1GctEtMissOf_->Fill(met->overFlow());
      }
      if (met->overFlow() == 0 && met->et() > 0) l1GctEtMissOccBx_->Fill(met->bx(),met->et()); // for all BX
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1EtMiss label was " << gctEnergySumsSource_ ;    
  }

  // Missing HT
  if (l1HtMiss.isValid()) { 
    for (L1GctHtMissCollection::const_iterator mht = l1HtMiss->begin(); mht != l1HtMiss->end(); mht++) {
      // only plot central BX
      if (mht->bx()==0) {
        if (mht->overFlow() == 0 && mht->et() > 0) {
          //Avoid problems with mht=0 candidates affecting MHT_PHI plots
          l1GctHtMiss_->Fill(mht->et());
          l1GctHtMissPhi_->Fill(mht->phi());
        }
        l1GctHtMissOf_->Fill(mht->overFlow());
      }
      if (mht->overFlow() == 0 && mht->et() > 0) l1GctHtMissOccBx_->Fill(mht->bx(),mht->et()); // for all BX
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1HtMiss label was " << gctEnergySumsSource_ ;    
  }

  // Missing ET HT correlations
  if (l1HtMiss.isValid() && l1EtMiss.isValid()) { 
    if (l1HtMiss->size() == l1EtMiss->size()) {
      for (unsigned i=0; i<l1HtMiss->size(); i++) {
        if (l1HtMiss->at(i).overFlow() == 0 && l1EtMiss->at(i).overFlow() == 0 && 
            l1HtMiss->at(i).bx() == 0 && l1EtMiss->at(i).bx() == 0) {
          // Avoid problems overflows and only plot central BX
          l1GctEtMissHtMissCorr_->Fill(l1EtMiss->at(i).et(),l1HtMiss->at(i).et());
          if (l1HtMiss->at(i).et() && l1EtMiss->at(i).et()){ // Don't plot phi if one or other is zero
            l1GctEtMissHtMissCorrPhi_->Fill(l1EtMiss->at(i).phi(),l1HtMiss->at(i).phi());
          }
        }
      }
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1EtMiss or l1HtMiss label was " << gctEnergySumsSource_ ;    
  }

  // HT 
  if (l1EtHad.isValid()) {
    for (L1GctEtHadCollection::const_iterator ht = l1EtHad->begin(); ht != l1EtHad->end(); ht++) {
      // only plot central BX
      if (ht->bx()==0) {
        l1GctEtHad_->Fill(ht->et());
        l1GctEtHadOf_->Fill(ht->overFlow());
      }
      l1GctEtHadOccBx_->Fill(ht->bx(),ht->et()); // for all BX
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1EtHad label was " << gctEnergySumsSource_ ;    
  }

  // Total ET
  if (l1EtTotal.isValid()) {
    for (L1GctEtTotalCollection::const_iterator et = l1EtTotal->begin(); et != l1EtTotal->end(); et++) {
      // only plot central BX
      if (et->bx()==0) {
        l1GctEtTotal_->Fill(et->et());
        l1GctEtTotalOf_->Fill(et->overFlow());
      }
      l1GctEtTotalOccBx_->Fill(et->bx(),et->et()); // for all BX
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1EtTotal label was " << gctEnergySumsSource_ ;    
  }

  // Total ET HT correlations
  if (l1EtTotal.isValid() && l1EtHad.isValid()) { 
    if (l1EtTotal->size() == l1EtHad->size()) {
      for (unsigned i=0; i<l1EtHad->size(); i++) {
        if (l1EtHad->at(i).overFlow() == 0 && l1EtTotal->at(i).overFlow() == 0 && 
            l1EtHad->at(i).bx() == 0 && l1EtTotal->at(i).bx() == 0) {
          // Avoid problems overflows and only plot central BX
          l1GctEtTotalEtHadCorr_->Fill(l1EtTotal->at(i).et(),l1EtHad->at(i).et());
        }
      }
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1EtTotal or l1EtHad label was " << gctEnergySumsSource_ ;    
  }

  //HF Ring Et Sums
  if (l1HFSums.isValid()) {
    for (L1GctHFRingEtSumsCollection::const_iterator hfs=l1HFSums->begin(); hfs!=l1HFSums->end(); hfs++){ 
      // only plot central BX
      if (hfs->bx()==0) {
        // Individual ring Et sums
        l1GctHFRing1ETSumPosEta_->Fill(hfs->etSum(0));
        l1GctHFRing1ETSumNegEta_->Fill(hfs->etSum(1));
        l1GctHFRing2ETSumPosEta_->Fill(hfs->etSum(2));
        l1GctHFRing2ETSumNegEta_->Fill(hfs->etSum(3));
        // Ratio of ring Et sums
        if (hfs->etSum(2)!=0) l1GctHFRingRatioPosEta_->Fill((hfs->etSum(0))/(hfs->etSum(2)));
        if (hfs->etSum(3)!=0) l1GctHFRingRatioNegEta_->Fill((hfs->etSum(1))/(hfs->etSum(3)));
        // Correlate positive and neagative eta
        l1GctHFRing1PosEtaNegEta_->Fill(hfs->etSum(0),hfs->etSum(1));
        l1GctHFRing2PosEtaNegEta_->Fill(hfs->etSum(2),hfs->etSum(3));
      }
      // Occupancy vs BX
      for (unsigned i=0; i<4; i++){
        l1GctHFRingETSumOccBx_->Fill(hfs->bx(),hfs->etSum(i));
      }
    }
  } else {    
    edm::LogWarning("DataNotFound") << " Could not find l1HFSums label was " << gctEnergySumsSource_ ;
  }

  // HF Ring Counts
  if (l1HFCounts.isValid()) {
    for (L1GctHFBitCountsCollection::const_iterator hfc=l1HFCounts->begin(); hfc!=l1HFCounts->end(); hfc++){ 
      // only plot central BX
      if (hfc->bx()==0) {
        // Individual ring counts
        l1GctHFRing1TowerCountPosEta_->Fill(hfc->bitCount(0));
        l1GctHFRing1TowerCountNegEta_->Fill(hfc->bitCount(1));
        l1GctHFRing2TowerCountPosEta_->Fill(hfc->bitCount(2));
        l1GctHFRing2TowerCountNegEta_->Fill(hfc->bitCount(3));
        // Correlate positive and negative eta
        l1GctHFRing1TowerCountPosEtaNegEta_->Fill(hfc->bitCount(0),hfc->bitCount(1));
        l1GctHFRing2TowerCountPosEtaNegEta_->Fill(hfc->bitCount(2),hfc->bitCount(3));
      }
      // Occupancy vs BX
      for (unsigned i=0; i<4; i++){
        l1GctHFRingTowerCountOccBx_->Fill(hfc->bx(),hfc->bitCount(i));
      }
    }
  } else {    
    edm::LogWarning("DataNotFound") << " Could not find l1HFCounts label was " << gctEnergySumsSource_ ;
  }
  
  // Isolated EM
  if (l1IsoEm.isValid()) {
    for (L1GctEmCandCollection::const_iterator ie=l1IsoEm->begin(); ie!=l1IsoEm->end(); ie++) {
      // only plot central BX
      if (ie->bx()==0) {
        l1GctIsoEmRank_->Fill(ie->rank());
        // only plot eta and phi maps for non-zero candidates
        if (ie->rank()){ 
          l1GctIsoEmRankEtaPhi_->Fill(ie->regionId().ieta(),ie->regionId().iphi(),ie->rank());
          l1GctIsoEmOccEtaPhi_->Fill(ie->regionId().ieta(),ie->regionId().iphi());
        }
      }
      if (ie->rank()) l1GctAllEmOccRankBx_->Fill(ie->bx(),ie->rank()); // for all BX
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1IsoEm label was " << gctIsoEmSource_ ;
  } 

  // Non-isolated EM
  if (l1NonIsoEm.isValid()) { 
    for (L1GctEmCandCollection::const_iterator ne=l1NonIsoEm->begin(); ne!=l1NonIsoEm->end(); ne++) {
      // only plot central BX
      if (ne->bx()==0) {
        l1GctNonIsoEmRank_->Fill(ne->rank());
        // only plot eta and phi maps for non-zero candidates
        if (ne->rank()){ 
          l1GctNonIsoEmRankEtaPhi_->Fill(ne->regionId().ieta(),ne->regionId().iphi(),ne->rank());
          l1GctNonIsoEmOccEtaPhi_->Fill(ne->regionId().ieta(),ne->regionId().iphi());
        }
      }
      if (ne->rank()) l1GctAllEmOccRankBx_->Fill(ne->bx(),ne->rank()); // for all BX
    }
  } else {
    edm::LogWarning("DataNotFound") << " Could not find l1NonIsoEm label was " << gctNonIsoEmSource_ ;
  }
  edm::LogInfo("L1TGCT") << "L1TGCT: end job...." << std::endl;
  edm::LogInfo("EndJob") << "analyzed " << nev_ << " events";
}

void L1TGCT::beginLuminosityBlock ( const edm::LuminosityBlock &  iLumi, const edm::EventSetup &  evSetup  )

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 232 of file L1TGCT.cc.

                                                                                               {
  //lumisecId_->Fill(iLumi.id().luminosityBlock());
}

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 110 of file L1TGCT.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookInt(), BXBINS, BXMAX, BXMIN, EMETABINS, EMETAMAX, EMETAMIN, MonitorElement::Fill(), JETETABINS, JETETAMAX, JETETAMIN, l1GctAllEmOccRankBx_, l1GctAllJetsEtEtaPhi_, l1GctAllJetsOccEtaPhi_, l1GctAllJetsOccRankBx_, l1GctCenJetsEtEtaPhi_, l1GctCenJetsOccEtaPhi_, l1GctCenJetsRank_, l1GctEtHad_, l1GctEtHadOccBx_, l1GctEtHadOf_, l1GctEtMiss_, l1GctEtMissHtMissCorr_, l1GctEtMissHtMissCorrPhi_, l1GctEtMissOccBx_, l1GctEtMissOf_, l1GctEtMissPhi_, l1GctEtTotal_, l1GctEtTotalEtHadCorr_, l1GctEtTotalOccBx_, l1GctEtTotalOf_, l1GctForJetsEtEtaPhi_, l1GctForJetsOccEtaPhi_, l1GctForJetsRank_, l1GctHFRing1ETSumNegEta_, l1GctHFRing1ETSumPosEta_, l1GctHFRing1PosEtaNegEta_, l1GctHFRing1TowerCountNegEta_, l1GctHFRing1TowerCountPosEta_, l1GctHFRing1TowerCountPosEtaNegEta_, l1GctHFRing2ETSumNegEta_, l1GctHFRing2ETSumPosEta_, l1GctHFRing2PosEtaNegEta_, l1GctHFRing2TowerCountNegEta_, l1GctHFRing2TowerCountPosEta_, l1GctHFRing2TowerCountPosEtaNegEta_, l1GctHFRingETSumOccBx_, l1GctHFRingRatioNegEta_, l1GctHFRingRatioPosEta_, l1GctHFRingTowerCountOccBx_, l1GctHtMiss_, l1GctHtMissOccBx_, l1GctHtMissOf_, l1GctHtMissPhi_, l1GctIsoEmOccEtaPhi_, l1GctIsoEmRank_, l1GctIsoEmRankEtaPhi_, l1GctIsoTauJetsEtEtaPhi_, l1GctIsoTauJetsOccEtaPhi_, l1GctIsoTauJetsRank_, l1GctNonIsoEmOccEtaPhi_, l1GctNonIsoEmRank_, l1GctNonIsoEmRankEtaPhi_, l1GctTauJetsEtEtaPhi_, l1GctTauJetsOccEtaPhi_, l1GctTauJetsRank_, lumisecId_, m_stage1_layer2_, METPHIBINS, METPHIMAX, METPHIMIN, MHTPHIBINS, MHTPHIMAX, MHTPHIMIN, nev_, OFBINS, OFMAX, OFMIN, PHIBINS, PHIMAX, PHIMIN, R12BINS, R12MAX, R12MIN, R3BINS, R3MAX, R3MIN, R5BINS, R5MAX, R5MIN, R6BINS, R6MAX, R6MIN, R7BINS, R7MAX, R7MIN, runId_, DQMStore::IBooker::setCurrentFolder(), and triggerType_.

{

  nev_ = 0;

  if (m_stage1_layer2_ == true){  
    ibooker.setCurrentFolder("L1T/stage1layer2");
  }  
  if (m_stage1_layer2_ == false){
    ibooker.setCurrentFolder("L1T/L1TGCT");
  }

  runId_     = ibooker.bookInt("iRun");
  runId_->Fill(-1);
  lumisecId_ = ibooker.bookInt("iLumiSection");
  lumisecId_->Fill(-1);
  
  triggerType_ =ibooker.book1D("TriggerType", "TriggerType", 17, -0.5, 16.5);

  l1GctAllJetsEtEtaPhi_ = ibooker.book2D("AllJetsEtEtaPhi", "CENTRAL AND FORWARD JET E_{T}",JETETABINS, JETETAMIN, JETETAMAX,PHIBINS, PHIMIN, PHIMAX);
  l1GctCenJetsEtEtaPhi_ = ibooker.book2D("CenJetsEtEtaPhi", "CENTRAL JET E_{T}",JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX); 
  l1GctForJetsEtEtaPhi_ = ibooker.book2D("ForJetsEtEtaPhi", "FORWARD JET E_{T}", JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX); 
  l1GctTauJetsEtEtaPhi_ = ibooker.book2D("TauJetsEtEtaPhi", "TAU JET E_{T}", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
  if (m_stage1_layer2_ == true){
    l1GctIsoTauJetsEtEtaPhi_ = ibooker.book2D("IsoTauJetsEtEtaPhi", "ISOTAU JET E_{T}", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
  }
  l1GctIsoEmRankEtaPhi_ = ibooker.book2D("IsoEmRankEtaPhi", "ISO EM E_{T}", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);            
  l1GctNonIsoEmRankEtaPhi_ = ibooker.book2D("NonIsoEmRankEtaPhi", "NON-ISO EM E_{T}", EMETABINS, EMETAMIN, EMETAMAX,PHIBINS, PHIMIN, PHIMAX); 
  l1GctAllJetsOccEtaPhi_ = ibooker.book2D("AllJetsOccEtaPhi", "CENTRAL AND FORWARD JET OCCUPANCY", JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX);
  l1GctCenJetsOccEtaPhi_ = ibooker.book2D("CenJetsOccEtaPhi", "CENTRAL JET OCCUPANCY", JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX); 
  l1GctForJetsOccEtaPhi_ = ibooker.book2D("ForJetsOccEtaPhi", "FORWARD JET OCCUPANCY",JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX);
  l1GctTauJetsOccEtaPhi_ = ibooker.book2D("TauJetsOccEtaPhi", "TAU JET OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
  if (m_stage1_layer2_ == true){
    l1GctIsoTauJetsOccEtaPhi_ = ibooker.book2D("IsoTauJetsOccEtaPhi", "ISOTAU JET OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
  }
  l1GctIsoEmOccEtaPhi_ = ibooker.book2D("IsoEmOccEtaPhi", "ISO EM OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
  l1GctNonIsoEmOccEtaPhi_ = ibooker.book2D("NonIsoEmOccEtaPhi", "NON-ISO EM OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX); 
  

    //HF Ring stuff

  l1GctHFRing1TowerCountPosEtaNegEta_ = ibooker.book2D("HFRing1TowerCountCorr", "HF RING1 TOWER COUNT CORRELATION +/-  #eta", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
  l1GctHFRing2TowerCountPosEtaNegEta_ = ibooker.book2D("HFRing2TowerCountCorr", "HF RING2 TOWER COUNT CORRELATION +/-  #eta", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
  
  l1GctHFRing1TowerCountPosEta_ = ibooker.book1D("HFRing1TowerCountPosEta", "HF RING1 TOWER COUNT  #eta  +", R3BINS, R3MIN, R3MAX);
  l1GctHFRing1TowerCountNegEta_ = ibooker.book1D("HFRing1TowerCountNegEta", "HF RING1 TOWER COUNT  #eta  -", R3BINS, R3MIN, R3MAX);
  l1GctHFRing2TowerCountPosEta_ = ibooker.book1D("HFRing2TowerCountPosEta", "HF RING2 TOWER COUNT  #eta  +", R3BINS, R3MIN, R3MAX);
  l1GctHFRing2TowerCountNegEta_ = ibooker.book1D("HFRing2TowerCountNegEta", "HF RING2 TOWER COUNT  #eta  -", R3BINS, R3MIN, R3MAX);

  l1GctHFRingTowerCountOccBx_ = ibooker.book2D("HFRingTowerCountOccBx", "HF RING TOWER COUNT PER BX",BXBINS, BXMIN, BXMAX, R3BINS, R3MIN, R3MAX);
  
  if (m_stage1_layer2_ == false){
    l1GctHFRing1PosEtaNegEta_ = ibooker.book2D("HFRing1Corr", "HF RING1 E_{T} CORRELATION +/-  #eta",  R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX); 
    l1GctHFRing2PosEtaNegEta_ = ibooker.book2D("HFRing2Corr", "HF RING2 E_{T} CORRELATION +/-  #eta", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
    l1GctHFRing1ETSumPosEta_ = ibooker.book1D("HFRing1ETSumPosEta", "HF RING1 E_{T}  #eta  +", R3BINS, R3MIN, R3MAX);
    l1GctHFRing1ETSumNegEta_ = ibooker.book1D("HFRing1ETSumNegEta", "HF RING1 E_{T}  #eta  -", R3BINS, R3MIN, R3MAX);
    l1GctHFRing2ETSumPosEta_ = ibooker.book1D("HFRing2ETSumPosEta", "HF RING2 E_{T}  #eta  +", R3BINS, R3MIN, R3MAX);
    l1GctHFRing2ETSumNegEta_ = ibooker.book1D("HFRing2ETSumNegEta", "HF RING2 E_{T}  #eta  -", R3BINS, R3MIN, R3MAX);
    l1GctHFRingETSumOccBx_ = ibooker.book2D("HFRingETSumOccBx", "HF RING E_{T} PER BX",BXBINS, BXMIN, BXMAX, R3BINS, R3MIN, R3MAX);
    l1GctHFRingRatioPosEta_  = ibooker.book1D("HFRingRatioPosEta", "HF RING E_{T} RATIO  #eta  +", R5BINS, R5MIN, R5MAX);
    l1GctHFRingRatioNegEta_  = ibooker.book1D("HFRingRatioNegEta", "HF RING E_{T} RATIO  #eta  -", R5BINS, R5MIN, R5MAX);
  }

  if (m_stage1_layer2_ == true){
    l1GctHFRing1PosEtaNegEta_ = ibooker.book2D("IsoTau 1 2 Corr", "IsoTau 1 IsoTau 2 E_{T} CORRELATION",  R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX); 
    l1GctHFRing2PosEtaNegEta_ = ibooker.book2D("IsoTau 3 4 Corr", "IsoTau 3 IsoTau 4 CORRELATION", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
    l1GctHFRing1ETSumPosEta_ = ibooker.book1D("Iso Tau 1 Et", "Isolated Tau1 E_{T}", 9, -0.5, 8.5);
    l1GctHFRing1ETSumNegEta_ = ibooker.book1D("Iso Tau 2 Et", "Isolated Tau2 E_{T}", 9, -0.5, 8.5);
    l1GctHFRing2ETSumPosEta_ = ibooker.book1D("Iso Tau 3 Et", "Isolated Tau3 E_{T}", 9, -0.5, 8.5);
    l1GctHFRing2ETSumNegEta_ = ibooker.book1D("Iso Tau 4 Et", "Isolated Tau4 E_{T}", 9, -0.5, 8.5);
    l1GctHFRingETSumOccBx_ = ibooker.book2D("IsoTau HFRingSum OccBx", "Iso Tau PER BX",BXBINS, BXMIN, BXMAX, R3BINS, R3MIN, R3MAX);
    l1GctHFRingRatioPosEta_  = ibooker.book1D("IsoTau Ratio 1 2", "IsoTau E_{T} RATIO", 9, -0.5, 8.5);
    l1GctHFRingRatioNegEta_  = ibooker.book1D("IsoTau Ratio 1 2", "IsoTau E_{T} RATIO", 9, -0.5, 8.5);
  }
    
    // Rank histograms
  l1GctCenJetsRank_  = ibooker.book1D("CenJetsRank", "CENTRAL JET E_{T}", R6BINS, R6MIN, R6MAX);
  l1GctForJetsRank_  = ibooker.book1D("ForJetsRank", "FORWARD JET E_{T}", R6BINS, R6MIN, R6MAX);
  l1GctTauJetsRank_  = ibooker.book1D("TauJetsRank", "TAU JET E_{T}", R6BINS, R6MIN, R6MAX);
  if (m_stage1_layer2_ == true){
    l1GctIsoTauJetsRank_ = ibooker.book1D("IsoTauJetsRank", "ISOTAU JET E_{T}", R6BINS, R6MIN, R6MAX);
  }
  l1GctIsoEmRank_    = ibooker.book1D("IsoEmRank", "ISO EM E_{T}", R6BINS, R6MIN, R6MAX);
  l1GctNonIsoEmRank_ = ibooker.book1D("NonIsoEmRank", "NON-ISO EM E_{T}", R6BINS, R6MIN, R6MAX);

  l1GctAllJetsOccRankBx_ = ibooker.book2D("AllJetsOccRankBx","ALL JETS E_{T} PER BX",BXBINS,BXMIN,BXMAX,R6BINS,R6MIN,R6MAX);
  l1GctAllEmOccRankBx_   = ibooker.book2D("AllEmOccRankBx","ALL EM E_{T} PER BX",BXBINS,BXMIN,BXMAX,R6BINS,R6MIN,R6MAX);

    // Energy sums
  l1GctEtMiss_    = ibooker.book1D("EtMiss", "MET", R12BINS, R12MIN, R12MAX);
  l1GctEtMissPhi_ = ibooker.book1D("EtMissPhi", "MET  #phi", METPHIBINS, METPHIMIN, METPHIMAX);
  l1GctEtMissOf_  = ibooker.book1D("EtMissOf", "MET OVERFLOW", OFBINS, OFMIN, OFMAX);
  l1GctEtMissOccBx_ = ibooker.book2D("EtMissOccBx","MET PER BX",BXBINS,BXMIN,BXMAX,R12BINS,R12MIN,R12MAX);
  if (m_stage1_layer2_ == false) {
    l1GctHtMiss_    = ibooker.book1D("HtMiss", "MHT", R7BINS, R7MIN, R7MAX);
    l1GctHtMissPhi_ = ibooker.book1D("HtMissPhi", "MHT  #phi", MHTPHIBINS, MHTPHIMIN, MHTPHIMAX);
    l1GctHtMissOf_  = ibooker.book1D("HtMissOf", "MHT OVERFLOW", OFBINS, OFMIN, OFMAX);
    l1GctHtMissOccBx_ = ibooker.book2D("HtMissOccBx","MHT PER BX",BXBINS,BXMIN,BXMAX,R7BINS,R7MIN,R7MAX);
  }
  if (m_stage1_layer2_ == true) {   
    l1GctHtMiss_    = ibooker.book1D("HtMissHtTotal", "MHTHT", R7BINS, R7MIN, R7MAX);
    l1GctHtMissPhi_ = ibooker.book1D("HtMissHtTotal Phi", "MHTHT  #phi", MHTPHIBINS, MHTPHIMIN, MHTPHIMAX);
    l1GctHtMissOf_  = ibooker.book1D("HtMissHtTotal Of", "MHTHT OVERFLOW", OFBINS, OFMIN, OFMAX);
    l1GctHtMissOccBx_ = ibooker.book2D("HtMissHtTotal OccBx","MHTHT PER BX",BXBINS,BXMIN,BXMAX,R7BINS,R7MIN,R7MAX);
  }
  l1GctEtMissHtMissCorr_ = ibooker.book2D("EtMissHtMissCorr", "MET MHT CORRELATION", R6BINS, R12MIN, R12MAX, R6BINS, R7MIN, R7MAX); 
  l1GctEtMissHtMissCorrPhi_ = ibooker.book2D("EtMissHtMissPhiCorr", "MET MHT  #phi  CORRELATION", METPHIBINS, METPHIMIN, METPHIMAX, MHTPHIBINS, MHTPHIMIN, MHTPHIMAX);
  l1GctEtTotal_   = ibooker.book1D("EtTotal", "SUM E_{T}", R12BINS, R12MIN, R12MAX);
  l1GctEtTotalOf_ = ibooker.book1D("EtTotalOf", "SUM E_{T} OVERFLOW", OFBINS, OFMIN, OFMAX);
  l1GctEtTotalOccBx_ = ibooker.book2D("EtTotalOccBx","SUM E_{T} PER BX",BXBINS,BXMIN,BXMAX,R12BINS,R12MIN,R12MAX);
  l1GctEtHad_     = ibooker.book1D("EtHad", "H_{T}", R12BINS, R12MIN, R12MAX);
  l1GctEtHadOf_   = ibooker.book1D("EtHadOf", "H_{T} OVERFLOW", OFBINS, OFMIN, OFMAX);
  l1GctEtHadOccBx_ = ibooker.book2D("EtHadOccBx","H_{T} PER BX",BXBINS,BXMIN,BXMAX,R12BINS,R12MIN,R12MAX);
  l1GctEtTotalEtHadCorr_ = ibooker.book2D("EtTotalEtHadCorr", "Sum E_{T} H_{T} CORRELATION", R6BINS, R12MIN, R12MAX, R6BINS, R12MIN, R12MAX); 
  //}
}

void L1TGCT::dqmBeginRun ( const edm::Run &  iRrun, const edm::EventSetup &  evSetup  )

protectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 228 of file L1TGCT.cc.

                                                                          {
  //runId_->Fill(iRrun.id().run());
}

Member Data Documentation

int L1TGCT::filterTriggerType_

private

Definition at line 230 of file L1TGCT.h.

Referenced by analyze().

edm::InputTag L1TGCT::gctCenJetsSource_

private

Definition at line 220 of file L1TGCT.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctJetCandCollection> L1TGCT::gctCenJetsSourceToken_

private

Definition at line 235 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::InputTag L1TGCT::gctEnergySumsSource_

private

Definition at line 224 of file L1TGCT.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctHFRingEtSumsCollection> L1TGCT::gctEnergySumsSourceToken_

private

Definition at line 239 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::InputTag L1TGCT::gctForJetsSource_

private

Definition at line 221 of file L1TGCT.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctJetCandCollection> L1TGCT::gctForJetsSourceToken_

private

Definition at line 236 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::InputTag L1TGCT::gctIsoEmSource_

private

Definition at line 225 of file L1TGCT.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctEmCandCollection> L1TGCT::gctIsoEmSourceToken_

private

Definition at line 233 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::InputTag L1TGCT::gctIsoTauJetsSource_

private

Definition at line 223 of file L1TGCT.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctJetCandCollection> L1TGCT::gctIsoTauJetsSourceToken_

private

Definition at line 238 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::InputTag L1TGCT::gctNonIsoEmSource_

private

Definition at line 226 of file L1TGCT.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctEmCandCollection> L1TGCT::gctNonIsoEmSourceToken_

private

Definition at line 234 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::InputTag L1TGCT::gctTauJetsSource_

private

Definition at line 222 of file L1TGCT.h.

Referenced by analyze().

edm::EDGetTokenT<L1GctJetCandCollection> L1TGCT::gctTauJetsSourceToken_

private

Definition at line 237 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::EDGetTokenT<L1GctEtHadCollection> L1TGCT::l1EtHadToken_

private

Definition at line 243 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::EDGetTokenT<L1GctEtMissCollection> L1TGCT::l1EtMissToken_

private

Definition at line 241 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::EDGetTokenT<L1GctEtTotalCollection> L1TGCT::l1EtTotalToken_

private

Definition at line 244 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

MonitorElement* L1TGCT::l1GctAllEmOccRankBx_

private

Definition at line 173 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctAllJetsEtEtaPhi_

private

Definition at line 149 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctAllJetsOccEtaPhi_

private

Definition at line 157 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctAllJetsOccRankBx_

private

Definition at line 172 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctCenJetsEtEtaPhi_

private

Definition at line 150 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctCenJetsOccEtaPhi_

private

Definition at line 158 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctCenJetsRank_

private

Definition at line 165 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtHad_

private

Definition at line 189 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtHadOccBx_

private

Definition at line 191 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtHadOf_

private

Definition at line 190 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtMiss_

private

Definition at line 176 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtMissHtMissCorr_

private

Definition at line 184 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtMissHtMissCorrPhi_

private

Definition at line 185 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtMissOccBx_

private

Definition at line 179 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtMissOf_

private

Definition at line 178 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtMissPhi_

private

Definition at line 177 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtTotal_

private

Definition at line 186 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtTotalEtHadCorr_

private

Definition at line 192 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtTotalOccBx_

private

Definition at line 188 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctEtTotalOf_

private

Definition at line 187 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctForJetsEtEtaPhi_

private

Definition at line 151 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctForJetsOccEtaPhi_

private

Definition at line 159 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctForJetsRank_

private

Definition at line 166 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing1ETSumNegEta_

private

Definition at line 204 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing1ETSumPosEta_

private

Definition at line 203 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing1PosEtaNegEta_

private

Definition at line 195 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing1TowerCountNegEta_

private

Definition at line 200 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing1TowerCountPosEta_

private

Definition at line 199 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing1TowerCountPosEtaNegEta_

private

Definition at line 197 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing2ETSumNegEta_

private

Definition at line 206 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing2ETSumPosEta_

private

Definition at line 205 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing2PosEtaNegEta_

private

Definition at line 196 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing2TowerCountNegEta_

private

Definition at line 202 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing2TowerCountPosEta_

private

Definition at line 201 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRing2TowerCountPosEtaNegEta_

private

Definition at line 198 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRingETSumOccBx_

private

Definition at line 209 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRingRatioNegEta_

private

Definition at line 208 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRingRatioPosEta_

private

Definition at line 207 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHFRingTowerCountOccBx_

private

Definition at line 210 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHtMiss_

private

Definition at line 180 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHtMissOccBx_

private

Definition at line 183 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHtMissOf_

private

Definition at line 182 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctHtMissPhi_

private

Definition at line 181 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctIsoEmOccEtaPhi_

private

Definition at line 162 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctIsoEmRank_

private

Definition at line 169 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctIsoEmRankEtaPhi_

private

Definition at line 154 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctIsoTauJetsEtEtaPhi_

private

Definition at line 153 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctIsoTauJetsOccEtaPhi_

private

Definition at line 161 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctIsoTauJetsRank_

private

Definition at line 168 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctNonIsoEmOccEtaPhi_

private

Definition at line 163 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctNonIsoEmRank_

private

Definition at line 170 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctNonIsoEmRankEtaPhi_

private

Definition at line 155 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctTauJetsEtEtaPhi_

private

Definition at line 152 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctTauJetsOccEtaPhi_

private

Definition at line 160 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* L1TGCT::l1GctTauJetsRank_

private

Definition at line 167 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<L1GctHFBitCountsCollection> L1TGCT::l1HFCountsToken_

private

Definition at line 240 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

edm::EDGetTokenT<L1GctHtMissCollection> L1TGCT::l1HtMissToken_

private

Definition at line 242 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().

std::ofstream L1TGCT::logFile_

private

Definition at line 218 of file L1TGCT.h.

MonitorElement* L1TGCT::lumisecId_

private

Definition at line 212 of file L1TGCT.h.

Referenced by bookHistograms().

bool L1TGCT::m_stage1_layer2_

private

filter TriggerType

Definition at line 229 of file L1TGCT.h.

Referenced by analyze(), bookHistograms(), and L1TGCT().

bool L1TGCT::monitorDaemon_

private

Definition at line 217 of file L1TGCT.h.

int L1TGCT::nev_

private

Definition at line 214 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

std::string L1TGCT::outputFile_

private

Definition at line 215 of file L1TGCT.h.

Referenced by L1TGCT().

MonitorElement* L1TGCT::runId_

private

Definition at line 211 of file L1TGCT.h.

Referenced by bookHistograms().

MonitorElement* L1TGCT::triggerType_

private

Definition at line 146 of file L1TGCT.h.

Referenced by analyze(), and bookHistograms().

bool L1TGCT::verbose_

private

Definition at line 216 of file L1TGCT.h.

Referenced by analyze(), and L1TGCT().