#include <HcalTrigPrimMonitor.h>

Inheritance diagram for HcalTrigPrimMonitor:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)

void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)

void	beginRun (const edm::Run &run, const edm::EventSetup &c)

void	cleanup ()

void	endJob ()

void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)

	HcalTrigPrimMonitor (const edm::ParameterSet &ps)

void	processEvent (const edm::Handle< HcalTrigPrimDigiCollection > &data_tp_col, const edm::Handle< HcalTrigPrimDigiCollection > &emul_tp_col)

void	reset ()

void	setup ()

	~HcalTrigPrimMonitor ()

Public Member Functions inherited from HcalBaseDQMonitor
	HcalBaseDQMonitor (const edm::ParameterSet &ps)

	HcalBaseDQMonitor ()

virtual	~HcalBaseDQMonitor ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
enum	ErrorFlag { kZeroTP =-1, kMatched = 0, kMismatchedEt = 1, kMismatchedFG = 2, kMissingData = 3, kMissingEmul = 4, kNErrorFlag = 5, kUnknown = kNErrorFlag }

Private Member Functions
MonitorElement *	create_errorflag (const std::string &folder, const std::string &name)

MonitorElement *	create_et_histogram (const std::string &folder, const std::string &name)

MonitorElement *	create_fg_correlation (const std::string &folder, const std::string &name)

MonitorElement *	create_map (const std::string &folder, const std::string &name)

MonitorElement *	create_summary (const std::string &folder, const std::string &name)

MonitorElement *	create_tp_correlation (const std::string &folder, const std::string &name)

Private Attributes
MonitorElement *	bad_tps [2]

MonitorElement *	bad_tps_oot [2]

edm::InputTag	dataLabel_

edm::InputTag	emulLabel_

MonitorElement *	errorflag [2]

MonitorElement *	errorflag_oot [2]

MonitorElement *	fg_corr [2][2]

MonitorElement *	fg_corr_oot [2][2]

MonitorElement *	good_tps [2]

MonitorElement *	good_tps_oot [2]

int	nBad_TP_per_LS_HB_

int	nBad_TP_per_LS_HE_

int	nBad_TP_per_LS_HF_

std::map< ErrorFlag, MonitorElement * >	problem_et [2][2]

std::map< ErrorFlag, MonitorElement * >	problem_et_oot [2][2]

std::map< ErrorFlag, MonitorElement * >	problem_map [2]

std::map< ErrorFlag, MonitorElement * >	problem_map_oot [2]

edm::EDGetTokenT < HcalTrigPrimDigiCollection >	tok_data_

edm::EDGetTokenT < HcalTrigPrimDigiCollection >	tok_emu_

MonitorElement *	tp_corr [2][2]

MonitorElement *	tp_corr_oot [2][2]

MonitorElement *	TPOccupancy_

MonitorElement *	TPOccupancyEta_

MonitorElement *	TPOccupancyPhi_

MonitorElement *	TPOccupancyPhiHFM_

MonitorElement *	TPOccupancyPhiHFP_

std::vector< int >	ZSAlarmThreshold_

std::vector< int >	ZSBadTPThreshold_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from HcalBaseDQMonitor
virtual void	beginJob ()

virtual void	endRun (const edm::Run &run, const edm::EventSetup &c)

void	getLogicalMap (const edm::EventSetup &c)

bool	IsAllowedCalibType ()

bool	LumiInOrder (int lumisec)

void	SetupEtaPhiHists (EtaPhiHists &hh, std::string Name, std::string Units)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from HcalBaseDQMonitor
std::vector< int >	AllowedCalibTypes_

int	badChannelStatusMask_

int	currentLS

int	currenttype_

DQMStore *	dbe_

int	debug_

bool	enableCleanup_

bool	eventAllowed_

bool	HBpresent_

bool	HEpresent_

bool	HFpresent_

bool	HOpresent_

int	ievt_

std::map< unsigned int, int >	KnownBadCells_

int	levt_

HcalLogicalMap *	logicalMap_

bool	makeDiagnostics_

MonitorElement *	meIevt_

MonitorElement *	meLevt_

bool	mergeRuns_

MonitorElement *	meTevt_

MonitorElement *	meTevtHist_

bool	needLogicalMap_

int	NLumiBlocks_

bool	Online_

std::string	prefixME_

MonitorElement *	ProblemsCurrentLB

MonitorElement *	ProblemsVsLB

MonitorElement *	ProblemsVsLB_HB

MonitorElement *	ProblemsVsLB_HBHEHF

MonitorElement *	ProblemsVsLB_HE

MonitorElement *	ProblemsVsLB_HF

MonitorElement *	ProblemsVsLB_HO

bool	skipOutOfOrderLS_

std::string	subdir_

int	tevt_

Detailed Description

Author: W. Fisher - FNAL

Definition at line 13 of file HcalTrigPrimMonitor.h.

Member Enumeration Documentation

enum HcalTrigPrimMonitor::ErrorFlag

private

Enumerator
kZeroTP
kMatched
kMismatchedEt
kMismatchedFG
kMissingData
kMissingEmul
kNErrorFlag
kUnknown

Definition at line 46 of file HcalTrigPrimMonitor.h.

                     {
          kZeroTP=-1,
          kMatched = 0,
          kMismatchedEt = 1,
          kMismatchedFG = 2,
          kMissingData = 3,
          kMissingEmul = 4,
          kNErrorFlag = 5,
          kUnknown = kNErrorFlag
       };

Constructor & Destructor Documentation

HcalTrigPrimMonitor::HcalTrigPrimMonitor ( const edm::ParameterSet & ps )

Definition at line 4 of file HcalTrigPrimMonitor.cc.

References HcalBaseDQMonitor::AllowedCalibTypes_, dataLabel_, HcalBaseDQMonitor::debug_, emulLabel_, HcalBaseDQMonitor::enableCleanup_, edm::ParameterSet::getUntrackedParameter(), HcalBaseDQMonitor::makeDiagnostics_, HcalBaseDQMonitor::mergeRuns_, HcalBaseDQMonitor::NLumiBlocks_, HcalBaseDQMonitor::Online_, HcalBaseDQMonitor::prefixME_, HcalBaseDQMonitor::skipOutOfOrderLS_, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, tok_data_, and tok_emu_.

                                                                    :
    dataLabel_(ps.getParameter<edm::InputTag>("dataLabel")),
    emulLabel_(ps.getParameter<edm::InputTag>("emulLabel")),
    ZSBadTPThreshold_(ps.getParameter< std::vector<int> >("ZSBadTPThreshold")),
    ZSAlarmThreshold_(ps.getParameter< std::vector<int> >("ZSAlarmThreshold"))
 {
    Online_                = ps.getUntrackedParameter<bool>("online",false);
    mergeRuns_             = ps.getUntrackedParameter<bool>("mergeRuns",false);
    enableCleanup_         = ps.getUntrackedParameter<bool>("enableCleanup",false);
    debug_                 = ps.getUntrackedParameter<int>("debug",false);
    prefixME_              = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
    if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
     prefixME_.append("/");
    subdir_                = ps.getUntrackedParameter<std::string>("TaskFolder","TrigPrimMonitor_Hcal"); 
    if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
     subdir_.append("/");
    subdir_=prefixME_+subdir_;
    AllowedCalibTypes_     = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
    skipOutOfOrderLS_      = ps.getUntrackedParameter<bool>("skipOutOfOrderLS",true);
    NLumiBlocks_           = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
    makeDiagnostics_       = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
 
   // register for data access
   tok_data_ = consumes<HcalTrigPrimDigiCollection>(dataLabel_);
   tok_emu_ = consumes<HcalTrigPrimDigiCollection>(emulLabel_);
 }

HcalTrigPrimMonitor::~HcalTrigPrimMonitor ( )

Definition at line 32 of file HcalTrigPrimMonitor.cc.

32 {

33 }

Member Function Documentation

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 175 of file HcalTrigPrimMonitor.cc.

References HcalBaseDQMonitor::analyze(), dataLabel_, emulLabel_, edm::Event::getByToken(), HcalBaseDQMonitor::IsAllowedCalibType(), HcalBaseDQMonitor::LumiInOrder(), edm::EventBase::luminosityBlock(), processEvent(), tok_data_, and tok_emu_.

                                                                        {
    if (!IsAllowedCalibType()) return;
    if (LumiInOrder(e.luminosityBlock())==false) return;
 
    edm::Handle<HcalTrigPrimDigiCollection> data_tp_col;
    if (!e.getByToken(tok_data_, data_tp_col)) {
       edm::LogWarning("HcalTrigPrimMonitor")<< dataLabel_<< " data TP not available";
       return;
    }
 
    edm::Handle<HcalTrigPrimDigiCollection> emul_tp_col;
    if (!e.getByToken(tok_emu_, emul_tp_col)) {
       edm::LogWarning("HcalTrigPrimMonitor")<< emulLabel_<< " emul TP not available";
       return;
    }
 
    HcalBaseDQMonitor::analyze(e,s); // base class increments ievt_, etc. counters
    processEvent(data_tp_col, emul_tp_col);
 }

void HcalTrigPrimMonitor::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	c
	)

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 570 of file HcalTrigPrimMonitor.cc.

References HcalBaseDQMonitor::beginLuminosityBlock(), HcalBaseDQMonitor::LumiInOrder(), edm::LuminosityBlockBase::luminosityBlock(), nBad_TP_per_LS_HB_, nBad_TP_per_LS_HE_, nBad_TP_per_LS_HF_, HcalBaseDQMonitor::ProblemsCurrentLB, and MonitorElement::Reset().

                                                                                                         {
    if (LumiInOrder(lumiSeg.luminosityBlock())==false) return;
    HcalBaseDQMonitor::beginLuminosityBlock(lumiSeg,c);
    ProblemsCurrentLB->Reset();
    // Rest counter
    nBad_TP_per_LS_HB_ = 0;
    nBad_TP_per_LS_HE_ = 0;
    nBad_TP_per_LS_HF_ = 0;
 }

void HcalTrigPrimMonitor::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	c
	)

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 166 of file HcalTrigPrimMonitor.cc.

References HcalBaseDQMonitor::beginRun(), HcalBaseDQMonitor::mergeRuns_, reset(), setup(), and HcalBaseDQMonitor::tevt_.

 {
   HcalBaseDQMonitor::beginRun(run,c);
   if (mergeRuns_ && tevt_>0) return; // don't reset counters if merging runs
   if (tevt_==0) this->setup(); // create all histograms; not necessary if merging runs together
   if (mergeRuns_==false) this->reset(); // call reset at start of all runs
 } // void HcalTrigPrimMonitor::beginRun()

void HcalTrigPrimMonitor::cleanup ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 541 of file HcalTrigPrimMonitor.cc.

References HcalBaseDQMonitor::dbe_, HcalBaseDQMonitor::enableCleanup_, DQMStore::removeContents(), DQMStore::setCurrentFolder(), and HcalBaseDQMonitor::subdir_.

Referenced by endJob().

                              {
    if (!enableCleanup_) return;
    if (dbe_) {
       dbe_->setCurrentFolder(subdir_);
       dbe_->removeContents();
 
       dbe_->setCurrentFolder(subdir_ + "noZS/Problem TPs/TP Values");
       dbe_->removeContents();
       dbe_->setCurrentFolder(subdir_ + "noZS/Problem TPs");
       dbe_->removeContents();
       dbe_->setCurrentFolder(subdir_ + "noZS");
 
       dbe_->setCurrentFolder(subdir_ + "Problem TPs/TP Values");
       dbe_->removeContents();
       dbe_->setCurrentFolder(subdir_ + "Problem TPs");
       dbe_->removeContents();
 
       dbe_->setCurrentFolder(subdir_ + "Problem OOT TPs/TP Values");
       dbe_->removeContents();
       dbe_->setCurrentFolder(subdir_ + "Problem OOT TPs");
       dbe_->removeContents();
    }
 }

MonitorElement * HcalTrigPrimMonitor::create_errorflag	(	const std::string &	folder,
		const std::string &	name
	)

private

Definition at line 604 of file HcalTrigPrimMonitor.cc.

References DQMStore::book2D(), HcalBaseDQMonitor::dbe_, MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().

Referenced by setup().

                                                                                     {
    edm::LogInfo("HcalTrigPrimMonitor") << "Creating MonitorElement " << name << " in folder " << folder << "\n";
 
    dbe_->setCurrentFolder(folder);
    MonitorElement* element = dbe_->book2D(name, name, 4, 1, 5, 2, 0, 2);
    element->setBinLabel(1, "Mismatched E");
    element->setBinLabel(2, "Mismatched FG");
    element->setBinLabel(3, "Missing Data");
    element->setBinLabel(4, "Missing Emul");
    element->setBinLabel(1, "HBHE", 2);
    element->setBinLabel(2, "HF", 2);
    return element;
 }

MonitorElement * HcalTrigPrimMonitor::create_et_histogram	(	const std::string &	folder,
		const std::string &	name
	)

private

Definition at line 650 of file HcalTrigPrimMonitor.cc.

References DQMStore::book1D(), HcalBaseDQMonitor::dbe_, and DQMStore::setCurrentFolder().

Referenced by setup().

                                                                                        {
    edm::LogInfo("HcalTrigPrimMonitor") << "Creating MonitorElement " << name << " in folder " << folder << "\n";
 
    dbe_->setCurrentFolder(folder);
    return dbe_->book1D(name, name, 256, 0, 256);
 }

MonitorElement * HcalTrigPrimMonitor::create_fg_correlation	(	const std::string &	folder,
		const std::string &	name
	)

private

Definition at line 630 of file HcalTrigPrimMonitor.cc.

References DQMStore::book2D(), HcalBaseDQMonitor::dbe_, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().

Referenced by setup().

                                                                                          {
    edm::LogInfo("HcalTrigPrimMonitor") << "Creating MonitorElement " << name << " in folder " << folder << "\n";
 
    dbe_->setCurrentFolder(folder);
    MonitorElement* element = dbe_->book2D(name, name, 2, 0, 2, 2, 0, 2);
    element->setAxisTitle("data FG", 1);
    element->setAxisTitle("emul FG", 2);
    return element;
 }

MonitorElement * HcalTrigPrimMonitor::create_map	(	const std::string &	folder,
		const std::string &	name
	)

private

Definition at line 641 of file HcalTrigPrimMonitor.cc.

References DQMStore::book2D(), HcalBaseDQMonitor::dbe_, DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and indexGen::title.

Referenced by setup().

                                                                               {
    edm::LogInfo("HcalTrigPrimMonitor") << "Creating MonitorElement " << name << " in folder " << folder << "\n";
 
    dbe_->setCurrentFolder(folder);
    std::string title = name +";ieta;iphi";
    return dbe_->book2D(name, title, 65, -32.5, 32.5, 72, 0.5, 72.5);
 }

MonitorElement * HcalTrigPrimMonitor::create_summary	(	const std::string &	folder,
		const std::string &	name
	)

private

Definition at line 596 of file HcalTrigPrimMonitor.cc.

References DQMStore::book2D(), HcalBaseDQMonitor::dbe_, and DQMStore::setCurrentFolder().

                                                                                   {
    edm::LogInfo("HcalTrigPrimMonitor") << "Creating MonitorElement " << name << " in folder " << folder << "\n";
 
    dbe_->setCurrentFolder(folder);
    return dbe_->book2D(name, name, 65, -32.5, 32.5, 72, 0.5, 72.5);
 }

MonitorElement * HcalTrigPrimMonitor::create_tp_correlation	(	const std::string &	folder,
		const std::string &	name
	)

private

Definition at line 619 of file HcalTrigPrimMonitor.cc.

References DQMStore::book2D(), HcalBaseDQMonitor::dbe_, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().

Referenced by setup().

                                                                                          {
    edm::LogInfo("HcalTrigPrimMonitor") << "Creating MonitorElement " << name << " in folder " << folder << "\n";
 
    dbe_->setCurrentFolder(folder);
    MonitorElement* element = dbe_->book2D(name, name, 50, 0, 256, 50, 0, 256);
    element->setAxisTitle("data TP", 1);
    element->setAxisTitle("emul TP", 2);
    return element;
 }

void HcalTrigPrimMonitor::endJob ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 565 of file HcalTrigPrimMonitor.cc.

References cleanup(), and HcalBaseDQMonitor::enableCleanup_.

 {
   if (enableCleanup_) cleanup(); 
 }

void HcalTrigPrimMonitor::endLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	c
	)

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 580 of file HcalTrigPrimMonitor.cc.

References HcalBaseDQMonitor::currentLS, MonitorElement::Fill(), HcalBaseDQMonitor::levt_, HcalBaseDQMonitor::LumiInOrder(), edm::LuminosityBlockBase::luminosityBlock(), nBad_TP_per_LS_HB_, nBad_TP_per_LS_HE_, nBad_TP_per_LS_HF_, HcalBaseDQMonitor::ProblemsCurrentLB, HcalBaseDQMonitor::ProblemsVsLB, HcalBaseDQMonitor::ProblemsVsLB_HB, HcalBaseDQMonitor::ProblemsVsLB_HE, and HcalBaseDQMonitor::ProblemsVsLB_HF.

                                                                                                       {
    if (LumiInOrder(lumiSeg.luminosityBlock())==false) return;
    // Fill histograms for this LS
    ProblemsVsLB_HB->Fill(currentLS, nBad_TP_per_LS_HB_);
    ProblemsVsLB_HE->Fill(currentLS, nBad_TP_per_LS_HE_);
    ProblemsVsLB_HF->Fill(currentLS, nBad_TP_per_LS_HF_);
    ProblemsVsLB->Fill(currentLS, nBad_TP_per_LS_HB_ + nBad_TP_per_LS_HE_ + nBad_TP_per_LS_HF_);
 
    ProblemsCurrentLB->Fill(-1,-1,levt_);
    ProblemsCurrentLB->Fill(0,0, nBad_TP_per_LS_HB_);
    ProblemsCurrentLB->Fill(1,0, nBad_TP_per_LS_HE_);
    ProblemsCurrentLB->Fill(3,0, nBad_TP_per_LS_HF_);
 }

void HcalTrigPrimMonitor::processEvent	(	const edm::Handle< HcalTrigPrimDigiCollection > &	data_tp_col,
		const edm::Handle< HcalTrigPrimDigiCollection > &	emul_tp_col
	)

Definition at line 197 of file HcalTrigPrimMonitor.cc.

References funct::abs(), bad_tps, bad_tps_oot, HcalBaseDQMonitor::dbe_, diffTreeTool::diff, errorflag, errorflag_oot, fg_corr, fg_corr_oot, MonitorElement::Fill(), good_tps, good_tps_oot, i, HcalBaseDQMonitor::ievt_, isHF(), kMismatchedEt, kMismatchedFG, kMissingData, kMissingEmul, kNErrorFlag, max(), nBad_TP_per_LS_HB_, nBad_TP_per_LS_HE_, nBad_TP_per_LS_HF_, problem_et, problem_et_oot, problem_map, problem_map_oot, MonitorElement::setBinContent(), tp_corr, tp_corr_oot, TPOccupancy_, TPOccupancyEta_, TPOccupancyPhi_, TPOccupancyPhiHFM_, TPOccupancyPhiHFP_, ZSAlarmThreshold_, and ZSBadTPThreshold_.

Referenced by analyze().

                                                                  {
 
    if(dbe_ == 0) 
       return;
 
    std::vector<int> errorflag_per_event[2][2];
    std::vector<int> errorflag_per_event_oot[2][2];
    for (int isZS = 0; isZS <= 1; ++isZS) {
       for (int isHF = 0; isHF <= 1; ++isHF) {
          errorflag_per_event[isZS][isHF] = std::vector<int>(kNErrorFlag, 0);
          errorflag_per_event_oot[isZS][isHF] = std::vector<int>(kNErrorFlag, 0);
       }//for isHF
    }//for isZS
    
    for (int isZS = 0; isZS <= 1; ++isZS) {
      good_tps[isZS]->setBinContent(-1,-1,ievt_);
      bad_tps[isZS]->setBinContent(-1,-1,ievt_);
      good_tps_oot[isZS]->setBinContent(-1,-1,ievt_);
      bad_tps_oot[isZS]->setBinContent(-1,-1,ievt_);
    }
 
    for (HcalTrigPrimDigiCollection::const_iterator data_tp = data_tp_col->begin();
                                                   data_tp != data_tp_col->end();
                                                   ++data_tp) {
       int ieta = data_tp->id().ieta();
       int iphi = data_tp->id().iphi();
       int isHF = data_tp->id().ietaAbs() >= 29 ? 1 : 0;
 
 
       //
       if (data_tp->SOI_compressedEt() > 0) {
          TPOccupancy_->Fill(ieta, iphi);
          TPOccupancyEta_->Fill(ieta);
          TPOccupancyPhi_->Fill(iphi);
 
          if (isHF) {
             if (ieta > 0) {
                TPOccupancyPhiHFP_->Fill(iphi);
             }
             else {
                TPOccupancyPhiHFM_->Fill(iphi);
             }
          }
       }
 
       //check missing from emulator
       HcalTrigPrimDigiCollection::const_iterator emul_tp = emul_tp_col->find(data_tp->id());
       if (emul_tp == emul_tp_col->end()) {
          bool pass_ZS = true;    
          bool pass_ZS_OOT = true;    
      
          for (int i=0; i<data_tp->size(); ++i) {
        if(i==2) {
          int dataEt(data_tp->sample(i).compressedEt());
          problem_et[0][isHF][kMissingEmul]->Fill(dataEt);
 
          if (dataEt > ZSAlarmThreshold_[abs(ieta)]) {
                problem_et[1][isHF][kMissingEmul]->Fill(dataEt);
                pass_ZS = false;
          }
          
          problem_map[0][kMissingEmul]->Fill(ieta, iphi);
          ++errorflag_per_event[0][isHF][kMissingEmul];
          bad_tps[0]->Fill(ieta, iphi);
          
          if (!pass_ZS) {
            problem_map[1][kMissingEmul]->Fill(ieta, iphi);
            ++errorflag_per_event[1][isHF][kMissingEmul];
            bad_tps[1]->Fill(ieta, iphi);
            
            // counts per LS, only for inTime TP
            if (abs(ieta) <= 16)
          ++nBad_TP_per_LS_HB_;
            else if(abs(ieta) <= 28)
          ++nBad_TP_per_LS_HE_;
            else
          ++nBad_TP_per_LS_HF_;
          }
        }
        
        if(i!=2) {
          int dataEt(data_tp->sample(i).compressedEt());
          problem_et_oot[0][isHF][kMissingEmul]->Fill(dataEt);
          
          if (dataEt > ZSAlarmThreshold_[abs(ieta)]) {
            problem_et_oot[1][isHF][kMissingEmul]->Fill(dataEt);
            pass_ZS_OOT = false;
          }
          
          problem_map_oot[0][kMissingEmul]->Fill(ieta, iphi);
          ++errorflag_per_event_oot[0][isHF][kMissingEmul];
          bad_tps_oot[0]->Fill(ieta, iphi);
          
          if (!pass_ZS_OOT) {
            problem_map_oot[1][kMissingEmul]->Fill(ieta, iphi);
            ++errorflag_per_event_oot[1][isHF][kMissingEmul];
            bad_tps_oot[1]->Fill(ieta, iphi);
          }       
        }       
      }
       } //emul tp not found
       else {
          bool mismatchedEt_noZS = false;
          bool mismatchedEt_ZS = false;
          bool mismatchedFG_noZS = false;
          bool mismatchedFG_ZS = false;
 
          bool mismatchedEt_OOT_noZS = false;
          bool mismatchedEt_OOT_ZS = false;
          bool mismatchedFG_OOT_noZS = false;
          bool mismatchedFG_OOT_ZS = false;
 
          for (int i=0; i<data_tp->size(); ++i) {
             int dataEt(data_tp->sample(i).compressedEt());
             int dataFG(data_tp->sample(i).fineGrain());
             int emulEt(emul_tp->sample(i).compressedEt());
             int emulFG(emul_tp->sample(i).fineGrain());
         
             int diff = abs(dataEt - emulEt);
             bool fill_corr_ZS = true;
             bool fill_corr_OOT_ZS = true;
 
         if (std::max(dataEt, emulEt) < ZSAlarmThreshold_.at(abs(ieta)))
           continue;
             
             if (diff == 0) {
                if (dataFG != emulFG) {
          if(i==2) {
            mismatchedFG_noZS = true;
            problem_et[0][isHF][kMismatchedFG]->Fill(dataEt);
            
            // exclude mismatched FG when HF TP < ZS_AlarmThreshold
            if (isHF == 1 && dataEt <= ZSAlarmThreshold_.at(abs(ieta))) {
              // Do not fill ZS correlation plots.
                      fill_corr_ZS = false;
            }
            else {
                      mismatchedFG_ZS = true;
                      problem_et[1][isHF][kMismatchedFG]->Fill(dataEt);
            }
          }
          if(i!=2){
            mismatchedFG_OOT_noZS = true;
            problem_et_oot[0][isHF][kMismatchedFG]->Fill(dataEt);
            
            // exclude mismatched FG when HF TP < ZS_AlarmThreshold
            if (isHF == 1 && dataEt <= ZSAlarmThreshold_.at(abs(ieta))) {
              // Do not fill ZS correlation plots.
              fill_corr_OOT_ZS = false;
            }
            else {
              mismatchedFG_OOT_ZS = true;
              problem_et_oot[1][isHF][kMismatchedFG]->Fill(dataEt);
            }            
          }       
                } // matched et but not fg
             }
             else {
           if(i==2) {
         mismatchedEt_noZS = true;
         //if (diff > ZSAlarmThreshold_.at(abs(ieta))) {
         if (diff > ZSBadTPThreshold_.at(abs(ieta))) {
           mismatchedEt_ZS = true;
           fill_corr_ZS = false;
         }
           }
           if(i!=2) {
         mismatchedEt_OOT_noZS = true;
         //if (diff > ZSAlarmThreshold_.at(abs(ieta))) {
         if (diff > ZSBadTPThreshold_.at(abs(ieta))) {
           mismatchedEt_OOT_ZS = true;
           fill_corr_OOT_ZS = false;       
         }
           }     
         } // mismatche et
 
             // Correlation plots
         if(i==2) {        
           tp_corr[0][isHF]->Fill(dataEt, emulEt);
           fg_corr[0][isHF]->Fill(dataFG, emulFG);
           
           if (fill_corr_ZS) {
         tp_corr[1][isHF]->Fill(dataEt, emulEt);
         fg_corr[1][isHF]->Fill(dataFG, emulFG);
           }
         }
         if(i!=2) {        
           tp_corr_oot[0][isHF]->Fill(dataEt, emulEt);
           fg_corr_oot[0][isHF]->Fill(dataFG, emulFG);
           
           if (fill_corr_OOT_ZS) {
         tp_corr_oot[1][isHF]->Fill(dataEt, emulEt);
         fg_corr_oot[1][isHF]->Fill(dataFG, emulFG);
           }
         }       
          }//for tp sample
 
          // Fill Problem Map and error counts
          if (mismatchedEt_noZS) {
             problem_map[0][kMismatchedEt]->Fill(ieta, iphi);
             ++errorflag_per_event[0][isHF][kMismatchedEt];
          }
          if (mismatchedEt_ZS) {
             problem_map[1][kMismatchedEt]->Fill(ieta, iphi);
             ++errorflag_per_event[1][isHF][kMismatchedEt];
          }
          if (mismatchedFG_noZS) {
             problem_map[0][kMismatchedFG]->Fill(ieta, iphi);
             ++errorflag_per_event[0][isHF][kMismatchedFG];
          }
          if (mismatchedFG_ZS) {
             problem_map[1][kMismatchedFG]->Fill(ieta, iphi);
             ++errorflag_per_event[1][isHF][kMismatchedFG];
          }
          if (mismatchedEt_noZS || mismatchedFG_noZS)
             bad_tps[0]->Fill(ieta, iphi);
          else
             good_tps[0]->Fill(ieta, iphi);
          if (mismatchedEt_ZS || mismatchedFG_ZS) {
        
        bad_tps[1]->Fill(ieta, iphi);
        
        // counts per LS
        if (abs(ieta) <= 16)
          ++nBad_TP_per_LS_HB_;
        else if(abs(ieta) <= 28)
          ++nBad_TP_per_LS_HE_;
        else
          ++nBad_TP_per_LS_HF_;
          }
          else
             good_tps[1]->Fill(ieta, iphi);
 
      // OOT Copy
      // Fill Problem Map and error counts
          if (mismatchedEt_OOT_noZS) {
             problem_map_oot[0][kMismatchedEt]->Fill(ieta, iphi);
             ++errorflag_per_event_oot[0][isHF][kMismatchedEt];
          }
          if (mismatchedEt_OOT_ZS) {
             problem_map_oot[1][kMismatchedEt]->Fill(ieta, iphi);
             ++errorflag_per_event_oot[1][isHF][kMismatchedEt];
          }
          if (mismatchedFG_noZS) {
             problem_map_oot[0][kMismatchedFG]->Fill(ieta, iphi);
             ++errorflag_per_event_oot[0][isHF][kMismatchedFG];
          }
          if (mismatchedFG_OOT_ZS) {
             problem_map_oot[1][kMismatchedFG]->Fill(ieta, iphi);
             ++errorflag_per_event_oot[1][isHF][kMismatchedFG];
          }
          if (mismatchedEt_OOT_noZS || mismatchedFG_OOT_noZS)
             bad_tps_oot[0]->Fill(ieta, iphi);
          else
             good_tps_oot[0]->Fill(ieta, iphi);
          if (mismatchedEt_OOT_ZS || mismatchedFG_OOT_ZS) {
        
        bad_tps_oot[1]->Fill(ieta, iphi);
          }
          else
             good_tps_oot[1]->Fill(ieta, iphi);
       }//emul tp found
    }//for data_tp_col
 
 
    //check missing from data
    for (HcalTrigPrimDigiCollection::const_iterator emul_tp = emul_tp_col->begin();
                                                    emul_tp != emul_tp_col->end();
                                                    ++emul_tp) {
       int ieta(emul_tp->id().ieta());
       int iphi(emul_tp->id().iphi());
       int isHF = emul_tp->id().ietaAbs() >= 29 ? 1 : 0;
 
       HcalTrigPrimDigiCollection::const_iterator data_tp = data_tp_col->find(emul_tp->id());
       if (data_tp == data_tp_col->end()) {
          bool pass_ZS = true;
          bool pass_OOT_ZS = true;
 
          for (int i=0; i<emul_tp->size(); ++i) {
             int emulEt(emul_tp->sample(i).compressedEt());
         if(i==2) {
           problem_et[0][isHF][kMissingData]->Fill(emulEt);
           
           if (emulEt > ZSAlarmThreshold_[abs(ieta)]) {
         problem_et[1][isHF][kMissingData]->Fill(emulEt);
         pass_ZS = false;
           }
 
           problem_map[0][kMissingData]->Fill(ieta, iphi);
           ++errorflag_per_event[0][isHF][kMissingData];
           bad_tps[0]->Fill(ieta, iphi);
           
           if (!pass_ZS) {
         problem_map[1][kMissingData]->Fill(ieta, iphi);
         ++errorflag_per_event[1][isHF][kMissingData];
         bad_tps[1]->Fill(ieta, iphi);
         
         // counts per LS
         if (abs(ieta) <= 16)
           ++nBad_TP_per_LS_HB_;
         else if(abs(ieta) <= 28)
           ++nBad_TP_per_LS_HE_;
         else
           ++nBad_TP_per_LS_HF_;
           }
         }
         else {
           problem_et_oot[0][isHF][kMissingData]->Fill(emulEt);
           
           if (emulEt > ZSAlarmThreshold_[abs(ieta)]) {
         problem_et_oot[1][isHF][kMissingData]->Fill(emulEt);
         pass_OOT_ZS = false;
           }   
           
           problem_map_oot[0][kMissingData]->Fill(ieta, iphi);
           ++errorflag_per_event_oot[0][isHF][kMissingData];
           bad_tps_oot[0]->Fill(ieta, iphi);
           
           if (!pass_OOT_ZS) {
         problem_map_oot[1][kMissingData]->Fill(ieta, iphi);
         ++errorflag_per_event_oot[1][isHF][kMissingData];
         bad_tps_oot[1]->Fill(ieta, iphi);
           }
         }
      }//for tp sample
       } //data tp not found
    } //for emul_tp_col
 
    // Fill error flag per event
    for (int isZS = 0; isZS <= 1; ++isZS) {
       for (int isHF = 0; isHF <= 1; ++isHF) {
          for (int i=0; i<kNErrorFlag; ++i) {
             if (errorflag_per_event[isZS][isHF][i] > 0)
                errorflag[isZS]->Fill(i, isHF);
             if (errorflag_per_event_oot[isZS][isHF][i] > 0)
                errorflag_oot[isZS]->Fill(i, isHF);
          }//for i
       }//for isHF
    }//for isZS
 }

void HcalTrigPrimMonitor::reset ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 37 of file HcalTrigPrimMonitor.cc.

Referenced by beginRun().

37 {

38 }

void HcalTrigPrimMonitor::setup ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 42 of file HcalTrigPrimMonitor.cc.

References bad_tps, bad_tps_oot, DQMStore::book1D(), DQMStore::bookProfile(), create_errorflag(), create_et_histogram(), create_fg_correlation(), create_map(), create_tp_correlation(), HcalBaseDQMonitor::dbe_, errorflag, errorflag_oot, fg_corr, fg_corr_oot, MonitorElement::getTProfile(), good_tps, good_tps_oot, isHF(), kMismatchedEt, kMismatchedFG, kMissingData, kMissingEmul, HcalBaseDQMonitor::NLumiBlocks_, problem_et, problem_et_oot, problem_map, problem_map_oot, HcalBaseDQMonitor::ProblemsVsLB, HcalBaseDQMonitor::ProblemsVsLB_HB, HcalBaseDQMonitor::ProblemsVsLB_HE, HcalBaseDQMonitor::ProblemsVsLB_HF, HcalBaseDQMonitor::ProblemsVsLB_HO, DQMStore::setCurrentFolder(), HcalBaseDQMonitor::setup(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, tp_corr, tp_corr_oot, TPOccupancy_, TPOccupancyEta_, TPOccupancyPhi_, TPOccupancyPhiHFM_, and TPOccupancyPhiHFP_.

Referenced by beginRun().

                            {
    HcalBaseDQMonitor::setup();
    
    if (dbe_ == 0)
       return;
 
    dbe_->setCurrentFolder(subdir_ + "TP Occupancy");
    TPOccupancyEta_ = dbe_->book1D("TPOccupancyVsEta", "TPOccupancyVsEta", 65, -32.5, 32.5);
    TPOccupancyPhi_ = dbe_->book1D("TPOccupancyVsPhi", "TPOccupancyVsPhi", 72, 0.5, 72.5);
    TPOccupancyPhiHFP_ = dbe_->book1D("TPOccupancyHFPVsPhi", "TPOccupancyHFPVsPhi", 72, 0.5, 72.5);
    TPOccupancyPhiHFM_ = dbe_->book1D("TPOccupancyHFMVsPhi", "TPOccupancyHFMVsPhi", 72, 0.5, 72.5);
    TPOccupancy_ = create_map(subdir_ + "TP Occupancy", "TPOccupancy");
 
    for (int isZS = 0; isZS <= 1; ++isZS) {
 
       std::string folder(subdir_);
       std::string zsname="_ZS";
       if (isZS == 0)
     {
       folder += "noZS/";
       zsname="_noZS";
     }
 
       std::string problem_folder(folder);
       problem_folder += "Problem TPs/";
       
       good_tps[isZS] = create_map(folder, "Good TPs"+zsname);
       bad_tps[isZS] = create_map(folder, "Bad TPs"+zsname);
 
       errorflag[isZS] = create_errorflag(folder, "Error Flag"+zsname);
       problem_map[isZS][kMismatchedEt] = create_map(problem_folder, "Mismatched Et"+zsname);
       problem_map[isZS][kMismatchedFG] = create_map(problem_folder, "Mismatched FG"+zsname);
       problem_map[isZS][kMissingData] = create_map(problem_folder, "Missing Data"+zsname);
       problem_map[isZS][kMissingEmul] = create_map(problem_folder, "Missing Emul"+zsname);
 
       for (int isHF = 0; isHF <= 1; ++isHF) {
          std::string subdet = (isHF == 0 ? "HBHE " : "HF ");
          tp_corr[isZS][isHF] = create_tp_correlation(folder, subdet + "TP Correlation"+zsname);
          fg_corr[isZS][isHF] = create_fg_correlation(folder, subdet + "FG Correlation"+zsname);
 
          problem_et[isZS][isHF][kMismatchedFG]
             = create_et_histogram(problem_folder + "TP Values/", subdet + "Mismatched FG"+zsname);
 
          problem_et[isZS][isHF][kMissingData]
             = create_et_histogram(problem_folder + "TP Values/", subdet + "Missing Data"+zsname);
 
          problem_et[isZS][isHF][kMissingEmul]
             = create_et_histogram(problem_folder + "TP Values/", subdet + "Missing Emul"+zsname);
       }//isHF
    }//isZS
 
    // Copy for OOT TPs
    for (int isZS = 0; isZS <= 1; ++isZS) {
 
       std::string folder(subdir_);
       std::string zsname="_ZS";
       if (isZS == 0)
     {
       folder += "noZS/";
       zsname="_noZS";
     }
 
       std::string problem_folder(folder);
       problem_folder += "Problem OOT TPs/";
       
       good_tps_oot[isZS] = create_map(folder, "Good OOT TPs"+zsname);
       bad_tps_oot[isZS] = create_map(folder, "Bad OOT TPs"+zsname);
 
       errorflag_oot[isZS] = create_errorflag(folder, "Error Flag OOT"+zsname);
       problem_map_oot[isZS][kMismatchedEt] = create_map(problem_folder, "Mismatched OOT Et"+zsname);
       problem_map_oot[isZS][kMismatchedFG] = create_map(problem_folder, "Mismatched OOT FG"+zsname);
       problem_map_oot[isZS][kMissingData] = create_map(problem_folder, "Missing OOT Data"+zsname);
       problem_map_oot[isZS][kMissingEmul] = create_map(problem_folder, "Missing OOT Emul"+zsname);
 
       for (int isHF = 0; isHF <= 1; ++isHF) {
          std::string subdet = (isHF == 0 ? "HBHE " : "HF ");
          tp_corr_oot[isZS][isHF] = create_tp_correlation(folder, subdet + "OOT TP Correlation"+zsname);
          fg_corr_oot[isZS][isHF] = create_fg_correlation(folder, subdet + "OOT FG Correlation"+zsname);
 
          problem_et_oot[isZS][isHF][kMismatchedFG]
             = create_et_histogram(problem_folder + "TP Values/", subdet + "OOT Mismatched FG"+zsname);
 
          problem_et_oot[isZS][isHF][kMissingData]
             = create_et_histogram(problem_folder + "TP Values/", subdet + "OOT Missing Data"+zsname);
 
          problem_et_oot[isZS][isHF][kMissingEmul]
             = create_et_histogram(problem_folder + "TP Values/", subdet + "OOT Missing Emul"+zsname);
       }//isHF
    }//isZS
 
    // Number of bad cells vs. luminosity block
    ProblemsVsLB = dbe_->bookProfile(
          "TotalBadTPs_HCAL_vs_LS",
          "Total Number of Bad HCAL TPs vs lumi section",
          NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,10000);
 
    ProblemsVsLB_HB = dbe_->bookProfile(
          "TotalBadTPs_HB_vs_LS",
          "Total Number of Bad HB TPs vs lumi section",
          NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,3000);
 
    ProblemsVsLB_HE = dbe_->bookProfile(
          "TotalBadTPs_HE_vs_LS",
          "Total Number of Bad HE TPs vs lumi section",
          NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,3000);
 
    ProblemsVsLB_HF = dbe_->bookProfile(
          "TotalBadTPs_HF_vs_LS",
          "Total Number of Bad HF TPs vs lumi section",
          NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,3000);
 
    // No TPs for HO, DO NOT fill this histogram
    ProblemsVsLB_HO = dbe_->bookProfile(
          "TotalBadTPs_HO_vs_LS",
          "Total Number of Bad HO TPs vs lumi section",
          NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,3000);
 
    ProblemsVsLB->getTProfile()->SetMarkerStyle(20);
    ProblemsVsLB_HB->getTProfile()->SetMarkerStyle(20);
    ProblemsVsLB_HE->getTProfile()->SetMarkerStyle(20);
    ProblemsVsLB_HO->getTProfile()->SetMarkerStyle(20);
    ProblemsVsLB_HF->getTProfile()->SetMarkerStyle(20);
 }