HcalDigisValidation.h

Go to the documentation of this file.

#ifndef HCALDIGISVALIDATION_H
#define HCALDIGISVALIDATION_H
 
#include <memory>
 
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
 
#include "FWCore/Framework/interface/Event.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DQMServices/Core/interface/DQMStore.h"
 
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
 
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
 
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"
#include "Geometry/Records/interface/HcalGeometryRecord.h"
#include "Geometry/HcalTowerAlgo/interface/HcalGeometry.h"
 
#include "CalibFormats/HcalObjects/interface/HcalDbRecord.h"
#include "CalibFormats/HcalObjects/interface/HcalCoderDb.h"
#include "CalibFormats/HcalObjects/interface/HcalCalibrations.h"
#include "CondFormats/HcalObjects/interface/HcalQIEShape.h"
 
#include "DataFormats/HcalDigi/interface/HBHEDataFrame.h"
 
#include "CalibFormats/HcalObjects/interface/HcalDbService.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
 
#include "Geometry/CaloTopology/interface/HcalTopology.h"
 
/*TP Code*/
#include "Geometry/CaloTopology/interface/HcalTopology.h"
#include "CalibFormats/CaloTPG/interface/CaloTPGTranscoder.h"
#include "CalibFormats/CaloTPG/interface/CaloTPGRecord.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
/*~TP Code*/
 
#include <map>
#include <vector>
#include <utility>
#include <ostream>
#include <string>
#include <algorithm>
#include <cmath>
#include <iostream>
 
class HcalDigisValidation : public DQMEDAnalyzer {
public:
  explicit HcalDigisValidation(const edm::ParameterSet &);
 
  ~HcalDigisValidation() override;
 
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  void dqmBeginRun(const edm::Run &run, const edm::EventSetup &c) override;
 
private:
  struct HistLim {
    HistLim(int nbin, double mini, double maxi) : n(nbin), min(mini), max(maxi) {}
    int n;
    double min;
    double max;
  };
 
  void analyze(const edm::Event &, const edm::EventSetup &) override;
 
  std::map<std::string, MonitorElement *> *msm_;
 
  void book1D(DQMStore::IBooker &ib, std::string name, int n, double min, double max);
 
  void book1D(DQMStore::IBooker &ib, std::string name, const HistLim &limX);
 
  void fill1D(std::string name, double X, double weight = 1);
 
  void book2D(DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY);
 
  void fill2D(std::string name, double X, double Y, double weight = 1);
 
  void bookPf(DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY);
 
  void bookPf(DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY, const char *option);
 
  void fillPf(std::string name, double X, double Y);
 
  MonitorElement *monitor(std::string name);
 
  void booking(DQMStore::IBooker &ib, std::string subdetopt, int bnoise, int bmc);
 
  std::string str(int x);
 
  template <class Digi>
  void reco(const edm::Event &iEvent,
            const edm::EventSetup &iSetup,
            const edm::EDGetTokenT<edm::SortedCollection<Digi> > &tok);
  template <class dataFrameType>
  void reco(const edm::Event &iEvent,
            const edm::EventSetup &iSetup,
            const edm::EDGetTokenT<HcalDataFrameContainer<dataFrameType> > &tok);
 
  std::string outputFile_;
  std::string subdet_;
  std::string zside_;
  std::string dirName_;
  //    std::string inputLabel_;
  edm::InputTag inputTag_;
  edm::InputTag QIE10inputTag_;
  edm::InputTag QIE11inputTag_;
  edm::InputTag emulTPsTag_;
  edm::InputTag dataTPsTag_;
  std::string mode_;
  std::string mc_;
  int noise_;
  bool testNumber_;
  bool hep17_;
  bool HEPhase1_;
  bool HBPhase1_;
  bool Plot_TP_ver_;
 
  edm::EDGetTokenT<edm::PCaloHitContainer> tok_mc_;
  edm::EDGetTokenT<HBHEDigiCollection> tok_hbhe_;
  edm::EDGetTokenT<HODigiCollection> tok_ho_;
  edm::EDGetTokenT<HFDigiCollection> tok_hf_;
  edm::EDGetTokenT<HcalTrigPrimDigiCollection> tok_emulTPs_;
  edm::EDGetTokenT<HcalTrigPrimDigiCollection> tok_dataTPs_;
 
  edm::EDGetTokenT<QIE10DigiCollection> tok_qie10_hf_;
  edm::EDGetTokenT<QIE11DigiCollection> tok_qie11_hbhe_;
 
  edm::ESHandle<CaloGeometry> geometry;
 
  edm::ESHandle<HcalDbService> conditions;
 
  //TP Code
  edm::ESHandle<HcalTopology> htopo;
  //~TP Code
 
  int nevent1;
  int nevent2;
  int nevent3;
  int nevent4;
  int nevtot;
 
  const HcalDDDRecConstants *hcons;
  const HcalTopology *htopology;
 
  int maxDepth_[5];   // 0:any, 1:HB, 2:HE, 3:HF
  int nChannels_[5];  // 0:any, 1:HB, 2:HE,
 
  bool skipDataTPs;
};
 
#endif