|
|
|
#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 edm::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) |
| DEFINE_FWK_MODULE (ZMuMu_MCanalyzer) | |
| 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) |
| typedef edm::ValueMap<float> IsolationCollection |
Definition at line 169 of file ZMuMu_MCanalyzer.cc.
| 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 137 of file ZMuMu_MCanalyzer.cc.
References Exception, edm::errors::InvalidReference, and reco::Candidate::masterClone().
Referenced by ZMuMu_MCanalyzer::analyze().
{
const pat::Muon * mu = dynamic_cast<const pat::Muon *>(&*c->masterClone());
if(mu != 0) 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 != 0) 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;
}
| DEFINE_FWK_MODULE | ( | ZMuMu_MCanalyzer | ) |
| 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 101 of file ZMuMu_MCanalyzer.cc.
References beta, dir, pat::EcalIso, pat::HcalIso, reco::IsoDeposit::sumWithin(), and pat::TrackIso.
{
// 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;
}
1.7.3