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 <math.h>

class HeavyChHiggsToTauNuSkim : public edm::EDFilter {

    public:
        explicit HeavyChHiggsToTauNuSkim(const edm::ParameterSet&);
        ~HeavyChHiggsToTauNuSkim();

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

   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