CastorMonitorModule.h

Go to the documentation of this file.

#ifndef CastorMonitorModule_H
#define CastorMonitorModule_H

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "SimG4CMS/Calo/interface/CaloHit.h"
#include "SimDataFormats/CaloHit/interface/PCaloHit.h"
#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"

#include "FWCore/Utilities/interface/CPUTimer.h"
#include "DataFormats/Provenance/interface/EventID.h"  

#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetup.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetupFwd.h"


#include "DataFormats/CastorReco/interface/CastorTower.h"
#include "DataFormats/CastorReco/interface/CastorCluster.h"
#include "DataFormats/CastorReco/interface/CastorJet.h"
#include "DataFormats/JetReco/interface/CastorJetID.h"
#include "RecoJets/JetProducers/interface/CastorJetIDHelper.h"
#include "RecoJets/JetProducers/plugins/CastorJetIDProducer.h"
#include "DataFormats/JetReco/interface/BasicJet.h"
#include "DataFormats/JetReco/interface/BasicJetCollection.h"
#include "DataFormats/JetReco/interface/Jet.h"

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/HcalDigi/interface/HcalUnpackerReport.h" //-- no CastorUnpackerReport at the moment !
#include "DataFormats/HcalDetId/interface/HcalCastorDetId.h" //-- HcalCastorDetId

#include "DQM/CastorMonitor/interface/CastorMonitorSelector.h"
#include "DQM/CastorMonitor/interface/CastorDigiMonitor.h"
#include "DQM/CastorMonitor/interface/CastorRecHitMonitor.h"
#include "DQM/CastorMonitor/interface/CastorChannelQualityMonitor.h"
#include "DQM/CastorMonitor/interface/CastorLEDMonitor.h"
#include "DQM/CastorMonitor/interface/CastorPSMonitor.h"
#include "DQM/CastorMonitor/interface/CastorHIMonitor.h"
#include "DQM/CastorMonitor/interface/CastorDataIntegrityMonitor.h"
#include "DQM/CastorMonitor/interface/CastorTowerJetMonitor.h"

#include "CalibCalorimetry/CastorCalib/interface/CastorDbASCIIIO.h" //-- use to get/dump Calib to DB 
#include "CondFormats/CastorObjects/interface/CastorChannelQuality.h" //-- use to get/hold channel status
#include "CondFormats/DataRecord/interface/CastorChannelQualityRcd.h"



#include <memory>
#include <iostream>
#include <fstream>
#include <vector>
#include <string>
#include <sys/time.h>



class CastorMonitorModule : public edm::EDAnalyzer{

public:
  
  CastorMonitorModule(const edm::ParameterSet& ps);

  ~CastorMonitorModule();
  
 protected:
  
  void analyze(const edm::Event& iEvent, const edm::EventSetup& eventSetup);
  
  void beginJob();
  
  void beginRun(const edm::Run& iRun, const edm::EventSetup& eventSetup);

  void beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg, 
                            const edm::EventSetup& eventSetup) ;

  void endLuminosityBlock(const edm::LuminosityBlock& lumiSeg, 
                          const edm::EventSetup& eventSetup);

  void endJob(void);
  
  void endRun(const edm::Run& run, const edm::EventSetup& eventSetup);

  void reset(void);

  bool prescale();

  void CheckCastorStatus     (const FEDRawDataCollection& rawraw, 
                  const HcalUnpackerReport& report, 
                  const CastorElectronicsMap& emap,
                  const CastorDigiCollection& castordigi
                  );
    
 private:
 
  int prescaleEvt_;    //-- units of events
  int prescaleLS_;     //-- units of lumi sections
  int prescaleTime_;   //-- units of minutes
  int prescaleUpdate_; //-- units of "updates", TBD

  std::string monitorName_;

  int fVerbosity;  

  int nevt_;
  int nlumisecs_;
  bool saved_;

  bool rawOK_    ;
  bool reportOK_ ;
  bool digiOK_   ;
  bool rechitOK_ ;
  bool towerOK_  ;
  bool jetOK_    ;
  bool jetIdOK_  ;
  
  int nRaw;
  int nDigi;
  int nRechit;
  int nTower;
  int nJet;
  int nJetId;
   
  
  bool showTiming_; 
  edm::CPUTimer cpu_timer; 

  struct{
    timeval startTV,updateTV;
    double elapsedTime; 
    double vetoTime; 
    double updateTime;
  } psTime_;    


  DQMStore* dbe_;  
  
  int irun_,ilumisec_,ievent_,itime_,ibunch_;
  bool actonLS_ ;
  std::string rootFolder_;

  int ievt_;
  int ievt_pre_; //-- copy of counter used for prescale purposes
  bool fedsListed_;
  
  edm::InputTag inputLabelRaw_;
  edm::InputTag inputLabelReport_;
  edm::InputTag inputLabelDigi_;
  edm::InputTag inputLabelRecHitCASTOR_;
  edm::InputTag inputLabelTowerCASTOR_;  
  edm::InputTag inputLabelBasicJetCASTOR_;  
  edm::InputTag inputLabelJetIdCASTOR_ ;
  edm::InputTag inputLabelCastorTowers_    ; 
  edm::InputTag inputLabelCastorBasicJets_ ; 
  edm::InputTag inputLabelCastorJetIDs_ ; 

  // typedef std::vector<reco::CastorTower> CastorTowerCollection;
  // edm::InputTag inputLabelCaloTower_;
  // edm::InputTag inputLabelLaser_;

  std::map<uint32_t, std::vector<HcalCastorDetId> > DCCtoCell;
  std::map<uint32_t, std::vector<HcalCastorDetId> > ::iterator thisDCC;
  std::map<std::pair <int,int> , std::vector<HcalCastorDetId> > HTRtoCell;
  std::map<std::pair <int,int> , std::vector<HcalCastorDetId> > ::iterator thisHTR;

  MonitorElement* meFEDS_;
  MonitorElement* meStatus_;
  MonitorElement* meRunType_;
  MonitorElement* meEvtMask_;
  MonitorElement* meTrigger_;
  MonitorElement* meLatency_;
  MonitorElement* meQuality_;
  MonitorElement* CastorEventProduct;
 

  CastorMonitorSelector*    evtSel_;
  CastorRecHitMonitor*      RecHitMon_;
  CastorChannelQualityMonitor*  CQMon_;
  CastorDigiMonitor*        DigiMon_;
  CastorLEDMonitor*         LedMon_;
  CastorPSMonitor*          PSMon_;
  CastorHIMonitor*          HIMon_;
  CastorDataIntegrityMonitor* DataIntMon_;
  CastorTowerJetMonitor*     TowerJetMon_;

  MonitorElement* meEVT_;

  edm::ESHandle<CastorDbService> conditions_;
  const CastorElectronicsMap*     CastorReadoutMap_;

  edm::ESHandle<CastorPedestals> dbPedestals;

  // pedestal width averaged over capIDs, calculated in beginRun
  // aware of the difference between index[0..15][0..13] 
  // and sector/module numeration[1..16][1..14]
  float        fPedestalNSigmaAverage[14][16];

  std::vector<HcalGenericDetId> listEMap; //electronics Emap


  std::ofstream m_logFile;

  bool checkCASTOR_;

  MonitorElement* meCASTOR_;
   
  bool EDMonOn_;

  bool dump2database_;
  std::map<HcalCastorDetId, unsigned int> myquality_;
  CastorChannelQuality* chanquality_;

  //define tokens
  edm::EDGetTokenT<FEDRawDataCollection> inputLabelRawToken_;
  edm::EDGetTokenT<HcalUnpackerReport> inputLabelReportToken_;
  edm::EDGetTokenT<CastorDigiCollection> inputLabelDigiToken_;
  edm::EDGetTokenT<CastorRecHitCollection> inputLabelRecHitCASTORToken_;
};

#endif