HeavyChHiggsToTauNuSkim.h

Go to the documentation of this file.

#ifndef HeavyChHiggsToTauNuSkim_h
#define HeavyChHiggsToTauNuSkim_h

// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDFilter.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/BTauReco/interface/IsolatedTauTagInfo.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"

#include <cmath>

class HeavyChHiggsToTauNuSkim : public edm::EDFilter {
public:
  explicit HeavyChHiggsToTauNuSkim(const edm::ParameterSet&);
  ~HeavyChHiggsToTauNuSkim() override;

  bool filter(edm::Event&, const edm::EventSetup&) override;

private:
  double deltaPhi(double phi1, double phi2) {
    const double PI = 3.1415926535;
    // in ORCA phi = [0,2pi], in TLorentzVector phi = [-pi,pi].
    // With the conversion below deltaPhi works ok despite the
    // 2*pi difference in phi definitions.
    if (phi1 < 0)
      phi1 += 2 * PI;
    if (phi2 < 0)
      phi2 += 2 * PI;

    double dphi = fabs(phi1 - phi2);

    if (dphi > PI)
      dphi = 2 * PI - dphi;
    return dphi;
  }

  double deltaR(double eta1, double eta2, double phi1, double phi2) {
    double dphi = deltaPhi(phi1, phi2);
    double deta = fabs(eta1 - eta2);
    return sqrt(dphi * dphi + deta * deta);
  }

  bool debug;

  edm::EDGetTokenT<reco::IsolatedTauTagInfoCollection> hltTauToken;
  edm::EDGetTokenT<reco::CaloJetCollection> jetToken;
  int minNumberOfjets;
  double jetEtMin;
  double jetEtaMin;
  double jetEtaMax;
  double minDRFromTau;

  int nEvents, nSelectedEvents;
};
#endif