NoPileUpMEtUtilities Class Reference

#include <NoPileUpMEtUtilities.h>

Public Types
enum	{ kOutsideJet = 0, kWithin, kAll }
enum	{   kChHSMET = 0, kChPUMET, kNeutralUncMET, kHadronicHSMET,   kHadronicPUMET }

Public Member Functions
reco::PUSubMETCandInfoCollection	cleanJets (const reco::PUSubMETCandInfoCollection &, const std::vector< reco::Candidate::LorentzVector > &, double, bool)
reco::PUSubMETCandInfoCollection	cleanPFCandidates (const reco::PUSubMETCandInfoCollection &, const std::vector< reco::Candidate::LorentzVector > &, double, bool)
void	computeAllSums (const reco::PUSubMETCandInfoCollection &jets, const reco::PUSubMETCandInfoCollection &pfCandidates)
CommonMETData	computeRecoil (int metType, double &sumAbsPx, double &sumAbsPy)
void	finalizeMEtData (CommonMETData &)
	NoPileUpMEtUtilities ()
	~NoPileUpMEtUtilities ()

Private Member Functions
CommonMETData	computeCandidateSum (const reco::PUSubMETCandInfoCollection &cands, bool neutralFracOnly, double &sumAbsPx, double &sumAbsPy)
reco::PUSubMETCandInfoCollection	selectCandidates (const reco::PUSubMETCandInfoCollection &cands, double minPt, double maxPt, int type, bool isCharged, int isWithinJet)

Private Attributes
CommonMETData	chHSPfcSum_
double	chHSPfcSumAbsPx_
double	chHSPfcSumAbsPy_
CommonMETData	chPUPfcSum_
double	chPUPfcSumAbsPx_
double	chPUPfcSumAbsPy_
double	maxPtDef_
double	minPtDef_
CommonMETData	nHSJetSum_
double	nHSJetSumAbsPx_
double	nHSJetSumAbsPy_
CommonMETData	nPUJetSum_
double	nPUJetSumAbsPx_
double	nPUJetSumAbsPy_
CommonMETData	nUncPfcSum_
double	nUncPfcSumAbsPx_
double	nUncPfcSumAbsPy_

Detailed Description

Definition at line 13 of file NoPileUpMEtUtilities.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
kOutsideJet
kWithin
kAll

Definition at line 15 of file NoPileUpMEtUtilities.h.

{ kOutsideJet = 0, kWithin, kAll };

anonymous enum

anonymous enum

Enumerator
kChHSMET
kChPUMET
kNeutralUncMET
kHadronicHSMET
kHadronicPUMET

Definition at line 16 of file NoPileUpMEtUtilities.h.

{ kChHSMET = 0, kChPUMET, kNeutralUncMET, kHadronicHSMET, kHadronicPUMET };

Constructor & Destructor Documentation

NoPileUpMEtUtilities()

NoPileUpMEtUtilities::NoPileUpMEtUtilities

(

)

Definition at line 8 of file NoPileUpMEtUtilities.cc.

                                           {
  minPtDef_ = -1;
  maxPtDef_ = 1000000;
}

References maxPtDef_, and minPtDef_.

~NoPileUpMEtUtilities()

NoPileUpMEtUtilities::~NoPileUpMEtUtilities

(

)

Definition at line 13 of file NoPileUpMEtUtilities.cc.

{}

Member Function Documentation

cleanJets()

reco::PUSubMETCandInfoCollection NoPileUpMEtUtilities::cleanJets	(	const reco::PUSubMETCandInfoCollection &	jets,
		const std::vector< reco::Candidate::LorentzVector > &	leptons,
		double	dRoverlap,
		bool	invert
	)

Definition at line 28 of file NoPileUpMEtUtilities.cc.

                 {
  reco::PUSubMETCandInfoCollection retVal;
  for (reco::PUSubMETCandInfoCollection::const_iterator jet = jets.begin(); jet != jets.end(); ++jet) {
    bool isOverlap = false;
    for (std::vector<reco::Candidate::LorentzVector>::const_iterator lepton = leptons.begin(); lepton != leptons.end();
         ++lepton) {
      if (deltaR2(jet->p4(), *lepton) < dRoverlap * dRoverlap) {
        isOverlap = true;
        break;
      }
    }
    if ((!isOverlap && !invert) || (isOverlap && invert))
      retVal.push_back(*jet);
  }
  return retVal;
}

References HLTMuonOfflineAnalyzer_cfi::deltaR2, l1tstage2_dqm_sourceclient-live_cfg::invert, metsig::jet, singleTopDQM_cfi::jets, and HLT_FULL_cff::leptons.

Referenced by NoPileUpPFMEtProducer::produce().

cleanPFCandidates()

reco::PUSubMETCandInfoCollection NoPileUpMEtUtilities::cleanPFCandidates	(	const reco::PUSubMETCandInfoCollection &	pfCandidates,
		const std::vector< reco::Candidate::LorentzVector > &	leptons,
		double	dRoverlap,
		bool	invert
	)

Definition at line 81 of file NoPileUpMEtUtilities.cc.

                 {
  // invert: false = PFCandidates are required not to overlap with leptons
  //         true  = PFCandidates are required to overlap with leptons
 
  reco::PUSubMETCandInfoCollection retVal;
  for (reco::PUSubMETCandInfoCollection::const_iterator pfCandidate = pfCandidates.begin();
       pfCandidate != pfCandidates.end();
       ++pfCandidate) {
    bool isOverlap = false;
    for (std::vector<reco::Candidate::LorentzVector>::const_iterator lepton = leptons.begin(); lepton != leptons.end();
         ++lepton) {
      if (deltaR2(pfCandidate->p4(), *lepton) < dRoverlap * dRoverlap) {
        isOverlap = true;
        break;
      }
    }
    if ((!isOverlap && !invert) || (isOverlap && invert))
      retVal.push_back(*pfCandidate);
  }
  return retVal;
}

References HLTMuonOfflineAnalyzer_cfi::deltaR2, l1tstage2_dqm_sourceclient-live_cfg::invert, HLT_FULL_cff::leptons, and zmumugammaAnalyzer_cfi::pfCandidates.

Referenced by NoPileUpPFMEtProducer::produce().

computeAllSums()

void NoPileUpMEtUtilities::computeAllSums	(	const reco::PUSubMETCandInfoCollection &	jets,
		const reco::PUSubMETCandInfoCollection &	pfCandidates
	)

Definition at line 140 of file NoPileUpMEtUtilities.cc.

                                                                                              {
  reco::PUSubMETCandInfoCollection pfcsCh = selectCandidates(
      pfCandidates, minPtDef_, maxPtDef_, reco::PUSubMETCandInfo::kUndefined, false, NoPileUpMEtUtilities::kAll);
 
  reco::PUSubMETCandInfoCollection pfcsChHS = selectCandidates(
      pfCandidates, minPtDef_, maxPtDef_, reco::PUSubMETCandInfo::kChHS, true, NoPileUpMEtUtilities::kAll);
 
  reco::PUSubMETCandInfoCollection pfcsChPU = selectCandidates(
      pfCandidates, minPtDef_, maxPtDef_, reco::PUSubMETCandInfo::kChPU, true, NoPileUpMEtUtilities::kAll);
 
  reco::PUSubMETCandInfoCollection pfcsNUnclustered = selectCandidates(
      pfCandidates, minPtDef_, maxPtDef_, reco::PUSubMETCandInfo::kNeutral, false, NoPileUpMEtUtilities::kOutsideJet);
 
  reco::PUSubMETCandInfoCollection jetsHS =
      selectCandidates(jets, 10.0, maxPtDef_, reco::PUSubMETCandInfo::kHS, false, NoPileUpMEtUtilities::kAll);
  reco::PUSubMETCandInfoCollection jetsPU =
      selectCandidates(jets, 10.0, maxPtDef_, reco::PUSubMETCandInfo::kPU, false, NoPileUpMEtUtilities::kAll);
 
  //not used so far
  //_chPfcSum = computeCandidateSum(pfcsCh, false, &_chPfcSumAbsPx, &_chPfcSumAbsPy);
  chHSPfcSum_ = computeCandidateSum(pfcsChHS, false, chHSPfcSumAbsPx_, chHSPfcSumAbsPy_);
  chPUPfcSum_ = computeCandidateSum(pfcsChPU, false, chPUPfcSumAbsPx_, chPUPfcSumAbsPy_);
  nUncPfcSum_ = computeCandidateSum(pfcsNUnclustered, false, nUncPfcSumAbsPx_, nUncPfcSumAbsPy_);
 
  nHSJetSum_ = computeCandidateSum(jetsHS, true, nHSJetSumAbsPx_, nHSJetSumAbsPy_);
  nPUJetSum_ = computeCandidateSum(jetsPU, true, nPUJetSumAbsPx_, nPUJetSumAbsPy_);
}

References chHSPfcSum_, chHSPfcSumAbsPx_, chHSPfcSumAbsPy_, chPUPfcSum_, chPUPfcSumAbsPx_, chPUPfcSumAbsPy_, computeCandidateSum(), singleTopDQM_cfi::jets, kAll, reco::PUSubMETCandInfo::kChHS, reco::PUSubMETCandInfo::kChPU, reco::PUSubMETCandInfo::kHS, reco::PUSubMETCandInfo::kNeutral, kOutsideJet, reco::PUSubMETCandInfo::kPU, reco::PUSubMETCandInfo::kUndefined, maxPtDef_, minPtDef_, nHSJetSum_, nHSJetSumAbsPx_, nHSJetSumAbsPy_, nPUJetSum_, nPUJetSumAbsPx_, nPUJetSumAbsPy_, nUncPfcSum_, nUncPfcSumAbsPx_, nUncPfcSumAbsPy_, zmumugammaAnalyzer_cfi::pfCandidates, and selectCandidates().

computeCandidateSum()

CommonMETData NoPileUpMEtUtilities::computeCandidateSum ( const reco::PUSubMETCandInfoCollection &  cands, bool  neutralFracOnly, double &  sumAbsPx, double &  sumAbsPy  )

private

Definition at line 49 of file NoPileUpMEtUtilities.cc.

                                                                          {
  CommonMETData retVal;
  retVal.mex = 0.;
  retVal.mey = 0.;
  retVal.sumet = 0.;
  double retVal_sumAbsPx = 0.;
  double retVal_sumAbsPy = 0.;
  double pFrac = 1;
  for (reco::PUSubMETCandInfoCollection::const_iterator cand = cands.begin(); cand != cands.end(); ++cand) {
    pFrac = 1;
    if (neutralFracOnly)
      pFrac = (1 - cand->chargedEnFrac());
 
    retVal.mex += cand->p4().px() * pFrac;
    retVal.mey += cand->p4().py() * pFrac;
    retVal.sumet += cand->p4().pt() * pFrac;
    retVal_sumAbsPx += std::abs(cand->p4().px());
    retVal_sumAbsPy += std::abs(cand->p4().py());
  }
  finalizeMEtData(retVal);
  sumAbsPx = retVal_sumAbsPx;
  sumAbsPy = retVal_sumAbsPy;
  return retVal;
}

References funct::abs(), HLT_FULL_cff::cands, finalizeMEtData(), CommonMETData::mex, CommonMETData::mey, and CommonMETData::sumet.

Referenced by computeAllSums().

computeRecoil()

CommonMETData NoPileUpMEtUtilities::computeRecoil	(	int	metType,
		double &	sumAbsPx,
		double &	sumAbsPy
	)

Definition at line 169 of file NoPileUpMEtUtilities.cc.

                                                                                                 {
  CommonMETData retVal;
  double retSumAbsPx = 0.;
  double retSumAbsPy = 0.;
 
  //never used....
  // if(metType==NoPileUpMEtUtilities::kChMET) {
  //   CommonMETData chPfcSum = computeCandidateSum(_pfcsCh, false, &trackSumAbsPx, &trackSumAbsPy);
  //   retVal.mex   = -chPfcSum.mex;
  //   retVal.mey   = -chPfcSum.mey;
  //   retVal.sumet = chPfcSum.sumet;
  //   retSumAbsPx = ;
  //   retSumAbsPy = ;
  // }
  if (metType == NoPileUpMEtUtilities::kChHSMET) {
    retVal.mex = chHSPfcSum_.mex;
    retVal.mey = chHSPfcSum_.mey;
    retVal.sumet = chHSPfcSum_.sumet;
    retSumAbsPx = chHSPfcSumAbsPx_;
    retSumAbsPy = chHSPfcSumAbsPy_;
  }
  if (metType == NoPileUpMEtUtilities::kChPUMET) {
    retVal.mex = chPUPfcSum_.mex;
    retVal.mey = chPUPfcSum_.mey;
    retVal.sumet = chPUPfcSum_.sumet;
    retSumAbsPx = chPUPfcSumAbsPx_;
    retSumAbsPy = chPUPfcSumAbsPy_;
  }
  if (metType == NoPileUpMEtUtilities::kNeutralUncMET) {
    retVal.mex = nUncPfcSum_.mex;
    retVal.mey = nUncPfcSum_.mey;
    retVal.sumet = nUncPfcSum_.sumet;
    retSumAbsPx = nUncPfcSumAbsPx_;
    retSumAbsPy = nUncPfcSumAbsPy_;
  }
  if (metType == NoPileUpMEtUtilities::kHadronicHSMET) {
    retVal.mex = chHSPfcSum_.mex + nHSJetSum_.mex;
    retVal.mey = chHSPfcSum_.mey + nHSJetSum_.mey;
    retVal.sumet = chHSPfcSum_.sumet + nHSJetSum_.sumet;
    retSumAbsPx = chHSPfcSumAbsPx_ + nHSJetSumAbsPx_;
    retSumAbsPy = chHSPfcSumAbsPy_ + nHSJetSumAbsPy_;
  }
  if (metType == NoPileUpMEtUtilities::kHadronicPUMET) {
    retVal.mex = chPUPfcSum_.mex + nHSJetSum_.mex;
    retVal.mey = chPUPfcSum_.mey + nHSJetSum_.mey;
    retVal.sumet = chPUPfcSum_.sumet + nHSJetSum_.sumet;
    retSumAbsPx = chPUPfcSumAbsPx_ + nHSJetSumAbsPx_;
    retSumAbsPy = chPUPfcSumAbsPy_ + nHSJetSumAbsPy_;
  }
 
  sumAbsPx = retSumAbsPx;
  sumAbsPy = retSumAbsPy;
 
  return retVal;
}

References chHSPfcSum_, chHSPfcSumAbsPx_, chHSPfcSumAbsPy_, chPUPfcSum_, chPUPfcSumAbsPx_, chPUPfcSumAbsPy_, kChHSMET, kChPUMET, kHadronicHSMET, kHadronicPUMET, kNeutralUncMET, objects.autophobj::metType, CommonMETData::mex, CommonMETData::mey, nHSJetSum_, nHSJetSumAbsPx_, nHSJetSumAbsPy_, nUncPfcSum_, nUncPfcSumAbsPx_, nUncPfcSumAbsPy_, and CommonMETData::sumet.

finalizeMEtData()

void NoPileUpMEtUtilities::finalizeMEtData

(

CommonMETData &

metData

)

Definition at line 19 of file NoPileUpMEtUtilities.cc.

                                                                 {
  metData.met = sqrt(metData.mex * metData.mex + metData.mey * metData.mey);
  metData.mez = 0.;
  metData.phi = atan2(metData.mey, metData.mex);
}

References CommonMETData::met, CommonMETData::mex, CommonMETData::mey, CommonMETData::mez, CommonMETData::phi, and mathSSE::sqrt().

Referenced by computeCandidateSum().

selectCandidates()

reco::PUSubMETCandInfoCollection NoPileUpMEtUtilities::selectCandidates ( const reco::PUSubMETCandInfoCollection &  cands, double  minPt, double  maxPt, int  type, bool  isCharged, int  isWithinJet  )

private

Definition at line 107 of file NoPileUpMEtUtilities.cc.

                                                                                         {
  reco::PUSubMETCandInfoCollection retVal;
  for (reco::PUSubMETCandInfoCollection::const_iterator cand = cands.begin(); cand != cands.end(); ++cand) {
    if (isCharged && cand->charge() == 0)
      continue;
    double jetPt = cand->p4().pt();
    if (jetPt < minPt || jetPt > maxPt)
      continue;
    if (type != reco::PUSubMETCandInfo::kUndefined && cand->type() != type)
      continue;
 
    //for pf candidates
    if (isWithinJet != NoPileUpMEtUtilities::kAll && (cand->isWithinJet() != isWithinJet))
      continue;
 
    retVal.push_back(*cand);
  }
  return retVal;
}

References HLT_FULL_cff::cands, reco::btau::jetPt, kAll, reco::PUSubMETCandInfo::kUndefined, and MuonErrorMatrixAnalyzer_cfi::maxPt.

Referenced by computeAllSums().

Member Data Documentation

chHSPfcSum_

CommonMETData NoPileUpMEtUtilities::chHSPfcSum_

private

Definition at line 62 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

chHSPfcSumAbsPx_

double NoPileUpMEtUtilities::chHSPfcSumAbsPx_

private

Definition at line 68 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

chHSPfcSumAbsPy_

double NoPileUpMEtUtilities::chHSPfcSumAbsPy_

private

Definition at line 74 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

chPUPfcSum_

CommonMETData NoPileUpMEtUtilities::chPUPfcSum_

private

Definition at line 63 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

chPUPfcSumAbsPx_

double NoPileUpMEtUtilities::chPUPfcSumAbsPx_

private

Definition at line 69 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

chPUPfcSumAbsPy_

double NoPileUpMEtUtilities::chPUPfcSumAbsPy_

private

Definition at line 75 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

maxPtDef_

double NoPileUpMEtUtilities::maxPtDef_

private

Definition at line 60 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and NoPileUpMEtUtilities().

minPtDef_

double NoPileUpMEtUtilities::minPtDef_

private

Definition at line 59 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and NoPileUpMEtUtilities().

nHSJetSum_

CommonMETData NoPileUpMEtUtilities::nHSJetSum_

private

Definition at line 65 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

nHSJetSumAbsPx_

double NoPileUpMEtUtilities::nHSJetSumAbsPx_

private

Definition at line 71 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

nHSJetSumAbsPy_

double NoPileUpMEtUtilities::nHSJetSumAbsPy_

private

Definition at line 77 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

nPUJetSum_

CommonMETData NoPileUpMEtUtilities::nPUJetSum_

private

Definition at line 66 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums().

nPUJetSumAbsPx_

double NoPileUpMEtUtilities::nPUJetSumAbsPx_

private

Definition at line 72 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums().

nPUJetSumAbsPy_

double NoPileUpMEtUtilities::nPUJetSumAbsPy_

private

Definition at line 78 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums().

nUncPfcSum_

CommonMETData NoPileUpMEtUtilities::nUncPfcSum_

private

Definition at line 64 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

nUncPfcSumAbsPx_

double NoPileUpMEtUtilities::nUncPfcSumAbsPx_

private

Definition at line 70 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

nUncPfcSumAbsPy_

double NoPileUpMEtUtilities::nUncPfcSumAbsPy_

private

Definition at line 76 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().