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