#include <METSignificance.h>

Public Member Functions
reco::METCovMatrix	getCovariance (const edm::View< reco::Jet > &jets, const std::vector< edm::Handle< reco::CandidateView > > &leptons, const edm::Handle< edm::View< reco::Candidate > > &pfCandidates, double rho, JME::JetResolution &resPtObj, JME::JetResolution &resPhiObj, JME::JetResolutionScaleFactor &resSFObj, bool isRealData)

	METSignificance (const edm::ParameterSet &iConfig)

	~METSignificance ()

Static Public Member Functions
static double	getSignificance (const reco::METCovMatrix &cov, const reco::MET &met)

Private Member Functions
bool	cleanJet (const reco::Jet &jet, const std::vector< edm::Handle< reco::CandidateView > > &leptons)

Private Attributes
double	dR2match_

std::vector< double >	jetEtas_

std::vector< double >	jetParams_

double	jetThreshold_

std::vector< double >	pjetParams_

Detailed Description

Definition at line 35 of file METSignificance.h.

Constructor & Destructor Documentation

METSignificance::METSignificance ( const edm::ParameterSet & iConfig )

Definition at line 21 of file METSignificance.cc.

References dR2match_, CaloRecoTauDiscriminationAgainstMuon_cfi::dRmatch, edm::ParameterSet::getParameter(), jetEtas_, jetParams_, jetThreshold_, and pjetParams_.

                                                                      {
 
   edm::ParameterSet cfgParams = iConfig.getParameter<edm::ParameterSet>("parameters");
 
   double dRmatch = cfgParams.getParameter<double>("dRMatch");
   dR2match_ = dRmatch*dRmatch;
   
   jetThreshold_ = cfgParams.getParameter<double>("jetThreshold");
   jetEtas_ = cfgParams.getParameter<std::vector<double> >("jeta");
   jetParams_ = cfgParams.getParameter<std::vector<double> >("jpar");
   pjetParams_ = cfgParams.getParameter<std::vector<double> >("pjpar");
   
 }

METSignificance::~METSignificance ( )

Definition at line 35 of file METSignificance.cc.

35 {

36 }

Member Function Documentation

bool METSignificance::cleanJet	(	const reco::Jet &	jet,
		const std::vector< edm::Handle< reco::CandidateView > > &	leptons
	)

private

Definition at line 192 of file METSignificance.cc.

References reco::deltaR2(), dR2match_, and AK4GenJetFlavourInfos_cfi::leptons.

Referenced by getCovariance().

                                                                 {
 
   for ( std::vector< edm::Handle<reco::CandidateView> >::const_iterator lep_i = leptons.begin();
         lep_i != leptons.end(); ++lep_i ) {
      for( reco::CandidateView::const_iterator lep = (*lep_i)->begin(); lep != (*lep_i)->end(); lep++ ){
         if ( reco::deltaR2(*lep, jet) < dR2match_ ) {
            return false;
         }
      }
   }
   return true;
 }

reco::METCovMatrix METSignificance::getCovariance	(	const edm::View< reco::Jet > &	jets,
		const std::vector< edm::Handle< reco::CandidateView > > &	leptons,
		const edm::Handle< edm::View< reco::Candidate > > &	pfCandidates,
		double	rho,
		JME::JetResolution &	resPtObj,
		JME::JetResolution &	resPhiObj,
		JME::JetResolutionScaleFactor &	resSFObj,
		bool	isRealData
	)

Definition at line 40 of file METSignificance.cc.

References funct::abs(), edm::View< T >::begin(), EnergyCorrector::c, cleanJet(), edm::View< T >::end(), JME::JetResolution::getResolution(), mps_fire::i, metsig::jet, METSignificanceParams_cfi::jeta, jetEtas_, jetParams_, jetThreshold_, AK4GenJetFlavourInfos_cfi::leptons, gen::n, slimmedMuons_cfi::pfCandidates, pjetParams_, edm::View< T >::ptrAt(), alignCSCRings::s, Scenarios_cff::scale, JME::JetParameters::setJetEta(), JME::JetParameters::setJetPt(), JME::JetParameters::setRho(), edm::View< T >::size(), and TtFullHadEvtBuilder_cfi::sumPt.

Referenced by cms::PFMETProducer::getMETCovMatrix(), pat::PATMETProducer::getMETCovMatrix(), and cms::METSignificanceProducer::produce().

                                         {
 
   //pfcandidates
   const edm::View<reco::Candidate>* pfCandidates=pfCandidatesH.product();
 
    // metsig covariance
    double cov_xx = 0;
    double cov_xy = 0;
    double cov_yy = 0;
  
    // for lepton and jet subtraction
    std::set<reco::CandidatePtr> footprint;
 
    // subtract leptons out of sumPt
    for ( std::vector< edm::Handle<reco::CandidateView> >::const_iterator lep_i = leptons.begin();
          lep_i != leptons.end(); ++lep_i ) {
      for( reco::CandidateView::const_iterator lep = (*lep_i)->begin(); lep != (*lep_i)->end(); lep++ ){
        if( lep->pt() > 10 ){
      for( unsigned int n=0; n < lep->numberOfSourceCandidatePtrs(); n++ ){
        if( lep->sourceCandidatePtr(n).isNonnull() and lep->sourceCandidatePtr(n).isAvailable() ){
          footprint.insert(lep->sourceCandidatePtr(n));
        } 
      }
        }
      }
    }
    // subtract jets out of sumPt
    for(edm::View<reco::Jet>::const_iterator jet = jets.begin(); jet != jets.end(); ++jet) {
 
      // disambiguate jets and leptons
      if(!cleanJet(*jet, leptons) ) continue;
      for( unsigned int n=0; n < jet->numberOfSourceCandidatePtrs(); n++){
        if( jet->sourceCandidatePtr(n).isNonnull() and jet->sourceCandidatePtr(n).isAvailable() ){
 
      footprint.insert(jet->sourceCandidatePtr(n));
        }
      }
 
    }
 
    // calculate sumPt
    double sumPt = 0;
    for(size_t i = 0; i< pfCandidates->size();  ++i) {
      
      // check if candidate exists in a lepton or jet
      bool cleancand = true;
      if(footprint.find( pfCandidates->ptrAt(i) )==footprint.end()) {
 
        //dP4 recovery
        for( std::set<reco::CandidatePtr>::const_iterator it=footprint.begin();it!=footprint.end();it++) {
      if( ((*it)->p4()-(*pfCandidates)[i].p4()).Et2()<0.000025 ){
        cleancand = false;
        break;
      }
        }
        // if not, add to sumPt
        if( cleancand ){
      sumPt += (*pfCandidates)[i].pt();
        }
      }
    }
    
    // add jets to metsig covariance matrix and subtract them from sumPt
    for(edm::View<reco::Jet>::const_iterator jet = jets.begin(); jet != jets.end(); ++jet) {
      
      // disambiguate jets and leptons
      if(!cleanJet(*jet, leptons) ) continue;
 
       double jpt  = jet->pt();
       double jeta = jet->eta();
       double feta = std::abs(jeta);
       double c = jet->px()/jet->pt();
       double s = jet->py()/jet->pt();
 
       JME::JetParameters parameters;
       parameters.setJetPt(jpt).setJetEta(jeta).setRho(rho);
 
       // jet energy resolutions
       double sigmapt = resPtObj.getResolution(parameters);
       double sigmaphi = resPhiObj.getResolution(parameters);
       // SF not needed since is already embedded in the sigma in the dataGlobalTag
       //      double sigmaSF = isRealData ? resSFObj.getScaleFactor(parameters) : 1.0;
 
       // split into high-pt and low-pt sector
       if( jpt > jetThreshold_ ){
          // high-pt jets enter into the covariance matrix via JER
 
          double scale = 0;
          if(feta<jetEtas_[0]) scale = jetParams_[0];
          else if(feta<jetEtas_[1]) scale = jetParams_[1];
          else if(feta<jetEtas_[2]) scale = jetParams_[2];
          else if(feta<jetEtas_[3]) scale = jetParams_[3];
          else scale = jetParams_[4];
 
      //         double dpt = sigmaSF*scale*jpt*sigmapt;
          double dpt = scale*jpt*sigmapt;
          double dph = jpt*sigmaphi;
 
          cov_xx += dpt*dpt*c*c + dph*dph*s*s;
          cov_xy += (dpt*dpt-dph*dph)*c*s;
          cov_yy += dph*dph*c*c + dpt*dpt*s*s;
 
       } else {
 
          // add the (corrected) jet to the sumPt
          sumPt += jpt;
 
       }
 
    }
 
    //protection against unphysical events
    if(sumPt<0) sumPt=0;
  
    // add pseudo-jet to metsig covariance matrix
    cov_xx += pjetParams_[0]*pjetParams_[0] + pjetParams_[1]*pjetParams_[1]*sumPt;
    cov_yy += pjetParams_[0]*pjetParams_[0] + pjetParams_[1]*pjetParams_[1]*sumPt;
 
    reco::METCovMatrix cov;
    cov(0,0) = cov_xx;
    cov(1,0) = cov_xy;
    cov(0,1) = cov_xy;
    cov(1,1) = cov_yy;
 
    return cov;
 }

double METSignificance::getSignificance	(	const reco::METCovMatrix &	cov,
		const reco::MET &	met
	)

static

Definition at line 175 of file METSignificance.cc.

References reco::LeafCandidate::px(), and reco::LeafCandidate::py().

Referenced by pat::PATMETProducer::produce(), cms::METSignificanceProducer::produce(), and CorrectedPATMETProducer::produce().

                                                                                       {
 
    // covariance matrix determinant
    double det = cov(0,0)*cov(1,1) - cov(0,1)*cov(1,0);
 
    // invert matrix
    double ncov_xx = cov(1,1) / det;
    double ncov_xy = -cov(0,1) / det;
    double ncov_yy = cov(0,0) / det;
 
    // product of met and inverse of covariance
    double sig = met.px()*met.px()*ncov_xx + 2*met.px()*met.py()*ncov_xy + met.py()*met.py()*ncov_yy;
 
    return sig;
 }