HcalZDCMonitor Class Reference

#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 *dbe)
	~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 (MonitorElement *&h, std::vector< MonitorElement * > &hh, std::string Name, std::string Units, int lowbound, int highbound, int Nbins)
void	setupDepthHists1D (std::vector< MonitorElement * > &hh, std::string Name, std::string Units, int lowbound, int highbound, int Nbins)
void	setupDepthHists2D (MonitorElement *&h, std::vector< MonitorElement * > &hh, std::string Name, std::string Units)
void	setupDepthHists2D (std::vector< MonitorElement * > &hh, std::string Name, std::string Units)
void	setupDepthHists2D (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 (std::vector< MonitorElement * > &hh, std::string Name, std::string Units, int nbinsx, int lowboundx, int highboundx, int nbinsy, int lowboundy, int highboundy)
void	SetupEtaPhiHists (MonitorElement *&h, EtaPhiHists &hh, std::string Name, std::string Units)
void	SetupEtaPhiHists (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
DQMStore *	m_dbe
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.

                               {
}

HcalZDCMonitor::~HcalZDCMonitor

(

)

Definition at line 5 of file HcalZDCMonitor.cc.

                                {
}

Member Function Documentation

void HcalZDCMonitor::endLuminosityBlock

(

void

)

Definition at line 441 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 417 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 245 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(), and edm::CPUTimer::start().

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)
      {     
    HcalZDCDetId id(rechit_iter->id());
    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().

                           {
}

void HcalZDCMonitor::setup ( const edm::ParameterSet &  ps, DQMStore *  dbe  )

virtual

Reimplemented from HcalBaseMonitor.

Definition at line 10 of file HcalZDCMonitor.cc.

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

Referenced by ZDCMonitorModule::ZDCMonitorModule().

                                                                     {
    HcalBaseMonitor::setup(ps, dbe);

    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 (m_dbe) {
        if (fVerbosity > 1)
            std::cout << "<HcalZDCMonitor::setup>  Setting up Histograms" << std::endl;

        m_dbe->setCurrentFolder(baseFolder_);
        meEVT_ = m_dbe->bookInt("ZDC Event Number");
        meEVT_->Fill(ievt_);
        char name[128];
        char title[128];
    
        h_2D_charge = m_dbe->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 = m_dbe->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 = m_dbe->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 = m_dbe->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 = m_dbe->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);
    
        m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->book1D(title, name, 100, -0.5, 9.5);
        h_ZDCM_HAD_TSMean[i]->setAxisTitle("Timing",1);
        h_ZDCM_HAD_TSMean[i]->setAxisTitle("Events",2);
        }

        m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->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] = m_dbe->book1D(title, name, 100, -100., 100.);  
        h_ZDCM_HAD_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
        h_ZDCM_HAD_RecHitTiming[i]->setAxisTitle("Events",2);
    }

    }
    return;
}

Member Data Documentation

HcalCalibrations HcalZDCMonitor::calibs_

private

Definition at line 44 of file HcalZDCMonitor.h.

const HcalQIECoder* HcalZDCMonitor::channelCoder_

private

Definition at line 43 of file HcalZDCMonitor.h.

bool HcalZDCMonitor::checkZDC_

private

Definition at line 38 of file HcalZDCMonitor.h.

MonitorElement* HcalZDCMonitor::h_2D_charge

private

Definition at line 49 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

MonitorElement* HcalZDCMonitor::h_2D_RecHitEnergy

private

Definition at line 51 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

MonitorElement* HcalZDCMonitor::h_2D_RecHitTime

private

Definition at line 52 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

MonitorElement* HcalZDCMonitor::h_2D_saturation

private

Definition at line 48 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_2D_TSMean

private

Definition at line 50 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_EM_Charge[5]

private

Definition at line 56 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_EM_Pulse[5]

private

Definition at line 54 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_EM_RecHitEnergy[5]

private

Definition at line 66 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_EM_RecHitTiming[5]

private

Definition at line 68 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_EM_TSMean[5]

private

Definition at line 58 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_HAD_Charge[4]

private

Definition at line 62 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_HAD_Pulse[4]

private

Definition at line 60 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_HAD_RecHitEnergy[4]

private

Definition at line 70 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_HAD_RecHitTiming[4]

private

Definition at line 72 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCM_HAD_TSMean[4]

private

Definition at line 64 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_EM_Charge[5]

private

Definition at line 55 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_EM_Pulse[5]

private

Definition at line 53 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_EM_RecHitEnergy[5]

private

Definition at line 65 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_EM_RecHitTiming[5]

private

Definition at line 67 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_EM_TSMean[5]

private

Definition at line 57 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_HAD_Charge[4]

private

Definition at line 61 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_HAD_Pulse[4]

private

Definition at line 59 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_HAD_RecHitEnergy[4]

private

Definition at line 69 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_HAD_RecHitTiming[4]

private

Definition at line 71 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

MonitorElement* HcalZDCMonitor::h_ZDCP_HAD_TSMean[4]

private

Definition at line 63 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), processEvent(), and setup().

int HcalZDCMonitor::ievt_

private

Definition at line 45 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalZDCMonitor::meEVT_

private

Definition at line 47 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

int HcalZDCMonitor::NumBadZDC

private

Definition at line 39 of file HcalZDCMonitor.h.

MonitorElement* HcalZDCMonitor::ProblemsVsLB_ZDC

private

Definition at line 40 of file HcalZDCMonitor.h.

const HcalQIEShape* HcalZDCMonitor::shape_

private

Definition at line 42 of file HcalZDCMonitor.h.