ZMuMu_MCanalyzer.cc File Reference

#include "DataFormats/Common/interface/AssociationVector.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/Candidate/interface/CandMatchMap.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "DataFormats/Candidate/interface/Particle.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Candidate/interface/OverlapChecker.h"
#include "DataFormats/Math/interface/deltaR.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "DataFormats/PatCandidates/interface/GenericParticle.h"
#include "DataFormats/PatCandidates/interface/TriggerObjectStandAlone.h"
#include "DataFormats/PatCandidates/interface/PATObject.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 "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include <vector>
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Candidate/interface/CandAssociation.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include <iostream>
#include <iterator>
#include <cmath>
#include "FWCore/Framework/interface/MakerMacros.h"

Go to the source code of this file.

Classes
class	ZMuMu_MCanalyzer

Typedefs
typedef ValueMap< float >	IsolationCollection

Functions
double	candidateIsolation (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< ZMuMu_MCanalyzer > >	s_filler__LINE__ ("ZMuMu_MCanalyzer")
static const edm::MakerPluginFactory::PMaker< edm::WorkerMaker< ZMuMu_MCanalyzer > >	s_maker__LINE__ ("ZMuMu_MCanalyzer")

Typedef Documentation

typedef ValueMap<float> IsolationCollection

Definition at line 47 of file ZMuMu_MCanalyzer.cc.

Function Documentation

double candidateIsolation	(	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 151 of file ZMuMu_MCanalyzer.cc.

References Exception, edm::errors::InvalidReference, isolation(), reco::Candidate::masterClone(), and amptDefaultParameters_cff::mu.

Referenced by ZMuMu_MCanalyzer::analyze().

                                                  {
  const pat::Muon* mu = dynamic_cast<const pat::Muon*>(&*c->masterClone());
  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->masterClone());
  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 105 of file ZMuMu_MCanalyzer.cc.

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

Referenced by candidateIsolation().

                                         {
  // 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< ZMuMu_MCanalyzer > > s_filler__LINE__ ( "ZMuMu_MCanalyzer"  )

static

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

static