ElectronConversionRejectionValidator.h

Go to the documentation of this file.

#ifndef ElectronConversionRejectionValidator_H
#define ElectronConversionRejectionValidator_H
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/Common/interface/RefToBaseVector.h"
#include "DataFormats/Common/interface/RefVector.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/EgammaCandidates/interface/ConversionFwd.h"
//
//DQM services
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

//
#include <map>
#include <vector>

// forward declarations
namespace reco {
  class BeamSpot;
}
class TFile;
class TH1F;
class TH2F;
class TProfile;
class TTree;
class SimVertex;
class SimTrack;
class ElectronConversionRejectionValidator : public DQMEDAnalyzer {
public:
  explicit ElectronConversionRejectionValidator(const edm::ParameterSet&);
  ~ElectronConversionRejectionValidator() override;

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

private:
  std::string fName_;

  int verbosity_;
  int nEvt_;
  int nEntry_;

  edm::ParameterSet parameters_;
  std::string conversionCollectionProducer_;
  std::string conversionCollection_;

  std::string gsfElectronCollectionProducer_;
  std::string gsfElectronCollection_;

  std::string dqmpath_;

  edm::EDGetTokenT<reco::GsfElectronCollection> gsfElecToken_;
  edm::EDGetTokenT<reco::ConversionCollection> convToken_;
  edm::EDGetTokenT<reco::VertexCollection> offline_pvToken_;
  edm::EDGetTokenT<reco::BeamSpot> beamspotToken_;

  bool isRunCentrally_;

  float elePtMin_;
  int eleExpectedHitsInnerMax_;
  float eleD0Max_;

  //
  MonitorElement* h_elePtAll_;
  MonitorElement* h_eleEtaAll_;
  MonitorElement* h_elePhiAll_;

  MonitorElement* h_elePtPass_;
  MonitorElement* h_eleEtaPass_;
  MonitorElement* h_elePhiPass_;

  MonitorElement* h_elePtFail_;
  MonitorElement* h_eleEtaFail_;
  MonitorElement* h_elePhiFail_;

  MonitorElement* h_elePtEff_;
  MonitorElement* h_eleEtaEff_;
  MonitorElement* h_elePhiEff_;

  MonitorElement* h_convPt_;
  MonitorElement* h_convEta_;
  MonitorElement* h_convPhi_;
  MonitorElement* h_convRho_;
  MonitorElement* h_convZ_;

  MonitorElement* h_convProb_;

  MonitorElement* h_convLeadTrackpt_;
  MonitorElement* h_convTrailTrackpt_;
  MonitorElement* h_convLog10TrailTrackpt_;
  MonitorElement* h_convLeadTrackAlgo_;
  MonitorElement* h_convTrailTrackAlgo_;
};

#endif