CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Private Member Functions | Static Private Member Functions | Private Attributes
GctErrorAnalyzer Class Reference

#include <L1Trigger/L1GctAnalyzer/src/GctErrorAnalyzer.cc>

Inheritance diagram for GctErrorAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 GctErrorAnalyzer ()=delete
 
 GctErrorAnalyzer (const GctErrorAnalyzer &)=delete
 
 GctErrorAnalyzer (const edm::ParameterSet &)
 
GctErrorAnalyzer operator= (const GctErrorAnalyzer &)=delete
 
 ~GctErrorAnalyzer () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex >
const & 
esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void beginJob () override
 
template<class T >
bool checkCollections (const T &collection, const unsigned int &constraint, const std::string &label)
 
void endJob () override
 
void plotCenJets (const edm::Handle< L1GctJetCandCollection > &data, const edm::Handle< L1GctJetCandCollection > &emu)
 
void plotEGErrors (const edm::Handle< L1GctEmCandCollection > &dataiso, const edm::Handle< L1GctEmCandCollection > &emuiso, const edm::Handle< L1GctEmCandCollection > &datanoniso, const edm::Handle< L1GctEmCandCollection > &emunoniso, const edm::Handle< L1CaloEmCollection > &regions)
 
void plotForJets (const edm::Handle< L1GctJetCandCollection > &data, const edm::Handle< L1GctJetCandCollection > &emu)
 
void plotHFBitCounts (const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsD, const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsE)
 
void plotHFErrors (const edm::Handle< L1GctHFRingEtSumsCollection > &hfRingSumsD, const edm::Handle< L1GctHFRingEtSumsCollection > &hfRingSumsE, const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsD, const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsE, const edm::Handle< L1CaloRegionCollection > &caloRegions)
 
void plotHFRingSums (const edm::Handle< L1GctHFRingEtSumsCollection > &data, const edm::Handle< L1GctHFRingEtSumsCollection > &emu)
 
void plotIntJets (const edm::Handle< L1GctInternJetDataCollection > &emu)
 
void plotIsoEm (const edm::Handle< L1GctEmCandCollection > &data, const edm::Handle< L1GctEmCandCollection > &emu)
 
void plotJetErrors (const edm::Handle< L1GctJetCandCollection > &cendata, const edm::Handle< L1GctJetCandCollection > &cenemu, const edm::Handle< L1GctJetCandCollection > &taudata, const edm::Handle< L1GctJetCandCollection > &tauemu, const edm::Handle< L1GctJetCandCollection > &fordata, const edm::Handle< L1GctJetCandCollection > &foremu, const edm::Handle< L1CaloRegionCollection > &regions)
 
void plotMissingEErrors (const edm::Handle< L1GctEtMissCollection > &missingEtD, const edm::Handle< L1GctEtMissCollection > &missingEtE, const edm::Handle< L1GctHtMissCollection > &missingHtD, const edm::Handle< L1GctHtMissCollection > &missingHtE, edm::Handle< L1CaloRegionCollection > &caloRegions, const edm::Handle< L1GctInternJetDataCollection > &intjetsemu, const edm::Handle< L1GctInternHtMissCollection > intMissingHtD)
 
void plotMissingEt (const edm::Handle< L1GctEtMissCollection > &missingEtD, const edm::Handle< L1GctEtMissCollection > &missingEtE)
 
void plotMissingHt (const edm::Handle< L1GctHtMissCollection > &missingHtD, const edm::Handle< L1GctHtMissCollection > &missingHtE)
 
void plotNonIsoEm (const edm::Handle< L1GctEmCandCollection > &data, const edm::Handle< L1GctEmCandCollection > &emu)
 
void plotRCTRegions (const edm::Handle< L1CaloRegionCollection > &caloRegions)
 
void plotTauJets (const edm::Handle< L1GctJetCandCollection > &data, const edm::Handle< L1GctJetCandCollection > &emu)
 
void plotTotalE (const edm::Handle< L1GctEtTotalCollection > &totalEtD, const edm::Handle< L1GctEtTotalCollection > &totalEtE)
 
void plotTotalEErrors (const edm::Handle< L1GctEtTotalCollection > &totalEtD, const edm::Handle< L1GctEtTotalCollection > &totalEtE, const edm::Handle< L1GctEtHadCollection > &totalHtD, const edm::Handle< L1GctEtHadCollection > &totalHtE, const edm::Handle< L1CaloRegionCollection > &caloRegions)
 
void plotTotalH (const edm::Handle< L1GctEtHadCollection > &totalHtD, const edm::Handle< L1GctEtHadCollection > &totalHtE)
 

Static Private Member Functions

static bool sortJets (const jetData &jet1, const jetData &jet2)
 

Private Attributes

TH1I * cenJet_errorFlag_
 
TH2I * cenJetD_EtEtaPhi_
 
TH2I * cenJetD_GlobalError_EtEtaPhi_
 
TH1I * cenJetD_GlobalError_Rank_
 
TH2I * cenJetD_OccEtaPhi_
 
TH1I * cenJetD_Rank_
 
TH2I * cenJetE_EtEtaPhi_
 
TH2I * cenJetE_GlobalError_EtEtaPhi_
 
TH1I * cenJetE_GlobalError_Rank_
 
TH2I * cenJetE_OccEtaPhi_
 
TH1I * cenJetE_Rank_
 
edm::InputTag dataTag_
 
bool doBitCountDebug_
 
bool doCenJetsDebug_
 
bool doEg_
 
bool doEmuMBx_
 
bool doExtraMissingHTDebug_
 
bool doForJetsDebug_
 
bool doGCTMBx_
 
bool doHF_
 
bool doIsoDebug_
 
bool doJets_
 
bool doMissingEnergySums_
 
bool doMissingETDebug_
 
bool doMissingHTDebug_
 
bool doNonIsoDebug_
 
bool doRCT_
 
bool doRCTMBx_
 
bool doRingSumDebug_
 
bool doTauJetsDebug_
 
bool doTotalEnergySums_
 
bool doTotalEtDebug_
 
bool doTotalHtDebug_
 
edm::InputTag emuTag_
 
int EmuTrigBx_
 
std::vector< TFileDirectoryerrorHistCat
 
unsigned int eventNumber
 
TH1I * forJet_errorFlag_
 
TH2I * forJetD_EtEtaPhi_
 
TH2I * forJetD_GlobalError_EtEtaPhi_
 
TH1I * forJetD_GlobalError_Rank_
 
TH2I * forJetD_OccEtaPhi_
 
TH1I * forJetD_Rank_
 
TH2I * forJetE_EtEtaPhi_
 
TH2I * forJetE_GlobalError_EtEtaPhi_
 
TH1I * forJetE_GlobalError_Rank_
 
TH2I * forJetE_OccEtaPhi_
 
TH1I * forJetE_Rank_
 
int GCTTrigBx_
 
TH1I * hfBitCount_errorFlag_
 
TH1I * hfBitCountD_1neg_
 
TH1I * hfBitCountD_1pos_
 
TH1I * hfBitCountD_2neg_
 
TH1I * hfBitCountD_2pos_
 
TH1I * hfBitCountE_1neg_
 
TH1I * hfBitCountE_1pos_
 
TH1I * hfBitCountE_2neg_
 
TH1I * hfBitCountE_2pos_
 
TH1I * hfRingSum_errorFlag_
 
TH1I * hfRingSumD_1neg_
 
TH1I * hfRingSumD_1pos_
 
TH1I * hfRingSumD_2neg_
 
TH1I * hfRingSumD_2pos_
 
TH1I * hfRingSumE_1neg_
 
TH1I * hfRingSumE_1pos_
 
TH1I * hfRingSumE_2neg_
 
TH1I * hfRingSumE_2pos_
 
TH1I * intJetE_Et_
 
TH1I * intJetE_Jet1Et_
 
TH1I * intJetE_Jet2Et_
 
TH1I * intJetE_Jet3Et_
 
TH1I * intJetE_Jet4Et_
 
TH1I * intJetE_Of_
 
TH2I * intJetEtEtaPhiE_
 
bool isBitCountError
 
bool isCenJetError
 
bool isForJetError
 
bool isIsoError
 
bool isMissingEError
 
bool isMissingHError
 
bool isNonIsoError
 
TH1I * isoEg_errorFlag_
 
TH2I * isoEgD_EtEtaPhi_
 
TH2I * isoEgD_GlobalError_EtEtaPhi_
 
TH1I * isoEgD_GlobalError_Rank_
 
TH2I * isoEgD_OccEtaPhi_
 
TH1I * isoEgD_Rank_
 
TH2I * isoEgE_EtEtaPhi_
 
TH2I * isoEgE_GlobalError_EtEtaPhi_
 
TH1I * isoEgE_GlobalError_Rank_
 
TH2I * isoEgE_OccEtaPhi_
 
TH1I * isoEgE_Rank_
 
bool isRingSumError
 
bool isTauJetError
 
bool isTotalEError
 
bool isTotalHError
 
GctErrorAnalyzerMBxInfo MBxInfo
 
TH1I * missingEt_errorFlag_
 
TH1I * missingEtD_
 
TH1I * missingEtD_Of_
 
TH1I * missingEtD_Phi_
 
TH1I * missingEtE_
 
TH1I * missingEtE_Of_
 
TH1I * missingEtE_Phi_
 
TH1I * missingHt_errorFlag_
 
TH1I * missingHtD_
 
TH1I * missingHtD_HtXNegLeaf1
 
TH1I * missingHtD_HtXNegLeaf2
 
TH1I * missingHtD_HtXNegLeaf3
 
TH1I * missingHtD_HtXPosLeaf1
 
TH1I * missingHtD_HtXPosLeaf2
 
TH1I * missingHtD_HtXPosLeaf3
 
TH1I * missingHtD_HtYNegLeaf1
 
TH1I * missingHtD_HtYNegLeaf2
 
TH1I * missingHtD_HtYNegLeaf3
 
TH1I * missingHtD_HtYPosLeaf1
 
TH1I * missingHtD_HtYPosLeaf2
 
TH1I * missingHtD_HtYPosLeaf3
 
TH1I * missingHtD_Of_
 
TH1I * missingHtD_Phi_
 
TH1I * missingHtE_
 
TH1I * missingHtE_Of_
 
TH1I * missingHtE_Phi_
 
TH1I * nonIsoEg_errorFlag_
 
TH2I * nonIsoEgD_EtEtaPhi_
 
TH2I * nonIsoEgD_GlobalError_EtEtaPhi_
 
TH1I * nonIsoEgD_GlobalError_Rank_
 
TH2I * nonIsoEgD_OccEtaPhi_
 
TH1I * nonIsoEgD_Rank_
 
TH2I * nonIsoEgE_EtEtaPhi_
 
TH2I * nonIsoEgE_GlobalError_EtEtaPhi_
 
TH1I * nonIsoEgE_GlobalError_Rank_
 
TH2I * nonIsoEgE_OccEtaPhi_
 
TH1I * nonIsoEgE_Rank_
 
TH2I * RCT_EtEtaPhi_
 
TH2I * RCT_FgEtaPhi_
 
TH2I * RCT_OfEtaPhi_
 
const unsigned int * RCT_REGION_QUANTA
 
TH2I * RCT_TvEtaPhi_
 
int RCTTrigBx_
 
TH1I * tauJet_errorFlag_
 
TH2I * tauJetD_EtEtaPhi_
 
TH2I * tauJetD_GlobalError_EtEtaPhi_
 
TH1I * tauJetD_GlobalError_Rank_
 
TH2I * tauJetD_OccEtaPhi_
 
TH1I * tauJetD_Rank_
 
TH2I * tauJetE_EtEtaPhi_
 
TH2I * tauJetE_GlobalError_EtEtaPhi_
 
TH1I * tauJetE_GlobalError_Rank_
 
TH2I * tauJetE_OccEtaPhi_
 
TH1I * tauJetE_Rank_
 
TH1I * totalEt_errorFlag_
 
TH1I * totalEtD_
 
TH1I * totalEtD_Of_
 
TH1I * totalEtE_
 
TH1I * totalEtE_Of_
 
TH1I * totalHt_errorFlag_
 
TH1I * totalHtD_
 
TH1I * totalHtD_Of_
 
TH1I * totalHtE_
 
TH1I * totalHtE_Of_
 
std::string useSys_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsInputProcessBlocks ()
 
static bool wantsProcessBlocks ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< Bconsumes (edm::InputTag tag) noexcept
 
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<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
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)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Description: Tool to debug the GCT with useful output

Implementation: <Notes on="" implementation>="">

Definition at line 57 of file GctErrorAnalyzer.cc.

Constructor & Destructor Documentation

GctErrorAnalyzer::GctErrorAnalyzer ( )
delete
GctErrorAnalyzer::GctErrorAnalyzer ( const GctErrorAnalyzer )
delete
GctErrorAnalyzer::GctErrorAnalyzer ( const edm::ParameterSet iConfig)
explicit

Definition at line 303 of file GctErrorAnalyzer.cc.

References cenJet_errorFlag_, cenJetD_EtEtaPhi_, cenJetD_GlobalError_EtEtaPhi_, cenJetD_GlobalError_Rank_, cenJetD_OccEtaPhi_, cenJetD_Rank_, cenJetE_EtEtaPhi_, cenJetE_GlobalError_EtEtaPhi_, cenJetE_GlobalError_Rank_, cenJetE_OccEtaPhi_, cenJetE_Rank_, GctErrorAnalyzerMBxInfo::EmuTrigBx, EmuTrigBx_, errorHistCat, forJet_errorFlag_, forJetD_EtEtaPhi_, forJetD_GlobalError_EtEtaPhi_, forJetD_GlobalError_Rank_, forJetD_OccEtaPhi_, forJetD_Rank_, forJetE_EtEtaPhi_, forJetE_GlobalError_EtEtaPhi_, forJetE_GlobalError_Rank_, forJetE_OccEtaPhi_, forJetE_Rank_, compareTotals::fs, GctErrorAnalyzerMBxInfo::GCTTrigBx, GCTTrigBx_, hfBitCount_errorFlag_, hfBitCountD_1neg_, hfBitCountD_1pos_, hfBitCountD_2neg_, hfBitCountD_2pos_, hfBitCountE_1neg_, hfBitCountE_1pos_, hfBitCountE_2neg_, hfBitCountE_2pos_, hfRingSum_errorFlag_, hfRingSumD_1neg_, hfRingSumD_1pos_, hfRingSumD_2neg_, hfRingSumD_2pos_, hfRingSumE_1neg_, hfRingSumE_1pos_, hfRingSumE_2neg_, hfRingSumE_2pos_, mps_fire::i, intJetE_Et_, intJetE_Jet1Et_, intJetE_Jet2Et_, intJetE_Jet3Et_, intJetE_Jet4Et_, intJetE_Of_, intJetEtEtaPhiE_, isBitCountError, isCenJetError, isForJetError, isIsoError, isMissingEError, isMissingHError, isNonIsoError, isoEg_errorFlag_, isoEgD_EtEtaPhi_, isoEgD_GlobalError_EtEtaPhi_, isoEgD_GlobalError_Rank_, isoEgD_OccEtaPhi_, isoEgD_Rank_, isoEgE_EtEtaPhi_, isoEgE_GlobalError_EtEtaPhi_, isoEgE_GlobalError_Rank_, isoEgE_OccEtaPhi_, isoEgE_Rank_, isRingSumError, isTauJetError, isTotalEError, isTotalHError, TFileDirectory::make(), MBxInfo, missingEt_errorFlag_, missingEtD_, missingEtD_Of_, missingEtD_Phi_, missingEtE_, missingEtE_Of_, missingEtE_Phi_, missingHt_errorFlag_, missingHtD_, missingHtD_HtXNegLeaf1, missingHtD_HtXNegLeaf2, missingHtD_HtXNegLeaf3, missingHtD_HtXPosLeaf1, missingHtD_HtXPosLeaf2, missingHtD_HtXPosLeaf3, missingHtD_HtYNegLeaf1, missingHtD_HtYNegLeaf2, missingHtD_HtYNegLeaf3, missingHtD_HtYPosLeaf1, missingHtD_HtYPosLeaf2, missingHtD_HtYPosLeaf3, missingHtD_Of_, missingHtD_Phi_, missingHtE_, missingHtE_Of_, missingHtE_Phi_, TFileDirectory::mkdir(), TFileService::mkdir(), nonIsoEg_errorFlag_, nonIsoEgD_EtEtaPhi_, nonIsoEgD_GlobalError_EtEtaPhi_, nonIsoEgD_GlobalError_Rank_, nonIsoEgD_OccEtaPhi_, nonIsoEgD_Rank_, nonIsoEgE_EtEtaPhi_, nonIsoEgE_GlobalError_EtEtaPhi_, nonIsoEgE_GlobalError_Rank_, nonIsoEgE_OccEtaPhi_, nonIsoEgE_Rank_, RCT_EtEtaPhi_, RCT_FgEtaPhi_, RCT_OfEtaPhi_, RCT_REGION_QUANTA, RCT_REGION_QUANTA_LAB, RCT_REGION_QUANTA_P5, RCT_TvEtaPhi_, GctErrorAnalyzerMBxInfo::RCTTrigBx, RCTTrigBx_, tauJet_errorFlag_, tauJetD_EtEtaPhi_, tauJetD_GlobalError_EtEtaPhi_, tauJetD_GlobalError_Rank_, tauJetD_OccEtaPhi_, tauJetD_Rank_, tauJetE_EtEtaPhi_, tauJetE_GlobalError_EtEtaPhi_, tauJetE_GlobalError_Rank_, tauJetE_OccEtaPhi_, tauJetE_Rank_, totalEt_errorFlag_, totalEtD_, totalEtD_Of_, totalEtE_, totalEtE_Of_, totalHt_errorFlag_, totalHtD_, totalHtD_Of_, totalHtE_, totalHtE_Of_, and useSys_.

304  : doRCT_(iConfig.getUntrackedParameter<bool>("doRCT", true)),
305  doEg_(iConfig.getUntrackedParameter<bool>("doEg", true)),
306  doIsoDebug_(iConfig.getUntrackedParameter<bool>("doIsoDebug", true)),
307  doNonIsoDebug_(iConfig.getUntrackedParameter<bool>("doNonIsoDebug", true)),
308  doJets_(iConfig.getUntrackedParameter<bool>("doJets", true)),
309  doCenJetsDebug_(iConfig.getUntrackedParameter<bool>("doCenJetsDebug", true)),
310  doTauJetsDebug_(iConfig.getUntrackedParameter<bool>("doTauJetsDebug", true)),
311  doForJetsDebug_(iConfig.getUntrackedParameter<bool>("doForJetsDebug", true)),
312  doHF_(iConfig.getUntrackedParameter<bool>("doHF", true)),
313  doRingSumDebug_(iConfig.getUntrackedParameter<bool>("doRingSumDebug", true)),
314  doBitCountDebug_(iConfig.getUntrackedParameter<bool>("doBitCountDebug", true)),
315  doTotalEnergySums_(iConfig.getUntrackedParameter<bool>("doTotalEnergySums", true)),
316  doTotalEtDebug_(iConfig.getUntrackedParameter<bool>("doTotalEtDebug", true)),
317  doTotalHtDebug_(iConfig.getUntrackedParameter<bool>("doTotalHtDebug", true)),
318  doMissingEnergySums_(iConfig.getUntrackedParameter<bool>("doMissingEnergySums", true)),
319  doMissingETDebug_(iConfig.getUntrackedParameter<bool>("doMissingETDebug", true)),
320  doMissingHTDebug_(iConfig.getUntrackedParameter<bool>("doMissingHTDebug", true)),
321  doExtraMissingHTDebug_(iConfig.getUntrackedParameter<bool>("doExtraMissingHTDebug", false)),
322  doRCTMBx_(iConfig.getUntrackedParameter<bool>("doRCTMBx", false)),
323  doEmuMBx_(iConfig.getUntrackedParameter<bool>("doEmuMBx", false)),
324  doGCTMBx_(iConfig.getUntrackedParameter<bool>("doGCTMBx", false)),
325  RCTTrigBx_(iConfig.getUntrackedParameter<int>("RCTTrigBx", 0)),
326  EmuTrigBx_(iConfig.getUntrackedParameter<int>("EmuTrigBx", 0)),
327  GCTTrigBx_(iConfig.getUntrackedParameter<int>("GCTTrigBx", 0)),
328  dataTag_(iConfig.getUntrackedParameter<edm::InputTag>("dataTag", edm::InputTag("gctDigis"))),
329  emuTag_(iConfig.getUntrackedParameter<edm::InputTag>("emuTag", edm::InputTag("gctEmuDigis"))),
330  useSys_(iConfig.getUntrackedParameter<std::string>("useSys", "P5")) {
331  //now do what ever initialization is needed
332  //make the root file
334 
335  //to try to make this look more elegant
336  //make a string for each folder we'd like for the Data and Emulator Histograms
337  std::vector<std::string> quantities;
338  quantities.push_back("IsoEm");
339  quantities.push_back("NonIsoEM");
340  quantities.push_back("CenJets");
341  quantities.push_back("TauJets");
342  quantities.push_back("ForJets");
343  quantities.push_back("HFRingSums");
344  quantities.push_back("HFBitCounts");
345  quantities.push_back("TotalESums");
346  quantities.push_back("MissingESums");
347 
348  //make the Emulator Histogram directory
349  TFileDirectory emuHist = fs->mkdir("EmulatorHistograms");
350  std::vector<TFileDirectory> emuHistCat;
351 
352  //make the Data Histogram directory
353  TFileDirectory dataHist = fs->mkdir("DataHistograms");
354  std::vector<TFileDirectory> dataHistCat;
355 
356  //make the ErrorFlags directory
357  TFileDirectory errorHistFlags = fs->mkdir("ErrorHistograms_Flags");
358 
359  //make the ErrorDebug directory
360  TFileDirectory errorHistDetails = fs->mkdir("ErrorHistograms_Details");
361 
362  for (unsigned int i = 0; i < quantities.size(); i++) {
363  //fill the data and emulator folders with the directories
364  emuHistCat.push_back(emuHist.mkdir(quantities.at(i)));
365  dataHistCat.push_back(dataHist.mkdir(quantities.at(i)));
366  }
367 
368  //add a folder for RCT Regions - which only exist in data
369  dataHistCat.push_back(dataHist.mkdir("RCTRegions"));
370  //and add a folder for the Intermediate Jets - which only exist in emulator
371  emuHistCat.push_back(emuHist.mkdir("IntJets"));
372 
373  //Fill the ErrorDebug folder with the directories
374  errorHistCat.push_back(errorHistDetails.mkdir("EM"));
375  errorHistCat.push_back(errorHistDetails.mkdir("Jets"));
376  errorHistCat.push_back(errorHistDetails.mkdir("HF"));
377  errorHistCat.push_back(errorHistDetails.mkdir("TotalE"));
378  errorHistCat.push_back(errorHistDetails.mkdir("MissingE"));
379 
380  //BOOK HISTOGRAMS
381  RCT_EtEtaPhi_ = dataHistCat.at(9).make<TH2I>(
382  "RCT_EtEtaPhi", "RCT_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
383  RCT_TvEtaPhi_ = dataHistCat.at(9).make<TH2I>(
384  "RCT_TvEtaPhi", "RCT_TvEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
385  RCT_FgEtaPhi_ = dataHistCat.at(9).make<TH2I>(
386  "RCT_FgEtaPhi", "RCT_FgEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
387  RCT_OfEtaPhi_ = dataHistCat.at(9).make<TH2I>(
388  "RCT_OfEtEtaPhi", "RCT_OfEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
389  //isoEg
390  isoEgD_Rank_ = dataHistCat.at(0).make<TH1I>("isoEgD_Rank", "isoEgD_Rank;Rank;Number of Events", 64, -0.5, 63.5);
391  isoEgE_Rank_ = emuHistCat.at(0).make<TH1I>("isoEgE_Rank", "isoEgE_Rank;Rank;Number of Events", 64, -0.5, 63.5);
392  isoEgD_EtEtaPhi_ = dataHistCat.at(0).make<TH2I>(
393  "isoEgD_EtEtaPhi", "isoEgD_EtEtaPhi;#eta (GCT Units);#phi(GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
394  isoEgE_EtEtaPhi_ = emuHistCat.at(0).make<TH2I>(
395  "isoEgE_EtEtaPhi", "isoEgE_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
396  isoEgD_OccEtaPhi_ = dataHistCat.at(0).make<TH2I>(
397  "isoEgD_OccEtaPhi", "isoEgD_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
398  isoEgE_OccEtaPhi_ = emuHistCat.at(0).make<TH2I>(
399  "isoEgE_OccEtaPhi", "isoEgE_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
401  errorHistFlags.make<TH1I>("isoEg_errorFlag", "isoEg_errorFlag;Status;Number of Candidates", 3, -0.5, 2.5);
402  //Global isoEg Error
403  isoEgD_GlobalError_Rank_ = errorHistCat.at(0).make<TH1I>(
404  "isoEgD_GlobalError_Rank", "isoEgD_GlobalError_Rank;Rank;Number of Events", 64, -0.5, 63.5);
405  isoEgE_GlobalError_Rank_ = errorHistCat.at(0).make<TH1I>(
406  "isoEgE_GlobalError_Rank", "isoEgE_GlobalError_Rank;Rank;Number of Events", 64, -0.5, 63.5);
407  isoEgD_GlobalError_EtEtaPhi_ = errorHistCat.at(0).make<TH2I>(
408  "isoEgD_GlobalError_EtEtaPhi", "isoEgD_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
409  isoEgE_GlobalError_EtEtaPhi_ = errorHistCat.at(0).make<TH2I>(
410  "isoEgE_GlobalError_EtEtaPhi", "isoEgE_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
411  //nonIsoEg
413  dataHistCat.at(1).make<TH1I>("nonIsoEgD_Rank", "nonIsoEgD_Rank;Rank;Number of Events", 64, -0.5, 63.5);
415  emuHistCat.at(1).make<TH1I>("nonIsoEgE_Rank", "nonIsoEgE_Rank;Rank;Number of Events", 64, -0.5, 63.5);
416  nonIsoEgD_EtEtaPhi_ = dataHistCat.at(1).make<TH2I>(
417  "nonIsoEgD_EtEtaPhi", "nonIsoEgD_EtEtaPhi;#eta (GCT Units);#phi(GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
418  nonIsoEgE_EtEtaPhi_ = emuHistCat.at(1).make<TH2I>(
419  "nonIsoEgE_EtEtaPhi", "nonIsoEgE_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
420  nonIsoEgD_OccEtaPhi_ = dataHistCat.at(1).make<TH2I>(
421  "nonIsoEgD_OccEtaPhi", "nonIsoEgD_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
422  nonIsoEgE_OccEtaPhi_ = emuHistCat.at(1).make<TH2I>(
423  "nonIsoEgE_OccEtaPhi", "nonIsoEgE_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
425  errorHistFlags.make<TH1I>("nonIsoEg_errorFlag", "nonIsoEg_errorFlag;Status;Number of Candidates", 3, -0.5, 2.5);
426  //Global nonIsoEg Error
427  nonIsoEgD_GlobalError_Rank_ = errorHistCat.at(0).make<TH1I>(
428  "nonIsoEgD_GlobalError_Rank", "nonIsoEgD_GlobalError_Rank;Rank;Number of Events", 64, -0.5, 63.5);
429  nonIsoEgE_GlobalError_Rank_ = errorHistCat.at(0).make<TH1I>(
430  "nonIsoEgE_GlobalError_Rank", "nonIsoEgE_GlobalError_Rank;Rank;Number of Events", 64, -0.5, 63.5);
432  "nonIsoEgD_GlobalError_EtEtaPhi", "nonIsoEgD_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
434  "nonIsoEgE_GlobalError_EtEtaPhi", "nonIsoEgE_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
435  //CenJets
436  cenJetD_Rank_ = dataHistCat.at(2).make<TH1I>("cenJetD_Rank", "cenJetD_Rank;Rank;Number of Events", 64, -0.5, 63.5);
437  cenJetE_Rank_ = emuHistCat.at(2).make<TH1I>("cenJetE_Rank", "cenJetE_Rank;Rank;Number of Events", 64, -0.5, 63.5);
438  cenJetD_EtEtaPhi_ = dataHistCat.at(2).make<TH2I>(
439  "cenJetD_EtEtaPhi", "cenJetD_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
440  cenJetE_EtEtaPhi_ = emuHistCat.at(2).make<TH2I>(
441  "cenJetE_EtEtaPhi", "cenJetE_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
442  cenJetD_OccEtaPhi_ = dataHistCat.at(2).make<TH2I>(
443  "cenJetD_OccEtaPhi", "cenJetD_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
444  cenJetE_OccEtaPhi_ = emuHistCat.at(2).make<TH2I>(
445  "cenJetE_OccEtaPhi", "cenJetE_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
447  errorHistFlags.make<TH1I>("cenJet_errorFlag", "cenJet_errorFlag;Status;Number of Candidates", 3, -0.5, 2.5);
448  //Global CenJet Error
450  errorHistCat.at(1).make<TH1I>("cenJetD_GlobalError_Rank", "cenJetD_GlobalError_Rank", 64, -0.5, 63.5);
452  errorHistCat.at(1).make<TH1I>("cenJetE_GlobalError_Rank", "cenJetE_GlobalError_Rank", 64, -0.5, 63.5);
453  cenJetD_GlobalError_EtEtaPhi_ = errorHistCat.at(1).make<TH2I>(
454  "cenJetD_GlobalError_EtEtaPhi", "cenJetD_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
455  cenJetE_GlobalError_EtEtaPhi_ = errorHistCat.at(1).make<TH2I>(
456  "cenJetE_GlobalError_EtEtaPhi", "cenJetE_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
457  //TauJets
458  tauJetD_Rank_ = dataHistCat.at(3).make<TH1I>("tauJetD_Rank", "tauJetD_Rank;Rank;Number of Events", 64, -0.5, 63.5);
459  tauJetE_Rank_ = emuHistCat.at(3).make<TH1I>("tauJetE_Rank", "tauJetE_Rank;Rank;Number of Events", 64, -0.5, 63.5);
460  tauJetD_EtEtaPhi_ = dataHistCat.at(3).make<TH2I>(
461  "tauJetD_EtEtaPhi", "tauJetD_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
462  tauJetE_EtEtaPhi_ = emuHistCat.at(3).make<TH2I>(
463  "tauJetE_EtEtaPhi", "tauJetE_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
464  tauJetD_OccEtaPhi_ = dataHistCat.at(3).make<TH2I>(
465  "tauJetD_OccEtaPhi", "tauJetD_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
466  tauJetE_OccEtaPhi_ = emuHistCat.at(3).make<TH2I>(
467  "tauJetE_OccEtaPhi", "tauJetE_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
469  errorHistFlags.make<TH1I>("tauJet_errorFlag", "tauJet_errorFlag;Status;Number of Candidates", 3, -0.5, 2.5);
470  //Global TauJet Error
472  errorHistCat.at(1).make<TH1I>("tauJetD_GlobalError_Rank", "tauJetD_GlobalError_Rank", 64, -0.5, 63.5);
474  errorHistCat.at(1).make<TH1I>("tauJetE_GlobalError_Rank", "tauJetE_GlobalError_Rank", 64, -0.5, 63.5);
475  tauJetD_GlobalError_EtEtaPhi_ = errorHistCat.at(1).make<TH2I>(
476  "tauJetD_GlobalError_EtEtaPhi", "tauJetD_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
477  tauJetE_GlobalError_EtEtaPhi_ = errorHistCat.at(1).make<TH2I>(
478  "tauJetE_GlobalError_EtEtaPhi", "tauJetE_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
479  //ForJets
480  forJetD_Rank_ = dataHistCat.at(4).make<TH1I>("forJetD_Rank", "forJetD_Rank;Rank;Number of Events", 64, -0.5, 63.5);
481  forJetE_Rank_ = emuHistCat.at(4).make<TH1I>("forJetE_Rank", "forJetE_Rank;Rank;Number of Events", 64, -0.5, 63.5);
482  forJetD_EtEtaPhi_ = dataHistCat.at(4).make<TH2I>(
483  "forJetD_EtEtaPhi", "forJetD_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
484  forJetE_EtEtaPhi_ = emuHistCat.at(4).make<TH2I>(
485  "forJetE_EtEtaPhi", "forJetE_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
486  forJetD_OccEtaPhi_ = dataHistCat.at(4).make<TH2I>(
487  "forJetD_OccEtaPhi", "forJetD_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
488  forJetE_OccEtaPhi_ = emuHistCat.at(4).make<TH2I>(
489  "forJetE_OccEtaPhi", "forJetE_OccEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
491  errorHistFlags.make<TH1I>("forJet_errorFlag", "forJet_errorFlag;Status;Number of Candidates", 3, -0.5, 2.5);
492  //Global ForJet Error
494  errorHistCat.at(1).make<TH1I>("forJetD_GlobalError_Rank", "forJetD_GlobalError_Rank", 64, -0.5, 63.5);
496  errorHistCat.at(1).make<TH1I>("forJetE_GlobalError_Rank", "forJetE_GlobalError_Rank", 64, -0.5, 63.5);
497  forJetD_GlobalError_EtEtaPhi_ = errorHistCat.at(1).make<TH2I>(
498  "forJetD_GlobalError_EtEtaPhi", "forJetD_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
499  forJetE_GlobalError_EtEtaPhi_ = errorHistCat.at(1).make<TH2I>(
500  "forJetE_GlobalError_EtEtaPhi", "forJetE_GlobalError_EtEtaPhi", 22, -0.5, 21.5, 18, -0.5, 17.5);
501  //IntJets
502  intJetEtEtaPhiE_ = emuHistCat.at(9).make<TH2I>(
503  "intJetEtEtaPhiE_", "intJetEtEtaPhiE_;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
504  intJetE_Et_ = emuHistCat.at(9).make<TH1I>("intJetE_Et", "intJetE_Et;E_{T};Number of Events", 1024, -0.5, 1023.5);
505  intJetE_Of_ =
506  emuHistCat.at(9).make<TH1I>("intJetE_Of", "intJetE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
508  emuHistCat.at(9).make<TH1I>("intJetE_Jet1Et", "intJetE_Jet1Et;E_{T};Number of Events", 1024, -0.5, 1023.5);
510  emuHistCat.at(9).make<TH1I>("intJetE_Jet2Et", "intJetE_Jet2Et;E_{T};Number of Events", 1024, -0.5, 1023.5);
512  emuHistCat.at(9).make<TH1I>("intJetE_Jet3Et", "intJetE_Jet3Et;E_{T};Number of Events", 1024, -0.5, 1023.5);
514  emuHistCat.at(9).make<TH1I>("intJetE_Jet4Et", "intJetE_Jet4Et;E_{T};Number of Events", 1024, -0.5, 1023.5);
515  //HFRing Sums
516  hfRingSumD_1pos_ = dataHistCat.at(5).make<TH1I>("hfRingSumD_1+", "hfRingSumD_1+;Rank;Number of Events", 8, -0.5, 7.5);
517  hfRingSumD_1neg_ = dataHistCat.at(5).make<TH1I>("hfRingSumD_1-", "hfRingSumD_1-;Rank;Number of Events", 8, -0.5, 7.5);
518  hfRingSumD_2pos_ = dataHistCat.at(5).make<TH1I>("hfRingSumD_2+", "hfRingSumD_2+;Rank;Number of Events", 8, -0.5, 7.5);
519  hfRingSumD_2neg_ = dataHistCat.at(5).make<TH1I>("hfRingSumD_2-", "hfRingSumD_2-;Rank;Number of Events", 8, -0.5, 7.5);
520  hfRingSumE_1pos_ = emuHistCat.at(5).make<TH1I>("hfRingSumE_1+", "hfRingSumE_1+;Rank;Number of Events", 8, -0.5, 7.5);
521  hfRingSumE_1neg_ = emuHistCat.at(5).make<TH1I>("hfRingSumE_1-", "hfRingSumE_1-;Rank;Number of Events", 8, -0.5, 7.5);
522  hfRingSumE_2pos_ = emuHistCat.at(5).make<TH1I>("hfRingSumE_2+", "hfRingSumE_2+;Rank;Number of Events", 8, -0.5, 7.5);
523  hfRingSumE_2neg_ = emuHistCat.at(5).make<TH1I>("hfRingSumE_2-", "hfRingSumE_2-;Rank;Number of Events", 8, -0.5, 7.5);
525  errorHistFlags.make<TH1I>("hfRingSum_errorFlag", "hfRingSum_errorFlag;Status;Number of Candidates", 2, -0.5, 1.5);
526  //HFRing BitCounts
528  dataHistCat.at(6).make<TH1I>("hfBitCountD_1+", "hfBitCountD_1+;Rank;Number of Events", 8, -0.5, 7.5);
530  dataHistCat.at(6).make<TH1I>("hfBitCountD_1-", "hfBitCountD_1-;Rank;Number of Events", 8, -0.5, 7.5);
532  dataHistCat.at(6).make<TH1I>("hfBitCountD_2+", "hfBitCountD_2+;Rank;Number of Events", 8, -0.5, 7.5);
534  dataHistCat.at(6).make<TH1I>("hfBitCountD_2-", "hfBitCountD_2-;Rank;Number of Events", 8, -0.5, 7.5);
536  emuHistCat.at(6).make<TH1I>("hfBitCountE_1+", "hfBitCountE_1+;Rank;Number of Events", 8, -0.5, 7.5);
538  emuHistCat.at(6).make<TH1I>("hfBitCountE_1-", "hfBitCountE_1-;Rank;Number of Events", 8, -0.5, 7.5);
540  emuHistCat.at(6).make<TH1I>("hfBitCountE_2+", "hfBitCountE_2+;Rank;Number of Events", 8, -0.5, 7.5);
542  emuHistCat.at(6).make<TH1I>("hfBitCountE_2-", "hfBitCountE_2-;Rank;Number of Events", 8, -0.5, 7.5);
543  hfBitCount_errorFlag_ = errorHistFlags.make<TH1I>(
544  "hfBitCount_errorFlag", "hfBitCount_errorFlag;Status;Number of Candidates", 2, -0.5, 1.5);
545  //Total ET
546  totalEtD_ = dataHistCat.at(7).make<TH1I>("totalEtD", "totalEtD;E_{T};Number of Events", 2048, -0.5, 2047.5);
547  totalEtD_Of_ =
548  dataHistCat.at(7).make<TH1I>("totalEtD_Of", "totalEtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
549  totalEtE_ = emuHistCat.at(7).make<TH1I>("totalEtE", "totalEtE;E_{T};Number of Events", 2048, -0.5, 2047.5);
550  totalEtE_Of_ =
551  emuHistCat.at(7).make<TH1I>("totalEtE_Of", "totalEtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
553  errorHistFlags.make<TH1I>("totalEt_errorFlag", "totalEt_errorFlag;Status;Number of Candidates", 2, -0.5, 1.5);
554  //Book the Global ET Error histograms in the errorHistCat
555  //totalEtD_GlobalError_ = errorHistCat.at(3).make<TH1I>("totalEtD_GlobalError", "totalEtD_GlobalError;E_{T};Number of Events", 1024, -0.5, 1023.5);
556  //totalEtE_GlobalError_ = errorHistCat.at(3).make<TH1I>("totalEtE_GlobalError", "totalEtE_GlobalError;E_{T};Number of Events", 1024, -0.5, 1023.5);
557  //totalEtD_GlobalError_Of_ = errorHistCat.at(3).make<TH1I>("totalEtD_GlobalError_Of", "totalEtD_GlobalError_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
558  //totalEtE_GlobalError_Of_ = errorHistCat.at(3).make<TH1I>("totalEtE_GlobalError_Of", "totalEtE_GlobalError_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
559  //Total HT
560  totalHtD_ = dataHistCat.at(7).make<TH1I>("totalHtD", "totalHtD;H_{T};Number of Events", 2048, -0.5, 2047.5);
561  totalHtD_Of_ =
562  dataHistCat.at(7).make<TH1I>("totalHtD_Of", "totalHtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
563  totalHtE_ = emuHistCat.at(7).make<TH1I>("totalHtE", "totalHtE;H_{T};Number of Events", 2048, -0.5, 2047.5);
564  totalHtE_Of_ =
565  emuHistCat.at(7).make<TH1I>("totalHtE_Of", "totalHtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
567  errorHistFlags.make<TH1I>("totalHt_errorFlag", "totalHt_errorFlag;Status;Number of Candidates", 2, -0.5, 1.5);
568  //Book the Global HT Error histograms in the errorHistCat
569  //totalHtD_GlobalError_ = errorHistCat.at(3).make<TH1I>("totalHtD_GlobalError", "totalHtD_GlobalError;E_{T};Number of Events", 1024, -0.5, 1023.5);
570  //totalHtE_GlobalError_ = errorHistCat.at(3).make<TH1I>("totalHtE_GlobalError", "totalHtE_GlobalError;E_{T};Number of Events", 1024, -0.5, 1023.5);
571  //totalHtD_GlobalError_Of_ = errorHistCat.at(3).make<TH1I>("totalHtD_GlobalError_Of", "totalHtD_GlobalError_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
572  //totalHtE_GlobalError_Of_ = errorHistCat.at(3).make<TH1I>("totalHtE_GlobalError_Of", "totalHtE_GlobalError_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
573  //MissingEt
574  missingEtD_ =
575  dataHistCat.at(8).make<TH1I>("missingEtD", "missingEtD;Missing E_{T};Number of Events", 1024, -0.5, 1023.5);
577  dataHistCat.at(8).make<TH1I>("missingEtD_Of", "missingEtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
578  missingEtD_Phi_ = dataHistCat.at(8).make<TH1I>(
579  "missingEtD_Phi", "missingEtD_Phi;Missing E_{T} #phi;Number of Events", 72, -0.5, 71.5);
580  missingEtE_ =
581  emuHistCat.at(8).make<TH1I>("missingEtE", "missingEtE;Missing E_{T};Number of Events", 1024, -0.5, 1023.5);
583  emuHistCat.at(8).make<TH1I>("missingEtE_Of", "missingEtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
584  missingEtE_Phi_ = emuHistCat.at(8).make<TH1I>(
585  "missingEtE_Phi", "missingEtE_Phi;Missing E_{T} #phi;Number of Events", 72, -0.5, 71.5);
587  errorHistFlags.make<TH1I>("missingEt_errorFlag", "missingEt_errorFlag;Status;Number of Candidates", 4, -0.5, 3.5);
588  //MissingHt
589  missingHtD_ =
590  dataHistCat.at(8).make<TH1I>("missingHtD", "missingHtD;Missing H_{T};Number of Events", 1024, -0.5, 1023.5);
592  dataHistCat.at(8).make<TH1I>("missingHtD_Of", "missingHtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
593  missingHtD_Phi_ = dataHistCat.at(8).make<TH1I>(
594  "missingHtD_Phi", "missingHtD_Phi;Missing H_{T} #phi;Number of Events", 72, -0.5, 71.5);
595  missingHtE_ =
596  emuHistCat.at(8).make<TH1I>("missingHtE", "missingHtE;Missing H_{T};Number of Events", 1024, -0.5, 1023.5);
598  emuHistCat.at(8).make<TH1I>("missingHtE_Of", "missingHtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
599  missingHtE_Phi_ = emuHistCat.at(8).make<TH1I>(
600  "missingHtE_Phi", "missingHtE_Phi;Missing H_{T} #phi;Number of Events", 72, -0.5, 71.5);
602  errorHistFlags.make<TH1I>("missingHt_errorFlag", "missingHt_errorFlag;Status;Number of Candidates", 4, -0.5, 3.5);
603  //Additional MissingHt Debug histograms
604  missingHtD_HtXPosLeaf1 = dataHistCat.at(8).make<TH1I>(
605  "missingHtD_HtXPosLeaf1", "missingHtD;Missing H_{T} X PosLeaf1;Number of Events", 4096, -2048.5, 2047.5);
606  missingHtD_HtXPosLeaf2 = dataHistCat.at(8).make<TH1I>(
607  "missingHtD_HtXPosLeaf2", "missingHtD;Missing H_{T} X PosLeaf2;Number of Events", 4096, -2048.5, 2047.5);
608  missingHtD_HtXPosLeaf3 = dataHistCat.at(8).make<TH1I>(
609  "missingHtD_HtXPosLeaf3", "missingHtD;Missing H_{T} X PosLeaf3;Number of Events", 4096, -2048.5, 2047.5);
610  missingHtD_HtXNegLeaf1 = dataHistCat.at(8).make<TH1I>(
611  "missingHtD_HtXNegLeaf1", "missingHtD;Missing H_{T} X NegLeaf1;Number of Events", 4096, -2048.5, 2047.5);
612  missingHtD_HtXNegLeaf2 = dataHistCat.at(8).make<TH1I>(
613  "missingHtD_HtXNegLeaf2", "missingHtD;Missing H_{T} X NegLeaf2;Number of Events", 4096, -2048.5, 2047.5);
614  missingHtD_HtXNegLeaf3 = dataHistCat.at(8).make<TH1I>(
615  "missingHtD_HtXNegLeaf3", "missingHtD;Missing H_{T} X NegLeaf3;Number of Events", 4096, -2048.5, 2047.5);
616 
617  missingHtD_HtYPosLeaf1 = dataHistCat.at(8).make<TH1I>(
618  "missingHtD_HtYPosLeaf1", "missingHtD;Missing H_{T} Y PosLeaf1;Number of Events", 4096, -2048.5, 2047.5);
619  missingHtD_HtYPosLeaf2 = dataHistCat.at(8).make<TH1I>(
620  "missingHtD_HtYPosLeaf2", "missingHtD;Missing H_{T} Y PosLeaf2;Number of Events", 4096, -2048.5, 2047.5);
621  missingHtD_HtYPosLeaf3 = dataHistCat.at(8).make<TH1I>(
622  "missingHtD_HtYPosLeaf3", "missingHtD;Missing H_{T} Y PosLeaf3;Number of Events", 4096, -2048.5, 2047.5);
623  missingHtD_HtYNegLeaf1 = dataHistCat.at(8).make<TH1I>(
624  "missingHtD_HtYNegLeaf1", "missingHtD;Missing H_{T} Y NegLeaf1;Number of Events", 4096, -2048.5, 2047.5);
625  missingHtD_HtYNegLeaf2 = dataHistCat.at(8).make<TH1I>(
626  "missingHtD_HtYNegLeaf2", "missingHtD;Missing H_{T} Y NegLeaf2;Number of Events", 4096, -2048.5, 2047.5);
627  missingHtD_HtYNegLeaf3 = dataHistCat.at(8).make<TH1I>(
628  "missingHtD_HtYNegLeaf3", "missingHtD;Missing H_{T} Y NegLeaf3;Number of Events", 4096, -2048.5, 2047.5);
629 
630  //Annotate the labels of the error flags
631  //For the electrons and jets
632  std::vector<std::string> errorFlagLabels;
633  errorFlagLabels.push_back("Matched");
634  errorFlagLabels.push_back("Unmatched Data Cand");
635  errorFlagLabels.push_back("Unmatched Emul Cand");
636 
637  for (unsigned int i = 0; i < errorFlagLabels.size(); i++) {
638  isoEg_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
639  nonIsoEg_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
640  cenJet_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
641  tauJet_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
642  forJet_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
643  }
644  errorFlagLabels.clear();
645 
646  //For the Total Energy Sums and HF
647  errorFlagLabels.push_back("Matched");
648  errorFlagLabels.push_back("Unmatched");
649 
650  for (unsigned int i = 0; i < errorFlagLabels.size(); i++) {
651  hfRingSum_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
652  hfBitCount_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
653  totalEt_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
654  totalHt_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
655  }
656  errorFlagLabels.clear();
657 
658  //For the Missing Energy Sums
659  errorFlagLabels.push_back("Matched");
660  errorFlagLabels.push_back("Matched Mag");
661  errorFlagLabels.push_back("Matched Phi");
662  errorFlagLabels.push_back("Unmatched");
663 
664  for (unsigned int i = 0; i < errorFlagLabels.size(); i++) {
665  missingEt_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
666  missingHt_errorFlag_->GetXaxis()->SetBinLabel(i + 1, errorFlagLabels.at(i).c_str());
667  }
668 
669  //initialise - set all flags to false as they will be set on an event-by-event basis
670  isIsoError = false;
671  isNonIsoError = false;
672  isCenJetError = false;
673  isTauJetError = false;
674  isForJetError = false;
675  isRingSumError = false;
676  isBitCountError = false;
677  isTotalEError = false;
678  isTotalHError = false;
679  isMissingEError = false;
680  isMissingHError = false;
681 
682  //fill the struct of MBXinformation. It is easier to pass this information to the respective functions as used below this way
686 
687  //set the parameters according to the system chosen
688  if (useSys_ == "P5") {
690  } else if (useSys_ == "Lab") {
692  } else {
693  edm::LogWarning("ChosenSystem") << " "
694  << "The system you chose to use (" << useSys_
695  << ") was not recognised. Defaulting to the full system geometry";
697  }
698 }
T getUntrackedParameter(std::string const &, T const &) const
TH1I * tauJetD_GlobalError_Rank_
GctErrorAnalyzerMBxInfo MBxInfo
TH1I * nonIsoEgE_GlobalError_Rank_
TH2I * isoEgD_GlobalError_EtEtaPhi_
TH2I * forJetE_GlobalError_EtEtaPhi_
TH2I * cenJetD_GlobalError_EtEtaPhi_
TH1I * nonIsoEgD_GlobalError_Rank_
TH2I * cenJetE_GlobalError_EtEtaPhi_
TH2I * forJetD_GlobalError_EtEtaPhi_
const unsigned int RCT_REGION_QUANTA_LAB
TH1I * cenJetE_GlobalError_Rank_
T * make(const Args &...args) const
make new ROOT object
edm::InputTag emuTag_
TH2I * tauJetE_GlobalError_EtEtaPhi_
TH2I * nonIsoEgE_GlobalError_EtEtaPhi_
TH1I * forJetE_GlobalError_Rank_
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
TH1I * cenJetD_GlobalError_Rank_
TH2I * isoEgE_GlobalError_EtEtaPhi_
edm::InputTag dataTag_
TH1I * forJetD_GlobalError_Rank_
std::vector< TFileDirectory > errorHistCat
const unsigned int RCT_REGION_QUANTA_P5
TH1I * tauJetE_GlobalError_Rank_
TH2I * nonIsoEgD_GlobalError_EtEtaPhi_
Log< level::Warning, false > LogWarning
TH2I * tauJetD_GlobalError_EtEtaPhi_
const unsigned int * RCT_REGION_QUANTA
GctErrorAnalyzer::~GctErrorAnalyzer ( )
override

Definition at line 700 of file GctErrorAnalyzer.cc.

700  {
701  // do anything here that needs to be done at desctruction time
702  // (e.g. close files, deallocate resources etc.)
703 }

Member Function Documentation

void GctErrorAnalyzer::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 710 of file GctErrorAnalyzer.cc.

References cenJet_errorFlag_, cenJetD_GlobalError_EtEtaPhi_, cenJetD_GlobalError_Rank_, cenJetE_GlobalError_EtEtaPhi_, cenJetE_GlobalError_Rank_, checkCollections(), dataTag_, doBitCountDebug_, doCenJetsDebug_, compareCands< T >::doCompare(), compareTotalEnergySums< T >::doCompare(), compareMissingEnergySums< T >::doCompare(), compareRingSums::doCompare(), compareBitCounts::doCompare(), doEg_, doExtraMissingHTDebug_, doForJetsDebug_, doGCTMBx_, doHF_, doIsoDebug_, doJets_, doMissingEnergySums_, doMissingETDebug_, doMissingHTDebug_, doNonIsoDebug_, doRCT_, doRingSumDebug_, doTauJetsDebug_, doTotalEnergySums_, doTotalEtDebug_, doTotalHtDebug_, emuTag_, edm::EventID::event(), eventNumber, forJet_errorFlag_, forJetD_GlobalError_EtEtaPhi_, forJetD_GlobalError_Rank_, forJetE_GlobalError_EtEtaPhi_, forJetE_GlobalError_Rank_, GCT_INT_HTMISS_QUANTA, GCT_OBJECT_QUANTA, GCT_SUMS_QUANTA, GCTTrigBx_, edm::Event::getByLabel(), hfBitCount_errorFlag_, hfRingSum_errorFlag_, mps_fire::i, edm::EventBase::id(), isBitCountError, isCenJetError, isForJetError, isIsoError, isMissingEError, isMissingHError, isNonIsoError, isoEg_errorFlag_, isoEgD_GlobalError_EtEtaPhi_, isoEgD_GlobalError_Rank_, isoEgE_GlobalError_EtEtaPhi_, isoEgE_GlobalError_Rank_, isRingSumError, isTauJetError, isTotalEError, isTotalHError, edm::InputTag::label(), MBxInfo, missingEt_errorFlag_, missingHt_errorFlag_, missingHtD_HtXNegLeaf1, missingHtD_HtXNegLeaf2, missingHtD_HtXNegLeaf3, missingHtD_HtXPosLeaf1, missingHtD_HtXPosLeaf2, missingHtD_HtXPosLeaf3, missingHtD_HtYNegLeaf1, missingHtD_HtYNegLeaf2, missingHtD_HtYNegLeaf3, missingHtD_HtYPosLeaf1, missingHtD_HtYPosLeaf2, missingHtD_HtYPosLeaf3, nonIsoEg_errorFlag_, nonIsoEgD_GlobalError_EtEtaPhi_, nonIsoEgD_GlobalError_Rank_, nonIsoEgE_GlobalError_EtEtaPhi_, nonIsoEgE_GlobalError_Rank_, NUM_INT_JETS, plotCenJets(), plotEGErrors(), plotForJets(), plotHFBitCounts(), plotHFErrors(), plotHFRingSums(), plotIntJets(), plotIsoEm(), plotJetErrors(), plotMissingEErrors(), plotMissingEt(), plotMissingHt(), plotNonIsoEm(), plotRCTRegions(), plotTauJets(), plotTotalE(), plotTotalEErrors(), plotTotalH(), RCT_EM_OBJECT_QUANTA, RCT_REGION_QUANTA, tauJet_errorFlag_, tauJetD_GlobalError_EtEtaPhi_, tauJetD_GlobalError_Rank_, tauJetE_GlobalError_EtEtaPhi_, tauJetE_GlobalError_Rank_, totalEt_errorFlag_, and totalHt_errorFlag_.

710  {
711  using namespace edm;
712  using namespace std;
713 
714  Handle<L1CaloRegionCollection> caloRegions;
715  Handle<L1CaloEmCollection> emRegions;
716 
721 
728 
730 
733 
736 
739 
742 
745 
748 
750 
751  //we need this for all user cases...
752  iEvent.getByLabel(dataTag_.label(), caloRegions);
753 
754  //in order to allow the debug folders to have a unique name (so that when jobs are split in crab, we can merge)
755  //use the eventnum in the folder name
756  eventNumber = iEvent.id().event();
757 
758  if (doRCT_) {
759  if (checkCollections(caloRegions, *RCT_REGION_QUANTA, "RCT CaloRegions"))
760  plotRCTRegions(caloRegions);
761  }
762 
763  if (doEg_) {
764  iEvent.getByLabel(dataTag_.label(), "nonIsoEm", nonIsoEgD);
765  iEvent.getByLabel(emuTag_.label(), "nonIsoEm", nonIsoEgE);
766 
767  iEvent.getByLabel(dataTag_.label(), "isoEm", isoEgD);
768  iEvent.getByLabel(emuTag_.label(), "isoEm", isoEgE);
769 
770  isIsoError = false;
771  isNonIsoError = false;
772 
773  if (checkCollections(isoEgD, GCT_OBJECT_QUANTA, "Iso e/g Data") &&
774  checkCollections(isoEgE, GCT_OBJECT_QUANTA, "Iso e/g Emulator")) {
775  plotIsoEm(isoEgD, isoEgE);
776  compareEG isoCompare(isoEgD, isoEgE, MBxInfo);
777  isIsoError = isoCompare.doCompare(isoEg_errorFlag_,
782  }
783 
784  if (checkCollections(nonIsoEgD, GCT_OBJECT_QUANTA, "NonIso e/g Data") &&
785  checkCollections(nonIsoEgE, GCT_OBJECT_QUANTA, "NonIso e/g Emulator")) {
786  plotNonIsoEm(nonIsoEgD, nonIsoEgE);
787  compareEG nonIsoCompare(nonIsoEgD, nonIsoEgE, MBxInfo);
788  isNonIsoError = nonIsoCompare.doCompare(nonIsoEg_errorFlag_,
793  }
794 
796  iEvent.getByLabel(dataTag_.label(), emRegions);
797  if (checkCollections(emRegions, RCT_EM_OBJECT_QUANTA, "RCT EMRegions"))
798  plotEGErrors(isoEgD, isoEgE, nonIsoEgD, nonIsoEgE, emRegions);
799  }
800  }
801 
802  if (doJets_) {
803  iEvent.getByLabel(emuTag_.label(), "cenJets", cenJetsE);
804  iEvent.getByLabel(dataTag_.label(), "cenJets", cenJetsD);
805 
806  iEvent.getByLabel(emuTag_.label(), "forJets", forJetsE);
807  iEvent.getByLabel(dataTag_.label(), "forJets", forJetsD);
808 
809  iEvent.getByLabel(emuTag_.label(), "tauJets", tauJetsE);
810  iEvent.getByLabel(dataTag_.label(), "tauJets", tauJetsD);
811 
812  iEvent.getByLabel(emuTag_.label(), intJetsE);
813 
814  isCenJetError = false;
815  isTauJetError = false;
816  isForJetError = false;
817 
818  //Central Jets
819  if (checkCollections(cenJetsD, GCT_OBJECT_QUANTA, "Central Jets Data") &&
820  checkCollections(cenJetsE, GCT_OBJECT_QUANTA, "Central Jets Emulator")) {
821  plotCenJets(cenJetsD, cenJetsE);
822  compareJets cenJetsCompare(cenJetsD, cenJetsE, MBxInfo);
823  isCenJetError = cenJetsCompare.doCompare(cenJet_errorFlag_,
828  }
829 
830  //Tau Jets
831  if (checkCollections(tauJetsD, GCT_OBJECT_QUANTA, "Tau Jets Data") &&
832  checkCollections(tauJetsE, GCT_OBJECT_QUANTA, "Tau Jets Emulator")) {
833  plotTauJets(tauJetsD, tauJetsE);
834  compareJets tauJetsCompare(tauJetsD, tauJetsE, MBxInfo);
835  isTauJetError = tauJetsCompare.doCompare(tauJet_errorFlag_,
840  }
841 
842  //For Jets
843  if (checkCollections(forJetsD, GCT_OBJECT_QUANTA, "Forward Jets Data") &&
844  checkCollections(forJetsE, GCT_OBJECT_QUANTA, "Forward Jets Emulator")) {
845  plotForJets(forJetsD, forJetsE);
846  compareJets forJetsCompare(forJetsD, forJetsE, MBxInfo);
847  isForJetError = forJetsCompare.doCompare(forJet_errorFlag_,
852  }
853 
854  //Emulator Intermediate Jets
855  if (checkCollections(intJetsE, NUM_INT_JETS, "Intermediate Jets Emulator"))
856  plotIntJets(intJetsE);
857 
860  plotJetErrors(cenJetsD, cenJetsE, tauJetsD, tauJetsE, forJetsD, forJetsE, caloRegions);
861  }
862  }
863 
864  if (doHF_) {
865  iEvent.getByLabel(dataTag_.label(), hfRingSumsD);
866  iEvent.getByLabel(emuTag_.label(), hfRingSumsE);
867 
868  iEvent.getByLabel(dataTag_.label(), hfBitCountsD);
869  iEvent.getByLabel(emuTag_.label(), hfBitCountsE);
870 
871  isRingSumError = false;
872  isBitCountError = false;
873 
874  if (checkCollections(hfRingSumsD, GCT_SUMS_QUANTA, "HF Ring Sums Data") &&
875  checkCollections(hfRingSumsE, GCT_SUMS_QUANTA, "HF Ring Sums Emulator")) {
876  plotHFRingSums(hfRingSumsD, hfRingSumsE);
877  compareRingSums HFRingSums(hfRingSumsD, hfRingSumsE, MBxInfo);
878  isRingSumError = HFRingSums.doCompare(hfRingSum_errorFlag_);
879  }
880 
881  if (checkCollections(hfBitCountsD, GCT_SUMS_QUANTA, "HF Bit Counts Data") &&
882  checkCollections(hfBitCountsE, GCT_SUMS_QUANTA, "HF Bit Counts Emulator")) {
883  plotHFBitCounts(hfBitCountsD, hfBitCountsE);
884  compareBitCounts HFBitCounts(hfBitCountsD, hfBitCountsE, MBxInfo);
885  isBitCountError = HFBitCounts.doCompare(hfBitCount_errorFlag_);
886  }
887 
889  plotHFErrors(hfRingSumsD, hfRingSumsE, hfBitCountsD, hfBitCountsE, caloRegions);
890  }
891  }
892 
893  if (doTotalEnergySums_) {
894  iEvent.getByLabel(dataTag_.label(), totalEtD);
895  iEvent.getByLabel(emuTag_.label(), totalEtE);
896 
897  iEvent.getByLabel(dataTag_.label(), totalHtD);
898  iEvent.getByLabel(emuTag_.label(), totalHtE);
899 
900  isTotalEError = false;
901  isTotalHError = false;
902 
903  if (checkCollections(totalEtD, GCT_SUMS_QUANTA, "Total Et Data") &&
904  checkCollections(totalEtE, GCT_SUMS_QUANTA, "Total Et Emulator")) {
905  plotTotalE(totalEtD, totalEtE);
906  compareTotalE compareET(totalEtD, totalEtE, MBxInfo);
907  isTotalEError = compareET.doCompare(totalEt_errorFlag_);
908  }
909 
910  if (checkCollections(totalHtD, GCT_SUMS_QUANTA, "Total Ht Data") &&
911  checkCollections(totalHtE, GCT_SUMS_QUANTA, "Total Ht Emulator")) {
912  plotTotalH(totalHtD, totalHtE);
913  compareTotalH compareHT(totalHtD, totalHtE, MBxInfo);
914  isTotalHError = compareHT.doCompare(totalHt_errorFlag_);
915  }
916 
918  plotTotalEErrors(totalEtD, totalEtE, totalHtD, totalHtE, caloRegions);
919  }
920  }
921 
922  if (doMissingEnergySums_) {
923  iEvent.getByLabel(dataTag_.label(), missingEtD);
924  iEvent.getByLabel(emuTag_.label(), missingEtE);
925 
926  iEvent.getByLabel(dataTag_.label(), missingHtD);
927  iEvent.getByLabel(emuTag_.label(), missingHtE);
928 
929  isMissingEError = false;
930  isMissingHError = false;
931 
932  if (checkCollections(missingEtD, GCT_SUMS_QUANTA, "Missing Et Data") &&
933  checkCollections(missingEtE, GCT_SUMS_QUANTA, "Missing Et Emulator")) {
934  plotMissingEt(missingEtD, missingEtE);
935  compareMissingE compareMET(missingEtD, missingEtE, MBxInfo);
936  isMissingEError = compareMET.doCompare(missingEt_errorFlag_);
937  }
938 
939  if (checkCollections(missingHtD, GCT_SUMS_QUANTA, "Missing Ht Data") &&
940  checkCollections(missingHtE, GCT_SUMS_QUANTA, "Missing Ht Emulator")) {
941  plotMissingHt(missingHtD, missingHtE);
942  compareMissingH compareMHT(missingHtD, missingHtE, MBxInfo);
943  isMissingHError = compareMHT.doCompare(missingHt_errorFlag_);
944 
945  //added 19/03/2010 for intermediate information on MissingHt quantities in the data
947  iEvent.getByLabel(dataTag_.label(), "", intHtMissD);
948  if (checkCollections(intHtMissD, GCT_INT_HTMISS_QUANTA, "Internal Missing Ht Data")) {
949  for (unsigned int i = 0; i < intHtMissD->size(); i++) {
950  if (doGCTMBx_ || intHtMissD->at(i).bx() == GCTTrigBx_) {
951  if (!intHtMissD->at(i).overflow()) {
952  //the capBlock 0x301 is the input pipeline at the wheel for positive eta, whereas 0x701 is for negative eta
953  if (intHtMissD->at(i).capBlock() == 0x301 && intHtMissD->at(i).capIndex() == 0 &&
954  intHtMissD->at(i).isThereHtx())
955  missingHtD_HtXPosLeaf1->Fill(intHtMissD->at(i).htx());
956  if (intHtMissD->at(i).capBlock() == 0x301 && intHtMissD->at(i).capIndex() == 1 &&
957  intHtMissD->at(i).isThereHtx())
958  missingHtD_HtXPosLeaf2->Fill(intHtMissD->at(i).htx());
959  if (intHtMissD->at(i).capBlock() == 0x301 && intHtMissD->at(i).capIndex() == 2 &&
960  intHtMissD->at(i).isThereHtx())
961  missingHtD_HtXPosLeaf3->Fill(intHtMissD->at(i).htx());
962  if (intHtMissD->at(i).capBlock() == 0x701 && intHtMissD->at(i).capIndex() == 0 &&
963  intHtMissD->at(i).isThereHtx())
964  missingHtD_HtXNegLeaf1->Fill(intHtMissD->at(i).htx());
965  if (intHtMissD->at(i).capBlock() == 0x701 && intHtMissD->at(i).capIndex() == 1 &&
966  intHtMissD->at(i).isThereHtx())
967  missingHtD_HtXNegLeaf2->Fill(intHtMissD->at(i).htx());
968  if (intHtMissD->at(i).capBlock() == 0x701 && intHtMissD->at(i).capIndex() == 2 &&
969  intHtMissD->at(i).isThereHtx())
970  missingHtD_HtXNegLeaf3->Fill(intHtMissD->at(i).htx());
971 
972  if (intHtMissD->at(i).capBlock() == 0x301 && intHtMissD->at(i).capIndex() == 0 &&
973  intHtMissD->at(i).isThereHty())
974  missingHtD_HtYPosLeaf1->Fill(intHtMissD->at(i).hty());
975  if (intHtMissD->at(i).capBlock() == 0x301 && intHtMissD->at(i).capIndex() == 1 &&
976  intHtMissD->at(i).isThereHty())
977  missingHtD_HtYPosLeaf2->Fill(intHtMissD->at(i).hty());
978  if (intHtMissD->at(i).capBlock() == 0x301 && intHtMissD->at(i).capIndex() == 2 &&
979  intHtMissD->at(i).isThereHty())
980  missingHtD_HtYPosLeaf3->Fill(intHtMissD->at(i).hty());
981  if (intHtMissD->at(i).capBlock() == 0x701 && intHtMissD->at(i).capIndex() == 0 &&
982  intHtMissD->at(i).isThereHty())
983  missingHtD_HtYNegLeaf1->Fill(intHtMissD->at(i).hty());
984  if (intHtMissD->at(i).capBlock() == 0x701 && intHtMissD->at(i).capIndex() == 1 &&
985  intHtMissD->at(i).isThereHty())
986  missingHtD_HtYNegLeaf2->Fill(intHtMissD->at(i).hty());
987  if (intHtMissD->at(i).capBlock() == 0x701 && intHtMissD->at(i).capIndex() == 2 &&
988  intHtMissD->at(i).isThereHty())
989  missingHtD_HtYNegLeaf3->Fill(intHtMissD->at(i).hty());
990  }
991  }
992  }
993  }
994  }
995  }
996 
998  plotMissingEErrors(missingEtD, missingEtE, missingHtD, missingHtE, caloRegions, intJetsE, intHtMissD);
999  }
1000  }
1001 }
void plotHFErrors(const edm::Handle< L1GctHFRingEtSumsCollection > &hfRingSumsD, const edm::Handle< L1GctHFRingEtSumsCollection > &hfRingSumsE, const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsD, const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsE, const edm::Handle< L1CaloRegionCollection > &caloRegions)
unsigned int eventNumber
EventNumber_t event() const
Definition: EventID.h:40
TH1I * tauJetD_GlobalError_Rank_
void plotMissingHt(const edm::Handle< L1GctHtMissCollection > &missingHtD, const edm::Handle< L1GctHtMissCollection > &missingHtE)
GctErrorAnalyzerMBxInfo MBxInfo
TH1I * nonIsoEgE_GlobalError_Rank_
TH2I * isoEgD_GlobalError_EtEtaPhi_
TH2I * forJetE_GlobalError_EtEtaPhi_
void plotHFBitCounts(const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsD, const edm::Handle< L1GctHFBitCountsCollection > &hfBitCountsE)
TH2I * cenJetD_GlobalError_EtEtaPhi_
const unsigned int GCT_SUMS_QUANTA
void plotForJets(const edm::Handle< L1GctJetCandCollection > &data, const edm::Handle< L1GctJetCandCollection > &emu)
void plotJetErrors(const edm::Handle< L1GctJetCandCollection > &cendata, const edm::Handle< L1GctJetCandCollection > &cenemu, const edm::Handle< L1GctJetCandCollection > &taudata, const edm::Handle< L1GctJetCandCollection > &tauemu, const edm::Handle< L1GctJetCandCollection > &fordata, const edm::Handle< L1GctJetCandCollection > &foremu, const edm::Handle< L1CaloRegionCollection > &regions)
TH1I * nonIsoEgD_GlobalError_Rank_
TH2I * cenJetE_GlobalError_EtEtaPhi_
void plotCenJets(const edm::Handle< L1GctJetCandCollection > &data, const edm::Handle< L1GctJetCandCollection > &emu)
void plotTotalH(const edm::Handle< L1GctEtHadCollection > &totalHtD, const edm::Handle< L1GctEtHadCollection > &totalHtE)
TH2I * forJetD_GlobalError_EtEtaPhi_
void plotEGErrors(const edm::Handle< L1GctEmCandCollection > &dataiso, const edm::Handle< L1GctEmCandCollection > &emuiso, const edm::Handle< L1GctEmCandCollection > &datanoniso, const edm::Handle< L1GctEmCandCollection > &emunoniso, const edm::Handle< L1CaloEmCollection > &regions)
void plotTotalE(const edm::Handle< L1GctEtTotalCollection > &totalEtD, const edm::Handle< L1GctEtTotalCollection > &totalEtE)
void plotTotalEErrors(const edm::Handle< L1GctEtTotalCollection > &totalEtD, const edm::Handle< L1GctEtTotalCollection > &totalEtE, const edm::Handle< L1GctEtHadCollection > &totalHtD, const edm::Handle< L1GctEtHadCollection > &totalHtE, const edm::Handle< L1CaloRegionCollection > &caloRegions)
void plotIsoEm(const edm::Handle< L1GctEmCandCollection > &data, const edm::Handle< L1GctEmCandCollection > &emu)
void plotNonIsoEm(const edm::Handle< L1GctEmCandCollection > &data, const edm::Handle< L1GctEmCandCollection > &emu)
void plotTauJets(const edm::Handle< L1GctJetCandCollection > &data, const edm::Handle< L1GctJetCandCollection > &emu)
void plotHFRingSums(const edm::Handle< L1GctHFRingEtSumsCollection > &data, const edm::Handle< L1GctHFRingEtSumsCollection > &emu)
void plotRCTRegions(const edm::Handle< L1CaloRegionCollection > &caloRegions)
TH1I * cenJetE_GlobalError_Rank_
void plotIntJets(const edm::Handle< L1GctInternJetDataCollection > &emu)
edm::InputTag emuTag_
const unsigned int GCT_OBJECT_QUANTA
const unsigned int RCT_EM_OBJECT_QUANTA
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:500
TH2I * tauJetE_GlobalError_EtEtaPhi_
TH2I * nonIsoEgE_GlobalError_EtEtaPhi_
TH1I * forJetE_GlobalError_Rank_
void plotMissingEt(const edm::Handle< L1GctEtMissCollection > &missingEtD, const edm::Handle< L1GctEtMissCollection > &missingEtE)
void plotMissingEErrors(const edm::Handle< L1GctEtMissCollection > &missingEtD, const edm::Handle< L1GctEtMissCollection > &missingEtE, const edm::Handle< L1GctHtMissCollection > &missingHtD, const edm::Handle< L1GctHtMissCollection > &missingHtE, edm::Handle< L1CaloRegionCollection > &caloRegions, const edm::Handle< L1GctInternJetDataCollection > &intjetsemu, const edm::Handle< L1GctInternHtMissCollection > intMissingHtD)
TH1I * cenJetD_GlobalError_Rank_
TH2I * isoEgE_GlobalError_EtEtaPhi_
edm::InputTag dataTag_
bool checkCollections(const T &collection, const unsigned int &constraint, const std::string &label)
std::string const & label() const
Definition: InputTag.h:36
edm::EventID id() const
Definition: EventBase.h:59
TH1I * forJetD_GlobalError_Rank_
TH1I * tauJetE_GlobalError_Rank_
const unsigned int NUM_INT_JETS
TH2I * nonIsoEgD_GlobalError_EtEtaPhi_
const unsigned int GCT_INT_HTMISS_QUANTA
TH2I * tauJetD_GlobalError_EtEtaPhi_
const unsigned int * RCT_REGION_QUANTA
void GctErrorAnalyzer::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1004 of file GctErrorAnalyzer.cc.

1004 {}
template<class T >
bool GctErrorAnalyzer::checkCollections ( const T collection,
const unsigned int &  constraint,
const std::string &  label 
)
private

Definition at line 1316 of file GctErrorAnalyzer.cc.

References HLT_FULL_cff::constraint.

Referenced by analyze(), and plotMissingEErrors().

1316  {
1317  //unfortunately, the dataformats are not consistent with the name() method (i.e. some have it, others don't)
1318  //and a typeof() function doesn't exist in ANSI C++, so to identify the templated type, we pass a std::string
1319 
1320  if (!collection.isValid()) {
1321  edm::LogWarning("DataNotFound") << " Could not find " << label << " label";
1322  return false;
1323  }
1324  if (collection->size() % constraint != 0 || collection->empty()) {
1325  edm::LogWarning("CollectionSizeError")
1326  << " " << label << " collection size is " << collection->size() << ", expected multiple of " << constraint;
1327  return false;
1328  }
1329 
1330  return true;
1331 }
char const * label
Log< level::Warning, false > LogWarning
void GctErrorAnalyzer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1007 of file GctErrorAnalyzer.cc.

1007 {}
GctErrorAnalyzer GctErrorAnalyzer::operator= ( const GctErrorAnalyzer )
delete
void GctErrorAnalyzer::plotCenJets ( const edm::Handle< L1GctJetCandCollection > &  data,
const edm::Handle< L1GctJetCandCollection > &  emu 
)
private

Definition at line 1171 of file GctErrorAnalyzer.cc.

References cenJetD_EtEtaPhi_, cenJetD_OccEtaPhi_, cenJetD_Rank_, cenJetE_EtEtaPhi_, cenJetE_OccEtaPhi_, cenJetE_Rank_, doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, and mps_fire::i.

Referenced by analyze().

1172  {
1173  for (unsigned int i = 0; i < cenJetsD->size(); i++) {
1174  if (doGCTMBx_ || cenJetsD->at(i).bx() == GCTTrigBx_) {
1175  cenJetD_Rank_->Fill(cenJetsD->at(i).rank());
1176  if (cenJetsD->at(i).rank() > 0) {
1177  cenJetD_EtEtaPhi_->Fill(
1178  cenJetsD->at(i).regionId().ieta(), cenJetsD->at(i).regionId().iphi(), cenJetsD->at(i).rank());
1179  cenJetD_OccEtaPhi_->Fill(cenJetsD->at(i).regionId().ieta(), cenJetsD->at(i).regionId().iphi());
1180  }
1181  }
1182  }
1183 
1184  for (unsigned int i = 0; i < cenJetsE->size(); i++) {
1185  if (doEmuMBx_ || cenJetsE->at(i).bx() == EmuTrigBx_) {
1186  cenJetE_Rank_->Fill(cenJetsE->at(i).rank());
1187  if (cenJetsE->at(i).rank() > 0) {
1188  cenJetE_EtEtaPhi_->Fill(
1189  cenJetsE->at(i).regionId().ieta(), cenJetsE->at(i).regionId().iphi(), cenJetsE->at(i).rank());
1190  cenJetE_OccEtaPhi_->Fill(cenJetsE->at(i).regionId().ieta(), cenJetsE->at(i).regionId().iphi());
1191  }
1192  }
1193  }
1194 }
void GctErrorAnalyzer::plotEGErrors ( const edm::Handle< L1GctEmCandCollection > &  dataiso,
const edm::Handle< L1GctEmCandCollection > &  emuiso,
const edm::Handle< L1GctEmCandCollection > &  datanoniso,
const edm::Handle< L1GctEmCandCollection > &  emunoniso,
const edm::Handle< L1CaloEmCollection > &  regions 
)
private

Definition at line 1079 of file GctErrorAnalyzer.cc.

References EmuTrigBx_, errorHistCat, eventNumber, GCTTrigBx_, mps_fire::i, isIsoError, isNonIsoError, TFileDirectory::make(), RCTTrigBx_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

1083  {
1084  std::string errorDirName = "err_";
1085  if (isIsoError)
1086  errorDirName.append("I");
1087  if (isNonIsoError)
1088  errorDirName.append("N");
1089  std::stringstream caseNumber;
1090  caseNumber << eventNumber;
1091  errorDirName.append(caseNumber.str());
1092  TFileDirectory errorDir = errorHistCat.at(0).mkdir(errorDirName);
1093 
1094  TH2I *errorEmRegionIsoEtEtaPhi_ = errorDir.make<TH2I>("errorEmRegionIsoEtEtaPhi",
1095  "errorEmRegionIsoEtEtaPhi;#eta (GCT Units);#phi (GCT Units)",
1096  22,
1097  -0.5,
1098  21.5,
1099  18,
1100  -0.5,
1101  17.5);
1102  TH2I *errorEmRegionNonIsoEtEtaPhi_ =
1103  errorDir.make<TH2I>("errorEmRegionNonIsoEtEtaPhi",
1104  "errorEmRegionNonIsoEtEtaPhi;#eta (GCT Units);#phi (GCT Units)",
1105  22,
1106  -0.5,
1107  21.5,
1108  18,
1109  -0.5,
1110  17.5);
1111  TH2I *errorIsoEtEtaPhiD_ = errorDir.make<TH2I>(
1112  "errorIsoEtEtaPhiD", "errorIsoEtEtaPhiD;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1113  TH2I *errorNonIsoEtEtaPhiD_ = errorDir.make<TH2I>(
1114  "errorNonIsoEtEtaPhiD", "errorNonIsoEtEtaPhiD;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1115  TH2I *errorIsoEtEtaPhiE_ = errorDir.make<TH2I>(
1116  "errorIsoEtEtaPhiE", "errorIsoEtEtaPhiE;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1117  TH2I *errorNonIsoEtEtaPhiE_ = errorDir.make<TH2I>(
1118  "errorNonIsoEtEtaPhiE", "errorNonIsoEtEtaPhiE;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1119 
1120  //fill the EM input collection
1121  //should only fill the correct bx for emRegions - and since this is showing an error in the comparison, we should plot the input to this comparison i.e. Bx=RCTTrigBx
1122  //this assumes that comparison is done on the central Bx i.e. RctBx=0 corresponds to GctBx=0, and EmuBx=0 takes RctBx=0
1123  for (unsigned int i = 0; i < emRegions->size(); i++) {
1124  if (emRegions->at(i).bx() == RCTTrigBx_) {
1125  if (emRegions->at(i).isolated()) {
1126  if (emRegions->at(i).rank() > 0)
1127  errorEmRegionIsoEtEtaPhi_->Fill(
1128  emRegions->at(i).regionId().ieta(), emRegions->at(i).regionId().iphi(), emRegions->at(i).rank());
1129  } else {
1130  if (emRegions->at(i).rank() > 0)
1131  errorEmRegionNonIsoEtEtaPhi_->Fill(
1132  emRegions->at(i).regionId().ieta(), emRegions->at(i).regionId().iphi(), emRegions->at(i).rank());
1133  }
1134  }
1135  }
1136 
1137  //no need to have the rank plot, because you can't have two electrons in the same place (eta,phi), in the same event...
1138  //in this case, since we're actually comparing the GCTTrigBx_ with the EmuTrigBx_, we plot these individually
1139  for (unsigned int i = 0; i < isoEgD->size(); i++) {
1140  if (isoEgD->at(i).bx() == GCTTrigBx_) {
1141  if (isoEgD->at(i).rank() > 0)
1142  errorIsoEtEtaPhiD_->Fill(
1143  isoEgD->at(i).regionId().ieta(), isoEgD->at(i).regionId().iphi(), isoEgD->at(i).rank());
1144  }
1145  }
1146  for (unsigned int i = 0; i < nonIsoEgD->size(); i++) {
1147  if (nonIsoEgD->at(i).bx() == GCTTrigBx_) {
1148  if (nonIsoEgD->at(i).rank() > 0)
1149  errorNonIsoEtEtaPhiD_->Fill(
1150  nonIsoEgD->at(i).regionId().ieta(), nonIsoEgD->at(i).regionId().iphi(), nonIsoEgD->at(i).rank());
1151  }
1152  }
1153 
1154  //now for the emulator candidates
1155  for (unsigned int i = 0; i < isoEgE->size(); i++) {
1156  if (isoEgE->at(i).bx() == EmuTrigBx_) {
1157  if (isoEgE->at(i).rank() > 0)
1158  errorIsoEtEtaPhiE_->Fill(
1159  isoEgE->at(i).regionId().ieta(), isoEgE->at(i).regionId().iphi(), isoEgE->at(i).rank());
1160  }
1161  }
1162  for (unsigned int i = 0; i < nonIsoEgE->size(); i++) {
1163  if (nonIsoEgE->at(i).bx() == EmuTrigBx_) {
1164  if (nonIsoEgE->at(i).rank() > 0)
1165  errorNonIsoEtEtaPhiE_->Fill(
1166  nonIsoEgE->at(i).regionId().ieta(), nonIsoEgE->at(i).regionId().iphi(), nonIsoEgE->at(i).rank());
1167  }
1168  }
1169 }
unsigned int eventNumber
T * make(const Args &...args) const
make new ROOT object
std::vector< TFileDirectory > errorHistCat
void GctErrorAnalyzer::plotForJets ( const edm::Handle< L1GctJetCandCollection > &  data,
const edm::Handle< L1GctJetCandCollection > &  emu 
)
private

Definition at line 1221 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, forJetD_EtEtaPhi_, forJetD_OccEtaPhi_, forJetD_Rank_, forJetE_EtEtaPhi_, forJetE_OccEtaPhi_, forJetE_Rank_, GCTTrigBx_, and mps_fire::i.

Referenced by analyze().

1222  {
1223  for (unsigned int i = 0; i < forJetsD->size(); i++) {
1224  if (doGCTMBx_ || forJetsD->at(i).bx() == GCTTrigBx_) {
1225  forJetD_Rank_->Fill(forJetsD->at(i).rank());
1226  if (forJetsD->at(i).rank() > 0) {
1227  forJetD_EtEtaPhi_->Fill(
1228  forJetsD->at(i).regionId().ieta(), forJetsD->at(i).regionId().iphi(), forJetsD->at(i).rank());
1229  forJetD_OccEtaPhi_->Fill(forJetsD->at(i).regionId().ieta(), forJetsD->at(i).regionId().iphi());
1230  }
1231  }
1232  }
1233 
1234  for (unsigned int i = 0; i < forJetsE->size(); i++) {
1235  if (doEmuMBx_ || forJetsE->at(i).bx() == EmuTrigBx_) {
1236  forJetE_Rank_->Fill(forJetsE->at(i).rank());
1237  if (forJetsE->at(i).rank() > 0) {
1238  forJetE_EtEtaPhi_->Fill(
1239  forJetsE->at(i).regionId().ieta(), forJetsE->at(i).regionId().iphi(), forJetsE->at(i).rank());
1240  forJetE_OccEtaPhi_->Fill(forJetsE->at(i).regionId().ieta(), forJetsE->at(i).regionId().iphi());
1241  }
1242  }
1243  }
1244 }
void GctErrorAnalyzer::plotHFBitCounts ( const edm::Handle< L1GctHFBitCountsCollection > &  hfBitCountsD,
const edm::Handle< L1GctHFBitCountsCollection > &  hfBitCountsE 
)
private

Definition at line 1489 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, hfBitCountD_1neg_, hfBitCountD_1pos_, hfBitCountD_2neg_, hfBitCountD_2pos_, hfBitCountE_1neg_, hfBitCountE_1pos_, hfBitCountE_2neg_, hfBitCountE_2pos_, and mps_fire::i.

Referenced by analyze().

1490  {
1491  for (unsigned int i = 0; i < hfBitCountsD->size(); i++) {
1492  if (doGCTMBx_ || hfBitCountsD->at(i).bx() == GCTTrigBx_) {
1493  //there are 4 rings - just fill the histograms
1494  hfBitCountD_1pos_->Fill(hfBitCountsD->at(i).bitCount(0));
1495  hfBitCountD_1neg_->Fill(hfBitCountsD->at(i).bitCount(1));
1496  hfBitCountD_2pos_->Fill(hfBitCountsD->at(i).bitCount(2));
1497  hfBitCountD_2neg_->Fill(hfBitCountsD->at(i).bitCount(3));
1498  }
1499  }
1500  for (unsigned int i = 0; i < hfBitCountsE->size(); i++) {
1501  if (doEmuMBx_ || hfBitCountsE->at(i).bx() == EmuTrigBx_) {
1502  hfBitCountE_1pos_->Fill(hfBitCountsE->at(i).bitCount(0));
1503  hfBitCountE_1neg_->Fill(hfBitCountsE->at(i).bitCount(1));
1504  hfBitCountE_2pos_->Fill(hfBitCountsE->at(i).bitCount(2));
1505  hfBitCountE_2neg_->Fill(hfBitCountsE->at(i).bitCount(3));
1506  }
1507  }
1508 }
void GctErrorAnalyzer::plotHFErrors ( const edm::Handle< L1GctHFRingEtSumsCollection > &  hfRingSumsD,
const edm::Handle< L1GctHFRingEtSumsCollection > &  hfRingSumsE,
const edm::Handle< L1GctHFBitCountsCollection > &  hfBitCountsD,
const edm::Handle< L1GctHFBitCountsCollection > &  hfBitCountsE,
const edm::Handle< L1CaloRegionCollection > &  caloRegions 
)
private

Definition at line 1510 of file GctErrorAnalyzer.cc.

References EmuTrigBx_, errorHistCat, eventNumber, GCTTrigBx_, mps_fire::i, isBitCountError, isRingSumError, TFileDirectory::make(), RCTTrigBx_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

1514  {
1515  std::string errorDirName = "err_";
1516  if (isRingSumError)
1517  errorDirName.append("R");
1518  if (isBitCountError)
1519  errorDirName.append("B");
1520  std::stringstream caseNumber;
1521  caseNumber << eventNumber;
1522  errorDirName.append(caseNumber.str());
1523  TFileDirectory errorDir = errorHistCat.at(2).mkdir(errorDirName);
1524 
1525  TH2I *errorRegionEtEtaPhi_ = errorDir.make<TH2I>(
1526  "errorRegionEtEtaPhi", "errorRegionEtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1527  TH2I *errorRegionFgEtaPhi_ = errorDir.make<TH2I>(
1528  "errorRegionFgEtaPhi", "errorRegionFgEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1529  TH2I *errorRegionOfEtaPhi_ = errorDir.make<TH2I>(
1530  "errorRegionOfEtaPhi", "errorRegionOfEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1531 
1532  TH1I *errorHFRingSumD_1pos_ =
1533  errorDir.make<TH1I>("errorHFRingSumD_1+", "errorHFRingSumD_1+;Rank;Number of Events", 8, -0.5, 7.5);
1534  TH1I *errorHFRingSumD_2pos_ =
1535  errorDir.make<TH1I>("errorHFRingSumD_2+", "errorHFRingSumD_2+;Rank;Number of Events", 8, -0.5, 7.5);
1536  TH1I *errorHFRingSumD_1neg_ =
1537  errorDir.make<TH1I>("errorHFRingSumD_1-", "errorHFRingSumD_1-;Rank;Number of Events", 8, -0.5, 7.5);
1538  TH1I *errorHFRingSumD_2neg_ =
1539  errorDir.make<TH1I>("errorHFRingSumD_2-", "errorHFRingSumD_2-;Rank;Number of Events", 8, -0.5, 7.5);
1540  TH1I *errorHFRingSumE_1pos_ =
1541  errorDir.make<TH1I>("errorHFRingSumE_1+", "errorHFRingSumE_1+;Rank;Number of Events", 8, -0.5, 7.5);
1542  TH1I *errorHFRingSumE_2pos_ =
1543  errorDir.make<TH1I>("errorHFRingSumE_2+", "errorHFRingSumE_2+;Rank;Number of Events", 8, -0.5, 7.5);
1544  TH1I *errorHFRingSumE_1neg_ =
1545  errorDir.make<TH1I>("errorHFRingSumE_1-", "errorHFRingSumE_1-;Rank;Number of Events", 8, -0.5, 7.5);
1546  TH1I *errorHFRingSumE_2neg_ =
1547  errorDir.make<TH1I>("errorHFRingSumE_2-", "errorHFRingSumE_2-;Rank;Number of Events", 8, -0.5, 7.5);
1548 
1549  TH1I *errorHFBitCountD_1pos_ =
1550  errorDir.make<TH1I>("errorHFBitCountD_1+", "errorHFBitCountD_1+;Rank;Number of Events", 8, -0.5, 7.5);
1551  TH1I *errorHFBitCountD_2pos_ =
1552  errorDir.make<TH1I>("errorHFBitCountD_2+", "errorHFBitCountD_2+;Rank;Number of Events", 8, -0.5, 7.5);
1553  TH1I *errorHFBitCountD_1neg_ =
1554  errorDir.make<TH1I>("errorHFBitCountD_1-", "errorHFBitCountD_1-;Rank;Number of Events", 8, -0.5, 7.5);
1555  TH1I *errorHFBitCountD_2neg_ =
1556  errorDir.make<TH1I>("errorHFBitCountD_2-", "errorHFBitCountD_2-;Rank;Number of Events", 8, -0.5, 7.5);
1557  TH1I *errorHFBitCountE_1pos_ =
1558  errorDir.make<TH1I>("errorHFBitCountE_1+", "errorHFBitCountE_1+;Rank;Number of Events", 8, -0.5, 7.5);
1559  TH1I *errorHFBitCountE_2pos_ =
1560  errorDir.make<TH1I>("errorHFBitCountE_2+", "errorHFBitCountE_2+;Rank;Number of Events", 8, -0.5, 7.5);
1561  TH1I *errorHFBitCountE_1neg_ =
1562  errorDir.make<TH1I>("errorHFBitCountE_1-", "errorHFBitCountE_1-;Rank;Number of Events", 8, -0.5, 7.5);
1563  TH1I *errorHFBitCountE_2neg_ =
1564  errorDir.make<TH1I>("errorHFBitCountE_2-", "errorHFBitCountE_2-;Rank;Number of Events", 8, -0.5, 7.5);
1565 
1566  for (unsigned int i = 0; i < caloRegions->size(); i++) {
1567  if (caloRegions->at(i).bx() == RCTTrigBx_) {
1568  if (caloRegions->at(i).et() > 0)
1569  errorRegionEtEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi(), caloRegions->at(i).et());
1570  if (caloRegions->at(i).fineGrain())
1571  errorRegionFgEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1572  if (caloRegions->at(i).overFlow())
1573  errorRegionOfEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1574  }
1575  }
1576 
1577  for (unsigned int i = 0; i < hfRingSumsD->size(); i++) {
1578  if (hfRingSumsD->at(i).bx() == GCTTrigBx_) {
1579  errorHFRingSumD_1pos_->Fill(hfRingSumsD->at(i).etSum(0));
1580  errorHFRingSumD_1neg_->Fill(hfRingSumsD->at(i).etSum(1));
1581  errorHFRingSumD_2pos_->Fill(hfRingSumsD->at(i).etSum(2));
1582  errorHFRingSumD_2neg_->Fill(hfRingSumsD->at(i).etSum(3));
1583  }
1584  }
1585  for (unsigned int i = 0; i < hfRingSumsE->size(); i++) {
1586  if (hfRingSumsE->at(i).bx() == EmuTrigBx_) {
1587  errorHFRingSumE_1pos_->Fill(hfRingSumsE->at(i).etSum(0));
1588  errorHFRingSumE_1neg_->Fill(hfRingSumsE->at(i).etSum(1));
1589  errorHFRingSumE_2pos_->Fill(hfRingSumsE->at(i).etSum(2));
1590  errorHFRingSumE_2neg_->Fill(hfRingSumsE->at(i).etSum(3));
1591  }
1592  }
1593 
1594  for (unsigned int i = 0; i < hfBitCountsD->size(); i++) {
1595  if (hfBitCountsD->at(i).bx() == GCTTrigBx_) {
1596  errorHFBitCountD_1pos_->Fill(hfBitCountsD->at(i).bitCount(0));
1597  errorHFBitCountD_1neg_->Fill(hfBitCountsD->at(i).bitCount(1));
1598  errorHFBitCountD_2pos_->Fill(hfBitCountsD->at(i).bitCount(2));
1599  errorHFBitCountD_2neg_->Fill(hfBitCountsD->at(i).bitCount(3));
1600  }
1601  }
1602  for (unsigned int i = 0; i < hfBitCountsE->size(); i++) {
1603  if (hfBitCountsE->at(i).bx() == EmuTrigBx_) {
1604  errorHFBitCountE_1pos_->Fill(hfBitCountsE->at(i).bitCount(0));
1605  errorHFBitCountE_1neg_->Fill(hfBitCountsE->at(i).bitCount(1));
1606  errorHFBitCountE_2pos_->Fill(hfBitCountsE->at(i).bitCount(2));
1607  errorHFBitCountE_2neg_->Fill(hfBitCountsE->at(i).bitCount(3));
1608  }
1609  }
1610 }
unsigned int eventNumber
T * make(const Args &...args) const
make new ROOT object
std::vector< TFileDirectory > errorHistCat
void GctErrorAnalyzer::plotHFRingSums ( const edm::Handle< L1GctHFRingEtSumsCollection > &  data,
const edm::Handle< L1GctHFRingEtSumsCollection > &  emu 
)
private

Definition at line 1467 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, hfRingSumD_1neg_, hfRingSumD_1pos_, hfRingSumD_2neg_, hfRingSumD_2pos_, hfRingSumE_1neg_, hfRingSumE_1pos_, hfRingSumE_2neg_, hfRingSumE_2pos_, and mps_fire::i.

Referenced by analyze().

1468  {
1469  for (unsigned int i = 0; i < hfRingSumsD->size(); i++) {
1470  if (doGCTMBx_ || hfRingSumsD->at(i).bx() == GCTTrigBx_) {
1471  //there are 4 rings - just fill the histograms
1472  hfRingSumD_1pos_->Fill(hfRingSumsD->at(i).etSum(0));
1473  hfRingSumD_1neg_->Fill(hfRingSumsD->at(i).etSum(1));
1474  hfRingSumD_2pos_->Fill(hfRingSumsD->at(i).etSum(2));
1475  hfRingSumD_2neg_->Fill(hfRingSumsD->at(i).etSum(3));
1476  }
1477  }
1478 
1479  for (unsigned int i = 0; i < hfRingSumsE->size(); i++) {
1480  if (doEmuMBx_ || hfRingSumsE->at(i).bx() == EmuTrigBx_) {
1481  hfRingSumE_1pos_->Fill(hfRingSumsE->at(i).etSum(0));
1482  hfRingSumE_1neg_->Fill(hfRingSumsE->at(i).etSum(1));
1483  hfRingSumE_2pos_->Fill(hfRingSumsE->at(i).etSum(2));
1484  hfRingSumE_2neg_->Fill(hfRingSumsE->at(i).etSum(3));
1485  }
1486  }
1487 }
void GctErrorAnalyzer::plotIntJets ( const edm::Handle< L1GctInternJetDataCollection > &  emu)
private

Definition at line 1246 of file GctErrorAnalyzer.cc.

References doEmuMBx_, EmuTrigBx_, jetData::et, jetData::eta, mps_fire::i, intJetE_Et_, intJetE_Jet1Et_, intJetE_Jet2Et_, intJetE_Jet3Et_, intJetE_Jet4Et_, intJetE_Of_, intJetEtEtaPhiE_, dqmiolumiharvest::j, NUM_INT_JETS, jetData::phi, and sortJets().

Referenced by analyze().

1246  {
1247  jetData intJet;
1248  std::vector<jetData> intJetCollection(
1249  NUM_INT_JETS); //define fixed size for the vector to avoid reallocation (i.e. max size possible)
1250 
1251  //since we don't read out the intermediate (i.e. leaf card) jets, we can only plot the emulator distributions
1252  //the 1st-4th jet Et will prove useful in understanding and motivating cuts on individual jets in HT and MHT.
1253  for (unsigned int i = 0; i < intJetsE->size(); i++) {
1254  if (doEmuMBx_ || intJetsE->at(i).bx() == EmuTrigBx_) {
1255  //the intermediate jets are not sorted in terms of Et so
1256  //in order to do this independently of the data format,
1257  //copy to a user defined struct and sort that way
1258  intJet.et = intJetsE->at(i).et();
1259  intJet.phi = intJetsE->at(i).phi();
1260  intJet.eta = intJetsE->at(i).eta();
1261  intJetCollection.at(i % NUM_INT_JETS) = intJet;
1262 
1263  //remember, if the event has 1 overflowed jet, then we fill the internal jet dist overflow histogram
1264  //and skip the event - there is no point looking at the leading jet distributions etc for an event
1265  //with an overflowed jet - this will imply HT, ET, MET and MHT all overflow too.
1266  if (intJetsE->at(i).oflow()) {
1267  intJetE_Of_->Fill(intJetsE->at(i).oflow());
1268  return;
1269  }
1270 
1271  //plot the (et,eta,phi) distribution of the intermediate jets (for non-zero et)
1272  if (intJetsE->at(i).et())
1273  intJetEtEtaPhiE_->Fill(
1274  intJetsE->at(i).regionId().ieta(), intJetsE->at(i).regionId().iphi(), intJetsE->at(i).et());
1275  }
1276  }
1277 
1278  //if we get this far, there are no jets with an overflow bit so fill the overflow histogram and
1279  //sort the intJetCollection according to the rule defined in sortJets (i.e. largest et first)
1280  intJetE_Of_->Fill(0);
1281  std::sort(intJetCollection.begin(), intJetCollection.end(), sortJets);
1282 
1283  std::vector<TH1I *> leadingJetDist(4);
1284  leadingJetDist.at(0) = intJetE_Jet1Et_;
1285  leadingJetDist.at(1) = intJetE_Jet2Et_;
1286  leadingJetDist.at(2) = intJetE_Jet3Et_;
1287  leadingJetDist.at(3) = intJetE_Jet4Et_;
1288 
1289  unsigned int i = 0;
1290  unsigned int j = 0;
1291  unsigned int currentEt = 0;
1292  while (intJetCollection.at(i).et > 0) {
1293  if (j < leadingJetDist.size()) {
1294  if (i == 0) {
1295  leadingJetDist.at(j)->Fill(intJetCollection.at(i).et);
1296  currentEt = intJetCollection.at(i).et;
1297  j++;
1298  } else {
1299  if (intJetCollection.at(i).et < currentEt) {
1300  leadingJetDist.at(j)->Fill(intJetCollection.at(i).et);
1301  currentEt = intJetCollection.at(i).et;
1302  j++;
1303  }
1304  }
1305  }
1306 
1307  intJetE_Et_->Fill(intJetCollection.at(i).et);
1308  i++;
1309  }
1310  return;
1311 }
static bool sortJets(const jetData &jet1, const jetData &jet2)
const unsigned int NUM_INT_JETS
void GctErrorAnalyzer::plotIsoEm ( const edm::Handle< L1GctEmCandCollection > &  data,
const edm::Handle< L1GctEmCandCollection > &  emu 
)
private

Definition at line 1027 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, mps_fire::i, isoEgD_EtEtaPhi_, isoEgD_OccEtaPhi_, isoEgD_Rank_, isoEgE_EtEtaPhi_, isoEgE_OccEtaPhi_, and isoEgE_Rank_.

Referenced by analyze().

1028  {
1029  //loop over all the data candidates - if multiple bx, then this should be a multiple of GCT_OBJECT_QUANTA
1030  for (unsigned int i = 0; i < isoEgD->size(); i++) {
1031  //if the GCTMBx flag is set, plot all Bx for this quantity, otherwise only plot Bx = GCTTrigBx_
1032  if (doGCTMBx_ || isoEgD->at(i).bx() == GCTTrigBx_) {
1033  isoEgD_Rank_->Fill(isoEgD->at(i).rank());
1034  if (isoEgD->at(i).rank() > 0) {
1035  isoEgD_EtEtaPhi_->Fill(isoEgD->at(i).regionId().ieta(), isoEgD->at(i).regionId().iphi(), isoEgD->at(i).rank());
1036  isoEgD_OccEtaPhi_->Fill(isoEgD->at(i).regionId().ieta(), isoEgD->at(i).regionId().iphi());
1037  }
1038  }
1039  }
1040  //now repeat for the emulator candidates
1041  for (unsigned int i = 0; i < isoEgE->size(); i++) {
1042  if (doEmuMBx_ || isoEgE->at(i).bx() == EmuTrigBx_) {
1043  isoEgE_Rank_->Fill(isoEgE->at(i).rank());
1044  if (isoEgE->at(i).rank() > 0) {
1045  isoEgE_EtEtaPhi_->Fill(isoEgE->at(i).regionId().ieta(), isoEgE->at(i).regionId().iphi(), isoEgE->at(i).rank());
1046  isoEgE_OccEtaPhi_->Fill(isoEgE->at(i).regionId().ieta(), isoEgE->at(i).regionId().iphi());
1047  }
1048  }
1049  }
1050 }
void GctErrorAnalyzer::plotJetErrors ( const edm::Handle< L1GctJetCandCollection > &  cendata,
const edm::Handle< L1GctJetCandCollection > &  cenemu,
const edm::Handle< L1GctJetCandCollection > &  taudata,
const edm::Handle< L1GctJetCandCollection > &  tauemu,
const edm::Handle< L1GctJetCandCollection > &  fordata,
const edm::Handle< L1GctJetCandCollection > &  foremu,
const edm::Handle< L1CaloRegionCollection > &  regions 
)
private

Definition at line 1333 of file GctErrorAnalyzer.cc.

References EmuTrigBx_, errorHistCat, eventNumber, GCTTrigBx_, mps_fire::i, isCenJetError, isForJetError, isTauJetError, TFileDirectory::make(), RCTTrigBx_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

1339  {
1340  std::string errorDirName = "err_";
1341  if (isCenJetError)
1342  errorDirName.append("C");
1343  if (isTauJetError)
1344  errorDirName.append("T");
1345  if (isForJetError)
1346  errorDirName.append("F");
1347  std::stringstream caseNumber;
1348  caseNumber << eventNumber;
1349  errorDirName.append(caseNumber.str());
1350  TFileDirectory errorDir = errorHistCat.at(1).mkdir(errorDirName);
1351 
1352  TH2I *errorRegionEtEtaPhi_ = errorDir.make<TH2I>(
1353  "errorRegionEtEtaPhi", "errorRegionEtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1354  TH2I *errorRegionTvEtaPhi_ = errorDir.make<TH2I>(
1355  "errorRegionTvEtaPhi", "errorRegionTvEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1356  TH2I *errorRegionOfEtaPhi_ = errorDir.make<TH2I>(
1357  "errorRegionOfEtaPhi", "errorRegionOfEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1358 
1359  //make sure to only plot the caloRegion bx which corresponds to the data vs emulator comparison
1360  for (unsigned int i = 0; i < caloRegions->size(); i++) {
1361  if (caloRegions->at(i).bx() == RCTTrigBx_) {
1362  if (caloRegions->at(i).et() > 0)
1363  errorRegionEtEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi(), caloRegions->at(i).et());
1364  if (caloRegions->at(i).tauVeto())
1365  errorRegionTvEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1366  if (caloRegions->at(i).overFlow())
1367  errorRegionOfEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1368  }
1369  }
1370 
1371  TH2I *cenJet_errorEtEtaPhiData_ = errorDir.make<TH2I>("cenJet_errorEtEtaPhiData",
1372  "cenJet_errorEtEtaPhiData;#eta (GCT Units);#phi (GCT Units)",
1373  22,
1374  -0.5,
1375  21.5,
1376  18,
1377  -0.5,
1378  17.5);
1379  TH2I *cenJet_errorEtEtaPhiEmu_ = errorDir.make<TH2I>("cenJet_errorEtEtaPhiEmu",
1380  "cenJet_errorEtEtaPhiEmu;#eta (GCT Units);#phi (GCT Units)",
1381  22,
1382  -0.5,
1383  21.5,
1384  18,
1385  -0.5,
1386  17.5);
1387  TH2I *tauJet_errorEtEtaPhiData_ = errorDir.make<TH2I>("tauJet_errorEtEtaPhiData",
1388  "tauJet_errorEtEtaPhiData;#eta (GCT Units);#phi (GCT Units)",
1389  22,
1390  -0.5,
1391  21.5,
1392  18,
1393  -0.5,
1394  17.5);
1395  TH2I *tauJet_errorEtEtaPhiEmu_ = errorDir.make<TH2I>("tauJet_errorEtEtaPhiEmu",
1396  "tauJet_errorEtEtaPhiEmu;#eta (GCT Units);#phi (GCT Units)",
1397  22,
1398  -0.5,
1399  21.5,
1400  18,
1401  -0.5,
1402  17.5);
1403  TH2I *forJet_errorEtEtaPhiData_ = errorDir.make<TH2I>("forJet_errorEtEtaPhiData",
1404  "forJet_errorEtEtaPhiData;#eta (GCT Units);#phi (GCT Units)",
1405  22,
1406  -0.5,
1407  21.5,
1408  18,
1409  -0.5,
1410  17.5);
1411  TH2I *forJet_errorEtEtaPhiEmu_ = errorDir.make<TH2I>("forJet_errorEtEtaPhiEmu",
1412  "forJet_errorEtEtaPhiEmu;#eta (GCT Units);#phi (GCT Units)",
1413  22,
1414  -0.5,
1415  21.5,
1416  18,
1417  -0.5,
1418  17.5);
1419 
1420  //first plot the data candiates for the Trigger Bx that this error corresponds to
1421  for (unsigned int i = 0; i < cenJetsD->size(); i++) {
1422  if (cenJetsD->at(i).bx() == GCTTrigBx_) {
1423  if (cenJetsD->at(i).rank() > 0)
1424  cenJet_errorEtEtaPhiData_->Fill(
1425  cenJetsD->at(i).regionId().ieta(), cenJetsD->at(i).regionId().iphi(), cenJetsD->at(i).rank());
1426  }
1427  }
1428  for (unsigned int i = 0; i < tauJetsD->size(); i++) {
1429  if (tauJetsD->at(i).bx() == GCTTrigBx_) {
1430  if (tauJetsD->at(i).rank() > 0)
1431  tauJet_errorEtEtaPhiData_->Fill(
1432  tauJetsD->at(i).regionId().ieta(), tauJetsD->at(i).regionId().iphi(), tauJetsD->at(i).rank());
1433  }
1434  }
1435  for (unsigned int i = 0; i < forJetsD->size(); i++) {
1436  if (forJetsD->at(i).bx() == GCTTrigBx_) {
1437  if (forJetsD->at(i).rank() > 0)
1438  forJet_errorEtEtaPhiData_->Fill(
1439  forJetsD->at(i).regionId().ieta(), forJetsD->at(i).regionId().iphi(), forJetsD->at(i).rank());
1440  }
1441  }
1442 
1443  //now the emulator candidates
1444  for (unsigned int i = 0; i < cenJetsE->size(); i++) {
1445  if (cenJetsE->at(i).bx() == EmuTrigBx_) {
1446  if (cenJetsE->at(i).rank() > 0)
1447  cenJet_errorEtEtaPhiEmu_->Fill(
1448  cenJetsE->at(i).regionId().ieta(), cenJetsE->at(i).regionId().iphi(), cenJetsE->at(i).rank());
1449  }
1450  }
1451  for (unsigned int i = 0; i < tauJetsE->size(); i++) {
1452  if (tauJetsE->at(i).bx() == EmuTrigBx_) {
1453  if (tauJetsE->at(i).rank() > 0)
1454  tauJet_errorEtEtaPhiEmu_->Fill(
1455  tauJetsE->at(i).regionId().ieta(), tauJetsE->at(i).regionId().iphi(), tauJetsE->at(i).rank());
1456  }
1457  }
1458  for (unsigned int i = 0; i < forJetsE->size(); i++) {
1459  if (forJetsE->at(i).bx() == EmuTrigBx_) {
1460  if (forJetsE->at(i).rank() > 0)
1461  forJet_errorEtEtaPhiEmu_->Fill(
1462  forJetsE->at(i).regionId().ieta(), forJetsE->at(i).regionId().iphi(), forJetsE->at(i).rank());
1463  }
1464  }
1465 }
unsigned int eventNumber
T * make(const Args &...args) const
make new ROOT object
std::vector< TFileDirectory > errorHistCat
void GctErrorAnalyzer::plotMissingEErrors ( const edm::Handle< L1GctEtMissCollection > &  missingEtD,
const edm::Handle< L1GctEtMissCollection > &  missingEtE,
const edm::Handle< L1GctHtMissCollection > &  missingHtD,
const edm::Handle< L1GctHtMissCollection > &  missingHtE,
edm::Handle< L1CaloRegionCollection > &  caloRegions,
const edm::Handle< L1GctInternJetDataCollection > &  intjetsemu,
const edm::Handle< L1GctInternHtMissCollection intMissingHtD 
)
private

Definition at line 1771 of file GctErrorAnalyzer.cc.

References checkCollections(), doExtraMissingHTDebug_, EmuTrigBx_, errorHistCat, eventNumber, GCT_INT_HTMISS_QUANTA, GCTTrigBx_, mps_fire::i, isMissingEError, isMissingHError, dqmiolumiharvest::j, TFileDirectory::make(), NUM_INT_JETS, RCTTrigBx_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

1777  {
1778  std::string errorDirName = "err_";
1779  if (isMissingEError)
1780  errorDirName.append("E");
1781  if (isMissingHError)
1782  errorDirName.append("H");
1783 
1784  //added 05.03.2010 to highlight overflow errors in the missing energy sum calculation
1785  for (unsigned int i = 0; i < missingHtE->size(); i++) {
1786  if (missingHtE->at(i).bx() == EmuTrigBx_) {
1787  for (unsigned int j = 0; j < missingHtD->size(); j++) {
1788  if (missingHtD->at(j).bx() == GCTTrigBx_) {
1789  if (missingHtD->at(j).overFlow() != missingHtE->at(i).overFlow())
1790  errorDirName.append("O");
1791  }
1792  }
1793  }
1794  }
1795 
1796  std::stringstream caseNumber;
1797  caseNumber << eventNumber;
1798  errorDirName.append(caseNumber.str());
1799  TFileDirectory errorDir = errorHistCat.at(4).mkdir(errorDirName);
1800 
1801  TH2I *errorRegionEtEtaPhi_ = errorDir.make<TH2I>(
1802  "errorRegionEtEtaPhi", "errorRegionEtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1803  TH2I *errorRegionOfEtaPhi_ = errorDir.make<TH2I>(
1804  "errorRegionOfEtaPhi", "errorRegionOfEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1805  TH1I *errorMissingEtD_ =
1806  errorDir.make<TH1I>("errorMissingEtD", "errorMissingEtD;E_{T};Number of Events", 1024, -0.5, 1023.5);
1807  TH1I *errorMissingEtD_Of_ = errorDir.make<TH1I>(
1808  "errorMissingEtD_Of", "errorMissingEtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1809  TH1I *errorMissingEtD_Phi_ = errorDir.make<TH1I>(
1810  "errorMissingEtD_Phi", "errorMissingEtD_Phi;Missing E_{T} #phi;Number of Events", 72, -0.5, 71.5);
1811  TH1I *errorMissingEtE_ =
1812  errorDir.make<TH1I>("errorMissingEtE", "errorMissingEtE;E_{T};Number of Events", 1024, -0.5, 1023.5);
1813  TH1I *errorMissingEtE_Of_ = errorDir.make<TH1I>(
1814  "errorMissingEtE_Of", "errorMissingEtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1815  TH1I *errorMissingEtE_Phi_ = errorDir.make<TH1I>(
1816  "errorMissingEtE_Phi", "errorMissingEtE_Phi;Missing E_{T} #phi;Number of Events", 72, -0.5, 71.5);
1817  TH1I *errorMissingHtD_ =
1818  errorDir.make<TH1I>("errorMissingHtD", "errorMissingHtD;H_{T};Number of Events", 1024, -0.5, 1023.5);
1819  TH1I *errorMissingHtD_Of_ = errorDir.make<TH1I>(
1820  "errorMissingHtD_Of", "errorMissingHtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1821  TH1I *errorMissingHtD_Phi_ = errorDir.make<TH1I>(
1822  "errorMissingHtD_Phi", "errorMissingHtD_Phi;Missing H_{T} #phi;Number of Events", 72, -0.5, 71.5);
1823  TH1I *errorMissingHtE_ =
1824  errorDir.make<TH1I>("errorMissingHtE", "errorMissingHtE;H_{T};Number of Events", 1024, -0.5, 1023.5);
1825  TH1I *errorMissingHtE_Of_ = errorDir.make<TH1I>(
1826  "errorMissingHtE_Of", "errorMissingHtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1827  TH1I *errorMissingHtE_Phi_ = errorDir.make<TH1I>(
1828  "errorMissingHtE_Phi", "errorMissingHtE_Phi;Missing H_{T} #phi;Number of Events", 72, -0.5, 71.5);
1829 
1830  //Added 19.03.2010 to provide additional information in the case of missingHt failures
1831  //1. The MHT from both wheel inputs (i.e. the leaf cards)
1832  //2. The emulator jet et,eta,phi for all jets found in an event
1833  if (doExtraMissingHTDebug_) {
1834  if (checkCollections(intMissingHtD, GCT_INT_HTMISS_QUANTA, "Internal Missing Ht Data")) {
1835  TH1I *errorMissingHtD_HtXPosLeaf1 = errorDir.make<TH1I>(
1836  "errorMissingHtD_HtXPosLeaf1", "missingHtD;Missing H_{T} X PosLeaf1;Number of Events", 4096, -2048.5, 2047.5);
1837  TH1I *errorMissingHtD_HtXPosLeaf2 = errorDir.make<TH1I>(
1838  "errorMissingHtD_HtXPosLeaf2", "missingHtD;Missing H_{T} X PosLeaf2;Number of Events", 4096, -2048.5, 2047.5);
1839  TH1I *errorMissingHtD_HtXPosLeaf3 = errorDir.make<TH1I>(
1840  "errorMissingHtD_HtXPosLeaf3", "missingHtD;Missing H_{T} X PosLeaf3;Number of Events", 4096, -2048.5, 2047.5);
1841  TH1I *errorMissingHtD_HtXNegLeaf1 = errorDir.make<TH1I>(
1842  "errorMissingHtD_HtXNegLeaf1", "missingHtD;Missing H_{T} X NegLeaf1;Number of Events", 4096, -2048.5, 2047.5);
1843  TH1I *errorMissingHtD_HtXNegLeaf2 = errorDir.make<TH1I>(
1844  "errorMissingHtD_HtXNegLeaf2", "missingHtD;Missing H_{T} X NegLeaf2;Number of Events", 4096, -2048.5, 2047.5);
1845  TH1I *errorMissingHtD_HtXNegLeaf3 = errorDir.make<TH1I>(
1846  "errorMissingHtD_HtXNegLeaf3", "missingHtD;Missing H_{T} X NegLeaf3;Number of Events", 4096, -2048.5, 2047.5);
1847 
1848  TH1I *errorMissingHtD_HtYPosLeaf1 = errorDir.make<TH1I>(
1849  "errorMissingHtD_HtYPosLeaf1", "missingHtD;Missing H_{T} Y PosLeaf1;Number of Events", 4096, -2048.5, 2047.5);
1850  TH1I *errorMissingHtD_HtYPosLeaf2 = errorDir.make<TH1I>(
1851  "errorMissingHtD_HtYPosLeaf2", "missingHtD;Missing H_{T} Y PosLeaf2;Number of Events", 4096, -2048.5, 2047.5);
1852  TH1I *errorMissingHtD_HtYPosLeaf3 = errorDir.make<TH1I>(
1853  "errorMissingHtD_HtYPosLeaf3", "missingHtD;Missing H_{T} Y PosLeaf3;Number of Events", 4096, -2048.5, 2047.5);
1854  TH1I *errorMissingHtD_HtYNegLeaf1 = errorDir.make<TH1I>(
1855  "errorMissingHtD_HtYNegLeaf1", "missingHtD;Missing H_{T} Y NegLeaf1;Number of Events", 4096, -2048.5, 2047.5);
1856  TH1I *errorMissingHtD_HtYNegLeaf2 = errorDir.make<TH1I>(
1857  "errorMissingHtD_HtYNegLeaf2", "missingHtD;Missing H_{T} Y NegLeaf2;Number of Events", 4096, -2048.5, 2047.5);
1858  TH1I *errorMissingHtD_HtYNegLeaf3 = errorDir.make<TH1I>(
1859  "errorMissingHtD_HtYNegLeaf3", "missingHtD;Missing H_{T} Y NegLeaf3;Number of Events", 4096, -2048.5, 2047.5);
1860 
1861  for (unsigned int i = 0; i < intMissingHtD->size(); i++) {
1862  if (intMissingHtD->at(i).bx() == GCTTrigBx_) {
1863  if (!intMissingHtD->at(i).overflow()) {
1864  if (intMissingHtD->at(i).capBlock() == 0x301 && intMissingHtD->at(i).capIndex() == 0 &&
1865  intMissingHtD->at(i).isThereHtx())
1866  errorMissingHtD_HtXPosLeaf1->Fill(intMissingHtD->at(i).htx());
1867  if (intMissingHtD->at(i).capBlock() == 0x301 && intMissingHtD->at(i).capIndex() == 1 &&
1868  intMissingHtD->at(i).isThereHtx())
1869  errorMissingHtD_HtXPosLeaf2->Fill(intMissingHtD->at(i).htx());
1870  if (intMissingHtD->at(i).capBlock() == 0x301 && intMissingHtD->at(i).capIndex() == 2 &&
1871  intMissingHtD->at(i).isThereHtx())
1872  errorMissingHtD_HtXPosLeaf3->Fill(intMissingHtD->at(i).htx());
1873  if (intMissingHtD->at(i).capBlock() == 0x701 && intMissingHtD->at(i).capIndex() == 0 &&
1874  intMissingHtD->at(i).isThereHtx())
1875  errorMissingHtD_HtXNegLeaf1->Fill(intMissingHtD->at(i).htx());
1876  if (intMissingHtD->at(i).capBlock() == 0x701 && intMissingHtD->at(i).capIndex() == 1 &&
1877  intMissingHtD->at(i).isThereHtx())
1878  errorMissingHtD_HtXNegLeaf2->Fill(intMissingHtD->at(i).htx());
1879  if (intMissingHtD->at(i).capBlock() == 0x701 && intMissingHtD->at(i).capIndex() == 2 &&
1880  intMissingHtD->at(i).isThereHtx())
1881  errorMissingHtD_HtXNegLeaf3->Fill(intMissingHtD->at(i).htx());
1882 
1883  if (intMissingHtD->at(i).capBlock() == 0x301 && intMissingHtD->at(i).capIndex() == 0 &&
1884  intMissingHtD->at(i).isThereHty())
1885  errorMissingHtD_HtYPosLeaf1->Fill(intMissingHtD->at(i).hty());
1886  if (intMissingHtD->at(i).capBlock() == 0x301 && intMissingHtD->at(i).capIndex() == 1 &&
1887  intMissingHtD->at(i).isThereHty())
1888  errorMissingHtD_HtYPosLeaf2->Fill(intMissingHtD->at(i).hty());
1889  if (intMissingHtD->at(i).capBlock() == 0x301 && intMissingHtD->at(i).capIndex() == 2 &&
1890  intMissingHtD->at(i).isThereHty())
1891  errorMissingHtD_HtYPosLeaf3->Fill(intMissingHtD->at(i).hty());
1892  if (intMissingHtD->at(i).capBlock() == 0x701 && intMissingHtD->at(i).capIndex() == 0 &&
1893  intMissingHtD->at(i).isThereHty())
1894  errorMissingHtD_HtYNegLeaf1->Fill(intMissingHtD->at(i).hty());
1895  if (intMissingHtD->at(i).capBlock() == 0x701 && intMissingHtD->at(i).capIndex() == 1 &&
1896  intMissingHtD->at(i).isThereHty())
1897  errorMissingHtD_HtYNegLeaf2->Fill(intMissingHtD->at(i).hty());
1898  if (intMissingHtD->at(i).capBlock() == 0x701 && intMissingHtD->at(i).capIndex() == 2 &&
1899  intMissingHtD->at(i).isThereHty())
1900  errorMissingHtD_HtYNegLeaf3->Fill(intMissingHtD->at(i).hty());
1901  }
1902  }
1903  }
1904  }
1905  }
1906 
1907  if (checkCollections(intJetsE, NUM_INT_JETS, "Intermediate Jets Emulator")) {
1908  TH2I *errorIntJetsE_EtEtaPhi = errorDir.make<TH2I>("errorIntJetsE_EtEtaPhi",
1909  "errorIntJetsE_EtEtaPhi;#eta (GCT Units);#phi (GCT Units)",
1910  22,
1911  -0.5,
1912  21.5,
1913  18,
1914  -0.5,
1915  17.5);
1916 
1917  for (unsigned int i = 0; i < intJetsE->size(); i++) {
1918  if (intJetsE->at(i).bx() == EmuTrigBx_) {
1919  if (!intJetsE->at(i).oflow() && intJetsE->at(i).et())
1920  errorIntJetsE_EtEtaPhi->Fill(
1921  intJetsE->at(i).regionId().ieta(), intJetsE->at(i).regionId().iphi(), intJetsE->at(i).et());
1922  }
1923  }
1924  }
1925 
1926  for (unsigned int i = 0; i < caloRegions->size(); i++) {
1927  if (caloRegions->at(i).bx() == RCTTrigBx_) {
1928  if (caloRegions->at(i).et() > 0)
1929  errorRegionEtEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi(), caloRegions->at(i).et());
1930  if (caloRegions->at(i).overFlow())
1931  errorRegionOfEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1932  }
1933  }
1934 
1935  //plot the data candidates relating to this mismatch
1936  for (unsigned int i = 0; i < missingEtD->size(); i++) {
1937  if (missingEtD->at(i).bx() == GCTTrigBx_) {
1938  errorMissingEtD_Of_->Fill(missingEtD->at(i).overFlow());
1939  if (!missingEtD->at(i).overFlow()) {
1940  errorMissingEtD_->Fill(missingEtD->at(i).et());
1941  errorMissingEtD_Phi_->Fill(missingEtD->at(i).phi());
1942  }
1943  }
1944  }
1945  for (unsigned int i = 0; i < missingHtD->size(); i++) {
1946  if (missingHtD->at(i).bx() == GCTTrigBx_) {
1947  errorMissingHtD_Of_->Fill(missingHtD->at(i).overFlow());
1948  if (!missingHtD->at(i).overFlow()) {
1949  errorMissingHtD_->Fill(missingHtD->at(i).et());
1950  errorMissingHtD_Phi_->Fill(missingHtD->at(i).phi());
1951  }
1952  }
1953  }
1954  //and now for the emulator candidates
1955  for (unsigned int i = 0; i < missingEtE->size(); i++) {
1956  if (missingEtE->at(i).bx() == EmuTrigBx_) {
1957  errorMissingEtE_Of_->Fill(missingEtE->at(i).overFlow());
1958  if (!missingEtE->at(i).overFlow()) {
1959  errorMissingEtE_->Fill(missingEtE->at(i).et());
1960  errorMissingEtE_Phi_->Fill(missingEtE->at(i).phi());
1961  }
1962  }
1963  }
1964  for (unsigned int i = 0; i < missingHtE->size(); i++) {
1965  if (missingHtE->at(i).bx() == EmuTrigBx_) {
1966  errorMissingHtE_Of_->Fill(missingHtE->at(i).overFlow());
1967  if (!missingHtE->at(i)
1968  .overFlow()) { //to see what values the emulator outputs despite the o/f bit being set comment out this statement
1969  errorMissingHtE_->Fill(missingHtE->at(i).et());
1970  errorMissingHtE_Phi_->Fill(missingHtE->at(i).phi());
1971  }
1972  }
1973  }
1974 }
unsigned int eventNumber
T * make(const Args &...args) const
make new ROOT object
bool checkCollections(const T &collection, const unsigned int &constraint, const std::string &label)
std::vector< TFileDirectory > errorHistCat
const unsigned int NUM_INT_JETS
const unsigned int GCT_INT_HTMISS_QUANTA
void GctErrorAnalyzer::plotMissingEt ( const edm::Handle< L1GctEtMissCollection > &  missingEtD,
const edm::Handle< L1GctEtMissCollection > &  missingEtE 
)
private

Definition at line 1725 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, mps_fire::i, missingEtD_, missingEtD_Of_, missingEtD_Phi_, missingEtE_, missingEtE_Of_, and missingEtE_Phi_.

Referenced by analyze().

1726  {
1727  for (unsigned int i = 0; i < missingEtD->size(); i++) {
1728  if (doGCTMBx_ || missingEtD->at(i).bx() == GCTTrigBx_) {
1729  missingEtD_Of_->Fill(missingEtD->at(i).overFlow());
1730  if (!missingEtD->at(i).overFlow() && missingEtD->at(i).et() > 0) {
1731  missingEtD_->Fill(missingEtD->at(i).et());
1732  missingEtD_Phi_->Fill(missingEtD->at(i).phi());
1733  }
1734  }
1735  }
1736 
1737  for (unsigned int i = 0; i < missingEtE->size(); i++) {
1738  if (doEmuMBx_ || missingEtE->at(i).bx() == EmuTrigBx_) {
1739  missingEtE_Of_->Fill(missingEtE->at(i).overFlow());
1740  if (!missingEtE->at(i).overFlow() && missingEtE->at(i).et()) {
1741  missingEtE_->Fill(missingEtE->at(i).et());
1742  missingEtE_Phi_->Fill(missingEtE->at(i).phi());
1743  }
1744  }
1745  }
1746 }
void GctErrorAnalyzer::plotMissingHt ( const edm::Handle< L1GctHtMissCollection > &  missingHtD,
const edm::Handle< L1GctHtMissCollection > &  missingHtE 
)
private

Definition at line 1748 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, mps_fire::i, missingHtD_, missingHtD_Of_, missingHtD_Phi_, missingHtE_, missingHtE_Of_, and missingHtE_Phi_.

Referenced by analyze().

1749  {
1750  for (unsigned int i = 0; i < missingHtD->size(); i++) {
1751  if (doGCTMBx_ || missingHtD->at(i).bx() == GCTTrigBx_) {
1752  missingHtD_Of_->Fill(missingHtD->at(i).overFlow());
1753  if (!missingHtD->at(i).overFlow() && missingHtD->at(i).et() > 0) {
1754  missingHtD_->Fill(missingHtD->at(i).et());
1755  missingHtD_Phi_->Fill(missingHtD->at(i).phi());
1756  }
1757  }
1758  }
1759 
1760  for (unsigned int i = 0; i < missingHtE->size(); i++) {
1761  if (doEmuMBx_ || missingHtE->at(i).bx() == EmuTrigBx_) {
1762  missingHtE_Of_->Fill(missingHtE->at(i).overFlow());
1763  if (!missingHtE->at(i).overFlow() && missingHtE->at(i).et() > 0) {
1764  missingHtE_->Fill(missingHtE->at(i).et());
1765  missingHtE_Phi_->Fill(missingHtE->at(i).phi());
1766  }
1767  }
1768  }
1769 }
void GctErrorAnalyzer::plotNonIsoEm ( const edm::Handle< L1GctEmCandCollection > &  data,
const edm::Handle< L1GctEmCandCollection > &  emu 
)
private

Definition at line 1052 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, mps_fire::i, nonIsoEgD_EtEtaPhi_, nonIsoEgD_OccEtaPhi_, nonIsoEgD_Rank_, nonIsoEgE_EtEtaPhi_, nonIsoEgE_OccEtaPhi_, and nonIsoEgE_Rank_.

Referenced by analyze().

1053  {
1054  //loop over all the data candidates - if multiple bx, then this should be a multiple of GCT_OBJECT_QUANTA
1055  for (unsigned int i = 0; i < nonIsoEgD->size(); i++) {
1056  //if the GCTMBx flag is set, plot all Bx for this quantity, otherwise only plot Bx = GCTTrigBx_
1057  if (doGCTMBx_ || nonIsoEgD->at(i).bx() == GCTTrigBx_) {
1058  nonIsoEgD_Rank_->Fill(nonIsoEgD->at(i).rank());
1059  if (nonIsoEgD->at(i).rank() > 0) {
1060  nonIsoEgD_EtEtaPhi_->Fill(
1061  nonIsoEgD->at(i).regionId().ieta(), nonIsoEgD->at(i).regionId().iphi(), nonIsoEgD->at(i).rank());
1062  nonIsoEgD_OccEtaPhi_->Fill(nonIsoEgD->at(i).regionId().ieta(), nonIsoEgD->at(i).regionId().iphi());
1063  }
1064  }
1065  }
1066  //now repeat for the emulator candidates
1067  for (unsigned int i = 0; i < nonIsoEgE->size(); i++) {
1068  if (doEmuMBx_ || nonIsoEgE->at(i).bx() == EmuTrigBx_) {
1069  nonIsoEgE_Rank_->Fill(nonIsoEgE->at(i).rank());
1070  if (nonIsoEgE->at(i).rank() > 0) {
1071  nonIsoEgE_EtEtaPhi_->Fill(
1072  nonIsoEgE->at(i).regionId().ieta(), nonIsoEgE->at(i).regionId().iphi(), nonIsoEgE->at(i).rank());
1073  nonIsoEgE_OccEtaPhi_->Fill(nonIsoEgE->at(i).regionId().ieta(), nonIsoEgE->at(i).regionId().iphi());
1074  }
1075  }
1076  }
1077 }
void GctErrorAnalyzer::plotRCTRegions ( const edm::Handle< L1CaloRegionCollection > &  caloRegions)
private

Definition at line 1009 of file GctErrorAnalyzer.cc.

References doRCTMBx_, mps_fire::i, RCT_EtEtaPhi_, RCT_FgEtaPhi_, RCT_OfEtaPhi_, RCT_TvEtaPhi_, and RCTTrigBx_.

Referenced by analyze().

1009  {
1010  //if more than one Bx is readout per event, then caloRegions->size() will be some multiple of 396
1011  for (unsigned int i = 0; i < caloRegions->size(); i++) {
1012  //if the RCTMBx flag is set to true, write out all the info into the same histogram
1013  //otherwise only the RCTTrigBx will be written out - could skip (RCT_REGION_QUANTA-1) events here to speed things up...
1014  if (doRCTMBx_ || caloRegions->at(i).bx() == RCTTrigBx_) {
1015  if (caloRegions->at(i).et() > 0)
1016  RCT_EtEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi(), caloRegions->at(i).et());
1017  if (caloRegions->at(i).tauVeto())
1018  RCT_TvEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1019  if (caloRegions->at(i).fineGrain())
1020  RCT_FgEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1021  if (caloRegions->at(i).overFlow())
1022  RCT_OfEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1023  }
1024  }
1025 }
void GctErrorAnalyzer::plotTauJets ( const edm::Handle< L1GctJetCandCollection > &  data,
const edm::Handle< L1GctJetCandCollection > &  emu 
)
private

Definition at line 1196 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, mps_fire::i, tauJetD_EtEtaPhi_, tauJetD_OccEtaPhi_, tauJetD_Rank_, tauJetE_EtEtaPhi_, tauJetE_OccEtaPhi_, and tauJetE_Rank_.

Referenced by analyze().

1197  {
1198  for (unsigned int i = 0; i < tauJetsD->size(); i++) {
1199  if (doGCTMBx_ || tauJetsD->at(i).bx() == GCTTrigBx_) {
1200  tauJetD_Rank_->Fill(tauJetsD->at(i).rank());
1201  if (tauJetsD->at(i).rank() > 0) {
1202  tauJetD_EtEtaPhi_->Fill(
1203  tauJetsD->at(i).regionId().ieta(), tauJetsD->at(i).regionId().iphi(), tauJetsD->at(i).rank());
1204  tauJetD_OccEtaPhi_->Fill(tauJetsD->at(i).regionId().ieta(), tauJetsD->at(i).regionId().iphi());
1205  }
1206  }
1207  }
1208 
1209  for (unsigned int i = 0; i < tauJetsE->size(); i++) {
1210  if (doEmuMBx_ || tauJetsE->at(i).bx() == EmuTrigBx_) {
1211  tauJetE_Rank_->Fill(tauJetsE->at(i).rank());
1212  if (tauJetsE->at(i).rank() > 0) {
1213  tauJetE_EtEtaPhi_->Fill(
1214  tauJetsE->at(i).regionId().ieta(), tauJetsE->at(i).regionId().iphi(), tauJetsE->at(i).rank());
1215  tauJetE_OccEtaPhi_->Fill(tauJetsE->at(i).regionId().ieta(), tauJetsE->at(i).regionId().iphi());
1216  }
1217  }
1218  }
1219 }
void GctErrorAnalyzer::plotTotalE ( const edm::Handle< L1GctEtTotalCollection > &  totalEtD,
const edm::Handle< L1GctEtTotalCollection > &  totalEtE 
)
private

Definition at line 1612 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, mps_fire::i, totalEtD_, totalEtD_Of_, totalEtE_, and totalEtE_Of_.

Referenced by analyze().

1613  {
1614  for (unsigned int i = 0; i < totalEtD->size(); i++) {
1615  if (doGCTMBx_ || totalEtD->at(i).bx() == GCTTrigBx_) {
1616  totalEtD_Of_->Fill(totalEtD->at(i).overFlow());
1617  if (!totalEtD->at(i).overFlow())
1618  totalEtD_->Fill(totalEtD->at(i).et());
1619  }
1620  }
1621  for (unsigned int i = 0; i < totalEtE->size(); i++) {
1622  if (doEmuMBx_ || totalEtE->at(i).bx() == EmuTrigBx_) {
1623  totalEtE_Of_->Fill(totalEtE->at(i).overFlow());
1624  if (!totalEtE->at(i).overFlow())
1625  totalEtE_->Fill(totalEtE->at(i).et());
1626  }
1627  }
1628 }
void GctErrorAnalyzer::plotTotalEErrors ( const edm::Handle< L1GctEtTotalCollection > &  totalEtD,
const edm::Handle< L1GctEtTotalCollection > &  totalEtE,
const edm::Handle< L1GctEtHadCollection > &  totalHtD,
const edm::Handle< L1GctEtHadCollection > &  totalHtE,
const edm::Handle< L1CaloRegionCollection > &  caloRegions 
)
private

Definition at line 1648 of file GctErrorAnalyzer.cc.

References EmuTrigBx_, errorHistCat, eventNumber, GCTTrigBx_, mps_fire::i, isTotalEError, isTotalHError, TFileDirectory::make(), RCTTrigBx_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

1652  {
1653  std::string errorDirName = "err_";
1654  if (isTotalEError)
1655  errorDirName.append("E");
1656  if (isTotalHError)
1657  errorDirName.append("H");
1658  std::stringstream caseNumber;
1659  caseNumber << eventNumber;
1660  errorDirName.append(caseNumber.str());
1661  TFileDirectory errorDir = errorHistCat.at(3).mkdir(errorDirName);
1662 
1663  TH2I *errorRegionEtEtaPhi_ = errorDir.make<TH2I>(
1664  "errorRegionEtEtaPhi", "errorRegionEtEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1665  TH2I *errorRegionOfEtaPhi_ = errorDir.make<TH2I>(
1666  "errorRegionOfEtaPhi", "errorRegionOfEtaPhi;#eta (GCT Units);#phi (GCT Units)", 22, -0.5, 21.5, 18, -0.5, 17.5);
1667  TH1I *errorTotalEtD_ =
1668  errorDir.make<TH1I>("errorTotalEtD", "errorTotalEtD;E_{T};Number of Events", 1024, -0.5, 1023.5);
1669  TH1I *errorTotalEtD_Of_ =
1670  errorDir.make<TH1I>("errorTotalEtD_Of", "errorTotalEtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1671  TH1I *errorTotalEtE_ =
1672  errorDir.make<TH1I>("errorTotalEtE", "errorTotalEtE;E_{T};Number of Events", 1024, -0.5, 1023.5);
1673  TH1I *errorTotalEtE_Of_ =
1674  errorDir.make<TH1I>("errorTotalEtE_Of", "errorTotalEtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1675  TH1I *errorTotalHtD_ =
1676  errorDir.make<TH1I>("errorTotalHtD", "errorTotalHtD;E_{T};Number of Events", 1024, -0.5, 1023.5);
1677  TH1I *errorTotalHtD_Of_ =
1678  errorDir.make<TH1I>("errorTotalHtD_Of", "errorTotalHtD_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1679  TH1I *errorTotalHtE_ =
1680  errorDir.make<TH1I>("errorTotalHtE", "errorTotalHtE;E_{T};Number of Events", 1024, -0.5, 1023.5);
1681  TH1I *errorTotalHtE_Of_ =
1682  errorDir.make<TH1I>("errorTotalHtE_Of", "errorTotalHtE_Of;Overflow Bit Status;Number of Events", 2, -0.5, 1.5);
1683 
1684  //plot the region ET and OF bits
1685  for (unsigned int i = 0; i < caloRegions->size(); i++) {
1686  if (caloRegions->at(i).bx() == RCTTrigBx_) {
1687  if (caloRegions->at(i).et() > 0)
1688  errorRegionEtEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi(), caloRegions->at(i).et());
1689  if (caloRegions->at(i).overFlow())
1690  errorRegionOfEtaPhi_->Fill(caloRegions->at(i).gctEta(), caloRegions->at(i).gctPhi());
1691  }
1692  }
1693  //now plot the error ET
1694  for (unsigned int i = 0; i < totalEtD->size(); i++) {
1695  if (totalEtD->at(i).bx() == GCTTrigBx_) {
1696  errorTotalEtD_Of_->Fill(totalEtD->at(i).overFlow());
1697  if (!totalEtD->at(i).overFlow())
1698  errorTotalEtD_->Fill(totalEtD->at(i).et());
1699  }
1700  }
1701  for (unsigned int i = 0; i < totalEtE->size(); i++) {
1702  if (totalEtE->at(i).bx() == EmuTrigBx_) {
1703  errorTotalEtE_Of_->Fill(totalEtE->at(i).overFlow());
1704  if (!totalEtE->at(i).overFlow())
1705  errorTotalEtE_->Fill(totalEtE->at(i).et());
1706  }
1707  }
1708  //and now the error HT
1709  for (unsigned int i = 0; i < totalHtD->size(); i++) {
1710  if (totalHtD->at(i).bx() == GCTTrigBx_) {
1711  errorTotalHtD_Of_->Fill(totalHtD->at(i).overFlow());
1712  if (!totalHtD->at(i).overFlow())
1713  errorTotalHtD_->Fill(totalHtD->at(i).et());
1714  }
1715  }
1716  for (unsigned int i = 0; i < totalHtE->size(); i++) {
1717  if (totalHtE->at(i).bx() == EmuTrigBx_) {
1718  errorTotalHtE_Of_->Fill(totalHtE->at(i).overFlow());
1719  if (!totalHtE->at(i).overFlow())
1720  errorTotalHtE_->Fill(totalHtE->at(i).et());
1721  }
1722  }
1723 }
unsigned int eventNumber
T * make(const Args &...args) const
make new ROOT object
std::vector< TFileDirectory > errorHistCat
void GctErrorAnalyzer::plotTotalH ( const edm::Handle< L1GctEtHadCollection > &  totalHtD,
const edm::Handle< L1GctEtHadCollection > &  totalHtE 
)
private

Definition at line 1630 of file GctErrorAnalyzer.cc.

References doEmuMBx_, doGCTMBx_, EmuTrigBx_, GCTTrigBx_, mps_fire::i, totalHtD_, totalHtD_Of_, totalHtE_, and totalHtE_Of_.

Referenced by analyze().

1631  {
1632  for (unsigned int i = 0; i < totalHtD->size(); i++) {
1633  if (doGCTMBx_ || totalHtD->at(i).bx() == GCTTrigBx_) {
1634  totalHtD_Of_->Fill(totalHtD->at(i).overFlow());
1635  if (!totalHtD->at(i).overFlow())
1636  totalHtD_->Fill(totalHtD->at(i).et());
1637  }
1638  }
1639  for (unsigned int i = 0; i < totalHtE->size(); i++) {
1640  if (doEmuMBx_ || totalHtE->at(i).bx() == EmuTrigBx_) {
1641  totalHtE_Of_->Fill(totalHtE->at(i).overFlow());
1642  if (!totalHtE->at(i).overFlow())
1643  totalHtE_->Fill(totalHtE->at(i).et());
1644  }
1645  }
1646 }
bool GctErrorAnalyzer::sortJets ( const jetData jet1,
const jetData jet2 
)
staticprivate

Definition at line 1313 of file GctErrorAnalyzer.cc.

References jetData::et.

Referenced by plotIntJets().

1313 { return jet1.et > jet2.et; }

Member Data Documentation

TH1I* GctErrorAnalyzer::cenJet_errorFlag_
private

Definition at line 192 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::cenJetD_EtEtaPhi_
private

Definition at line 190 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotCenJets().

TH2I* GctErrorAnalyzer::cenJetD_GlobalError_EtEtaPhi_
private

Definition at line 196 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::cenJetD_GlobalError_Rank_
private

Definition at line 194 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::cenJetD_OccEtaPhi_
private

Definition at line 191 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotCenJets().

TH1I* GctErrorAnalyzer::cenJetD_Rank_
private

Definition at line 189 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotCenJets().

TH2I * GctErrorAnalyzer::cenJetE_EtEtaPhi_
private

Definition at line 190 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotCenJets().

TH2I* GctErrorAnalyzer::cenJetE_GlobalError_EtEtaPhi_
private

Definition at line 197 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::cenJetE_GlobalError_Rank_
private

Definition at line 195 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I * GctErrorAnalyzer::cenJetE_OccEtaPhi_
private

Definition at line 191 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotCenJets().

TH1I * GctErrorAnalyzer::cenJetE_Rank_
private

Definition at line 189 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotCenJets().

edm::InputTag GctErrorAnalyzer::dataTag_
private

Definition at line 157 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doBitCountDebug_
private

Definition at line 137 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doCenJetsDebug_
private

Definition at line 132 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doEg_
private

Definition at line 128 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doEmuMBx_
private
bool GctErrorAnalyzer::doExtraMissingHTDebug_
private

Definition at line 144 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and plotMissingEErrors().

bool GctErrorAnalyzer::doForJetsDebug_
private

Definition at line 134 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doGCTMBx_
private
bool GctErrorAnalyzer::doHF_
private

Definition at line 135 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doIsoDebug_
private

Definition at line 129 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doJets_
private

Definition at line 131 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doMissingEnergySums_
private

Definition at line 141 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doMissingETDebug_
private

Definition at line 142 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doMissingHTDebug_
private

Definition at line 143 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doNonIsoDebug_
private

Definition at line 130 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doRCT_
private

Definition at line 127 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doRCTMBx_
private

Definition at line 149 of file GctErrorAnalyzer.cc.

Referenced by plotRCTRegions().

bool GctErrorAnalyzer::doRingSumDebug_
private

Definition at line 136 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doTauJetsDebug_
private

Definition at line 133 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doTotalEnergySums_
private

Definition at line 138 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doTotalEtDebug_
private

Definition at line 139 of file GctErrorAnalyzer.cc.

Referenced by analyze().

bool GctErrorAnalyzer::doTotalHtDebug_
private

Definition at line 140 of file GctErrorAnalyzer.cc.

Referenced by analyze().

edm::InputTag GctErrorAnalyzer::emuTag_
private

Definition at line 158 of file GctErrorAnalyzer.cc.

Referenced by analyze().

int GctErrorAnalyzer::EmuTrigBx_
private
std::vector<TFileDirectory> GctErrorAnalyzer::errorHistCat
private
unsigned int GctErrorAnalyzer::eventNumber
private
TH1I* GctErrorAnalyzer::forJet_errorFlag_
private

Definition at line 212 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::forJetD_EtEtaPhi_
private

Definition at line 210 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotForJets().

TH2I* GctErrorAnalyzer::forJetD_GlobalError_EtEtaPhi_
private

Definition at line 216 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::forJetD_GlobalError_Rank_
private

Definition at line 214 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::forJetD_OccEtaPhi_
private

Definition at line 211 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotForJets().

TH1I* GctErrorAnalyzer::forJetD_Rank_
private

Definition at line 209 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotForJets().

TH2I * GctErrorAnalyzer::forJetE_EtEtaPhi_
private

Definition at line 210 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotForJets().

TH2I* GctErrorAnalyzer::forJetE_GlobalError_EtEtaPhi_
private

Definition at line 217 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::forJetE_GlobalError_Rank_
private

Definition at line 215 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I * GctErrorAnalyzer::forJetE_OccEtaPhi_
private

Definition at line 211 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotForJets().

TH1I * GctErrorAnalyzer::forJetE_Rank_
private

Definition at line 209 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotForJets().

int GctErrorAnalyzer::GCTTrigBx_
private
TH1I* GctErrorAnalyzer::hfBitCount_errorFlag_
private

Definition at line 233 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::hfBitCountD_1neg_
private

Definition at line 231 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I* GctErrorAnalyzer::hfBitCountD_1pos_
private

Definition at line 231 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I * GctErrorAnalyzer::hfBitCountD_2neg_
private

Definition at line 231 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I * GctErrorAnalyzer::hfBitCountD_2pos_
private

Definition at line 231 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I * GctErrorAnalyzer::hfBitCountE_1neg_
private

Definition at line 232 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I* GctErrorAnalyzer::hfBitCountE_1pos_
private

Definition at line 232 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I * GctErrorAnalyzer::hfBitCountE_2neg_
private

Definition at line 232 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I * GctErrorAnalyzer::hfBitCountE_2pos_
private

Definition at line 232 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFBitCounts().

TH1I* GctErrorAnalyzer::hfRingSum_errorFlag_
private

Definition at line 229 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::hfRingSumD_1neg_
private

Definition at line 227 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I* GctErrorAnalyzer::hfRingSumD_1pos_
private

Definition at line 227 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I * GctErrorAnalyzer::hfRingSumD_2neg_
private

Definition at line 227 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I * GctErrorAnalyzer::hfRingSumD_2pos_
private

Definition at line 227 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I * GctErrorAnalyzer::hfRingSumE_1neg_
private

Definition at line 228 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I* GctErrorAnalyzer::hfRingSumE_1pos_
private

Definition at line 228 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I * GctErrorAnalyzer::hfRingSumE_2neg_
private

Definition at line 228 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I * GctErrorAnalyzer::hfRingSumE_2pos_
private

Definition at line 228 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotHFRingSums().

TH1I* GctErrorAnalyzer::intJetE_Et_
private

Definition at line 220 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIntJets().

TH1I* GctErrorAnalyzer::intJetE_Jet1Et_
private

Definition at line 222 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIntJets().

TH1I* GctErrorAnalyzer::intJetE_Jet2Et_
private

Definition at line 223 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIntJets().

TH1I* GctErrorAnalyzer::intJetE_Jet3Et_
private

Definition at line 224 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIntJets().

TH1I* GctErrorAnalyzer::intJetE_Jet4Et_
private

Definition at line 225 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIntJets().

TH1I* GctErrorAnalyzer::intJetE_Of_
private

Definition at line 221 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIntJets().

TH2I* GctErrorAnalyzer::intJetEtEtaPhiE_
private

Definition at line 219 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIntJets().

bool GctErrorAnalyzer::isBitCountError
private

Definition at line 270 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotHFErrors().

bool GctErrorAnalyzer::isCenJetError
private

Definition at line 266 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotJetErrors().

bool GctErrorAnalyzer::isForJetError
private

Definition at line 268 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotJetErrors().

bool GctErrorAnalyzer::isIsoError
private

Definition at line 264 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotEGErrors().

bool GctErrorAnalyzer::isMissingEError
private

Definition at line 273 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotMissingEErrors().

bool GctErrorAnalyzer::isMissingHError
private

Definition at line 274 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotMissingEErrors().

bool GctErrorAnalyzer::isNonIsoError
private

Definition at line 265 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotEGErrors().

TH1I* GctErrorAnalyzer::isoEg_errorFlag_
private

Definition at line 172 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::isoEgD_EtEtaPhi_
private

Definition at line 170 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIsoEm().

TH2I* GctErrorAnalyzer::isoEgD_GlobalError_EtEtaPhi_
private

Definition at line 176 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::isoEgD_GlobalError_Rank_
private

Definition at line 174 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::isoEgD_OccEtaPhi_
private

Definition at line 171 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIsoEm().

TH1I* GctErrorAnalyzer::isoEgD_Rank_
private

Definition at line 169 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIsoEm().

TH2I * GctErrorAnalyzer::isoEgE_EtEtaPhi_
private

Definition at line 170 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIsoEm().

TH2I* GctErrorAnalyzer::isoEgE_GlobalError_EtEtaPhi_
private

Definition at line 177 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::isoEgE_GlobalError_Rank_
private

Definition at line 175 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I * GctErrorAnalyzer::isoEgE_OccEtaPhi_
private

Definition at line 171 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIsoEm().

TH1I * GctErrorAnalyzer::isoEgE_Rank_
private

Definition at line 169 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotIsoEm().

bool GctErrorAnalyzer::isRingSumError
private

Definition at line 269 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotHFErrors().

bool GctErrorAnalyzer::isTauJetError
private

Definition at line 267 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotJetErrors().

bool GctErrorAnalyzer::isTotalEError
private

Definition at line 271 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotTotalEErrors().

bool GctErrorAnalyzer::isTotalHError
private

Definition at line 272 of file GctErrorAnalyzer.cc.

Referenced by analyze(), GctErrorAnalyzer(), and plotTotalEErrors().

GctErrorAnalyzerMBxInfo GctErrorAnalyzer::MBxInfo
private

Definition at line 163 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::missingEt_errorFlag_
private

Definition at line 252 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::missingEtD_
private

Definition at line 249 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingEt().

TH1I* GctErrorAnalyzer::missingEtD_Of_
private

Definition at line 250 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingEt().

TH1I* GctErrorAnalyzer::missingEtD_Phi_
private

Definition at line 251 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingEt().

TH1I * GctErrorAnalyzer::missingEtE_
private

Definition at line 249 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingEt().

TH1I * GctErrorAnalyzer::missingEtE_Of_
private

Definition at line 250 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingEt().

TH1I * GctErrorAnalyzer::missingEtE_Phi_
private

Definition at line 251 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingEt().

TH1I* GctErrorAnalyzer::missingHt_errorFlag_
private

Definition at line 257 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::missingHtD_
private

Definition at line 254 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingHt().

TH1I * GctErrorAnalyzer::missingHtD_HtXNegLeaf1
private

Definition at line 258 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtXNegLeaf2
private

Definition at line 258 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtXNegLeaf3
private

Definition at line 258 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::missingHtD_HtXPosLeaf1
private

Definition at line 258 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtXPosLeaf2
private

Definition at line 258 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtXPosLeaf3
private

Definition at line 258 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtYNegLeaf1
private

Definition at line 260 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtYNegLeaf2
private

Definition at line 260 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtYNegLeaf3
private

Definition at line 260 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::missingHtD_HtYPosLeaf1
private

Definition at line 260 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtYPosLeaf2
private

Definition at line 260 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I * GctErrorAnalyzer::missingHtD_HtYPosLeaf3
private

Definition at line 260 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::missingHtD_Of_
private

Definition at line 255 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingHt().

TH1I* GctErrorAnalyzer::missingHtD_Phi_
private

Definition at line 256 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingHt().

TH1I * GctErrorAnalyzer::missingHtE_
private

Definition at line 254 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingHt().

TH1I * GctErrorAnalyzer::missingHtE_Of_
private

Definition at line 255 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingHt().

TH1I * GctErrorAnalyzer::missingHtE_Phi_
private

Definition at line 256 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotMissingHt().

TH1I* GctErrorAnalyzer::nonIsoEg_errorFlag_
private

Definition at line 182 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::nonIsoEgD_EtEtaPhi_
private

Definition at line 180 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotNonIsoEm().

TH2I* GctErrorAnalyzer::nonIsoEgD_GlobalError_EtEtaPhi_
private

Definition at line 186 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::nonIsoEgD_GlobalError_Rank_
private

Definition at line 184 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::nonIsoEgD_OccEtaPhi_
private

Definition at line 181 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotNonIsoEm().

TH1I* GctErrorAnalyzer::nonIsoEgD_Rank_
private

Definition at line 179 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotNonIsoEm().

TH2I * GctErrorAnalyzer::nonIsoEgE_EtEtaPhi_
private

Definition at line 180 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotNonIsoEm().

TH2I* GctErrorAnalyzer::nonIsoEgE_GlobalError_EtEtaPhi_
private

Definition at line 187 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::nonIsoEgE_GlobalError_Rank_
private

Definition at line 185 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I * GctErrorAnalyzer::nonIsoEgE_OccEtaPhi_
private

Definition at line 181 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotNonIsoEm().

TH1I * GctErrorAnalyzer::nonIsoEgE_Rank_
private

Definition at line 179 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotNonIsoEm().

TH2I* GctErrorAnalyzer::RCT_EtEtaPhi_
private

Definition at line 167 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotRCTRegions().

TH2I * GctErrorAnalyzer::RCT_FgEtaPhi_
private

Definition at line 167 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotRCTRegions().

TH2I * GctErrorAnalyzer::RCT_OfEtaPhi_
private

Definition at line 167 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotRCTRegions().

const unsigned int* GctErrorAnalyzer::RCT_REGION_QUANTA
private

Definition at line 284 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I * GctErrorAnalyzer::RCT_TvEtaPhi_
private

Definition at line 167 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotRCTRegions().

int GctErrorAnalyzer::RCTTrigBx_
private
TH1I* GctErrorAnalyzer::tauJet_errorFlag_
private

Definition at line 202 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::tauJetD_EtEtaPhi_
private

Definition at line 200 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTauJets().

TH2I* GctErrorAnalyzer::tauJetD_GlobalError_EtEtaPhi_
private

Definition at line 206 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::tauJetD_GlobalError_Rank_
private

Definition at line 204 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I* GctErrorAnalyzer::tauJetD_OccEtaPhi_
private

Definition at line 201 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTauJets().

TH1I* GctErrorAnalyzer::tauJetD_Rank_
private

Definition at line 199 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTauJets().

TH2I * GctErrorAnalyzer::tauJetE_EtEtaPhi_
private

Definition at line 200 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTauJets().

TH2I* GctErrorAnalyzer::tauJetE_GlobalError_EtEtaPhi_
private

Definition at line 207 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::tauJetE_GlobalError_Rank_
private

Definition at line 205 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH2I * GctErrorAnalyzer::tauJetE_OccEtaPhi_
private

Definition at line 201 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTauJets().

TH1I * GctErrorAnalyzer::tauJetE_Rank_
private

Definition at line 199 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTauJets().

TH1I* GctErrorAnalyzer::totalEt_errorFlag_
private

Definition at line 237 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::totalEtD_
private

Definition at line 235 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalE().

TH1I* GctErrorAnalyzer::totalEtD_Of_
private

Definition at line 236 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalE().

TH1I * GctErrorAnalyzer::totalEtE_
private

Definition at line 235 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalE().

TH1I * GctErrorAnalyzer::totalEtE_Of_
private

Definition at line 236 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalE().

TH1I* GctErrorAnalyzer::totalHt_errorFlag_
private

Definition at line 244 of file GctErrorAnalyzer.cc.

Referenced by analyze(), and GctErrorAnalyzer().

TH1I* GctErrorAnalyzer::totalHtD_
private

Definition at line 242 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalH().

TH1I* GctErrorAnalyzer::totalHtD_Of_
private

Definition at line 243 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalH().

TH1I * GctErrorAnalyzer::totalHtE_
private

Definition at line 242 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalH().

TH1I * GctErrorAnalyzer::totalHtE_Of_
private

Definition at line 243 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer(), and plotTotalH().

std::string GctErrorAnalyzer::useSys_
private

Definition at line 160 of file GctErrorAnalyzer.cc.

Referenced by GctErrorAnalyzer().