ZToLLEdmNtupleDumper.cc File Reference

#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Provenance/interface/RunLumiEventNumber.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/Candidate/interface/CandMatchMap.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "DataFormats/PatCandidates/interface/GenericParticle.h"
#include "FWCore/Utilities/interface/EDMException.h"
#include "DataFormats/PatCandidates/interface/TriggerObjectStandAlone.h"
#include "DataFormats/RecoCandidate/interface/IsoDeposit.h"
#include "DataFormats/RecoCandidate/interface/IsoDepositFwd.h"
#include "DataFormats/PatCandidates/interface/Isolation.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/RecoCandidate/interface/IsoDepositVetos.h"
#include "DataFormats/RecoCandidate/interface/IsoDepositDirection.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include <vector>
#include "FWCore/Framework/interface/MakerMacros.h"

Go to the source code of this file.

Classes
class	ZToLLEdmNtupleDumper

Functions
double	candIsolation (const reco::Candidate *c, double ptThreshold, double etEcalThreshold, double etHcalThreshold, double dRVetoTrk, double dRTrk, double dREcal, double dRHcal, double alpha, double beta, bool relativeIsolation)
template<typename T >
double	isolation (const T *t, double ptThreshold, double etEcalThreshold, double etHcalThreshold, double dRVetoTrk, double dRTrk, double dREcal, double dRHcal, double alpha, double beta, bool relativeIsolation)
static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker< edm::ParameterSetDescriptionFiller< ZToLLEdmNtupleDumper > >	s_filler__LINE__ ("ZToLLEdmNtupleDumper")
static const edm::MakerPluginFactory::PMaker< edm::WorkerMaker< ZToLLEdmNtupleDumper > >	s_maker__LINE__ ("ZToLLEdmNtupleDumper")

Function Documentation

double candIsolation	(	const reco::Candidate *	c,
		double	ptThreshold,
		double	etEcalThreshold,
		double	etHcalThreshold,
		double	dRVetoTrk,
		double	dRTrk,
		double	dREcal,
		double	dRHcal,
		double	alpha,
		double	beta,
		bool	relativeIsolation
	)

Definition at line 105 of file ZToLLEdmNtupleDumper.cc.

References HltBtagPostValidation_cff::c, Exception, edm::errors::InvalidReference, isolation(), and amptDefaultParameters_cff::mu.

Referenced by ZToMuMuIsoDepositSelector< Isolator >::operator()(), and ZToLLEdmNtupleDumper::produce().

                                             {
  const pat::Muon *mu = dynamic_cast<const pat::Muon *>(c);
  if (mu != nullptr)
    return isolation(mu,
                     ptThreshold,
                     etEcalThreshold,
                     etHcalThreshold,
                     dRVetoTrk,
                     dRTrk,
                     dREcal,
                     dRHcal,
                     alpha,
                     beta,
                     relativeIsolation);
  const pat::GenericParticle *trk = dynamic_cast<const pat::GenericParticle *>(c);
  if (trk != nullptr)
    return isolation(trk,
                     ptThreshold,
                     etEcalThreshold,
                     etHcalThreshold,
                     dRVetoTrk,
                     dRTrk,
                     dREcal,
                     dRHcal,
                     alpha,
                     beta,
                     relativeIsolation);
  throw edm::Exception(edm::errors::InvalidReference)
      << "Candidate daughter #0 is neither pat::Muons nor pat::GenericParticle\n";
  return -1;
}

template<typename T >

double isolation	(	const T *	t,
		double	ptThreshold,
		double	etEcalThreshold,
		double	etHcalThreshold,
		double	dRVetoTrk,
		double	dRTrk,
		double	dREcal,
		double	dRHcal,
		double	alpha,
		double	beta,
		bool	relativeIsolation
	)

Definition at line 59 of file ZToLLEdmNtupleDumper.cc.

References zMuMuMuonUserData::beta, DeadROC_duringRun::dir, pat::EcalIso, pat::HcalIso, reco::IsoDeposit::sumWithin(), pat::TrackIso, and ZMuMuCategoriesSequences_cff::trkIso.

Referenced by candIsolation().

                                         {
  // on 34X:
  const pat::IsoDeposit *trkIso = t->isoDeposit(pat::TrackIso);
  //  const pat::IsoDeposit * trkIso = t->trackerIsoDeposit();
  // on 34X
  const pat::IsoDeposit *ecalIso = t->isoDeposit(pat::EcalIso);
  //  const pat::IsoDeposit * ecalIso = t->ecalIsoDeposit();
  //    on 34X
  const pat::IsoDeposit *hcalIso = t->isoDeposit(pat::HcalIso);
  //    const pat::IsoDeposit * hcalIso = t->hcalIsoDeposit();

  Direction dir = Direction(t->eta(), t->phi());

  pat::IsoDeposit::AbsVetos vetosTrk;
  vetosTrk.push_back(new ConeVeto(dir, dRVetoTrk));
  vetosTrk.push_back(new ThresholdVeto(ptThreshold));

  pat::IsoDeposit::AbsVetos vetosEcal;
  vetosEcal.push_back(new ConeVeto(dir, 0.));
  vetosEcal.push_back(new ThresholdVeto(etEcalThreshold));

  pat::IsoDeposit::AbsVetos vetosHcal;
  vetosHcal.push_back(new ConeVeto(dir, 0.));
  vetosHcal.push_back(new ThresholdVeto(etHcalThreshold));

  double isovalueTrk = (trkIso->sumWithin(dRTrk, vetosTrk));
  double isovalueEcal = (ecalIso->sumWithin(dREcal, vetosEcal));
  double isovalueHcal = (hcalIso->sumWithin(dRHcal, vetosHcal));

  double iso =
      alpha * (((1 + beta) / 2 * isovalueEcal) + ((1 - beta) / 2 * isovalueHcal)) + ((1 - alpha) * isovalueTrk);
  if (relativeIsolation)
    iso /= t->pt();
  return iso;
}

static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker<edm::ParameterSetDescriptionFiller< ZToLLEdmNtupleDumper > > s_filler__LINE__ ( "ZToLLEdmNtupleDumper"  )

static

static const edm::MakerPluginFactory ::PMaker< edm::WorkerMaker< ZToLLEdmNtupleDumper > > s_maker__LINE__ ( "ZToLLEdmNtupleDumper"  )

static