|
|
|
#include <HLTFiltCand.h>
Public Member Functions | |
| HLTFiltCand (const edm::ParameterSet &) | |
| virtual bool | hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) |
| ~HLTFiltCand () | |
Private Attributes | |
| edm::InputTag | ecalTag_ |
| edm::InputTag | elecTag_ |
| edm::InputTag | jetsTag_ |
| edm::InputTag | metsTag_ |
| edm::InputTag | mhtsTag_ |
| double | min_Pt_ |
| edm::InputTag | muonTag_ |
| edm::InputTag | photTag_ |
| edm::InputTag | tausTag_ |
| edm::InputTag | trckTag_ |
This class is an HLTFilter (-> EDFilter) implementing a very basic HLT trigger acting on candidates, requiring a g/e/m/j tuple above pt cuts
See header file for documentation
This class is an HLTFilter (-> EDFilter) implementing a minimum-bias HLT trigger acting on candidates, requiring tracks in Pixel det
Definition at line 24 of file HLTFiltCand.h.
| HLTFiltCand::HLTFiltCand | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 45 of file HLTFiltCand.cc.
References ecalTag_, elecTag_, edm::InputTag::encode(), jetsTag_, LogDebug, metsTag_, mhtsTag_, min_Pt_, muonTag_, photTag_, tausTag_, and trckTag_.
: HLTFilter(iConfig), photTag_ (iConfig.getParameter<edm::InputTag>("photTag")), elecTag_ (iConfig.getParameter<edm::InputTag>("elecTag")), muonTag_ (iConfig.getParameter<edm::InputTag>("muonTag")), tausTag_ (iConfig.getParameter<edm::InputTag>("tausTag")), jetsTag_ (iConfig.getParameter<edm::InputTag>("jetsTag")), metsTag_ (iConfig.getParameter<edm::InputTag>("metsTag")), mhtsTag_ (iConfig.getParameter<edm::InputTag>("mhtsTag")), trckTag_ (iConfig.getParameter<edm::InputTag>("trckTag")), ecalTag_ (iConfig.getParameter<edm::InputTag>("ecalTag")), min_Pt_ (iConfig.getParameter<double>("MinPt")) { LogDebug("") << "MinPt cut " << min_Pt_ << " g: " << photTag_.encode() << " e: " << elecTag_.encode() << " m: " << muonTag_.encode() << " t: " << tausTag_.encode() << " j: " << jetsTag_.encode() << " M: " << metsTag_.encode() << " H: " << mhtsTag_.encode() <<" TR: " << trckTag_.encode() <<" SC: " << ecalTag_.encode() ; }
| HLTFiltCand::~HLTFiltCand | ( | ) |
Definition at line 70 of file HLTFiltCand.cc.
{
}
| bool HLTFiltCand::hltFilter | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| trigger::TriggerFilterObjectWithRefs & | filterproduct | ||
| ) | [virtual] |
Implements HLTFilter.
Definition at line 80 of file HLTFiltCand.cc.
References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), ecalTag_, elecTag_, HI_PhotonSkim_cff::electrons, edm::Event::getByLabel(), fwrapper::jets, jetsTag_, LogDebug, metsTag_, mhtsTag_, min_Pt_, patZpeak::muons, muonTag_, interactiveExample::photons, photTag_, dt_dqm_sourceclient_common_cff::reco, HLTFilter::saveTags(), tausTag_, trckTag_, trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerJet, trigger::TriggerMET, trigger::TriggerMHT, trigger::TriggerMuon, trigger::TriggerPhoton, trigger::TriggerTau, and trigger::TriggerTrack.
{
using namespace std;
using namespace edm;
using namespace reco;
using namespace trigger;
// All HLT filters must create and fill an HLT filter object,
// recording any reconstructed physics objects satisfying (or not)
// this HLT filter, and place it in the Event.
// The filter object
if (saveTags()) {
filterproduct.addCollectionTag(photTag_);
filterproduct.addCollectionTag(elecTag_);
filterproduct.addCollectionTag(muonTag_);
filterproduct.addCollectionTag(tausTag_);
filterproduct.addCollectionTag(jetsTag_);
filterproduct.addCollectionTag(metsTag_);
filterproduct.addCollectionTag(mhtsTag_);
filterproduct.addCollectionTag(trckTag_);
filterproduct.addCollectionTag(ecalTag_);
}
// Specific filter code
// get hold of products from Event
Handle<RecoEcalCandidateCollection> photons;
Handle<ElectronCollection> electrons;
Handle<RecoChargedCandidateCollection> muons;
Handle<CaloJetCollection> taus;
Handle<CaloJetCollection> jets;
Handle<CaloMETCollection> mets;
Handle<METCollection> mhts;
Handle<RecoChargedCandidateCollection> trcks;
Handle<RecoEcalCandidateCollection> ecals;
iEvent.getByLabel(photTag_,photons );
iEvent.getByLabel(elecTag_,electrons);
iEvent.getByLabel(muonTag_,muons );
iEvent.getByLabel(tausTag_,taus );
iEvent.getByLabel(jetsTag_,jets );
iEvent.getByLabel(metsTag_,mets );
iEvent.getByLabel(mhtsTag_,mhts );
iEvent.getByLabel(trckTag_,trcks );
iEvent.getByLabel(ecalTag_,ecals );
// look for at least one g,e,m,t,j,M,H,TR,SC above its pt cut
// photons
int nphot(0);
RecoEcalCandidateCollection::const_iterator aphot(photons->begin());
RecoEcalCandidateCollection::const_iterator ophot(photons->end());
RecoEcalCandidateCollection::const_iterator iphot;
for (iphot=aphot; iphot!=ophot; iphot++) {
if (iphot->pt() >= min_Pt_) {
nphot++;
RecoEcalCandidateRef ref(RecoEcalCandidateRef(photons,distance(aphot,iphot)));
filterproduct.addObject(TriggerPhoton,ref);
}
}
// electrons
int nelec(0);
ElectronCollection::const_iterator aelec(electrons->begin());
ElectronCollection::const_iterator oelec(electrons->end());
ElectronCollection::const_iterator ielec;
for (ielec=aelec; ielec!=oelec; ielec++) {
if (ielec->pt() >= min_Pt_) {
nelec++;
ElectronRef ref(ElectronRef(electrons,distance(aelec,ielec)));
filterproduct.addObject(-TriggerElectron,ref);
}
}
// muon
int nmuon(0);
RecoChargedCandidateCollection::const_iterator amuon(muons->begin());
RecoChargedCandidateCollection::const_iterator omuon(muons->end());
RecoChargedCandidateCollection::const_iterator imuon;
for (imuon=amuon; imuon!=omuon; imuon++) {
if (imuon->pt() >= min_Pt_) {
nmuon++;
RecoChargedCandidateRef ref(RecoChargedCandidateRef(muons,distance(amuon,imuon)));
filterproduct.addObject(TriggerMuon,ref);
}
}
// taus (are stored as jets)
int ntaus(0);
CaloJetCollection::const_iterator ataus(taus->begin());
CaloJetCollection::const_iterator otaus(taus->end());
CaloJetCollection::const_iterator itaus;
for (itaus=ataus; itaus!=otaus; itaus++) {
if (itaus->pt() >= min_Pt_) {
ntaus++;
CaloJetRef ref(CaloJetRef(taus,distance(ataus,itaus)));
filterproduct.addObject(-TriggerTau,ref);
}
}
// jets
int njets(0);
CaloJetCollection::const_iterator ajets(jets->begin());
CaloJetCollection::const_iterator ojets(jets->end());
CaloJetCollection::const_iterator ijets;
for (ijets=ajets; ijets!=ojets; ijets++) {
if (ijets->pt() >= min_Pt_) {
njets++;
CaloJetRef ref(CaloJetRef(jets,distance(ajets,ijets)));
filterproduct.addObject(TriggerJet,ref);
}
}
// mets
int nmets(0);
CaloMETCollection::const_iterator amets(mets->begin());
CaloMETCollection::const_iterator omets(mets->end());
CaloMETCollection::const_iterator imets;
for (imets=amets; imets!=omets; imets++) {
if (imets->pt() >= min_Pt_) {
nmets++;
CaloMETRef ref(CaloMETRef(mets,distance(amets,imets)));
filterproduct.addObject(TriggerMET,ref);
}
}
// mhts
int nmhts(0);
METCollection::const_iterator amhts(mhts->begin());
METCollection::const_iterator omhts(mhts->end());
METCollection::const_iterator imhts;
for (imhts=amhts; imhts!=omhts; imhts++) {
if (imhts->pt() >= min_Pt_) {
nmhts++;
METRef ref(METRef(mhts,distance(amhts,imhts)));
filterproduct.addObject(TriggerMHT,ref);
}
}
// trcks
int ntrck(0);
RecoChargedCandidateCollection::const_iterator atrcks(trcks->begin());
RecoChargedCandidateCollection::const_iterator otrcks(trcks->end());
RecoChargedCandidateCollection::const_iterator itrcks;
for (itrcks=atrcks; itrcks!=otrcks; itrcks++) {
if (itrcks->pt() >= min_Pt_) {
ntrck++;
RecoChargedCandidateRef ref(RecoChargedCandidateRef(trcks,distance(atrcks,itrcks)));
filterproduct.addObject(TriggerTrack,ref);
}
}
// ecals
int necal(0);
RecoEcalCandidateCollection::const_iterator aecals(ecals->begin());
RecoEcalCandidateCollection::const_iterator oecals(ecals->end());
RecoEcalCandidateCollection::const_iterator iecals;
for (iecals=aecals; iecals!=oecals; iecals++) {
if (iecals->pt() >= min_Pt_) {
necal++;
RecoEcalCandidateRef ref(RecoEcalCandidateRef(ecals,distance(aecals,iecals)));
filterproduct.addObject(TriggerCluster,ref);
}
}
// error case
// filterproduct.addObject(0,Ref<vector<int> >());
// final filter decision:
const bool accept ( (nphot>0) && (nelec>0) && (nmuon>0) && (ntaus>0) &&
// (njets>0) && (nmets>0) && (nmhts>=0) && (ntrck>0) && (necal>0) );
(njets>0) && (nmets>0) && (ntrck>0) && (necal>0) );
LogDebug("") << "Number of g/e/m/t/j/M/H/TR/SC objects accepted:"
<< " " << nphot
<< " " << nelec
<< " " << nmuon
<< " " << ntaus
<< " " << njets
<< " " << nmets
<< " " << nmhts
<< " " << ntrck
<< " " << necal
;
// return with final filter decision
return accept;
}
edm::InputTag HLTFiltCand::ecalTag_ [private] |
Definition at line 40 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::elecTag_ [private] |
Definition at line 33 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::jetsTag_ [private] |
Definition at line 36 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::metsTag_ [private] |
Definition at line 37 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::mhtsTag_ [private] |
Definition at line 38 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
double HLTFiltCand::min_Pt_ [private] |
Definition at line 42 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::muonTag_ [private] |
Definition at line 34 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::photTag_ [private] |
Definition at line 32 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::tausTag_ [private] |
Definition at line 35 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
edm::InputTag HLTFiltCand::trckTag_ [private] |
Definition at line 39 of file HLTFiltCand.h.
Referenced by HLTFiltCand(), and hltFilter().
1.7.3