EcalEndcapSimHitsValidation.h

Go to the documentation of this file.

#ifndef EcalEndcapSimHitsValidation_H
#define EcalEndcapSimHitsValidation_H
 
/*
 * \file EcalEndcapSimHitsValidation.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/ValidationFormats/interface/PValidationFormats.h"
 
#include <fstream>
#include <iostream>
#include <map>
#include <vector>
 
class EcalEndcapSimHitsValidation : public DQMEDAnalyzer {
  typedef std::map<uint32_t, float, std::less<uint32_t>> MapType;
 
public:
  EcalEndcapSimHitsValidation(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:
  uint32_t getUnitWithMaxEnergy(MapType &themap);
 
  virtual float energyInMatrixEE(int nCellInX, int nCellInY, int centralX, int centralY, int centralZ, MapType &themap);
 
  std::vector<uint32_t> getIdsAroundMax(
      int nCellInX, int nCellInY, int centralX, int centralY, int centralZ, MapType &themap);
 
  bool fillEEMatrix(
      int nCellInX, int nCellInY, int CentralX, int CentralY, int CentralZ, MapType &fillmap, MapType &themap);
 
  float eCluster2x2(MapType &themap);
  float eCluster4x4(float e33, MapType &themap);
 
  std::string g4InfoLabel;
  std::string EEHitsCollection;
  std::string ValidationCollection;
 
  edm::EDGetTokenT<edm::PCaloHitContainer> EEHitsToken;
  edm::EDGetTokenT<PEcalValidInfo> ValidationCollectionToken;
 
  bool verbose_;
 
  int myEntries;
  float eRLength[26];
 
  MonitorElement *meEEzpHits_;
  MonitorElement *meEEzmHits_;
 
  MonitorElement *meEEzpCrystals_;
  MonitorElement *meEEzmCrystals_;
 
  MonitorElement *meEEzpOccupancy_;
  MonitorElement *meEEzmOccupancy_;
 
  MonitorElement *meEELongitudinalShower_;
 
  MonitorElement *meEEHitEnergy_;
 
  MonitorElement *meEEhitLog10Energy_;
 
  MonitorElement *meEEhitLog10EnergyNorm_;
 
  MonitorElement *meEEhitLog10Energy25Norm_;
 
  MonitorElement *meEEHitEnergy2_;
 
  MonitorElement *meEEcrystalEnergy_;
  MonitorElement *meEEcrystalEnergy2_;
 
  MonitorElement *meEEe1_;
  MonitorElement *meEEe4_;
  MonitorElement *meEEe9_;
  MonitorElement *meEEe16_;
  MonitorElement *meEEe25_;
 
  MonitorElement *meEEe1oe4_;
  MonitorElement *meEEe1oe9_;
  MonitorElement *meEEe4oe9_;
  MonitorElement *meEEe9oe16_;
  MonitorElement *meEEe1oe25_;
  MonitorElement *meEEe9oe25_;
  MonitorElement *meEEe16oe25_;
};
 
#endif