LeptonRecoSkim

Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125
// Class:      LeptonRecoSkim
//
/**\class LeptonRecoSkim LeptonRecoSkim.cc Configuration/Skimming/src/LeptonRecoSkim.cc

   Description: [one line class summary]

   Implementation:
   [Notes on implementation]
*/
//
// Original Author:  Massimiliano Chiorboli,40 4-A01,+41227671535,
//         Created:  Wed Mar 31 21:49:08 CEST 2010
// $Id: LeptonRecoSkim.h,v 1.1 2010/11/05 18:37:50 torimoto Exp $
//
//
#include <memory>
#include <fstream>

// user include files
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/one/EDFilter.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/EDGetToken.h"

//includes for reco objects
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalClusterTools.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "Geometry/Records/interface/CaloTopologyRecord.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "Geometry/CaloTopology/interface/EcalBarrelTopology.h"
#include "Geometry/CaloTopology/interface/EcalEndcapTopology.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/BTauReco/interface/JetTag.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"

//includes for PF
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"

//includes for Trigger
#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Common/interface/TriggerNames.h"

class LeptonRecoSkim : public edm::one::EDFilter<> {
public:
  explicit LeptonRecoSkim(const edm::ParameterSet&);
  ~LeptonRecoSkim() override;

private:
  void beginJob() override;
  bool filter(edm::Event&, const edm::EventSetup&) override;
  void endJob() override;

  void handleObjects(const edm::Event&, const edm::EventSetup& iSetup);

  // ----------member data ---------------------------
  const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> m_CaloGeoToken;
  const edm::ESGetToken<CaloTopology, CaloTopologyRecord> m_CaloTopoToken;

  edm::InputTag hltLabel;
  std::string filterName;

  edm::EDGetTokenT<reco::GsfElectronCollection> gsfElectronCollectionToken_;
  edm::EDGetTokenT<reco::PFCandidateCollection> pfCandidateCollectionToken_;
  edm::EDGetTokenT<reco::MuonCollection> muonCollectionToken_;
  edm::EDGetTokenT<reco::CaloJetCollection> caloJetCollectionToken_;
  edm::EDGetTokenT<reco::PFJetCollection> pfJetCollectionToken_;
  edm::EDGetTokenT<EcalRecHitCollection> ebRecHitCollectionToken_;
  edm::EDGetTokenT<EcalRecHitCollection> eeRecHitCollectionToken_;

  const reco::GsfElectronCollection* theElectronCollection;
  const reco::PFCandidateCollection* thePfCandidateCollection;
  const reco::MuonCollection* theMuonCollection;
  const reco::CaloJetCollection* theCaloJetCollection;
  const reco::PFJetCollection* thePFJetCollection;
  const EcalRecHitCollection* theEcalBarrelCollection;
  const EcalRecHitCollection* theEcalEndcapCollection;
  const CaloTopology* theCaloTopology;
  const CaloGeometry* theCaloGeometry;

  bool firstEvent;

  bool useElectronSelection;
  bool usePfElectronSelection;
  bool useMuonSelection;
  bool useHtSelection;
  bool usePFHtSelection;

  bool ElectronCutPassed;
  bool PfElectronCutPassed;
  bool MuonCutPassed;
  bool HtCutPassed;
  bool PFHtCutPassed;

  double ptElecMin;
  double ptPfElecMin;
  int nSelectedElectrons;
  int nSelectedPfElectrons;
  double ptGlobalMuonMin;
  double ptTrackerMuonMin;
  int nSelectedMuons;
  double htMin;
  double pfHtMin;
  double htJetThreshold;
  double pfHtJetThreshold;

  int NeventsTotal;
  int NeventsFiltered;
  int NHltMu9;
  int NHltDiMu3;

  int NtotalElectrons;
  int NmvaElectrons;
};