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

Enumerator
kOutsideJet
kWithin
kAll

Definition at line 18 of file NoPileUpMEtUtilities.h.

{kOutsideJet=0,kWithin, kAll};

anonymous enum

Enumerator
kChHSMET
kChPUMET
kNeutralUncMET
kHadronicHSMET
kHadronicPUMET

Definition at line 19 of file NoPileUpMEtUtilities.h.

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

Constructor & Destructor Documentation

NoPileUpMEtUtilities::NoPileUpMEtUtilities

(

)

Definition at line 9 of file NoPileUpMEtUtilities.cc.

References maxPtDef_, and minPtDef_.

                                           {
  minPtDef_=-1;
  maxPtDef_=1000000;
}

NoPileUpMEtUtilities::~NoPileUpMEtUtilities

(

)

Definition at line 14 of file NoPileUpMEtUtilities.cc.

                                             {
}

Member Function Documentation

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

Definition at line 32 of file NoPileUpMEtUtilities.cc.

References reco::deltaR2(), TopAnalysis_sequences::invert, and metsig::jet.

Referenced by NoPileUpPFMEtProducer::produce().

                                               {
  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;
}

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

Definition at line 86 of file NoPileUpMEtUtilities.cc.

References reco::deltaR2(), and TopAnalysis_sequences::invert.

Referenced by NoPileUpPFMEtProducer::produce().

                                                   {
// 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;
}

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

Definition at line 139 of file NoPileUpMEtUtilities.cc.

References chHSPfcSum_, chHSPfcSumAbsPx_, chHSPfcSumAbsPy_, chPUPfcSum_, chPUPfcSumAbsPx_, chPUPfcSumAbsPy_, computeCandidateSum(), 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_, and selectCandidates().

                                                                          {
  
  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_);
  
}

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

private

Definition at line 53 of file NoPileUpMEtUtilities.cc.

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

Referenced by computeAllSums().

                                                                                {
  
  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;
}

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

Definition at line 175 of file NoPileUpMEtUtilities.cc.

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

                                                                                   {
  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;
}

void NoPileUpMEtUtilities::finalizeMEtData

(

CommonMETData &

metData

)

Definition at line 22 of file NoPileUpMEtUtilities.cc.

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

Referenced by computeCandidateSum().

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

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

private

Definition at line 110 of file NoPileUpMEtUtilities.cc.

References reco::btau::jetPt, kAll, and reco::PUSubMETCandInfo::kUndefined.

Referenced by computeAllSums().

                                                        {
  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;
}

Member Data Documentation

CommonMETData NoPileUpMEtUtilities::chHSPfcSum_

private

Definition at line 64 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::chHSPfcSumAbsPx_

private

Definition at line 70 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::chHSPfcSumAbsPy_

private

Definition at line 76 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

CommonMETData NoPileUpMEtUtilities::chPUPfcSum_

private

Definition at line 65 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::chPUPfcSumAbsPx_

private

Definition at line 71 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::chPUPfcSumAbsPy_

private

Definition at line 77 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::maxPtDef_

private

Definition at line 62 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and NoPileUpMEtUtilities().

double NoPileUpMEtUtilities::minPtDef_

private

Definition at line 61 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and NoPileUpMEtUtilities().

CommonMETData NoPileUpMEtUtilities::nHSJetSum_

private

Definition at line 67 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::nHSJetSumAbsPx_

private

Definition at line 73 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::nHSJetSumAbsPy_

private

Definition at line 79 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

CommonMETData NoPileUpMEtUtilities::nPUJetSum_

private

Definition at line 68 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums().

double NoPileUpMEtUtilities::nPUJetSumAbsPx_

private

Definition at line 74 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums().

double NoPileUpMEtUtilities::nPUJetSumAbsPy_

private

Definition at line 80 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums().

CommonMETData NoPileUpMEtUtilities::nUncPfcSum_

private

Definition at line 66 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::nUncPfcSumAbsPx_

private

Definition at line 72 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().

double NoPileUpMEtUtilities::nUncPfcSumAbsPy_

private

Definition at line 78 of file NoPileUpMEtUtilities.h.

Referenced by computeAllSums(), and computeRecoil().