DQMSourcePi0.h

Go to the documentation of this file.

#ifndef DQMSourcePi0_H
#define DQMSourcePi0_H

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"

// Geometry
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "Geometry/CaloTopology/interface/EcalBarrelTopology.h"
#include "Geometry/CaloTopology/interface/EcalEndcapTopology.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "RecoEcal/EgammaCoreTools/interface/PositionCalc.h"

// Less than operator for sorting EcalRecHits according to energy.
bool ecalRecHitGreater(EcalRecHit x, EcalRecHit y) { return (x.energy() > y.energy()); }

class DQMSourcePi0 : public DQMEDAnalyzer {
public:
  DQMSourcePi0(const edm::ParameterSet &);
  ~DQMSourcePi0() override;

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

  void convxtalid(int &, int &);
  int diff_neta_s(int, int);
  int diff_nphi_s(int, int);

private:
  int eventCounter_;
  PositionCalc posCalculator_;

  MonitorElement *hiPhiDistrEBpi0_;

  MonitorElement *hiXDistrEEpi0_;

  MonitorElement *hiPhiDistrEBeta_;

  MonitorElement *hiXDistrEEeta_;

  MonitorElement *hiEtaDistrEBpi0_;

  MonitorElement *hiYDistrEEpi0_;

  MonitorElement *hiEtaDistrEBeta_;

  MonitorElement *hiYDistrEEeta_;

  MonitorElement *hRechitEnergyEBpi0_;

  MonitorElement *hRechitEnergyEEpi0_;

  MonitorElement *hRechitEnergyEBeta_;

  MonitorElement *hRechitEnergyEEeta_;

  MonitorElement *hEventEnergyEBpi0_;

  MonitorElement *hEventEnergyEEpi0_;

  MonitorElement *hEventEnergyEBeta_;

  MonitorElement *hEventEnergyEEeta_;

  MonitorElement *hNRecHitsEBpi0_;

  MonitorElement *hNRecHitsEEpi0_;

  MonitorElement *hNRecHitsEBeta_;

  MonitorElement *hNRecHitsEEeta_;

  MonitorElement *hMeanRecHitEnergyEBpi0_;

  MonitorElement *hMeanRecHitEnergyEEpi0_;

  MonitorElement *hMeanRecHitEnergyEBeta_;

  MonitorElement *hMeanRecHitEnergyEEeta_;

  MonitorElement *hMinvPi0EB_;

  MonitorElement *hMinvPi0EE_;

  MonitorElement *hMinvEtaEB_;

  MonitorElement *hMinvEtaEE_;

  MonitorElement *hPt1Pi0EB_;

  MonitorElement *hPt1Pi0EE_;

  MonitorElement *hPt1EtaEB_;

  MonitorElement *hPt1EtaEE_;

  MonitorElement *hPt2Pi0EB_;

  MonitorElement *hPt2Pi0EE_;

  MonitorElement *hPt2EtaEB_;

  MonitorElement *hPt2EtaEE_;

  MonitorElement *hPtPi0EB_;

  MonitorElement *hPtPi0EE_;

  MonitorElement *hPtEtaEB_;

  MonitorElement *hPtEtaEE_;

  MonitorElement *hIsoPi0EB_;

  MonitorElement *hIsoPi0EE_;

  MonitorElement *hIsoEtaEB_;

  MonitorElement *hIsoEtaEE_;

  MonitorElement *hS4S91Pi0EB_;

  MonitorElement *hS4S91Pi0EE_;

  MonitorElement *hS4S91EtaEB_;

  MonitorElement *hS4S91EtaEE_;

  MonitorElement *hS4S92Pi0EB_;

  MonitorElement *hS4S92Pi0EE_;

  MonitorElement *hS4S92EtaEB_;

  MonitorElement *hS4S92EtaEE_;

  edm::EDGetTokenT<EcalRecHitCollection> productMonitoredEBpi0_;
  edm::EDGetTokenT<EcalRecHitCollection> productMonitoredEBeta_;

  edm::EDGetTokenT<EcalRecHitCollection> productMonitoredEEpi0_;
  edm::EDGetTokenT<EcalRecHitCollection> productMonitoredEEeta_;

  int gammaCandEtaSize_;
  int gammaCandPhiSize_;

  double seleXtalMinEnergy_;
  double seleXtalMinEnergyEndCap_;

  double clusSeedThr_;
  int clusEtaSize_;
  int clusPhiSize_;

  double clusSeedThrEndCap_;

  double selePtGamma_;
  double selePtPi0_;
  double seleMinvMaxPi0_;
  double seleMinvMinPi0_;
  double seleS4S9Gamma_;
  double selePi0BeltDR_;
  double selePi0BeltDeta_;
  double selePi0Iso_;
  double ptMinForIsolation_;

  double selePtGammaEndCap_;
  double selePtPi0EndCap_;
  double seleMinvMaxPi0EndCap_;
  double seleMinvMinPi0EndCap_;
  double seleS4S9GammaEndCap_;
  double selePi0IsoEndCap_;
  double selePi0BeltDREndCap_;
  double selePi0BeltDetaEndCap_;
  double ptMinForIsolationEndCap_;

  double selePtGammaEta_;
  double selePtEta_;
  double seleS4S9GammaEta_;
  double seleS9S25GammaEta_;
  double seleMinvMaxEta_;
  double seleMinvMinEta_;
  double ptMinForIsolationEta_;
  double seleEtaIso_;
  double seleEtaBeltDR_;
  double seleEtaBeltDeta_;

  double selePtGammaEtaEndCap_;
  double seleS4S9GammaEtaEndCap_;
  double seleS9S25GammaEtaEndCap_;
  double selePtEtaEndCap_;
  double seleMinvMaxEtaEndCap_;
  double seleMinvMinEtaEndCap_;
  double ptMinForIsolationEtaEndCap_;
  double seleEtaIsoEndCap_;
  double seleEtaBeltDREndCap_;
  double seleEtaBeltDetaEndCap_;

  bool ParameterLogWeighted_;
  double ParameterX0_;
  double ParameterT0_barl_;
  double ParameterT0_endc_;
  double ParameterT0_endcPresh_;
  double ParameterW0_;

  std::vector<EBDetId> detIdEBRecHits;
  std::vector<EcalRecHit> EBRecHits;

  std::vector<EEDetId> detIdEERecHits;
  std::vector<EcalRecHit> EERecHits;

  unsigned int prescaleFactor_;

  std::string folderName_;

  bool saveToFile_;

  bool isMonEBpi0_;
  bool isMonEBeta_;
  bool isMonEEpi0_;
  bool isMonEEeta_;

  std::string fileName_;
};

#endif