WMuNuCandidatePtr.h

Go to the documentation of this file.

#ifndef AnalysisDataFormat_EWK_WMuNuCandidatePtr_h
#define AnalysisDataFormat_EWK_WMuNuCandidatePtr_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/CompositePtrCandidate.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/Candidate/interface/ShallowClonePtrCandidate.h"
namespace reco{

class WMuNuCandidatePtr : public reco::CompositePtrCandidate {
public:
  WMuNuCandidatePtr();
  WMuNuCandidatePtr(const reco::CandidatePtr,const reco::CandidatePtr);
  ~WMuNuCandidatePtr() 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... WMuNuCandidatePtrs 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::WMuNuCandidatePtr> WMuNuCandidatePtrCollection;
}


#endif