HcalZDCMonitor.h

Go to the documentation of this file.

#ifndef DQM_HCALMONITORTASKS_HCALZDCMONITOR_H
#define DQM_HCALMONITORTASKS_HCALZDCMONITOR_H

#include "DQM/HcalMonitorTasks/interface/HcalBaseDQMonitor.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "EventFilter/HcalRawToDigi/interface/HcalUnpacker.h"
#include "CondFormats/HcalObjects/interface/HcalPedestal.h"
#include "CondFormats/HcalObjects/interface/HcalPedestalWidth.h"

#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "CalibFormats/HcalObjects/interface/HcalCoderDb.h"
#include "CalibFormats/HcalObjects/interface/HcalCalibrations.h"
#include "CalibFormats/HcalObjects/interface/HcalDbService.h"
#include "CalibFormats/HcalObjects/interface/HcalDbRecord.h"
#include "CondFormats/HcalObjects/interface/HcalQIECoder.h"
#include "CalibFormats/HcalObjects/interface/HcalCoder.h"
#include "DataFormats/HcalDetId/interface/HcalZDCDetId.h"
#include "DataFormats/DetId/interface/DetId.h"

#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "EventFilter/HcalRawToDigi/interface/HcalDCCHeader.h"
#include "CondFormats/HcalObjects/interface/HcalElectronicsMap.h"
#include "FWCore/Utilities/interface/CPUTimer.h"

#include <cmath>

/*  Revision!!!!
 *
 * $DATE: 2012/19/12
 * $Revision: 1.2
 * \author: J. Gomez
 */

class HcalZDCMonitor: public HcalBaseDQMonitor
{
 public:
  HcalZDCMonitor(const edm::ParameterSet& ps);
  ~HcalZDCMonitor();
  void setup();
  void processEvent(const ZDCDigiCollection& digi,
                    const ZDCRecHitCollection& rechit, const HcalUnpackerReport& report);
  //void done();
  void reset();
  void endLuminosityBlock(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& c);
  void beginRun(const edm::Run& run, const edm::EventSetup& c);
  void beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& c);
  void endRun(const edm::Run& run, const edm::EventSetup& c);
  void analyze(const edm::Event& e, const edm::EventSetup& c);
 private:
  //virtual void endJob();
  //void endLuminosityBlock(void);
  //void zeroCounters();
  //void fillHistos();
  //void setZDClabels(MonitorElement* h);
  double getTime(std::vector<double> fData, unsigned int ts_min, unsigned int ts_max, double& fSum);
  //int getTSMax(std::vector<double> fData);
  //bool isGood(std::vector<double> fData, double fCut, double fPercentage);


  bool checkZDC_;
  int NumBadZDC;
  MonitorElement* ProblemsVsLB_ZDC;
  edm::InputTag rechitLabel_;
  edm::InputTag digiLabel_;

  float ChannelRatio[18];//errorevents/total events in 1 LS
  int EventCounter;//events in lumi
  int TotalChannelErrors[18];//total events with an error per channel in a LS
  std::vector<double> ChannelWeighting_; //Quality index(QI) see description below
  std::vector<double> MaxErrorRates_; //the fractional error rate before a channel is called bad for a LS
  int OfflineColdThreshold_;
  int OnlineColdThreshold_;
  int OfflineDeadThreshold_;
  int OnlineDeadThreshold_;
  int ColdADCThreshold_;
  int DeadChannelCounter[18];
  int ColdChannelCounter[18];
  bool DeadChannelError[18];
  bool HotChannelError[18];
  bool DigiErrorCAPID[18];
  bool DigiErrorDVER[18];
  bool ChannelHasDigiError[18];

  const HcalQIEShape* shape_;
  const HcalQIECoder* channelCoder_;
  HcalCalibrations calibs_;
  int ievt_;
  //int zdc_checkNevents_;
  MonitorElement* meEVT_;
  MonitorElement* h_2D_saturation;
  MonitorElement* h_2D_charge;
  MonitorElement* h_2D_TSMean;
  MonitorElement* h_2D_RecHitEnergy;
  MonitorElement* h_2D_RecHitTime;
  MonitorElement* h_ZDCP_EM_Pulse[5];
  MonitorElement* h_ZDCM_EM_Pulse[5];
  MonitorElement* h_ZDCP_EM_Charge[5];
  MonitorElement* h_ZDCM_EM_Charge[5];
  MonitorElement* h_ZDCP_EM_TSMean[5];
  MonitorElement* h_ZDCM_EM_TSMean[5];
  MonitorElement* h_ZDCP_HAD_Pulse[4];
  MonitorElement* h_ZDCM_HAD_Pulse[4];
  MonitorElement* h_ZDCP_HAD_Charge[4];
  MonitorElement* h_ZDCM_HAD_Charge[4];
  MonitorElement* h_ZDCP_HAD_TSMean[4];
  MonitorElement* h_ZDCM_HAD_TSMean[4];
  MonitorElement* h_ZDCP_EM_RecHitEnergy[5];
  MonitorElement* h_ZDCM_EM_RecHitEnergy[5];
  MonitorElement* h_ZDCP_EM_RecHitTiming[5];
  MonitorElement* h_ZDCM_EM_RecHitTiming[5];
  MonitorElement* h_ZDCP_HAD_RecHitEnergy[4];
  MonitorElement* h_ZDCM_HAD_RecHitEnergy[4];
  MonitorElement* h_ZDCP_HAD_RecHitTiming[4];
  MonitorElement* h_ZDCM_HAD_RecHitTiming[4];

  MonitorElement* ZDC_Digi_Errors;
  MonitorElement* ZDC_DigiErrorsVsLS;
  MonitorElement* ZDC_DigiErrors_DVER;
  MonitorElement* ZDC_DigiErrors_CAPID;

  MonitorElement* ZDC_Hot_Channel_Errors;
  MonitorElement* ZDC_HotChannelErrorsVsLS;

  MonitorElement* ZDC_Cold_Channel_Errors;
  MonitorElement* ZDC_ColdChannelErrorsVsLS;

  MonitorElement* ZDC_Dead_Channel_Errors;
  MonitorElement* ZDC_DeadChannelErrorsVsLS;

  MonitorElement* ZDC_TotalChannelErrors;
  MonitorElement* EventsVsLS;
  // Quality index(QI) per LS...QI is a number between 0 and 1. Each ZDC channel is assigned a custom weighting factor, found in /python/HcalZDCMonitor_cfi.py. This weighting factor is representative of the "value" of that channel to data taking. For example, HAD1 is more valuable than HAD4, this is because HAD1 is closer to the IP and will absorb more of the energy of the neutron, so if it is broken, we care more about it being broken than if HAD4 is broken. If all channels are working properly in a given ZDC (ZDC+/-),the QI for that ZDC will be 1. If a channel is deemed malfunctional for a given LS, then the assigned weightfactor will be subtracted from 1.
  MonitorElement* PZDC_QualityIndexVsLB_; 
  MonitorElement* NZDC_QualityIndexVsLB_;
};

#endif