noPileUpMEtUtilities Namespace Reference

Functions
reco::MVAMEtJetInfoCollection	cleanJets (const reco::MVAMEtJetInfoCollection &, const std::vector< reco::Candidate::LorentzVector > &, double, bool)
reco::MVAMEtPFCandInfoCollection	cleanPFCandidates (const reco::MVAMEtPFCandInfoCollection &, const std::vector< reco::Candidate::LorentzVector > &, double, bool)
CommonMETData	computeHadRecoilNoPU (const reco::MVAMEtJetInfoCollection &, const reco::MVAMEtPFCandInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computeHadRecoilPU (const reco::MVAMEtJetInfoCollection &, const reco::MVAMEtPFCandInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computeJetSum (const reco::MVAMEtJetInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computeJetSum_neutral (const reco::MVAMEtJetInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computeNeutralRecoil_unclustered (const reco::MVAMEtPFCandInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computePFCandidateSum (const reco::MVAMEtPFCandInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computeTrackRecoil (const reco::MVAMEtPFCandInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computeTrackRecoilNoPU (const reco::MVAMEtPFCandInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
CommonMETData	computeTrackRecoilPU (const reco::MVAMEtPFCandInfoCollection &, double *sumAbsPx=0, double *sumAbsPy=0)
void	finalizeMEtData (CommonMETData &)
reco::MVAMEtJetInfo	jet (const reco::MVAMEtJetInfoCollection &, unsigned)
reco::MVAMEtJetInfo	leadJet (const reco::MVAMEtJetInfoCollection &)
reco::MVAMEtJetInfoCollection	selectJets (const reco::MVAMEtJetInfoCollection &, double, double, int)
reco::MVAMEtPFCandInfoCollection	selectPFCandidates (const reco::MVAMEtPFCandInfoCollection &, double, double, int, int)
reco::MVAMEtJetInfo	subleadJet (const reco::MVAMEtJetInfoCollection &)

Function Documentation

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

Definition at line 22 of file noPileUpMEtUtilities.cc.

References reco::deltaR2(), jet(), and reco::MVAMEtJetInfo::p4_.

Referenced by NoPileUpPFMEtProducer::produce().

  {
    reco::MVAMEtJetInfoCollection retVal;
    for ( reco::MVAMEtJetInfoCollection::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;
  }

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

Definition at line 124 of file noPileUpMEtUtilities.cc.

References reco::deltaR2().

Referenced by NoPileUpPFMEtProducer::produce().

  {
    reco::MVAMEtPFCandInfoCollection retVal;
    for ( reco::MVAMEtPFCandInfoCollection::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;
  }

CommonMETData noPileUpMEtUtilities::computeHadRecoilNoPU	(	const reco::MVAMEtJetInfoCollection &	jets,
		const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 232 of file noPileUpMEtUtilities.cc.

References computeJetSum_neutral(), computeTrackRecoilNoPU(), alignCSCRings::e, finalizeMEtData(), reco::MVAMEtJetInfo::kNoPileUp, CommonMETData::mex, CommonMETData::mey, selectJets(), and CommonMETData::sumet.

  {
    double trackSumAbsPx = 0.;
    double trackSumAbsPy = 0.;
    CommonMETData trackSumNoPU = computeTrackRecoilNoPU(pfCandidates, &trackSumAbsPx, &trackSumAbsPy);
    reco::MVAMEtJetInfoCollection jetsNoPU = selectJets(jets, 10.0, 1.e+6, reco::MVAMEtJetInfo::kNoPileUp);
    double jetSumNoPUabsPx_neutral = 0.;
    double jetSumNoPUabsPy_neutral = 0.;
    CommonMETData jetSumNoPU_neutral = computeJetSum_neutral(jetsNoPU, &jetSumNoPUabsPx_neutral, &jetSumNoPUabsPy_neutral);
    CommonMETData retVal;
    retVal.mex   = trackSumNoPU.mex   + jetSumNoPU_neutral.mex;
    retVal.mey   = trackSumNoPU.mey   + jetSumNoPU_neutral.mey;
    retVal.sumet = trackSumNoPU.sumet + jetSumNoPU_neutral.sumet;
    finalizeMEtData(retVal);
    if ( sumAbsPx ) (*sumAbsPx) = trackSumAbsPx + jetSumNoPUabsPx_neutral;
    if ( sumAbsPy ) (*sumAbsPy) = trackSumAbsPy + jetSumNoPUabsPy_neutral;
    return retVal;
  }

CommonMETData noPileUpMEtUtilities::computeHadRecoilPU	(	const reco::MVAMEtJetInfoCollection &	jets,
		const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 252 of file noPileUpMEtUtilities.cc.

References computeJetSum_neutral(), computeTrackRecoilPU(), alignCSCRings::e, finalizeMEtData(), reco::MVAMEtJetInfo::kPileUp, CommonMETData::mex, CommonMETData::mey, selectJets(), and CommonMETData::sumet.

  {
    double trackSumPUabsPx = 0.;
    double trackSumPUabsPy = 0.;
    CommonMETData trackSumPU = computeTrackRecoilPU(pfCandidates, &trackSumPUabsPx, &trackSumPUabsPy);
    reco::MVAMEtJetInfoCollection jetsPU = selectJets(jets, 10.0, 1.e+6, reco::MVAMEtJetInfo::kPileUp);
    double jetSumPUabsPx_neutral = 0.;
    double jetSumPUabsPy_neutral = 0.;
    CommonMETData jetSumPU_neutral = computeJetSum_neutral(jetsPU, &jetSumPUabsPx_neutral, &jetSumPUabsPy_neutral);
    CommonMETData retVal;
    retVal.mex   = trackSumPU.mex   + jetSumPU_neutral.mex;
    retVal.mey   = trackSumPU.mey   + jetSumPU_neutral.mey;
    retVal.sumet = trackSumPU.sumet + jetSumPU_neutral.sumet;
    finalizeMEtData(retVal);
    if ( sumAbsPx ) (*sumAbsPx) = trackSumPUabsPx + jetSumPUabsPx_neutral;
    if ( sumAbsPy ) (*sumAbsPy) = trackSumPUabsPy + jetSumPUabsPy_neutral;
    return retVal;
  }

CommonMETData noPileUpMEtUtilities::computeJetSum	(	const reco::MVAMEtJetInfoCollection &	jets,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 56 of file noPileUpMEtUtilities.cc.

References funct::abs(), finalizeMEtData(), jet(), CommonMETData::mex, CommonMETData::mey, reco::MVAMEtJetInfo::p4_, and CommonMETData::sumet.

  {
    CommonMETData retVal;
    retVal.mex   = 0.;
    retVal.mey   = 0.;
    retVal.sumet = 0.;
    double retVal_sumAbsPx = 0.;
    double retVal_sumAbsPy = 0.;
    for ( reco::MVAMEtJetInfoCollection::const_iterator jet = jets.begin();
      jet != jets.end(); ++jet ) {
      retVal.mex   += jet->p4_.px();
      retVal.mey   += jet->p4_.py();
      retVal.sumet += jet->p4_.pt();
      retVal_sumAbsPx += std::abs(jet->p4_.px());
      retVal_sumAbsPy += std::abs(jet->p4_.py());
    }
    finalizeMEtData(retVal);
    if ( sumAbsPx ) (*sumAbsPx) = retVal_sumAbsPx;
    if ( sumAbsPy ) (*sumAbsPy) = retVal_sumAbsPy;
    return retVal;
  }

CommonMETData noPileUpMEtUtilities::computeJetSum_neutral	(	const reco::MVAMEtJetInfoCollection &	jets,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 78 of file noPileUpMEtUtilities.cc.

References funct::abs(), finalizeMEtData(), jet(), CommonMETData::mex, CommonMETData::mey, reco::MVAMEtJetInfo::neutralEnFrac_, reco::MVAMEtJetInfo::p4_, and CommonMETData::sumet.

Referenced by computeHadRecoilNoPU(), and computeHadRecoilPU().

  {
    CommonMETData retVal;
    retVal.mex   = 0.;
    retVal.mey   = 0.;
    retVal.sumet = 0.;
    double retVal_sumAbsPx = 0.;
    double retVal_sumAbsPy = 0.;
    for ( reco::MVAMEtJetInfoCollection::const_iterator jet = jets.begin();
      jet != jets.end(); ++jet ) {
      retVal.mex   += (jet->p4_.px()*jet->neutralEnFrac_);
      retVal.mey   += (jet->p4_.py()*jet->neutralEnFrac_);
      retVal.sumet += (jet->p4_.pt()*jet->neutralEnFrac_);
      retVal_sumAbsPx += std::abs(jet->p4_.px());
      retVal_sumAbsPy += std::abs(jet->p4_.py());
    }
    finalizeMEtData(retVal);
    if ( sumAbsPx ) (*sumAbsPx) = retVal_sumAbsPx;
    if ( sumAbsPy ) (*sumAbsPy) = retVal_sumAbsPy;
    return retVal;
  }

CommonMETData noPileUpMEtUtilities::computeNeutralRecoil_unclustered	(	const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 225 of file noPileUpMEtUtilities.cc.

References computePFCandidateSum(), reco::MVAMEtPFCandInfo::kNeutral, and selectPFCandidates().

  {
    reco::MVAMEtPFCandInfoCollection neutralPFCandidates_unclustered = selectPFCandidates(
      pfCandidates, -0.5, +0.5, reco::MVAMEtPFCandInfo::kNeutral, 0);
    return computePFCandidateSum(neutralPFCandidates_unclustered, sumAbsPx, sumAbsPy);
  }

CommonMETData noPileUpMEtUtilities::computePFCandidateSum	(	const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 165 of file noPileUpMEtUtilities.cc.

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

Referenced by computeNeutralRecoil_unclustered(), computeTrackRecoil(), computeTrackRecoilNoPU(), and computeTrackRecoilPU().

  {
    CommonMETData retVal;
    retVal.mex   = 0.;
    retVal.mey   = 0.;
    retVal.sumet = 0.;
    double retVal_sumAbsPx = 0.;
    double retVal_sumAbsPy = 0.;
    for ( reco::MVAMEtPFCandInfoCollection::const_iterator pfCandidate = pfCandidates.begin();
      pfCandidate != pfCandidates.end(); ++pfCandidate ) {
      retVal.mex   += pfCandidate->p4_.px();
      retVal.mey   += pfCandidate->p4_.py();
      retVal.sumet += pfCandidate->p4_.pt();
      retVal_sumAbsPx += std::abs(pfCandidate->p4_.px());
      retVal_sumAbsPy += std::abs(pfCandidate->p4_.py());
    }
    finalizeMEtData(retVal);
    if ( sumAbsPx ) (*sumAbsPx) = retVal_sumAbsPx;
    if ( sumAbsPy ) (*sumAbsPy) = retVal_sumAbsPy;
    return retVal;
  }

CommonMETData noPileUpMEtUtilities::computeTrackRecoil	(	const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 194 of file noPileUpMEtUtilities.cc.

References PFRecoTauChargedHadronBuilderPlugins_cfi::chargedPFCandidates, computePFCandidateSum(), alignCSCRings::e, finalizeMEtData(), reco::MVAMEtPFCandInfo::kUndefined, CommonMETData::mex, CommonMETData::mey, selectPFCandidates(), and CommonMETData::sumet.

  {
    reco::MVAMEtPFCandInfoCollection chargedPFCandidates = selectPFCandidates(
      pfCandidates, 0.5, 1.e+3, reco::MVAMEtPFCandInfo::kUndefined, -1);
    double trackSumAbsPx = 0.;
    double trackSumAbsPy = 0.;
    CommonMETData trackSum = computePFCandidateSum(chargedPFCandidates, &trackSumAbsPx, &trackSumAbsPy);
    CommonMETData retVal;
    retVal.mex   = -trackSum.mex;
    retVal.mey   = -trackSum.mey;
    retVal.sumet =  trackSum.sumet;
    finalizeMEtData(retVal);
    if ( sumAbsPx ) (*sumAbsPx) = trackSumAbsPx;
    if ( sumAbsPy ) (*sumAbsPy) = trackSumAbsPy;
    return retVal;
  }

CommonMETData noPileUpMEtUtilities::computeTrackRecoilNoPU	(	const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 211 of file noPileUpMEtUtilities.cc.

References computePFCandidateSum(), alignCSCRings::e, reco::MVAMEtPFCandInfo::kNoPileUpCharged, and selectPFCandidates().

Referenced by computeHadRecoilNoPU().

  {
    reco::MVAMEtPFCandInfoCollection chargedPFCandidatesNoPU = selectPFCandidates(
      pfCandidates, 0.5, 1.e+3, reco::MVAMEtPFCandInfo::kNoPileUpCharged, -1);
    return computePFCandidateSum(chargedPFCandidatesNoPU, sumAbsPx, sumAbsPy);
  }

CommonMETData noPileUpMEtUtilities::computeTrackRecoilPU	(	const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double *	sumAbsPx = 0,
		double *	sumAbsPy = 0
	)

Definition at line 218 of file noPileUpMEtUtilities.cc.

References computePFCandidateSum(), alignCSCRings::e, reco::MVAMEtPFCandInfo::kPileUpCharged, and selectPFCandidates().

Referenced by computeHadRecoilPU().

  {
    reco::MVAMEtPFCandInfoCollection chargedPFCandidatesPU = selectPFCandidates(
      pfCandidates, 0.5, 1.e+3, reco::MVAMEtPFCandInfo::kPileUpCharged, -1);
    return computePFCandidateSum(chargedPFCandidatesPU, sumAbsPx, sumAbsPy);
  }

void noPileUpMEtUtilities::finalizeMEtData

(

CommonMETData &

metData

)

Definition at line 12 of file noPileUpMEtUtilities.cc.

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

Referenced by computeHadRecoilNoPU(), computeHadRecoilPU(), computeJetSum(), computeJetSum_neutral(), computePFCandidateSum(), and computeTrackRecoil().

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

reco::MVAMEtJetInfo noPileUpMEtUtilities::jet	(	const reco::MVAMEtJetInfoCollection &	jets,
		unsigned	idx
	)

Definition at line 100 of file noPileUpMEtUtilities.cc.

References customizeTrackingMonitorSeedNumber::idx, fwrapper::jets, and python.multivaluedict::sort().

Referenced by cleanJets(), computeJetSum(), computeJetSum_neutral(), leadJet(), selectJets(), and subleadJet().

  {
    reco::MVAMEtJetInfo retVal;
    if ( idx < jets.size() ) {
      reco::MVAMEtJetInfoCollection jets_sorted = jets;
      std::sort(jets_sorted.begin(), jets_sorted.end());
      retVal = jets_sorted[idx];
    }
    return retVal;
  }

reco::MVAMEtJetInfo noPileUpMEtUtilities::leadJet

(

const reco::MVAMEtJetInfoCollection &

jets

)

Definition at line 111 of file noPileUpMEtUtilities.cc.

References jet().

Referenced by EwkMuDQM::analyze().

  {
    return jet(jets, 0);
  }

reco::MVAMEtJetInfoCollection noPileUpMEtUtilities::selectJets	(	const reco::MVAMEtJetInfoCollection &	jets,
		double	minJetPt,
		double	maxJetPt,
		int	type
	)

Definition at line 42 of file noPileUpMEtUtilities.cc.

References jet(), reco::btau::jetPt, reco::MVAMEtJetInfo::kUndefined, reco::MVAMEtJetInfo::p4_, and reco::MVAMEtJetInfo::type_.

Referenced by computeHadRecoilNoPU(), and computeHadRecoilPU().

  {
    reco::MVAMEtJetInfoCollection retVal;
    for ( reco::MVAMEtJetInfoCollection::const_iterator jet = jets.begin();
      jet != jets.end(); ++jet ) {
      double jetPt = jet->p4_.pt();      
      if (  jetPt > minJetPt && 
        jetPt < maxJetPt && 
       (type == reco::MVAMEtJetInfo::kUndefined || jet->type_ == type) ) retVal.push_back(*jet);
    }
    return retVal;
  }

reco::MVAMEtPFCandInfoCollection noPileUpMEtUtilities::selectPFCandidates	(	const reco::MVAMEtPFCandInfoCollection &	pfCandidates,
		double	minCharge,
		double	maxCharge,
		int	type,
		int	isWithinJet
	)

Definition at line 146 of file noPileUpMEtUtilities.cc.

References funct::abs(), and reco::MVAMEtJetInfo::kUndefined.

Referenced by computeNeutralRecoil_unclustered(), computeTrackRecoil(), computeTrackRecoilNoPU(), and computeTrackRecoilPU().

  {
    reco::MVAMEtPFCandInfoCollection retVal;
    for ( reco::MVAMEtPFCandInfoCollection::const_iterator pfCandidate = pfCandidates.begin();
      pfCandidate != pfCandidates.end(); ++pfCandidate ) {
      int charge = abs(pfCandidate->charge_);
      if (  charge > minCharge &&
        charge < maxCharge &&
       (type == reco::MVAMEtJetInfo::kUndefined || pfCandidate->type_ == type) &&
       (isWithinJet == -1 || (isWithinJet == 1 && pfCandidate->isWithinJet_) || (isWithinJet == 0 && !pfCandidate->isWithinJet_)) )
    retVal.push_back(*pfCandidate);
    }
    return retVal;
  }

reco::MVAMEtJetInfo noPileUpMEtUtilities::subleadJet

(

const reco::MVAMEtJetInfoCollection &

jets

)

Definition at line 116 of file noPileUpMEtUtilities.cc.

References jet().

  {
    return jet(jets, 1);
  }