#include <HcalZDCMonitor.h>

Inheritance diagram for HcalZDCMonitor:

Public Member Functions
void	endLuminosityBlock (void)

	HcalZDCMonitor ()

void	processEvent (const ZDCDigiCollection &digi, const ZDCRecHitCollection &rechit, const HcalUnpackerReport &report)

void	reset ()

void	setup (const edm::ParameterSet &ps, DQMStore::IBooker &ib)

	~HcalZDCMonitor ()

Public Member Functions inherited from HcalBaseMonitor
void	beginLuminosityBlock (int lb)

virtual void	beginRun ()

virtual void	clearME ()

virtual void	done ()

void	endLuminosityBlock ()

bool	getDiagnostics () const

int	getVerbosity () const

	HcalBaseMonitor ()

void	hideKnownBadCells ()

virtual void	periodicReset ()

void	processEvent ()

void	setDiagnostics (bool myval)

void	setMinMaxHists1D (std::vector< MonitorElement * > &hh, double min, double max)

void	setMinMaxHists2D (std::vector< MonitorElement * > &hh, double min, double max)

void	setupDepthHists1D (DQMStore::IBooker &ib, MonitorElement &h, std::vector< MonitorElement > &hh, std::string Name, std::string Units, int lowbound, int highbound, int Nbins)

void	setupDepthHists1D (DQMStore::IBooker &ib, std::vector< MonitorElement * > &hh, std::string Name, std::string Units, int lowbound, int highbound, int Nbins)

void	setupDepthHists2D (DQMStore::IBooker &ib, MonitorElement &h, std::vector< MonitorElement > &hh, std::string Name, std::string Units)

void	setupDepthHists2D (DQMStore::IBooker &ib, std::vector< MonitorElement * > &hh, std::string Name, std::string Units)

void	setupDepthHists2D (DQMStore::IBooker &ib, MonitorElement &h, std::vector< MonitorElement > &hh, std::string Name, std::string Units, int nbinsx, int lowboundx, int highboundx, int nbinsy, int lowboundy, int highboundy)

void	setupDepthHists2D (DQMStore::IBooker &ib, std::vector< MonitorElement * > &hh, std::string Name, std::string Units, int nbinsx, int lowboundx, int highboundx, int nbinsy, int lowboundy, int highboundy)

void	SetupEtaPhiHists (DQMStore::IBooker &ib, MonitorElement *&h, EtaPhiHists &hh, std::string Name, std::string Units)

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

void	setVerbosity (int verb)

bool	vetoCell (HcalDetId &id)

virtual	~HcalBaseMonitor ()

Private Member Functions
double	getTime (const std::vector< double > &fData, unsigned int ts_min, unsigned int ts_max, double &fSum)

Private Attributes
bool	ChannelHasDigiError [18]

float	ChannelRatio [18]

std::vector< double >	ChannelWeighting_

bool	checkZDC_

int	ColdChannelCounter [18]

int	DeadChannelCounter [18]

bool	DeadChannelError [18]

bool	DigiErrorCAPID [18]

bool	DigiErrorDVER [18]

int	EventCounter

MonitorElement *	EventsVsLS

MonitorElement *	h_2D_charge

MonitorElement *	h_2D_RecHitEnergy

MonitorElement *	h_2D_RecHitTime

MonitorElement *	h_2D_saturation

MonitorElement *	h_2D_TSMean

MonitorElement *	h_ZDCM_EM_Charge [5]

MonitorElement *	h_ZDCM_EM_Pulse [5]

MonitorElement *	h_ZDCM_EM_RecHitEnergy [5]

MonitorElement *	h_ZDCM_EM_RecHitTiming [5]

MonitorElement *	h_ZDCM_EM_TSMean [5]

MonitorElement *	h_ZDCM_HAD_Charge [4]

MonitorElement *	h_ZDCM_HAD_Pulse [4]

MonitorElement *	h_ZDCM_HAD_RecHitEnergy [4]

MonitorElement *	h_ZDCM_HAD_RecHitTiming [4]

MonitorElement *	h_ZDCM_HAD_TSMean [4]

MonitorElement *	h_ZDCP_EM_Charge [5]

MonitorElement *	h_ZDCP_EM_Pulse [5]

MonitorElement *	h_ZDCP_EM_RecHitEnergy [5]

MonitorElement *	h_ZDCP_EM_RecHitTiming [5]

MonitorElement *	h_ZDCP_EM_TSMean [5]

MonitorElement *	h_ZDCP_HAD_Charge [4]

MonitorElement *	h_ZDCP_HAD_Pulse [4]

MonitorElement *	h_ZDCP_HAD_RecHitEnergy [4]

MonitorElement *	h_ZDCP_HAD_RecHitTiming [4]

MonitorElement *	h_ZDCP_HAD_TSMean [4]

bool	HotChannelError [18]

int	ievt_

std::vector< double >	MaxErrorRates_

MonitorElement *	meEVT_

int	NLumiBlocks_

int	NumBadZDC

MonitorElement *	NZDC_QualityIndexVsLB_

int	OfflineColdThreshold_

int	OfflineDeadThreshold_

int	OnlineColdThreshold_

int	OnlineDeadThreshold_

MonitorElement *	ProblemsVsLB_ZDC

MonitorElement *	PZDC_QualityIndexVsLB_

int	TotalChannelErrors [18]

MonitorElement *	ZDC_Cold_Channel_Errors

MonitorElement *	ZDC_ColdChannelErrorsVsLS

MonitorElement *	ZDC_Dead_Channel_Errors

MonitorElement *	ZDC_DeadChannelErrorsVsLS

MonitorElement *	ZDC_Digi_Errors

MonitorElement *	ZDC_DigiErrors_CAPID

MonitorElement *	ZDC_DigiErrors_DVER

MonitorElement *	ZDC_DigiErrorsVsLS

MonitorElement *	ZDC_Hot_Channel_Errors

MonitorElement *	ZDC_HotChannelErrorsVsLS

MonitorElement *	ZDC_TotalChannelErrors

Additional Inherited Members
Protected Member Functions inherited from HcalBaseMonitor
void	LumiBlockUpdate (int lb)

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

std::vector< std::string >	badCells_

std::string	baseFolder_

bool	checkHB_

bool	checkHE_

bool	checkHF_

bool	checkHO_

int	checkNevents_

edm::CPUTimer	cpu_timer

bool	dump2database

int	etaBins_

double	etaMax_

double	etaMin_

int	fVerbosity

int	ievt_

bool	LBprocessed_

int	levt_

int	lumiblock

bool	makeDiagnostics

MonitorElement *	meEVT_

MonitorElement *	meTOTALEVT_

double	minErrorFlag_

int	Nlumiblocks_

int	NumBadHB

int	NumBadHE

int	NumBadHF

int	NumBadHO

int	oldlumiblock

bool	Online_

int	phiBins_

double	phiMax_

double	phiMin_

MonitorElement *	ProblemCells

EtaPhiHists	ProblemCellsByDepth

MonitorElement *	ProblemsVsLB

MonitorElement *	ProblemsVsLB_HB

MonitorElement *	ProblemsVsLB_HBHEHF

MonitorElement *	ProblemsVsLB_HE

MonitorElement *	ProblemsVsLB_HF

MonitorElement *	ProblemsVsLB_HO

int	resetNevents_

std::string	rootFolder_

bool	showTiming

int	tevt_

Detailed Description

Author: S.Sen

Definition at line 22 of file HcalZDCMonitor.h.

Constructor & Destructor Documentation

HcalZDCMonitor::HcalZDCMonitor ( )

Definition at line 3 of file HcalZDCMonitor.cc.

                                : TotalChannelErrors{}, DeadChannelCounter{}, ColdChannelCounter{},
                    DeadChannelError{}, HotChannelError{}, DigiErrorCAPID{}, DigiErrorDVER{},
                    ChannelHasDigiError{}
 {
 }

HcalZDCMonitor::~HcalZDCMonitor ( )

Definition at line 8 of file HcalZDCMonitor.cc.

8 {

9 }

Member Function Documentation

void HcalZDCMonitor::endLuminosityBlock ( void )

Definition at line 743 of file HcalZDCMonitor.cc.

References ChannelRatio, ChannelWeighting_, ColdChannelCounter, DeadChannelCounter, HcalBaseMonitor::endLuminosityBlock(), EventCounter, EventsVsLS, MonitorElement::Fill(), h_2D_charge, h_2D_RecHitEnergy, h_2D_RecHitTime, h_2D_TSMean, h_ZDCM_EM_Charge, h_ZDCM_EM_RecHitEnergy, h_ZDCM_EM_RecHitTiming, h_ZDCM_EM_TSMean, h_ZDCM_HAD_Charge, h_ZDCM_HAD_RecHitEnergy, h_ZDCM_HAD_RecHitTiming, h_ZDCM_HAD_TSMean, h_ZDCP_EM_Charge, h_ZDCP_EM_RecHitEnergy, h_ZDCP_EM_RecHitTiming, h_ZDCP_EM_TSMean, h_ZDCP_HAD_Charge, h_ZDCP_HAD_RecHitEnergy, h_ZDCP_HAD_RecHitTiming, h_ZDCP_HAD_TSMean, i, HcalBaseMonitor::lumiblock, MaxErrorRates_, NZDC_QualityIndexVsLB_, OfflineColdThreshold_, OfflineDeadThreshold_, HcalBaseMonitor::Online_, OnlineColdThreshold_, OnlineDeadThreshold_, PZDC_QualityIndexVsLB_, MonitorElement::setBinContent(), TotalChannelErrors, ZDC_Cold_Channel_Errors, ZDC_ColdChannelErrorsVsLS, ZDC_Dead_Channel_Errors, ZDC_DeadChannelErrorsVsLS, and ZDC_TotalChannelErrors.

Referenced by ZDCMonitorModule::endLuminosityBlock().

 {
     bool HadLumiError[18]={false, false, false, false, false, false, false, false, false, false, false, false, false, false, false, false, false, false};
 
     EventsVsLS->Fill(lumiblock, EventCounter);
 
     if(Online_ == false)
     {//checks if DQM is in OFFLINE MODE
         for (int i=0;i<18;++i) {
             ChannelRatio[i]=(TotalChannelErrors[i])*1./EventCounter;
             if (ChannelRatio[i] <= MaxErrorRates_[i]) {
                 if (i<9) {
                     PZDC_QualityIndexVsLB_->Fill(lumiblock,ChannelWeighting_[i]);
                 } else {
                     NZDC_QualityIndexVsLB_->Fill(lumiblock,ChannelWeighting_[i]);
                 }
                 if (ColdChannelCounter[i] >= OfflineColdThreshold_)//Begin Cold Error Plots
                 {
                     ZDC_Cold_Channel_Errors->Fill(i/9,i%9,ColdChannelCounter[i]);
                     ZDC_ColdChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,ColdChannelCounter[i]);//Can change this between 1, or the amount of errors (Currently the latter)
                     ColdChannelCounter[i]=0;
                     HadLumiError[i]=true;
                 }//END OF Cold Error Plot
                 if (DeadChannelCounter[i] >= OfflineDeadThreshold_)
                 {//Begin Dead Error Plots
                     ZDC_Dead_Channel_Errors->Fill(i/9,i%9,DeadChannelCounter[i]);
                     ZDC_DeadChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,DeadChannelCounter[i]); //Could fill this with 1 or total dead errors (which is currently done)
                     DeadChannelCounter[i]=0;
                     HadLumiError[i]=true;
                 }//END OF Dead Channel Plots
                 if (HadLumiError[i]==true)
                 {//Removing the QI for Dead of Cold Channels
                     if (i<9) {
                         PZDC_QualityIndexVsLB_->Fill(lumiblock,(-1*ChannelWeighting_[i]));
                     } else {
                         NZDC_QualityIndexVsLB_->Fill(lumiblock,(-1*ChannelWeighting_[i]));
                     }
                 }//END OF QI Removal
             }//END OF ChannelRatio[i]<=MaxErrorRates_[i]
             else {
                 //This part only happens if ChannelRatio[i] > MaxErrorRates_[i]...
                 //Above you notice the QI plots become 'un-filled'. 
                 //If the plot was never filled to begin with, 
                 //then we do not want to remove from the plot  causing there to be a negative QI
                 if (ColdChannelCounter[i] >= OfflineColdThreshold_)
                 {//Begin Cold Error Plots
                     ZDC_Cold_Channel_Errors->Fill(i/9,i%9,ColdChannelCounter[i]);
                     ZDC_ColdChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,ColdChannelCounter[i]);
                     ColdChannelCounter[i]=0;
                 }//END OF Cold Error Plot
                 if (DeadChannelCounter[i] >= OfflineDeadThreshold_)
                 {//Begin Dead Error plots
                     ZDC_Dead_Channel_Errors->Fill(i/9,i%9,DeadChannelCounter[i]);
                     ZDC_DeadChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,DeadChannelCounter[i]);
                     DeadChannelCounter[i]=0;
                 }//END OF Dead Error Plots
             }
         }//END OF FOR LOOP
     }//END OF DQM OFFLINE PART
 
     if(Online_ == true)
     {//checks if DQM is in ONLINE MODE
         for (int i=0;i<18;++i) {
             ChannelRatio[i]=(TotalChannelErrors[i])*1./EventCounter;
             if (ChannelRatio[i] <= MaxErrorRates_[i]) {
                 if (i<9) {
                     PZDC_QualityIndexVsLB_->Fill(lumiblock,ChannelWeighting_[i]);
                 } else {
                     NZDC_QualityIndexVsLB_->Fill(lumiblock,ChannelWeighting_[i]);
                 }
                 if (ColdChannelCounter[i] >= OnlineColdThreshold_)//Begin Cold Error Plots
                 {
                     ZDC_Cold_Channel_Errors->Fill(i/9,i%9,ColdChannelCounter[i]);
                     ZDC_ColdChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,ColdChannelCounter[i]);//Can change this between 1, or the amount of errors (Currently the latter)
                     ColdChannelCounter[i]=0;
                     HadLumiError[i]=true;
                 }//END OF Cold Error Plot
                 if (DeadChannelCounter[i] >= OnlineDeadThreshold_)
                 {//Begin Dead Error Plots
                     ZDC_Dead_Channel_Errors->Fill(i/9,i%9,DeadChannelCounter[i]);
                     ZDC_DeadChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,DeadChannelCounter[i]); //Could fill this with 1 or total dead errors (which is currently done)
                     DeadChannelCounter[i]=0;
                     HadLumiError[i]=true;
                 }//END OF Dead Channel Plots
                 if (HadLumiError[i]==true)
                 {//Removing the QI for Dead of Cold Channels
                     if (i<9) {
                         PZDC_QualityIndexVsLB_->Fill(lumiblock,(-1*ChannelWeighting_[i]));
                     } else {
                         NZDC_QualityIndexVsLB_->Fill(lumiblock,(-1*ChannelWeighting_[i]));
                     }
                 }//END OF QI Removal
             }//END OF ChannelRatio[i]<=MaxErrorRates_[i]
             else {
                 //This part only happens if ChannelRatio[i] > MaxErrorRates_[i]...
                 //Above you notice the QI plots become 'un-filled'. 
                 //If the plot was never filled to begin with, 
                 //then we do not want to remove from the plot  causing there to be a negative QI
                 if (ColdChannelCounter[i] >= OnlineColdThreshold_)
                 {//Begin Cold Error Plots
                     ZDC_Cold_Channel_Errors->Fill(i/9,i%9,ColdChannelCounter[i]);
                     ZDC_ColdChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,ColdChannelCounter[i]);
                     ColdChannelCounter[i]=0;
                 }//END OF Cold Error Plot
                 if (DeadChannelCounter[i] >= OnlineDeadThreshold_)
                 {//Begin Dead Error plots
                     ZDC_Dead_Channel_Errors->Fill(i/9,i%9,DeadChannelCounter[i]);
                     ZDC_DeadChannelErrorsVsLS->Fill(lumiblock,1);
                     ZDC_TotalChannelErrors->Fill(i/9,i%9,DeadChannelCounter[i]);
                     DeadChannelCounter[i]=0;
                 }//END OF Dead Error Plots
             }//end of ChannelRatio[i] > MaxErrorRates_[i] part
         }//END OF FOR LOOP
     }//END OF DQM ONLINE PART
 
 
 
 
     HcalBaseMonitor::endLuminosityBlock();
 
     for (int i = 0; i < 5; ++i) {   // EM Channels
         // ZDC Plus
         h_2D_charge->setBinContent(1, i + 1, h_ZDCP_EM_Charge[i]->getMean());
         h_2D_TSMean->setBinContent(1, i + 1, h_ZDCP_EM_TSMean[i]->getMean());
         h_2D_RecHitEnergy->setBinContent(1, i + 1, h_ZDCP_EM_RecHitEnergy[i]->getMean());
         h_2D_RecHitTime->setBinContent(1, i + 1, h_ZDCP_EM_RecHitTiming[i]->getMean());
         // ZDC Minus
         h_2D_charge->setBinContent(2, i + 1, h_ZDCM_EM_Charge[i]->getMean());
         h_2D_TSMean->setBinContent(2, i + 1, h_ZDCM_EM_TSMean[i]->getMean());
         h_2D_RecHitEnergy->setBinContent(2, i + 1, h_ZDCM_EM_RecHitEnergy[i]->getMean());
         h_2D_RecHitTime->setBinContent(2, i + 1, h_ZDCM_EM_RecHitTiming[i]->getMean());
     }
 
     for (int i = 0; i < 4; ++i) {   // HAD channels
         // ZDC Plus
         h_2D_charge->setBinContent(1, i + 6, h_ZDCP_HAD_Charge[i]->getMean());
         h_2D_TSMean->setBinContent(1, i + 6, h_ZDCP_HAD_TSMean[i]->getMean());
         h_2D_RecHitEnergy->setBinContent(1, i + 6, h_ZDCP_HAD_RecHitEnergy[i]->getMean());
         h_2D_RecHitTime->setBinContent(1, i + 6, h_ZDCP_HAD_RecHitTiming[i]->getMean());
         // ZDC Minus
         //h_ZDCM_HAD_Pulse[i]->Scale(10. / h_ZDCM_HAD_Pulse[i]->getEntries());
         h_2D_charge->setBinContent(2, i + 6, h_ZDCM_HAD_Charge[i]->getMean());
         h_2D_TSMean->setBinContent(2, i + 6, h_ZDCM_HAD_TSMean[i]->getMean());
         h_2D_RecHitEnergy->setBinContent(2, i + 6, h_ZDCM_HAD_RecHitEnergy[i]->getMean());
         h_2D_RecHitTime->setBinContent(2, i + 6, h_ZDCM_HAD_RecHitTiming[i]->getMean());
     }
 } // void HcalZDCMonitor::endLuminosityBlock()

double HcalZDCMonitor::getTime	(	const std::vector< double > &	fData,
		unsigned int	ts_min,
		unsigned int	ts_max,
		double &	fSum
	)

private

Definition at line 719 of file HcalZDCMonitor.cc.

Referenced by processEvent().

                                                                                                                      {
     double weightedTime = 0.;
     double SumT = 0.;
     double Time = -999.;
     double digiThreshf = 99.5;
 
     for (unsigned int ts=ts_min; ts<=ts_max; ++ts) {
         if (fData[ts] > digiThreshf){
             weightedTime += ts * fData[ts];
             SumT += fData[ts];
         }
     }
 
     if (SumT > 0.) {
         Time = weightedTime / SumT;
     }
 
     fSum = SumT;
 
     return Time;
 
 } //double HcalZDCMonitor::getTime()

void HcalZDCMonitor::processEvent	(	const ZDCDigiCollection &	digi,
		const ZDCRecHitCollection &	rechit,
		const HcalUnpackerReport &	report
	)

DEAD CELL ERROR/////////////////////////// Right now we are simply checking that the digi exists

Definition at line 389 of file HcalZDCMonitor.cc.

References ecalMGPA::adc(), HcalUnpackerReport::bad_quality_begin(), HcalUnpackerReport::bad_quality_end(), edm::SortedCollection< T, SORT >::begin(), DetId::Calo, HcalQIESample::capid(), ChannelHasDigiError, ColdChannelCounter, gather_cfg::cout, HcalBaseMonitor::cpu_timer, DeadChannelCounter, DeadChannelError, DigiErrorCAPID, DigiErrorDVER, HcalQIESample::dv(), edm::SortedCollection< T, SORT >::end(), HcalQIESample::er(), EventCounter, MonitorElement::Fill(), HcalBaseMonitor::fVerbosity, getTime(), h_2D_saturation, h_ZDCM_EM_Charge, h_ZDCM_EM_Pulse, h_ZDCM_EM_RecHitEnergy, h_ZDCM_EM_RecHitTiming, h_ZDCM_EM_TSMean, h_ZDCM_HAD_Charge, h_ZDCM_HAD_Pulse, h_ZDCM_HAD_RecHitEnergy, h_ZDCM_HAD_RecHitTiming, h_ZDCM_HAD_TSMean, h_ZDCP_EM_Charge, h_ZDCP_EM_Pulse, h_ZDCP_EM_RecHitEnergy, h_ZDCP_EM_RecHitTiming, h_ZDCP_EM_TSMean, h_ZDCP_HAD_Charge, h_ZDCP_HAD_Pulse, h_ZDCP_HAD_RecHitEnergy, h_ZDCP_HAD_RecHitTiming, h_ZDCP_HAD_TSMean, HotChannelError, i, ievt_, HcalBaseMonitor::lumiblock, meEVT_, edm::CPUTimer::reset(), ZDCDataFrame::sample(), HcalBaseMonitor::showTiming, ZDCDataFrame::size(), edm::SortedCollection< T, SORT >::size(), edm::CPUTimer::start(), HcalZDCDetId::SubdetectorId, TotalChannelErrors, ZDC_Cold_Channel_Errors, ZDC_Dead_Channel_Errors, ZDC_Digi_Errors, ZDC_DigiErrors_CAPID, ZDC_DigiErrors_DVER, ZDC_DigiErrorsVsLS, ZDC_Hot_Channel_Errors, ZDC_HotChannelErrorsVsLS, ZDC_TotalChannelErrors, HcalZDCDetId::zside(), and ecaldqm::zside().

Referenced by ZDCMonitorModule::analyze().

                                                                                                                                     {
     if (fVerbosity > 0)
         std::cout << "<HcalZDCMonitor::processEvent> Processing Event..." << std::endl;
     if (showTiming)
     {
         cpu_timer.reset();
         cpu_timer.start();
     }
     ++ievt_;
     meEVT_->Fill(ievt_);
 
     //--------------------------------------
     // ZDC Digi part
     //--------------------------------------
     double fSum = 0.;
     std::vector<double> fData;
     double digiThresh = 99.5; //corresponds to 40 ADC counts
     //int digiThreshADC = 40;
     int digiSaturation = 127;
     //double ZDCQIEConst = 2.6;
 
     if (digi.size()>0) {
         ZDC_Digi_Errors->Fill(-1,-1,1);
         ZDC_Hot_Channel_Errors->Fill(-1,-1,1);
         ZDC_TotalChannelErrors->Fill(-1,-1,1);
         ZDC_Cold_Channel_Errors->Fill(-1,-1,1);
         ZDC_Dead_Channel_Errors->Fill(-1,-1,1);
         EventCounter+=1;
     }
 
     for (int i=0;i<18;++i) {
         ChannelHasDigiError[i]=false;
         DigiErrorDVER[i]=false;
         DigiErrorCAPID[i]=false;
         HotChannelError[i]=false;
         DeadChannelError[i]=true;
     }
 
     typedef std::vector<DetId> DetIdVector;
 
     for (DetIdVector::const_iterator baddigi_iter=report.bad_quality_begin();
             baddigi_iter != report.bad_quality_end();
             ++baddigi_iter)
     {
         DetId id(baddigi_iter->rawId());
         if (id.det()==DetId::Calo && id.subdetId()==HcalZDCDetId::SubdetectorId)
         {
             HcalZDCDetId id(baddigi_iter->rawId());
             int iSide = id.zside();
             int iSection = id.section();
             int iChannel = id.channel();
             if(iSection==1 || iSection==2){
                 ChannelHasDigiError[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=true;
                 DeadChannelError[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=false;
                 //do stuff
             }//end of if i(Section==1 || iSection==2)
         }
         else continue;
 
     }//end unpacker section
 
 
     // ChannelHasError[18]:  [0-8] are iSide=1, [9-17] are iSide=2
     //  First 5 bins ([0-4],[9-13]) are EM bins
     //  Last 4 bins are HAD bins
 
     for (ZDCDigiCollection::const_iterator digi_iter = digi.begin();
             digi_iter != digi.end(); ++digi_iter)
     {
         const ZDCDataFrame digi = (const ZDCDataFrame) (*digi_iter);
         //HcalZDCDetId id(digi_iter->id());
         int iSide = digi_iter->id().zside();
         int iSection = digi_iter->id().section();
         int iChannel = digi_iter->id().channel();
 
         unsigned int fTS = digi_iter->size();
         while (fData.size()<fTS)
             fData.push_back(-999);
         while (fData.size()>fTS)
             fData.pop_back(); // delete last elements
 
         if (iSection==1 || iSection==2)
         {
 
 
             DeadChannelError[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=false;
 
             int iCapID=27;
             int iCapIDPrevious=27;
             int HotCounter=0;
             int ColdCounter=0;
 
             for (int iTS=0; iTS<digi.size(); ++iTS) //looping over all ten timeslices
             {
                 if (digi[iTS].adc()==127) HotCounter+=1;
                 else HotCounter=0;//require 3 consecutive saturated Time Slices in a single channel in a single event
                 if (HotCounter >= 3) HotChannelError[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=true;
 
 
                 if (digi[iTS].adc()<=10) ColdCounter+=1;
                 if (ColdCounter==10)
                 {
                     ColdChannelCounter[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]+=1;
                 }
 
 
                 if ((ChannelHasDigiError[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=false))
                 {
                     iCapID=digi.sample(iTS).capid();
                     if (iTS>0) iCapIDPrevious=digi.sample(iTS-1).capid();
 
                     if (digi.sample(iTS).dv()==0 || digi.sample(iTS).er()==1)
                     {
                         ChannelHasDigiError[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=true;
                         DigiErrorDVER[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=true;
                         break;
                     }
                     else
                     {
                         if (iTS==0) continue;
                         else
                         {
                             if ((iCapID-iCapIDPrevious+4)%4!=1)
                             {
                                 ChannelHasDigiError[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=true;
                                 DigiErrorCAPID[(9*((1-iSide)/2))+(iChannel-1)+(5*((iSection-1)%2))]=true;
                                 break;
                             }//end of capid rotation check
                         }//checking if TS!=0
                     } // end of the check for dv/er
                 }//END of unpacker double check
             } // end of TS loop
 
 
 
             fSum = 0.;
             bool saturated = false;
             for (unsigned int i = 0; i < fTS; ++i)
             {
                 //fData[i]=digi[i].nominal_fC() * ZDCQIEConst;
                 fData[i]=digi[i].nominal_fC();
                 if (digi[i].adc()==digiSaturation){
                     saturated=true;
                 }
             }
 
             double fTSMean = 0;
             if (fData.size()>6)
                 fTSMean = getTime(fData, 4, 6, fSum); // tsmin = 4, tsmax = 6.
             //std::cout << "Side= " << iSide << " Section= " << iSection << " Channel= " << iChannel << "\tCharge\t" << fSum <<std::endl;
             if (saturated==true){
                 h_2D_saturation->Fill(iSide==1?0:1,iSection==1?iChannel-1:iChannel+4,1);
             }
 
             if (iSection == 1)
             {    // EM
                 if (iSide == 1) {   // Plus
                     for (unsigned int i = 0; i < fTS; ++i) {
                         if (fData[i] > digiThresh) h_ZDCP_EM_Pulse[iChannel - 1]->Fill(i, fData[i]);
                     }
                     if (fSum > digiThresh) {
                         h_ZDCP_EM_Charge[iChannel - 1]->Fill(fSum);
                         h_ZDCP_EM_TSMean[iChannel - 1]->Fill(fTSMean);
                         //std::cout<< "fSum " << fSum << " fTSMean " << fTSMean <<std::endl;
                     }
                 } // Plus
                 if (iSide == -1) {  // Minus
                     for (unsigned int i = 0; i < fTS; ++i) {
                         if (fData[i] > digiThresh) h_ZDCM_EM_Pulse[iChannel - 1]->Fill(i, fData[i]);
                     }
                     if (fSum > digiThresh) {
                         h_ZDCM_EM_Charge[iChannel - 1]->Fill(fSum);
                         h_ZDCM_EM_TSMean[iChannel - 1]->Fill(fTSMean);
                     }
                 } // Minus
             }// EM
 
             else if (iSection == 2)
             {    // HAD
                 if (iSide == 1) {   // Plus
                     for (unsigned int i = 0; i < fTS; ++i) {
                         if (fData[i] > digiThresh) h_ZDCP_HAD_Pulse[iChannel - 1]->Fill(i, fData[i]);
                     }
                     if (fSum > digiThresh) {
                         h_ZDCP_HAD_Charge[iChannel - 1]->Fill(fSum);
                         h_ZDCP_HAD_TSMean[iChannel - 1]->Fill(fTSMean);
                     }
                 } // Plus
                 if (iSide == -1) {  // Minus
                     for (unsigned int i = 0; i < fTS; ++i) {
                         if (fData[i] > digiThresh) h_ZDCM_HAD_Pulse[iChannel - 1]->Fill(i, fData[i]);
                     }
                     if (fSum > digiThresh) {
                         h_ZDCM_HAD_Charge[iChannel - 1]->Fill(fSum);
                         h_ZDCM_HAD_TSMean[iChannel - 1]->Fill(fTSMean);
                     }
                 }// Minus
             } // HAD
         }//end of if (iSection==1 || iSection==2)
     } // loop on zdc digi collection
 
 
 
     // Fill Fall 2012 histograms
     int numdigierrors=0;
     int numhoterrors=0;
 
     for (int i = 0; i<18; i++){
         if (ChannelHasDigiError[i]==true)
         {
             ++numdigierrors;
             ZDC_Digi_Errors->Fill(i/9,i%9,1);
         }
         if (DigiErrorDVER[i]==true)
         {
             ZDC_DigiErrors_DVER->Fill(i/9,i%9,1);
         }
         if (DigiErrorCAPID[i]==true)
         {
             ZDC_DigiErrors_CAPID->Fill(i/9,i%9,1);
         }
         if(HotChannelError[i]==true)
         {
             ++numhoterrors;
             ZDC_Hot_Channel_Errors->Fill(i/9,(i%9),1);
         }
         if(DeadChannelError[i]==true)
         {
             DeadChannelCounter[i]+=1;
         }
         // If any of the above is true, fill the total channel errors
         if (ChannelHasDigiError[i] || HotChannelError[i])
         {
             ZDC_TotalChannelErrors->Fill(i/9,i%9,1);
             TotalChannelErrors[i]+=1;
         }
     }//end the for i<18 loop
 
 
     if (numdigierrors>0)
         ZDC_DigiErrorsVsLS->Fill(lumiblock,numdigierrors); //could make this ->Fill(currentLS,1) if I want to count each event as only having one error even if multiple channels are in error
     if (numhoterrors>0)
         ZDC_HotChannelErrorsVsLS->Fill(lumiblock,numhoterrors);
     //End of  filling Fall 2012  histograms
 
 
     //--------------------------------------
     // ZDC RecHit part
     //--------------------------------------
     for (ZDCRecHitCollection::const_iterator rechit_iter = rechit.begin();
             rechit_iter != rechit.end(); ++rechit_iter)
     {
         int Side      = (rechit_iter->id()).zside();
         int Section   = (rechit_iter->id()).section();
         int Channel   = (rechit_iter->id()).channel();
         //std::cout << "RecHitEnergy  " << zhit->energy() << "  RecHitTime  " << zhit->time() << std::endl;
 
         if(Section==1)
         { //EM
             if (Side ==1 ){ // Plus
                 h_ZDCP_EM_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
                 h_ZDCP_EM_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
             }
             if (Side == -1 ){ //Minus
                 h_ZDCM_EM_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
                 h_ZDCM_EM_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
             }
         } //EM
         else if(Section==2)
         { //HAD
             if (Side ==1 ){ //Plus
                 h_ZDCP_HAD_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
                 h_ZDCP_HAD_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
             }
             if (Side == -1 ){ //Minus
                 h_ZDCM_HAD_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
                 h_ZDCM_HAD_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
             }
         } // HAD
     } // loop on rechits
 
 } // end of event processing

void HcalZDCMonitor::reset ( void )

Definition at line 10 of file HcalZDCMonitor.cc.

Referenced by ZDCMonitorModule::reset().

10 {

11 }

void HcalZDCMonitor::setup	(	const edm::ParameterSet &	ps,
		DQMStore::IBooker &	ib
	)

virtual

Reimplemented from HcalBaseMonitor.

Definition at line 13 of file HcalZDCMonitor.cc.

References HcalBaseMonitor::baseFolder_, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookInt(), ChannelWeighting_, ColdChannelCounter, gather_cfg::cout, HcalBaseMonitor::cpu_timer, DeadChannelCounter, EventCounter, EventsVsLS, MonitorElement::Fill(), HcalBaseMonitor::fVerbosity, edm::ParameterSet::getParameter(), MonitorElement::getTH2F(), edm::ParameterSet::getUntrackedParameter(), h_2D_charge, h_2D_RecHitEnergy, h_2D_RecHitTime, h_2D_saturation, h_2D_TSMean, h_ZDCM_EM_Charge, h_ZDCM_EM_Pulse, h_ZDCM_EM_RecHitEnergy, h_ZDCM_EM_RecHitTiming, h_ZDCM_EM_TSMean, h_ZDCM_HAD_Charge, h_ZDCM_HAD_Pulse, h_ZDCM_HAD_RecHitEnergy, h_ZDCM_HAD_RecHitTiming, h_ZDCM_HAD_TSMean, h_ZDCP_EM_Charge, h_ZDCP_EM_Pulse, h_ZDCP_EM_RecHitEnergy, h_ZDCP_EM_RecHitTiming, h_ZDCP_EM_TSMean, h_ZDCP_HAD_Charge, h_ZDCP_HAD_Pulse, h_ZDCP_HAD_RecHitEnergy, h_ZDCP_HAD_RecHitTiming, h_ZDCP_HAD_TSMean, i, ievt_, MaxErrorRates_, meEVT_, mergeVDriftHistosByStation::name, NLumiBlocks_, NZDC_QualityIndexVsLB_, OfflineColdThreshold_, OfflineDeadThreshold_, OnlineColdThreshold_, OnlineDeadThreshold_, PZDC_QualityIndexVsLB_, edm::CPUTimer::reset(), HcalBaseMonitor::rootFolder_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), HcalBaseMonitor::setup(), HcalBaseMonitor::showTiming, edm::CPUTimer::start(), indexGen::title, ZDC_Cold_Channel_Errors, ZDC_ColdChannelErrorsVsLS, ZDC_Dead_Channel_Errors, ZDC_DeadChannelErrorsVsLS, ZDC_Digi_Errors, ZDC_DigiErrors_CAPID, ZDC_DigiErrors_DVER, ZDC_DigiErrorsVsLS, ZDC_Hot_Channel_Errors, ZDC_HotChannelErrorsVsLS, and ZDC_TotalChannelErrors.

Referenced by ZDCMonitorModule::bookHistograms().

                                                                            {
     HcalBaseMonitor::setup(ps, ib);
 
     baseFolder_ = rootFolder_ + "ZDCMonitor_Hcal";
 
     const edm::ParameterSet& psZDC(ps.getParameter<edm::ParameterSet>("zdcMonitorTask"));
     NLumiBlocks_           = psZDC.getUntrackedParameter<int>("NLumiBlocks",4000);
     ChannelWeighting_      = psZDC.getUntrackedParameter<std::vector<double>> ("ZDC_ChannelWeighting");
     MaxErrorRates_         = psZDC.getUntrackedParameter<std::vector<double>> ("ZDC_AcceptableChannelErrorRates");
     OfflineColdThreshold_  = psZDC.getUntrackedParameter<int>("ZDC_OfflineColdThreshold");
     OfflineDeadThreshold_  = psZDC.getUntrackedParameter<int>("ZDC_OfflineDeadThreshold");
     OnlineDeadThreshold_   = psZDC.getUntrackedParameter<int>("ZDC_OnlineDeadThreshold");
     OnlineColdThreshold_   = psZDC.getUntrackedParameter<int>("ZDC_OnlineColdThreshold");
 
     for (int i=0;i<18;++i)
     {
         ColdChannelCounter[i]=0;
         DeadChannelCounter[i]=0;
     }
 
     EventCounter=0;
 
     if (showTiming) {
         cpu_timer.reset();
         cpu_timer.start();
     }
 
     if (fVerbosity > 0)
         std::cout << "<HcalZDCMonitor::setup>  Setting up histograms" << std::endl;
 
     if (fVerbosity > 1)
         std::cout << "<HcalZDCMonitor::setup> Getting variable values from cfg files" << std::endl;
 
     // Set initial event # to 0
     ievt_ = 0;
 
     //Histograms
     if (fVerbosity > 1)
         std::cout << "<HcalZDCMonitor::setup>  Setting up Histograms" << std::endl;
 
     ib.setCurrentFolder(baseFolder_);
     meEVT_ = ib.bookInt("ZDC Event Number");
     meEVT_->Fill(ievt_);
     char name[128];
     char title[128];
 
     PZDC_QualityIndexVsLB_ = ib.book1D("PZDC_QualityIndexVSLB","Quality Index for the ZDC+ vs LS; LS; Quality Index", NLumiBlocks_,0,NLumiBlocks_);
     NZDC_QualityIndexVsLB_ = ib.book1D("NZDC_QualityIndexVSLB","Quality Index for the ZDC- vs  LS; LS; Quality Index", NLumiBlocks_,0,NLumiBlocks_);
     EventsVsLS             = ib.book1D("EventsVsLS", "Total Number of Events per LS; LS; # of Events", NLumiBlocks_,0,NLumiBlocks_);
 
     //
     h_2D_charge = ib.book2D("2D_DigiCharge", "Digi Charge (fC)", 2, 0, 2, 9, 0, 9);
     h_2D_charge->setBinLabel(1,"ZDC+",1);
     h_2D_charge->setBinLabel(2,"ZDC-",1);
     h_2D_charge->setBinLabel(1,"EM1",2);
     h_2D_charge->setBinLabel(2,"EM2",2);
     h_2D_charge->setBinLabel(3,"EM3",2);
     h_2D_charge->setBinLabel(4,"EM4",2);
     h_2D_charge->setBinLabel(5,"EM5",2);
     h_2D_charge->setBinLabel(6,"HAD1",2);
     h_2D_charge->setBinLabel(7,"HAD2",2);
     h_2D_charge->setBinLabel(8,"HAD3",2);
     h_2D_charge->setBinLabel(9,"HAD4",2);
 
     h_2D_TSMean = ib.book2D("2D_DigiTiming", "Digi Timing", 2, 0, 2, 9, 0, 9);
     h_2D_TSMean->setBinLabel(1,"ZDC+",1);
     h_2D_TSMean->setBinLabel(2,"ZDC-",1);
     h_2D_TSMean->setBinLabel(1,"EM1",2);
     h_2D_TSMean->setBinLabel(2,"EM2",2);
     h_2D_TSMean->setBinLabel(3,"EM3",2);
     h_2D_TSMean->setBinLabel(4,"EM4",2);
     h_2D_TSMean->setBinLabel(5,"EM5",2);
     h_2D_TSMean->setBinLabel(6,"HAD1",2);
     h_2D_TSMean->setBinLabel(7,"HAD2",2);
     h_2D_TSMean->setBinLabel(8,"HAD3",2);
     h_2D_TSMean->setBinLabel(9,"HAD4",2);
 
     h_2D_RecHitEnergy = ib.book2D("2D_RecHitEnergy", "Rechit Energy", 2, 0, 2, 9, 0, 9);
     h_2D_RecHitEnergy->setBinLabel(1,"ZDC+",1);
     h_2D_RecHitEnergy->setBinLabel(2,"ZDC-",1);
     h_2D_RecHitEnergy->setBinLabel(1,"EM1",2);
     h_2D_RecHitEnergy->setBinLabel(2,"EM2",2);
     h_2D_RecHitEnergy->setBinLabel(3,"EM3",2);
     h_2D_RecHitEnergy->setBinLabel(4,"EM4",2);
     h_2D_RecHitEnergy->setBinLabel(5,"EM5",2);
     h_2D_RecHitEnergy->setBinLabel(6,"HAD1",2);
     h_2D_RecHitEnergy->setBinLabel(7,"HAD2",2);
     h_2D_RecHitEnergy->setBinLabel(8,"HAD3",2);
     h_2D_RecHitEnergy->setBinLabel(9,"HAD4",2);
 
     h_2D_RecHitTime = ib.book2D("2D_RecHitTime", "Rechit Timing", 2, 0, 2, 9, 0, 9);
     h_2D_RecHitTime->setBinLabel(1,"ZDC+",1);
     h_2D_RecHitTime->setBinLabel(2,"ZDC-",1);
     h_2D_RecHitTime->setBinLabel(1,"EM1",2);
     h_2D_RecHitTime->setBinLabel(2,"EM2",2);
     h_2D_RecHitTime->setBinLabel(3,"EM3",2);
     h_2D_RecHitTime->setBinLabel(4,"EM4",2);
     h_2D_RecHitTime->setBinLabel(5,"EM5",2);
     h_2D_RecHitTime->setBinLabel(6,"HAD1",2);
     h_2D_RecHitTime->setBinLabel(7,"HAD2",2);
     h_2D_RecHitTime->setBinLabel(8,"HAD3",2);
     h_2D_RecHitTime->setBinLabel(9,"HAD4",2);
 
     h_2D_saturation = ib.book2D("h_2D_QIE", "Saturation Check", 2, 0, 2, 9, 0, 9);
     h_2D_saturation->setBinLabel(1,"ZDC+",1);
     h_2D_saturation->setBinLabel(2,"ZDC-",1);
     h_2D_saturation->setBinLabel(1,"EM1",2);
     h_2D_saturation->setBinLabel(2,"EM2",2);
     h_2D_saturation->setBinLabel(3,"EM3",2);
     h_2D_saturation->setBinLabel(4,"EM4",2);
     h_2D_saturation->setBinLabel(5,"EM5",2);
     h_2D_saturation->setBinLabel(6,"HAD1",2);
     h_2D_saturation->setBinLabel(7,"HAD2",2);
     h_2D_saturation->setBinLabel(8,"HAD3",2);
     h_2D_saturation->setBinLabel(9,"HAD4",2);
 
 
     // digi errors
     ib.setCurrentFolder(baseFolder_ + "/Errors/Digis");
         ZDC_Digi_Errors = ib.book2D("ZDC_Digi_Errors", "Raw Number of Digi Errors Per ZDC Channel", 2, 0, 2, 9, 0, 9);
         ZDC_Digi_Errors->setBinLabel(1,"ZDC+",1);
         ZDC_Digi_Errors->setBinLabel(2,"ZDC-",1);
     ZDC_Digi_Errors->setBinLabel(1,"EM1",2);
         ZDC_Digi_Errors->setBinLabel(2,"EM2",2);
         ZDC_Digi_Errors->setBinLabel(3,"EM3",2);
         ZDC_Digi_Errors->setBinLabel(4,"EM4",2);
         ZDC_Digi_Errors->setBinLabel(5,"EM5",2);
         ZDC_Digi_Errors->setBinLabel(6,"HAD1",2);
         ZDC_Digi_Errors->setBinLabel(7,"HAD2",2);
         ZDC_Digi_Errors->setBinLabel(8,"HAD3",2);
         ZDC_Digi_Errors->setBinLabel(9,"HAD4",2);
         ZDC_Digi_Errors->getTH2F()->SetOption("coltext");
 
 
     ZDC_DigiErrorsVsLS = ib.book1D("ZDC_DigiErrorsVsLS","Total Number of (Digi) Errors found in the ZDCs vs. Lumi Section;LS;# errors",NLumiBlocks_,0,NLumiBlocks_);
 
     ib.setCurrentFolder(baseFolder_ + "/Errors/Digis/DigiErrorCauses");
     ZDC_DigiErrors_DVER = ib.book2D("ZDC_DigiErrors_DVER","Raw Number of Digi Errors Caused by Finding .dv()=0 or .er()=1",2,0,2,9,0,9);
     ZDC_DigiErrors_DVER->setBinLabel(1,"ZDC+",1);
         ZDC_DigiErrors_DVER->setBinLabel(2,"ZDC-",1);
     ZDC_DigiErrors_DVER->setBinLabel(1,"EM1",2);
         ZDC_DigiErrors_DVER->setBinLabel(2,"EM2",2);
         ZDC_DigiErrors_DVER->setBinLabel(3,"EM3",2);
         ZDC_DigiErrors_DVER->setBinLabel(4,"EM4",2);
         ZDC_DigiErrors_DVER->setBinLabel(5,"EM5",2);
         ZDC_DigiErrors_DVER->setBinLabel(6,"HAD1",2);
         ZDC_DigiErrors_DVER->setBinLabel(7,"HAD2",2);
         ZDC_DigiErrors_DVER->setBinLabel(8,"HAD3",2);
         ZDC_DigiErrors_DVER->setBinLabel(9,"HAD4",2);
         ZDC_DigiErrors_DVER->getTH2F()->SetOption("coltext");
 
     ZDC_DigiErrors_CAPID = ib.book2D("ZDC_DigiErrors_CAPID","Raw Number of Digi Errors Caused by the Caps not Alternating",2,0,2,9,0,9);
         ZDC_DigiErrors_CAPID->setBinLabel(1,"ZDC+",1);
         ZDC_DigiErrors_CAPID->setBinLabel(2,"ZDC-",1);
     ZDC_DigiErrors_CAPID->setBinLabel(1,"EM1",2);
         ZDC_DigiErrors_CAPID->setBinLabel(2,"EM2",2);
         ZDC_DigiErrors_CAPID->setBinLabel(3,"EM3",2);
         ZDC_DigiErrors_CAPID->setBinLabel(4,"EM4",2);
         ZDC_DigiErrors_CAPID->setBinLabel(5,"EM5",2);
         ZDC_DigiErrors_CAPID->setBinLabel(6,"HAD1",2);
         ZDC_DigiErrors_CAPID->setBinLabel(7,"HAD2",2);
         ZDC_DigiErrors_CAPID->setBinLabel(8,"HAD3",2);
         ZDC_DigiErrors_CAPID->setBinLabel(9,"HAD4",2);
         ZDC_DigiErrors_CAPID->getTH2F()->SetOption("coltext");
 
 
     // hot channels
     ib.setCurrentFolder(baseFolder_ + "/Errors/HotChannel");
         ZDC_Hot_Channel_Errors = ib.book2D("ZDC_Hot_Channel_Errors", "Raw Number of Times Each Channel Appeared Hot", 2, 0, 2, 9, 0, 9);
         ZDC_Hot_Channel_Errors->setBinLabel(1,"ZDC+",1);
         ZDC_Hot_Channel_Errors->setBinLabel(2,"ZDC-",1);
         ZDC_Hot_Channel_Errors->setBinLabel(1,"EM1",2);
         ZDC_Hot_Channel_Errors->setBinLabel(2,"EM2",2);
         ZDC_Hot_Channel_Errors->setBinLabel(3,"EM3",2);
         ZDC_Hot_Channel_Errors->setBinLabel(4,"EM4",2);
         ZDC_Hot_Channel_Errors->setBinLabel(5,"EM5",2);
         ZDC_Hot_Channel_Errors->setBinLabel(6,"HAD1",2);
         ZDC_Hot_Channel_Errors->setBinLabel(7,"HAD2",2);
         ZDC_Hot_Channel_Errors->setBinLabel(8,"HAD3",2);
         ZDC_Hot_Channel_Errors->setBinLabel(9,"HAD4",2);
         ZDC_Hot_Channel_Errors->getTH2F()->SetOption("coltext");
 
 
     ZDC_HotChannelErrorsVsLS = ib.book1D("ZDC_HotChannelErrorsVsLS","Total Number of Hot Channel Errors in the ZDCs vs. Lumi Section; LS; # Hot channels", NLumiBlocks_,0,NLumiBlocks_);
 
     // dead channels
     ib.setCurrentFolder(baseFolder_ + "/Errors/DeadChannel");
         ZDC_Dead_Channel_Errors = ib.book2D("ZDC_Dead_Channel_Errors", "Raw Number of Times Each Channel Appeared Dead", 2, 0, 2, 9, 0, 9);
         ZDC_Dead_Channel_Errors->setBinLabel(1,"ZDC+",1);
         ZDC_Dead_Channel_Errors->setBinLabel(2,"ZDC-",1);
         ZDC_Dead_Channel_Errors->setBinLabel(1,"EM1",2);
         ZDC_Dead_Channel_Errors->setBinLabel(2,"EM2",2);
         ZDC_Dead_Channel_Errors->setBinLabel(3,"EM3",2);
         ZDC_Dead_Channel_Errors->setBinLabel(4,"EM4",2);
         ZDC_Dead_Channel_Errors->setBinLabel(5,"EM5",2);
         ZDC_Dead_Channel_Errors->setBinLabel(6,"HAD1",2);
         ZDC_Dead_Channel_Errors->setBinLabel(7,"HAD2",2);
         ZDC_Dead_Channel_Errors->setBinLabel(8,"HAD3",2);
         ZDC_Dead_Channel_Errors->setBinLabel(9,"HAD4",2);
         ZDC_Dead_Channel_Errors->getTH2F()->SetOption("coltext");
 
     ZDC_DeadChannelErrorsVsLS = ib.book1D("ZDC_DeadChannelErrorsVsLS","Total Number of Dead Channel Errors in the ZDC vs. Lumi Section; LS; # of Dead Chanels", NLumiBlocks_, 0, NLumiBlocks_);
 
     // cold channels
     ib.setCurrentFolder(baseFolder_ + "/Errors/ColdChannel");
         ZDC_Cold_Channel_Errors = ib.book2D("ZDC_Cold_Channel_Errors", "Raw Number of Times Each Channel Appeared Cold", 2, 0, 2, 9, 0, 9);
         ZDC_Cold_Channel_Errors->setBinLabel(1,"ZDC+",1);
         ZDC_Cold_Channel_Errors->setBinLabel(2,"ZDC-",1);
         ZDC_Cold_Channel_Errors->setBinLabel(1,"EM1",2);
         ZDC_Cold_Channel_Errors->setBinLabel(2,"EM2",2);
         ZDC_Cold_Channel_Errors->setBinLabel(3,"EM3",2);
         ZDC_Cold_Channel_Errors->setBinLabel(4,"EM4",2);
         ZDC_Cold_Channel_Errors->setBinLabel(5,"EM5",2);
         ZDC_Cold_Channel_Errors->setBinLabel(6,"HAD1",2);
         ZDC_Cold_Channel_Errors->setBinLabel(7,"HAD2",2);
         ZDC_Cold_Channel_Errors->setBinLabel(8,"HAD3",2);
         ZDC_Cold_Channel_Errors->setBinLabel(9,"HAD4",2);
         ZDC_Cold_Channel_Errors->getTH2F()->SetOption("coltext");
 
     ZDC_ColdChannelErrorsVsLS=ib.book1D("ZDC_ColdChannelErrorsVsLS","Total Number of Cold Channels in the ZDC vs. Lumi Section; LS; # of Cold Chanels", NLumiBlocks_, 0, NLumiBlocks_);
 
     // total errors
     ib.setCurrentFolder(baseFolder_ + "/Errors");
         ZDC_TotalChannelErrors = ib.book2D("ZDC_TotalChannelErrors","Total Number of Errors(Digi Error, Hot Cell or Dead Cell) Per Channel in the ZDC" ,2,0,2,9,0,9);
         ZDC_TotalChannelErrors->setBinLabel(1,"ZDC+",1);
         ZDC_TotalChannelErrors->setBinLabel(2,"ZDC-",1);
         ZDC_TotalChannelErrors->setBinLabel(1,"EM1",2);
         ZDC_TotalChannelErrors->setBinLabel(2,"EM2",2);
         ZDC_TotalChannelErrors->setBinLabel(3,"EM3",2);
         ZDC_TotalChannelErrors->setBinLabel(4,"EM4",2);
         ZDC_TotalChannelErrors->setBinLabel(5,"EM5",2);
         ZDC_TotalChannelErrors->setBinLabel(6,"HAD1",2);
         ZDC_TotalChannelErrors->setBinLabel(7,"HAD2",2);
         ZDC_TotalChannelErrors->setBinLabel(8,"HAD3",2);
         ZDC_TotalChannelErrors->setBinLabel(9,"HAD4",2);
         ZDC_TotalChannelErrors->getTH2F()->SetOption("coltext");
 
     ib.setCurrentFolder(baseFolder_ + "/Digis");
 
     for (int i = 0; i < 5; ++i) {
         // pulse Plus Side
         sprintf(title, "h_ZDCP_EMChan_%i_Pulse", i + 1);
         sprintf(name, "ZDC Plus EM Section Pulse for channel %i", i + 1);
         h_ZDCP_EM_Pulse[i] = ib.book1D(title, name, 10, -0.5, 9.5);
         h_ZDCP_EM_Pulse[i]->setAxisTitle("Time Slice id",1);
         h_ZDCP_EM_Pulse[i]->setAxisTitle("Pulse Height",2);
         // pulse Minus Side
         sprintf(title, "h_ZDCM_EMChan_%i_Pulse", i + 1);
         sprintf(name, "ZDC Minus EM Section Pulse for channel %i", i + 1);
         h_ZDCM_EM_Pulse[i] = ib.book1D(title, name, 10, -0.5, 9.5);
         h_ZDCM_EM_Pulse[i]->setAxisTitle("Time Slice id",1);
         h_ZDCM_EM_Pulse[i]->setAxisTitle("Pulse Height",2);
         // integrated charge over 10 time samples
         sprintf(title, "h_ZDCP_EMChan_%i_Charge", i + 1);
         sprintf(name, "ZDC Plus EM Section Charge for channel %i", i + 1);
         h_ZDCP_EM_Charge[i] = ib.book1D(title, name, 1000, 0., 30000.);
         h_ZDCP_EM_Charge[i]->setAxisTitle("Charge (fC)",1);
         h_ZDCP_EM_Charge[i]->setAxisTitle("Events",2);
         // integrated charge over 10 time samples
         sprintf(title, "h_ZDCM_EMChan_%i_Charge", i + 1);
         sprintf(name, "ZDC Minus EM Section Charge for channel %i", i + 1);
         h_ZDCM_EM_Charge[i] = ib.book1D(title, name, 1000, 0., 30000.);
         h_ZDCM_EM_Charge[i]->setAxisTitle("Charge (fC)",1);
         h_ZDCM_EM_Charge[i]->setAxisTitle("Events",2);
         // charge weighted time slice
         sprintf(title, "h_ZDCP_EMChan_%i_TSMean", i + 1);
         sprintf(name, "ZDC Plus EM Section TSMean for channel %i", i + 1);
         h_ZDCP_EM_TSMean[i] = ib.book1D(title, name, 100, -0.5, 9.5);
         h_ZDCP_EM_TSMean[i]->setAxisTitle("Timing",1);
         h_ZDCP_EM_TSMean[i]->setAxisTitle("Events",2);
         // charge weighted time slice
         sprintf(title, "h_ZDCM_EMChan_%i_TSMean", i + 1);
         sprintf(name, "ZDC Minus EM Section TSMean for channel %i", i + 1);
         h_ZDCM_EM_TSMean[i] = ib.book1D(title, name, 100, -0.5, 9.5);
         h_ZDCM_EM_TSMean[i]->setAxisTitle("Timing",1);
         h_ZDCM_EM_TSMean[i]->setAxisTitle("Events",2);
     }
 
     for (int i = 0; i < 4; ++i) {
         // pulse Plus Side
         sprintf(title, "h_ZDCP_HADChan_%i_Pulse", i + 1);
         sprintf(name, "ZDC Plus HAD Section Pulse for channel %i", i + 1);
         h_ZDCP_HAD_Pulse[i] = ib.book1D(title, name, 10, -0.5, 9.5);
         h_ZDCP_HAD_Pulse[i]->setAxisTitle("Time Slice id",1);
         h_ZDCP_HAD_Pulse[i]->setAxisTitle("Pulse Height",2);
         // pulse Minus Side
         sprintf(title, "h_ZDCM_HADChan_%i_Pulse", i + 1);
         sprintf(name, "ZDC Minus HAD Section Pulse for channel %i", i + 1);
         h_ZDCM_HAD_Pulse[i] = ib.book1D(title, name, 10, -0.5, 9.5);
         h_ZDCM_HAD_Pulse[i]->setAxisTitle("Time Slice id",1);
         h_ZDCM_HAD_Pulse[i]->setAxisTitle("Pulse Height",2);
         // integrated charge over 10 time samples
         sprintf(title, "h_ZDCP_HADChan_%i_Charge", i + 1);
         sprintf(name, "ZDC Plus HAD Section Charge for channel %i", i + 1);
         h_ZDCP_HAD_Charge[i] = ib.book1D(title, name, 1000, 0., 30000.);
         h_ZDCP_HAD_Charge[i]->setAxisTitle("Charge (fC)",1);
         h_ZDCP_HAD_Charge[i]->setAxisTitle("Events",2);
         // integrated charge over 10 time samples
         sprintf(title, "h_ZDCM_HADChan_%i_Charge", i + 1);
         sprintf(name, "ZDC Minus HAD Section Charge for channel %i", i + 1);
         h_ZDCM_HAD_Charge[i] = ib.book1D(title, name, 1000, 0., 30000.);
         h_ZDCM_HAD_Charge[i]->setAxisTitle("Charge (fC)",1);
         h_ZDCM_HAD_Charge[i]->setAxisTitle("Events",2);
         // charge weighted time slice
         sprintf(title, "h_ZDCP_HADChan_%i_TSMean", i + 1);
         sprintf(name, "ZDC Plus HAD Section TSMean for channel %i", i + 1);
         h_ZDCP_HAD_TSMean[i] = ib.book1D(title, name, 100, -0.5, 9.5);
         h_ZDCP_HAD_TSMean[i]->setAxisTitle("Timing",1);
         h_ZDCP_HAD_TSMean[i]->setAxisTitle("Events",2);
         // charge weighted time slice
         sprintf(title, "h_ZDCM_HADChan_%i_TSMean", i + 1);
         sprintf(name, "ZDC Minus HAD Section TSMean for channel %i", i + 1);
         h_ZDCM_HAD_TSMean[i] = ib.book1D(title, name, 100, -0.5, 9.5);
         h_ZDCM_HAD_TSMean[i]->setAxisTitle("Timing",1);
         h_ZDCM_HAD_TSMean[i]->setAxisTitle("Events",2);
     }
 
     ib.setCurrentFolder(baseFolder_ + "/RecHits");
 
     for (int i = 0; i < 5; ++i) {
         //RecHitEnergy Plus Side
         sprintf(title,"h_ZDCP_EMChan_%i_RecHit_Energy",i+1);
         sprintf(name,"ZDC EM Section Rechit Energy for channel %i",i+1);
         h_ZDCP_EM_RecHitEnergy[i] = ib.book1D(title, name, 1010, -100., 10000.);
         h_ZDCP_EM_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
         h_ZDCP_EM_RecHitEnergy[i]->setAxisTitle("Events",2);
         //RecHitEnergy Minus Side
         sprintf(title,"h_ZDCM_EMChan_%i_RecHit_Energy",i+1);
         sprintf(name,"ZDC EM Section Rechit Energy for channel %i",i+1);
         h_ZDCM_EM_RecHitEnergy[i] = ib.book1D(title, name, 1010, -100., 10000.);
         h_ZDCM_EM_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
         h_ZDCM_EM_RecHitEnergy[i]->setAxisTitle("Events",2);
         //RecHit Timing Plus Side
         sprintf(title,"h_ZDCP_EMChan_%i_RecHit_Timing",i+1);
         sprintf(name,"ZDC EM Section Rechit Timing for channel %i",i+1);
         h_ZDCP_EM_RecHitTiming[i] = ib.book1D(title, name, 100, -100., 100.);
         h_ZDCP_EM_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
         h_ZDCP_EM_RecHitTiming[i]->setAxisTitle("Events",2);
         //RecHit Timing Minus Side
         sprintf(title,"h_ZDCM_EMChan_%i_RecHit_Timing",i+1);
         sprintf(name,"ZDC EM Section Rechit Timing for channel %i",i+1);
         h_ZDCM_EM_RecHitTiming[i] = ib.book1D(title, name, 100, -100., 100.);
         h_ZDCM_EM_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
         h_ZDCM_EM_RecHitTiming[i]->setAxisTitle("Events",2);
     }
 
     for (int i = 0; i < 4; ++i) {
         //RecHitEnergy Plus Side
         sprintf(title,"h_ZDCP_HADChan_%i_RecHit_Energy",i+1);
         sprintf(name,"ZDC HAD Section Rechit Energy for channel %i",i+1);
         h_ZDCP_HAD_RecHitEnergy[i] = ib.book1D(title, name, 1010, -100., 10000.);
         h_ZDCP_HAD_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
         h_ZDCP_HAD_RecHitEnergy[i]->setAxisTitle("Events",2);
         //RecHitEnergy Minus Side
         sprintf(title,"h_ZDCM_HADChan_%i_RecHit_Energy",i+1);
         sprintf(name,"ZDC HAD Section Rechit Energy for channel %i",i+1);
         h_ZDCM_HAD_RecHitEnergy[i] = ib.book1D(title, name, 1010, -100., 10000.);
         h_ZDCM_HAD_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
         h_ZDCM_HAD_RecHitEnergy[i]->setAxisTitle("Events",2);
         //RecHit Timing Plus Side
         sprintf(title,"h_ZDCP_HADChan_%i_RecHit_Timing",i+1);
         sprintf(name,"ZDC HAD Section Rechit Timing for channel %i",i+1);
         h_ZDCP_HAD_RecHitTiming[i] = ib.book1D(title, name, 100, -100., 100.);
         h_ZDCP_HAD_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
         h_ZDCP_HAD_RecHitTiming[i]->setAxisTitle("Events",2);
         //RecHit Timing Minus Side
         sprintf(title,"h_ZDCM_HADChan_%i_RecHit_Timing",i+1);
         sprintf(name,"ZDC HAD Section Rechit Timing for channel %i",i+1);
         h_ZDCM_HAD_RecHitTiming[i] = ib.book1D(title, name, 100, -100., 100.);
         h_ZDCM_HAD_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
         h_ZDCM_HAD_RecHitTiming[i]->setAxisTitle("Events",2);
     }
 
     return;
 }