/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/DataFormats/PatCandidates/interface/MET.h

Go to the documentation of this file.

//
// $Id: MET.h,v 1.20 2010/08/09 14:39:00 mbluj Exp $
//

#ifndef DataFormats_PatCandidates_MET_h
#define DataFormats_PatCandidates_MET_h

#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/GenMET.h"
#include "DataFormats/PatCandidates/interface/PATObject.h"


// Define typedefs for convenience
namespace pat {
  class MET;
  typedef std::vector<MET>              METCollection; 
  typedef edm::Ref<METCollection>       METRef; 
  typedef edm::RefVector<METCollection> METRefVector; 
}


// Class definition
namespace pat {


  class MET : public PATObject<reco::MET> {

    public:

      MET();
      MET(const reco::MET & aMET);
      MET(const edm::RefToBase<reco::MET> & aMETRef);
      MET(const edm::Ptr<reco::MET> & aMETRef);
      virtual ~MET();

      virtual MET * clone() const { return new MET(*this); }

      // ---- methods for generated MET link ----
      const reco::GenMET * genMET() const;
      void setGenMET(const reco::GenMET & gm);

      // ---- methods for MET corrections ----
      unsigned int nCorrections() const;
      enum UncorrectionType {
        uncorrNONE = -1, 
        uncorrALL = 0, 
        uncorrJES,     
        uncorrMUON,    
        uncorrTAU,    
        uncorrMAXN
      };
      float corEx(UncorrectionType ix = uncorrALL) const;
      float corEy(UncorrectionType ix = uncorrALL) const;
      float corSumEt(UncorrectionType ix = uncorrALL) const;
      float uncorrectedPt(UncorrectionType ix = uncorrALL) const;
      float uncorrectedPhi(UncorrectionType ix = uncorrALL) const;

      // ---- methods to know what the pat::MET was constructed from ----
      bool isCaloMET() const { return !caloMET_.empty(); }
      bool isPFMET() const { return !pfMET_.empty(); }
      bool isRecoMET() const { return  ( caloMET_.empty() && pfMET_.empty() ); }

      // ---- methods for CaloMET specific stuff ----
      double maxEtInEmTowers() const {return caloSpecific().MaxEtInEmTowers;}
      double maxEtInHadTowers() const {return caloSpecific().MaxEtInHadTowers;}
      double etFractionHadronic () const {return caloSpecific().EtFractionHadronic;}
      double emEtFraction() const {return caloSpecific().EtFractionEm;}
      double hadEtInHB() const {return caloSpecific().HadEtInHB;}
      double hadEtInHO() const {return caloSpecific().HadEtInHO;}
      double hadEtInHE() const {return caloSpecific().HadEtInHE;}
      double hadEtInHF() const {return caloSpecific().HadEtInHF;}
      double emEtInEB() const {return caloSpecific().EmEtInEB;}
      double emEtInEE() const {return caloSpecific().EmEtInEE;}
      double emEtInHF() const {return caloSpecific().EmEtInHF;}
      double metSignificance() const {return caloSpecific().METSignificance;}
      double CaloSETInpHF() const {return caloSpecific().CaloSETInpHF;}
      double CaloSETInmHF() const {return caloSpecific().CaloSETInmHF;}
      double CaloMETInpHF() const {return caloSpecific().CaloMETInpHF;}
      double CaloMETInmHF() const {return caloSpecific().CaloMETInmHF;}
      double CaloMETPhiInpHF() const {return caloSpecific().CaloMETPhiInpHF;}
      double CaloMETPhiInmHF() const {return caloSpecific().CaloMETPhiInmHF;}
      const SpecificCaloMETData & caloSpecific() const {
          if (!isCaloMET()) throw cms::Exception("pat::MET") << "This pat::MET has not been made from a reco::CaloMET\n";
          return caloMET_[0];
      }

      // ---- methods for PFMET specific stuff ----
      double NeutralEMFraction() const { return pfSpecific().NeutralEMFraction; }
      double NeutralHadEtFraction() const { return pfSpecific().NeutralHadFraction; }
      double ChargedEMEtFraction() const { return pfSpecific().ChargedEMFraction; }
      double ChargedHadEtFraction() const { return pfSpecific().ChargedHadFraction; }
      double MuonEtFraction() const { return pfSpecific().MuonFraction; }
      double Type6EtFraction() const { return pfSpecific().Type6Fraction; }
      double Type7EtFraction() const { return pfSpecific().Type7Fraction; }
      const SpecificPFMETData & pfSpecific() const {
          if (!isPFMET()) throw cms::Exception("pat::MET") << "This pat::MET has not been made from a reco::PFMET\n";
          return pfMET_[0];
      }

    protected:

      // ---- GenMET holder ----
      std::vector<reco::GenMET> genMET_;
      // ---- holder for CaloMET specific info ---
      std::vector<SpecificCaloMETData> caloMET_;
      // ---- holder for pfMET specific info ---
      std::vector<SpecificPFMETData> pfMET_;

      // ---- members for MET corrections ----
      struct UncorInfo {
        UncorInfo(): corEx(0), corEy(0), corSumEt(0), pt(0), phi(0) {}
        float corEx;
        float corEy;
        float corSumEt;
        float pt;
        float phi;
      };
      // uncorrection transients
      mutable std::vector<UncorInfo> uncorInfo_;
      mutable unsigned int nCorrections_;
      mutable float oldPt_;
      
    protected:

      // ---- non-public correction utilities ----
      void checkUncor_() const;
      void setPtPhi_(UncorInfo& uci) const;

  };


}

#endif