DQMSourcePi0.h

Go to the documentation of this file.

#ifndef DQMSourcePi0_H
#define DQMSourcePi0_H

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.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 "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloTopology/interface/EcalEndcapTopology.h"
#include "Geometry/CaloTopology/interface/EcalBarrelTopology.h"
#include "RecoEcal/EgammaCoreTools/interface/PositionCalc.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

typedef std::map<DetId, EcalRecHit> RecHitsMap;
// Less than operator for sorting EcalRecHits according to energy.
class ecalRecHitLess : public std::binary_function<EcalRecHit, EcalRecHit, bool> 
{
public:
  bool operator()(EcalRecHit x, EcalRecHit y) 
  { 
    return (x.energy() > y.energy()); 
  }
};




class DQMStore;
class MonitorElement;

class DQMSourcePi0 : public DQMEDAnalyzer {

public:

  DQMSourcePi0( const edm::ParameterSet& );
  ~DQMSourcePi0();

protected:
   
  void beginJob();

//  void beginRun(const edm::Run& r, const edm::EventSetup& c);
  virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  void analyze(const edm::Event& e, const edm::EventSetup& c) ;

  void beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg, 
                            const edm::EventSetup& context) ;

  void endLuminosityBlock(const edm::LuminosityBlock& lumiSeg, 
                          const edm::EventSetup& c);

  void endRun(const edm::Run& r, const edm::EventSetup& c);

  void endJob();

  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