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/Records/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.
inline 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_;
 
  edm::ESGetToken<CaloTopology, CaloTopologyRecord> caloTopoToken_;
  edm::ESGetToken<CaloGeometry, CaloGeometryRecord> caloGeomToken_;
 
  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