BoostedTopProducer.h

Go to the documentation of this file.

#ifndef TopQuarkAnalysis_TopPairBSM_interface_BoostedTopProducer_h
#define TopQuarkAnalysis_TopPairBSM_interface_BoostedTopProducer_h

// -*- C++ -*-
//
// Package:    BoostedTopProducer
// Class:      BoostedTopProducer
//
//
// Original Author:  "Salvatore Rappoccio"
//         Created:  Thu May  1 11:37:48 CDT 2008
//
//


// system include files
#include <memory>

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

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/Candidate/interface/CompositeCandidate.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "DataFormats/PatCandidates/interface/Electron.h"
#include "DataFormats/PatCandidates/interface/Jet.h"
#include "DataFormats/PatCandidates/interface/MET.h"
#include "DataFormats/PatCandidates/interface/Hemisphere.h"

#include "METzCalculator.h"
#include "AnalysisDataFormats/TopObjects/interface/TtSemiLeptonicEvent.h"
#include "TLorentzVector.h"
#include "TMath.h"



//
// class decleration
//

class BoostedTopProducer : public edm::EDProducer {
   public:
      explicit BoostedTopProducer(const edm::ParameterSet&);
      ~BoostedTopProducer();

   private:
      virtual void beginJob(const edm::EventSetup&) ;
      virtual void produce(edm::Event&, const edm::EventSetup&);
      virtual void endJob() ;

      // ----------member data ---------------------------

      // data labels
      edm::EDGetTokenT<std::vector<pat::Electron> > eleToken_;
      edm::EDGetTokenT<std::vector<pat::Muon> > muoToken_;
      edm::EDGetTokenT<std::vector<pat::Jet> > jetToken_;
      edm::EDGetTokenT<std::vector<pat::MET> > metToken_;
      edm::EDGetTokenT<TtSemiLeptonicEvent> solToken_;

      // Cut variables
      double        caloIsoCut_;     // isolation cut to consider a lepton isolated
      double        mTop_;           // input top mass

      // Rapidity-invariant deltaR
      double Psi(const TLorentzVector& p1, const TLorentzVector& p2, double mass);
};

#endif