EcalPreshowerSimHitsValidation.h

Go to the documentation of this file.

#ifndef EcalPreshowerSimHitsValidation_H
#define EcalPreshowerSimHitsValidation_H

/*
 * \file EcalPreshowerSimHitsValidation.h
 *
 * \author C.Rovelli
 *
 */

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/ESHandle.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/DQMEDAnalyzer.h"

#include "SimDataFormats/CaloHit/interface/PCaloHit.h"
#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/ValidationFormats/interface/PValidationFormats.h"

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

class EcalPreshowerSimHitsValidation : public DQMEDAnalyzer {
  typedef std::map<uint32_t, float, std::less<uint32_t>> MapType;

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

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

  void analyze(const edm::Event &e, const edm::EventSetup &c) override;

private:
  std::string g4InfoLabel;
  std::string EEHitsCollection;
  std::string ESHitsCollection;

  edm::EDGetTokenT<edm::HepMCProduct> HepMCToken;
  edm::EDGetTokenT<edm::PCaloHitContainer> EEHitsToken;
  edm::EDGetTokenT<edm::PCaloHitContainer> ESHitsToken;

  bool verbose_;

  MonitorElement *menESHits1zp_;
  MonitorElement *menESHits2zp_;

  MonitorElement *menESHits1zm_;
  MonitorElement *menESHits2zm_;

  MonitorElement *meEShitLog10Energy_;
  MonitorElement *meEShitLog10EnergyNorm_;

  MonitorElement *meESEnergyHits1zp_;
  MonitorElement *meESEnergyHits2zp_;

  MonitorElement *meESEnergyHits1zm_;
  MonitorElement *meESEnergyHits2zm_;

  MonitorElement *meE1alphaE2zp_;
  MonitorElement *meE1alphaE2zm_;

  MonitorElement *meEEoverESzp_;
  MonitorElement *meEEoverESzm_;

  MonitorElement *me2eszpOver1eszp_;
  MonitorElement *me2eszmOver1eszm_;
};

#endif