WMuNuCandidate.h

Go to the documentation of this file.

#ifndef AnalysisDataFormats_EWK_WMuNu_h
#define AnalysisDataFormats_EWK_WMuNu_h

#include <map>
#include <memory>

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/METReco/interface/MET.h"
#include "DataFormats/METReco/interface/METFwd.h"
#include "DataFormats/Candidate/interface/CompositeCandidate.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/Candidate/interface/ShallowClonePtrCandidate.h"
namespace reco {

  class WMuNuCandidate : public reco::CompositeCandidate {
  public:
    WMuNuCandidate();
    WMuNuCandidate(edm::Ptr<reco::Muon>, edm::Ptr<reco::MET>);
    ~WMuNuCandidate() override;

    //WARNING: W Candidates combine the information from a Muon with the (px,py) information of the MET as the Neutrino
    // --> There is no Pz information!!!!
    // Be very careful when using the default Candidate functions (.mass, .mt, .et, etc). They may not be what you are looking for :-).

    // Example: Candidates have a mt() function which computes the tranverse mass from E & pz.
    // As MET does not have pz information... WMuNuCandidates have an alternative function used in the WMuNu Inclusive Analysis
    // to compute mt just from px, py:
    // Transverse Mass from px, py:
    double massT() const;
    // Transverse Energy from px, py:
    double eT() const;

    // Acoplanarity between the muon and the MET
    double acop() const;

    const reco::Muon& getMuon() const { return *muon_; };
    const reco::MET& getNeutrino() const { return *neutrino_; };

  private:
    edm::Ptr<Muon> muon_;
    edm::Ptr<MET> neutrino_;
  };

  typedef std::vector<reco::WMuNuCandidate> WMuNuCandidateCollection;
}  // namespace reco

#endif