#include <HcalCoarsePedestalMonitor.h>

Inheritance diagram for HcalCoarsePedestalMonitor:

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

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

void	bookHistograms (DQMStore::IBooker &ib, const edm::Run &run, const edm::EventSetup &c)

void	endJob ()

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

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

	HcalCoarsePedestalMonitor (const edm::ParameterSet &ps)

void	processEvent (const HBHEDigiCollection &hbhe, const HODigiCollection &ho, const HFDigiCollection &hf, const HcalUnpackerReport &report)

void	reset ()

void	setup (DQMStore::IBooker &)

	~HcalCoarsePedestalMonitor ()

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

	HcalBaseDQMonitor ()

virtual	~HcalBaseDQMonitor ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	fill_Nevents (const HcalTopology &)
	Methods, variables accessible only within class code. More...

void	zeroCounters ()

Private Attributes
double	ADCDiffThresh_

EtaPhiHists	CoarsePedestalsOccByDepth

EtaPhiHists	CoarsePedestalsSumByDepth

edm::InputTag	digiLabel_

bool	excludeHORing2_

int	minEvents_

int	pedestalocc_ [85][72][4]

double	pedestalsum_ [85][72][4]

edm::EDGetTokenT < HBHEDigiCollection >	tok_hbhe_

edm::EDGetTokenT < HFDigiCollection >	tok_hf_

edm::EDGetTokenT < HODigiCollection >	tok_ho_

edm::EDGetTokenT < HcalUnpackerReport >	tok_report_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from HcalBaseDQMonitor
virtual void	cleanup (void)

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

void	getLogicalMap (const edm::EventSetup &c)

bool	IsAllowedCalibType ()

bool	LumiInOrder (int lumisec)

void	SetupEtaPhiHists (DQMStore::IBooker &ib, 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_

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: J. Temple - Univ. of Maryland

Definition at line 17 of file HcalCoarsePedestalMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 15 of file HcalCoarsePedestalMonitor.cc.

                                                                               : HcalBaseDQMonitor(ps) {
   Online_                = ps.getUntrackedParameter<bool>("online",false);
   mergeRuns_             = ps.getUntrackedParameter<bool>("mergeRuns",false);
   enableCleanup_         = ps.getUntrackedParameter<bool>("enableCleanup",false);
   debug_                 = ps.getUntrackedParameter<int>("debug",0);
   prefixME_              = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
   if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
     prefixME_.append("/");
   subdir_                = ps.getUntrackedParameter<std::string>("TaskFolder","CoarsePedestalMonitor_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);
   digiLabel_             = ps.getUntrackedParameter<edm::InputTag>("digiLabel");
   ADCDiffThresh_         = ps.getUntrackedParameter<double>("ADCDiffThresh",1.);
   minEvents_             = ps.getUntrackedParameter<int>("minEvents",100); // minimum number of events needed before histograms are filled
   excludeHORing2_       = ps.getUntrackedParameter<bool>("excludeHORing2",false);
 
   tok_hbhe_ = consumes<HBHEDigiCollection>(digiLabel_);
   tok_ho_ = consumes<HODigiCollection>(digiLabel_);
   tok_hf_ = consumes<HFDigiCollection>(digiLabel_);
   tok_report_ = consumes<HcalUnpackerReport>(digiLabel_);
 
 }

HcalCoarsePedestalMonitor::~HcalCoarsePedestalMonitor ( )

Definition at line 45 of file HcalCoarsePedestalMonitor.cc.

45 { }

Member Function Documentation

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 104 of file HcalCoarsePedestalMonitor.cc.

References HcalBaseDQMonitor::analyze(), gather_cfg::cout, HcalBaseDQMonitor::debug_, digiLabel_, edm::Event::getByToken(), HcalBaseDQMonitor::ievt_, HcalBaseDQMonitor::IsAllowedCalibType(), HcalBaseDQMonitor::LumiInOrder(), edm::EventBase::luminosityBlock(), processEvent(), zeeHLT_cff::report, tok_hbhe_, tok_hf_, tok_ho_, and tok_report_.

                                                                                {
   HcalBaseDQMonitor::analyze(e,s);
 
   if (!IsAllowedCalibType()) return;
   if (LumiInOrder(e.luminosityBlock())==false) return;
 
 
   // try to get digis
   edm::Handle<HBHEDigiCollection> hbhe_digi;
   edm::Handle<HODigiCollection> ho_digi;
   edm::Handle<HFDigiCollection> hf_digi;
   edm::Handle<HcalUnpackerReport> report;
 
   if (!(e.getByToken(tok_hbhe_,hbhe_digi)))
     {
       edm::LogWarning("HcalCoarsePedestalMonitor")<< digiLabel_<<" hbhe_digi not available";
       return;
     }
   
   if (!(e.getByToken(tok_hf_,hf_digi)))
     {
       edm::LogWarning("HcalCoarsePedestalMonitor")<< digiLabel_<<" hf_digi not available";
       return;
     }
   if (!(e.getByToken(tok_ho_,ho_digi)))
     {
       edm::LogWarning("HcalCoarsePedestalMonitor")<< digiLabel_<<" ho_digi not available";
       return;
     }
   if (!(e.getByToken(tok_report_,report)))
     {
       edm::LogWarning("HcalCoarsePedestalMonitor")<< digiLabel_<<" unpacker report not available";
       return;
     }
 
   // all objects grabbed; event is good
   if (debug_>1) std::cout <<"\t<HcalCoarsePedestalMonitor::analyze>  Processing good event! event # = "<<ievt_<<std::endl;
 
 //  HcalBaseDQMonitor::analyze(e,s); // base class increments ievt_, etc. counters
 
   // Digi collection was grabbed successfully; process the Event
   processEvent(*hbhe_digi, *ho_digi, *hf_digi, *report);
 
 } //void HcalCoarsePedestalMonitor::analyze(...)

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 265 of file HcalCoarsePedestalMonitor.cc.

References HcalBaseDQMonitor::beginLuminosityBlock(), HcalBaseDQMonitor::ProblemsCurrentLB, and MonitorElement::Reset().

 {
   HcalBaseDQMonitor::beginLuminosityBlock(lumiSeg,c);
   ProblemsCurrentLB->Reset();
 }

void HcalCoarsePedestalMonitor::bookHistograms	(	DQMStore::IBooker &	ib,
		const edm::Run &	run,
		const edm::EventSetup &	c
	)

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 94 of file HcalCoarsePedestalMonitor.cc.

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

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

void HcalCoarsePedestalMonitor::endJob ( void )

Definition at line 53 of file HcalCoarsePedestalMonitor.cc.

References HcalBaseDQMonitor::cleanup(), gather_cfg::cout, HcalBaseDQMonitor::debug_, and HcalBaseDQMonitor::enableCleanup_.

 {
   if (debug_>0) std::cout <<"HcalCoarsePedestalMonitor::endJob()"<<std::endl;
   if (enableCleanup_) cleanup(); // when do we force cleanup?
 }

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 272 of file HcalCoarsePedestalMonitor.cc.

References fill_Nevents(), edm::EventSetup::get(), HcalBaseDQMonitor::LumiInOrder(), and edm::LuminosityBlockBase::luminosityBlock().

                                                {
   if (LumiInOrder(lumiSeg.luminosityBlock())==false) return;
 
   edm::ESHandle<HcalTopology> topo;
   c.get<HcalRecNumberingRecord>().get(topo);
   fill_Nevents(*topo);
   return;
 }

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 48 of file HcalCoarsePedestalMonitor.cc.

 {
   // Anything to do here?
 }

void HcalCoarsePedestalMonitor::fill_Nevents ( const HcalTopology & topology )

private

Methods, variables accessible only within class code.

Definition at line 282 of file HcalCoarsePedestalMonitor.cc.

References funct::abs(), CalcEtaBin(), CoarsePedestalsOccByDepth, CoarsePedestalsSumByDepth, ztail::d, EtaPhiHists::depth, eta, HcalObjRepresent::FillUnphysicalHEHFBins(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, i, HcalBaseDQMonitor::ievt_, minEvents_, pedestalocc_, pedestalsum_, phi, HcalTopology::validDetId(), and ecaldqm::zside().

Referenced by endLuminosityBlock().

                                                                          {
 
   // require minimum number of events before processing
   if (ievt_<minEvents_)
     return;
 
 
   // Set underflow bin to the number of pedestal events processed
   // (Assume that number of ped events is the same for all channels/subdetectors)
   for (unsigned int i=0;i<CoarsePedestalsSumByDepth.depth.size();++i)
     CoarsePedestalsSumByDepth.depth[i]->setBinContent(0,0,ievt_);
   for (unsigned int i=0;i<CoarsePedestalsOccByDepth.depth.size();++i)
     CoarsePedestalsOccByDepth.depth[i]->setBinContent(0,0,ievt_);
 
 
   int iphi=-1, ieta=-99, idepth=0, calcEta=-99;
   // Loop over all depths, eta, phi
   for (int d=0;d<4;++d) {
     idepth=d+1;
     for (int phi=0;phi<72;++phi) {
       iphi=phi+1; // actual iphi value
       for (int eta=0;eta<83;++eta) {
     ieta=eta-41; // actual ieta value;
     if (topology.validDetId(HcalBarrel, ieta, iphi, idepth)) {
       calcEta = CalcEtaBin(HcalBarrel,ieta,idepth);
       CoarsePedestalsSumByDepth.depth[d]->Fill(ieta,iphi,pedestalsum_[calcEta][phi][d]);
       CoarsePedestalsOccByDepth.depth[d]->Fill(ieta,iphi,pedestalocc_[calcEta][phi][d]);
     }
     if (topology.validDetId(HcalEndcap, ieta, iphi, idepth)) {
       calcEta = CalcEtaBin(HcalEndcap,ieta,idepth);
       CoarsePedestalsSumByDepth.depth[d]->Fill(ieta,iphi,pedestalsum_[calcEta][phi][d]);
       CoarsePedestalsOccByDepth.depth[d]->Fill(ieta,iphi,pedestalocc_[calcEta][phi][d]);
     }
     if (topology.validDetId(HcalOuter, ieta, iphi, idepth)) {
       calcEta = CalcEtaBin(HcalOuter,ieta,idepth);
       CoarsePedestalsSumByDepth.depth[d]->Fill(ieta,iphi,pedestalsum_[calcEta][phi][d]);
       CoarsePedestalsOccByDepth.depth[d]->Fill(ieta,iphi,pedestalocc_[calcEta][phi][d]);
     }
     if (topology.validDetId(HcalForward, ieta, iphi, idepth))   {
       calcEta = CalcEtaBin(HcalBarrel,ieta,idepth);
       int zside=ieta/abs(ieta);
       CoarsePedestalsSumByDepth.depth[d]->Fill(ieta+zside,iphi,pedestalsum_[calcEta][phi][d]);
       CoarsePedestalsOccByDepth.depth[d]->Fill(ieta+zside,iphi,pedestalocc_[calcEta][phi][d]);
     }
       }
     }
   }
   // Fill unphysical bins
   FillUnphysicalHEHFBins(CoarsePedestalsSumByDepth);
   FillUnphysicalHEHFBins(CoarsePedestalsOccByDepth);
 
   return;
 } // void HcalCoarsePedestalMonitor::fill_Nevents(const HcalTopology&)

void HcalCoarsePedestalMonitor::processEvent	(	const HBHEDigiCollection &	hbhe,
		const HODigiCollection &	ho,
		const HFDigiCollection &	hf,
		const HcalUnpackerReport &	report
	)

Definition at line 150 of file HcalCoarsePedestalMonitor.cc.

References funct::abs(), HcalQIESample::adc(), HcalUnpackerReport::badQualityDigis(), edm::SortedCollection< T, SORT >::begin(), CalcEtaBin(), HcalDetId::depth(), HLT_25ns10e33_v2_cff::depth, edm::SortedCollection< T, SORT >::end(), excludeHORing2_, HcalBaseDQMonitor::HBpresent_, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HcalBaseDQMonitor::HEpresent_, HcalBaseDQMonitor::HFpresent_, HcalBaseDQMonitor::HOpresent_, i, HFDataFrame::id(), HBHEDataFrame::id(), HODataFrame::id(), HcalDetId::ieta(), HcalDetId::iphi(), isSiPM(), j, pedestalocc_, pedestalsum_, HFDataFrame::sample(), HBHEDataFrame::sample(), HODataFrame::sample(), HFDataFrame::size(), HBHEDataFrame::size(), HODataFrame::size(), edm::SortedCollection< T, SORT >::size(), HcalDetId::subdet(), and relativeConstraints::value.

Referenced by analyze().

 { 
   // Skip events in which minimal good digis found -- still getting some strange (calib?) events through DQM
   
   unsigned int allgooddigis= hbhe.size()+ho.size()+hf.size();
   // bad threshold:  ignore events in which bad outnumber good by more than 100:1
   // (one RBX in HBHE seems to send valid data occasionally even on QIE resets, which is why we can't just require allgooddigis==0 when looking for events to skip)
   if ((allgooddigis==0) ||
       (1.*report.badQualityDigis()>100*allgooddigis))
     {
       return;
     }
 
 
   unsigned int digisize=0;
   int depth=0, iphi=0, ieta=0, binEta=-9999;
 
   double value=0;
 
   for (HBHEDigiCollection::const_iterator j=hbhe.begin(); j!=hbhe.end(); ++j)
     {
     const HBHEDataFrame digi = (const HBHEDataFrame)(*j);
     if (digi.id().subdet()==HcalBarrel)
       {if (!HBpresent_) continue;}
     else if (digi.id().subdet()==HcalEndcap)
       {if (!HEpresent_) continue;}
     else continue;
     digisize=digi.size();
     if (digisize<8) 
       continue;
     
     depth=digi.id().depth();
     iphi=digi.id().iphi();
     ieta=digi.id().ieta();
     
     digi.id().subdet()==HcalBarrel ? 
       binEta=CalcEtaBin(HcalBarrel, ieta, depth) :
       binEta=CalcEtaBin(HcalEndcap, ieta, depth);
       
     // 'value' is the average pedestal over the 8 time slices.
     // In the CoarsePedestal client, we will divide the summed value by Nevents (in underflow bin)
     // in order to calculate average pedestals.
     value=0;
     for (unsigned int i=0;i<8;++i)
       value+=digi.sample(i).adc()/8.;
     pedestalsum_[binEta][iphi-1][depth-1]+=value;
     ++pedestalocc_[binEta][iphi-1][depth-1];
     }
 
 
   if (HOpresent_)
     {
       for (HODigiCollection::const_iterator j=ho.begin(); j!=ho.end(); ++j)
     {
       const HODataFrame digi = (const HODataFrame)(*j);
       digisize=digi.size();
       if (digisize<8) 
         continue;
       
       depth=digi.id().depth();
       iphi=digi.id().iphi();
       ieta=digi.id().ieta();
       binEta=CalcEtaBin(HcalOuter, ieta, depth);
 
       // Don't fill cells that are part of HO ring 2 if an exclusion is applied
       if (excludeHORing2_==true && abs(ieta)>10 && isSiPM(ieta,iphi,4)==false)
         continue;
 
       // 'value' is the average pedestal over the 8 time slices.
       // In the CoarsePedestal client, we will divide the summed value by Nevents (in underflow bin)
       // in order to calculate average pedestals.
       value=0;
       for (unsigned int i=0;i<8;++i)
         value+=digi.sample(i).adc()/8.;
       pedestalsum_[binEta][iphi-1][depth-1]+=value;
       ++pedestalocc_[binEta][iphi-1][depth-1];
     } // for (HODigiCollection)
     }
   
   if (HFpresent_)
     {
       for (HFDigiCollection::const_iterator j=hf.begin(); j!=hf.end(); ++j)
     {
       const HFDataFrame digi = (const HFDataFrame)(*j);
       digisize=digi.size();
       if (digisize<8) 
         continue;
       digisize=digi.size();
       if (digisize<8) 
         continue;
       
       depth=digi.id().depth();
       iphi=digi.id().iphi();
       ieta=digi.id().ieta();
       binEta=CalcEtaBin(HcalForward, ieta, depth);
       // 'value' is the average pedestal over the 8 time slices.
       // In the CoarsePedestal client, we will divide the summed value by Nevents (in underflow bin)
       // in order to calculate average pedestals.
       value=0;
       for (unsigned int i=0;i<8;++i)
         value+=digi.sample(i).adc()/8.;
       pedestalsum_[binEta][iphi-1][depth-1]+=value;
       ++pedestalocc_[binEta][iphi-1][depth-1];
     } // for (HFDigiCollection)
     } // if (HFpresent_)
 
   return;
 } // void HcalCoarsePedestalMonitor::processEvent(...)

void HcalCoarsePedestalMonitor::reset ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 337 of file HcalCoarsePedestalMonitor.cc.

References CoarsePedestalsOccByDepth, CoarsePedestalsSumByDepth, EtaPhiHists::Reset(), and zeroCounters().

Referenced by bookHistograms(), and setup().

 {
   // then reset the MonitorElements
   zeroCounters();
   CoarsePedestalsSumByDepth.Reset();
   CoarsePedestalsOccByDepth.Reset();
   return;
 }

void HcalCoarsePedestalMonitor::setup ( DQMStore::IBooker & ib )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 60 of file HcalCoarsePedestalMonitor.cc.

References ADCDiffThresh_, DQMStore::IBooker::bookFloat(), DQMStore::IBooker::bookInt(), CoarsePedestalsOccByDepth, CoarsePedestalsSumByDepth, gather_cfg::cout, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, excludeHORing2_, MonitorElement::Fill(), i, minEvents_, mergeVDriftHistosByStation::name, reset(), DQMStore::IBooker::setCurrentFolder(), HcalBaseDQMonitor::setup(), HcalBaseDQMonitor::SetupEtaPhiHists(), and HcalBaseDQMonitor::subdir_.

Referenced by bookHistograms().

 {
   // Call base class setup
   HcalBaseDQMonitor::setup(ib);
 
   /******* Set up all histograms  ********/
   if (debug_>1)
     std::cout <<"<HcalCoarsePedestalMonitor::setup>  Setting up histograms"<<std::endl;
 
   std::ostringstream name;
   ib.setCurrentFolder(subdir_ +"CoarsePedestalSumPlots");
   SetupEtaPhiHists(ib,CoarsePedestalsSumByDepth,"Coarse Pedestal Summed Map","");
   SetupEtaPhiHists(ib,CoarsePedestalsOccByDepth,"Coarse Pedestal Occupancy Map","");
   for (unsigned int i=0;i<CoarsePedestalsSumByDepth.depth.size();++i)
     (CoarsePedestalsSumByDepth.depth[i]->getTH2F())->SetOption("colz");
   for (unsigned int i=0;i<CoarsePedestalsOccByDepth.depth.size();++i)
     (CoarsePedestalsOccByDepth.depth[i]->getTH2F())->SetOption("colz");
 
 
   ib.setCurrentFolder(subdir_+"CoarsePedestal_parameters");
   MonitorElement* ADCDiffThresh = ib.bookFloat("ADCdiff_Problem_Threshold");
   ADCDiffThresh->Fill(ADCDiffThresh_);
   MonitorElement* minevents = ib.bookInt("minEventsNeededForPedestalCalculation");
   minevents->Fill(minEvents_);
   MonitorElement* excludeHORing2 = ib.bookInt("excludeHORing2");
   if (excludeHORing2_==true)
     excludeHORing2->Fill(1);
   else
     excludeHORing2->Fill(0);
 
   this->reset();
   return;
 } // void HcalCoarsePedestalMonitor::setup()

void HcalCoarsePedestalMonitor::zeroCounters ( )

private

Definition at line 346 of file HcalCoarsePedestalMonitor.cc.

References i, j, relval_2017::k, pedestalocc_, and pedestalsum_.

Referenced by reset().

 {
   for (int i=0;i<85;++i)
     for (int j=0;j<72;++j)
       for (int k=0;k<4;++k)
     {
       pedestalsum_[i][j][k]=0;
       pedestalocc_[i][j][k]=0;
     }
 }