SimHitsValidationHcal.h

Go to the documentation of this file.

#ifndef ValidationSimHitsValidationHcal_H
#define ValidationSimHitsValidationHcal_H

#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "Geometry/HcalCommonData/interface/HcalDDDRecConstants.h"
#include "Geometry/Records/interface/HcalRecNumberingRecord.h"
#include "SimDataFormats/CaloHit/interface/PCaloHit.h"
#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"

#include <fstream>
#include <iostream>
#include <map>
#include <string>
#include <vector>

class SimHitsValidationHcal : public DQMEDAnalyzer {
public:
  SimHitsValidationHcal(const edm::ParameterSet &ps);
  ~SimHitsValidationHcal() override = default;
  static void fillDescriptions(edm::ConfigurationDescriptions &descriptions);

protected:
  void analyze(const edm::Event &e, const edm::EventSetup &c) override;
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;

  void analyzeHits(std::vector<PCaloHit> &);

private:
  struct energysum {
    double e25, e50, e100, e250;
    energysum() { e25 = e50 = e100 = e250 = 0.0; }
  };

  struct idType {
    idType() {
      subdet = HcalEmpty;
      z = depth1 = depth2 = 0;
    }
    idType(HcalSubdetector det, int iz, int d1, int d2) {
      subdet = det;
      z = iz;
      depth1 = d1;
      depth2 = d2;
    }
    HcalSubdetector subdet;
    int z, depth1, depth2;
  };

  struct etaRange {
    etaRange() {
      bins = 0;
      low = high = 0;
    }
    etaRange(int bin, double min, double max) {
      bins = bin;
      low = min;
      high = max;
    }
    int bins;
    double low, high;
  };

  std::vector<std::pair<std::string, std::string>> getHistogramTypes();
  etaRange getLimits(idType);
  std::pair<int, int> histId(int subdet, int eta, int depth, unsigned int dep);

  bool initialized;
  const std::string g4Label_, hcalHits_;
  const bool verbose_, testNumber_;
  const edm::EDGetTokenT<edm::PCaloHitContainer> tok_hits_;
  const edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> tok_HRNDC_;
  const HcalDDDRecConstants *hcons;
  std::vector<idType> types;
  int maxDepthHB_, maxDepthHE_;
  int maxDepthHO_, maxDepthHF_;

  std::vector<MonitorElement *> meHcalHitEta_, meHcalHitTimeEta_;
  std::vector<MonitorElement *> meHcalEnergyl25_, meHcalEnergyl50_;
  std::vector<MonitorElement *> meHcalEnergyl100_, meHcalEnergyl250_;
  MonitorElement *meEnergy_HB, *metime_HB, *metime_enweighted_HB;
  MonitorElement *meEnergy_HE, *metime_HE, *metime_enweighted_HE;
  MonitorElement *meEnergy_HO, *metime_HO, *metime_enweighted_HO;
  MonitorElement *meEnergy_HF, *metime_HF, *metime_enweighted_HF;
};

#endif