CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

HcalZDCMonitor Class Reference

#include <HcalZDCMonitor.h>

Inheritance diagram for HcalZDCMonitor:
HcalBaseMonitor

List of all members.

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 ()

Private Member Functions

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

Private Attributes

HcalCalibrations calibs_
const HcalQIECoderchannelCoder_
bool checkZDC_
MonitorElementh_2D_charge
MonitorElementh_2D_RecHitEnergy
MonitorElementh_2D_RecHitTime
MonitorElementh_2D_saturation
MonitorElementh_2D_TSMean
MonitorElementh_ZDCM_EM_Charge [5]
MonitorElementh_ZDCM_EM_Pulse [5]
MonitorElementh_ZDCM_EM_RecHitEnergy [5]
MonitorElementh_ZDCM_EM_RecHitTiming [5]
MonitorElementh_ZDCM_EM_TSMean [5]
MonitorElementh_ZDCM_HAD_Charge [4]
MonitorElementh_ZDCM_HAD_Pulse [4]
MonitorElementh_ZDCM_HAD_RecHitEnergy [4]
MonitorElementh_ZDCM_HAD_RecHitTiming [4]
MonitorElementh_ZDCM_HAD_TSMean [4]
MonitorElementh_ZDCP_EM_Charge [5]
MonitorElementh_ZDCP_EM_Pulse [5]
MonitorElementh_ZDCP_EM_RecHitEnergy [5]
MonitorElementh_ZDCP_EM_RecHitTiming [5]
MonitorElementh_ZDCP_EM_TSMean [5]
MonitorElementh_ZDCP_HAD_Charge [4]
MonitorElementh_ZDCP_HAD_Pulse [4]
MonitorElementh_ZDCP_HAD_RecHitEnergy [4]
MonitorElementh_ZDCP_HAD_RecHitTiming [4]
MonitorElementh_ZDCP_HAD_TSMean [4]
int ievt_
MonitorElementmeEVT_
int NumBadZDC
MonitorElementProblemsVsLB_ZDC
const HcalQIEShapeshape_

Detailed Description

$DATE: 2010/02/04 $Revision:

Author:
S.Sen

Definition at line 19 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  )

Reimplemented from HcalBaseMonitor.

Definition at line 439 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 ( std::vector< double >  fData,
unsigned int  ts_min,
unsigned int  ts_max,
double &  fSum 
) [private]

Definition at line 415 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, ExpressReco_HICollisions_FallBack::id, 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 = 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_, AlCaRecoCosmics_cfg::name, edm::CPUTimer::reset(), HcalBaseMonitor::rootFolder_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), 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

Definition at line 46 of file HcalZDCMonitor.h.

Definition at line 45 of file HcalZDCMonitor.h.

bool HcalZDCMonitor::checkZDC_ [private]

Definition at line 40 of file HcalZDCMonitor.h.

Definition at line 51 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

Definition at line 53 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

Definition at line 54 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

Definition at line 50 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

Definition at line 52 of file HcalZDCMonitor.h.

Referenced by endLuminosityBlock(), and setup().

Definition at line 58 of file HcalZDCMonitor.h.

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

Definition at line 56 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

Definition at line 68 of file HcalZDCMonitor.h.

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

Definition at line 70 of file HcalZDCMonitor.h.

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

Definition at line 60 of file HcalZDCMonitor.h.

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

Definition at line 64 of file HcalZDCMonitor.h.

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

Definition at line 62 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

Definition at line 72 of file HcalZDCMonitor.h.

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

Definition at line 74 of file HcalZDCMonitor.h.

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

Definition at line 66 of file HcalZDCMonitor.h.

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

Definition at line 57 of file HcalZDCMonitor.h.

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

Definition at line 55 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

Definition at line 67 of file HcalZDCMonitor.h.

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

Definition at line 69 of file HcalZDCMonitor.h.

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

Definition at line 59 of file HcalZDCMonitor.h.

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

Definition at line 63 of file HcalZDCMonitor.h.

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

Definition at line 61 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

Definition at line 71 of file HcalZDCMonitor.h.

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

Definition at line 73 of file HcalZDCMonitor.h.

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

Definition at line 65 of file HcalZDCMonitor.h.

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

int HcalZDCMonitor::ievt_ [private]

Reimplemented from HcalBaseMonitor.

Definition at line 47 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

Reimplemented from HcalBaseMonitor.

Definition at line 49 of file HcalZDCMonitor.h.

Referenced by processEvent(), and setup().

Definition at line 41 of file HcalZDCMonitor.h.

Definition at line 42 of file HcalZDCMonitor.h.

Definition at line 44 of file HcalZDCMonitor.h.