da/dbf/DQMSourcePi0_8h_source.html

#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