#include <HcalZDCMonitor.h>

Inheritance diagram for HcalZDCMonitor:

Public Member Functions
void	endLuminosityBlock (void)

	HcalZDCMonitor ()

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

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
HcalCalibrations	calibs_

const HcalQIECoder *	channelCoder_

bool	checkZDC_

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]

int	ievt_

MonitorElement *	meEVT_

int	NumBadZDC

MonitorElement *	ProblemsVsLB_ZDC

const HcalQIEShape *	shape_

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 17 of file HcalZDCMonitor.h.

Constructor & Destructor Documentation

HcalZDCMonitor::HcalZDCMonitor ( )

Definition at line 3 of file HcalZDCMonitor.cc.

3 {

4 }

HcalZDCMonitor::~HcalZDCMonitor ( )

Definition at line 5 of file HcalZDCMonitor.cc.

5 {

6 }

Member Function Documentation

void HcalZDCMonitor::endLuminosityBlock ( void )

Definition at line 438 of file HcalZDCMonitor.cc.

References 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, and MonitorElement::setBinContent().

Referenced by ZDCMonitorModule::endLuminosityBlock(), and ZDCMonitorModule::endRun().

 {
 
     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 414 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
	)

Definition at line 243 of file HcalZDCMonitor.cc.

References ecalMGPA::adc(), edm::SortedCollection< T, SORT >::begin(), gather_cfg::cout, HcalBaseMonitor::cpu_timer, edm::SortedCollection< T, SORT >::end(), 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, i, ievt_, meEVT_, edm::CPUTimer::reset(), HcalBaseMonitor::showTiming, edm::SortedCollection< T, SORT >::size(), edm::CPUTimer::start(), 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; 
 
     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 
 
     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 
       } // loop on zdc digi collection
 
 
     //--------------------------------------
     // 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 7 of file HcalZDCMonitor.cc.

Referenced by ZDCMonitorModule::reset().

7 {

8 }

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

virtual

Reimplemented from HcalBaseMonitor.

Definition at line 10 of file HcalZDCMonitor.cc.

References HcalBaseMonitor::baseFolder_, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookInt(), gather_cfg::cout, HcalBaseMonitor::cpu_timer, MonitorElement::Fill(), HcalBaseMonitor::fVerbosity, 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_, meEVT_, mergeVDriftHistosByStation::name, edm::CPUTimer::reset(), HcalBaseMonitor::rootFolder_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), HcalBaseMonitor::setup(), HcalBaseMonitor::showTiming, edm::CPUTimer::start(), and indexGen::title.

Referenced by ZDCMonitorModule::bookHistograms().

                                                                            {
     HcalBaseMonitor::setup(ps, ib);
 
     baseFolder_ = rootFolder_ + "ZDCMonitor_Hcal";
 
     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];
     
         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);
     
         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_ZDCP_HAD_Pulse[i]->setAxisTitle("Time Slice id",1);
         h_ZDCP_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;
 }