Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

HLT dilepton conditional DZ filter, 81X #15683

Merged
merged 1 commit into from Sep 5, 2016
Merged

HLT dilepton conditional DZ filter, 81X #15683

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
1 change: 1 addition & 0 deletions HLTrigger/HLTfilters/interface/HLTDoubletDZ.h
View file
Open in desktop
Expand Up @@ -53,6 +53,7 @@ class HLTDoubletDZ : public HLTFilter {
const int triggerType2_;
const double minDR_; // minimum dR between two objects to be considered a pair
const double maxDZ_; // number of pairs passing cuts required
const int minPixHitsForDZ_; // minimum number of required pixel hits to check DZ
const int min_N_; // number of pairs passing cuts required
const bool checkSC_; // make sure SC constituents are different
const bool same_; // 1st and 2nd product are one and the same
Expand Down
94 changes: 89 additions & 5 deletions HLTrigger/HLTfilters/src/HLTDoubletDZ.cc
View file
Open in desktop
Expand Up @@ -39,6 +39,7 @@ HLTDoubletDZ<T1,T2>::HLTDoubletDZ(const edm::ParameterSet& iConfig) :
triggerType2_(iConfig.template getParameter<int>("triggerType2")),
minDR_ (iConfig.template getParameter<double>("MinDR")),
maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
min_N_ (iConfig.template getParameter<int>("MinN")),
checkSC_ (iConfig.template getParameter<bool>("checkSC")),
same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
Expand All @@ -58,6 +59,7 @@ HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>::HLTDoubletDZ(con
triggerType2_(iConfig.template getParameter<int>("triggerType2")),
minDR_ (iConfig.template getParameter<double>("MinDR")),
maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
min_N_ (iConfig.template getParameter<int>("MinN")),
checkSC_ (iConfig.template getParameter<bool>("checkSC")),
same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
Expand All @@ -77,6 +79,7 @@ HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>::HLTDoubletDZ(
triggerType2_(iConfig.template getParameter<int>("triggerType2")),
minDR_ (iConfig.template getParameter<double>("MinDR")),
maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
min_N_ (iConfig.template getParameter<int>("MinN")),
checkSC_ (iConfig.template getParameter<bool>("checkSC")),
same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
Expand All @@ -96,6 +99,27 @@ HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>::HLTDoubletDZ(
triggerType2_(iConfig.template getParameter<int>("triggerType2")),
minDR_ (iConfig.template getParameter<double>("MinDR")),
maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
min_N_ (iConfig.template getParameter<int>("MinN")),
checkSC_ (iConfig.template getParameter<bool>("checkSC")),
same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
{}

template<>
HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>::HLTDoubletDZ(const edm::ParameterSet& iConfig) :
HLTFilter(iConfig),
originTag1_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag1")),
originTag2_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag2")),
inputTag1_(iConfig.template getParameter<edm::InputTag>("inputTag1")),
inputTag2_(iConfig.template getParameter<edm::InputTag>("inputTag2")),
inputToken1_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag1_)),
inputToken2_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag2_)),
//electronToken_ (consumes<reco::ElectronCollection>(iConfig.template getParameter<edm::InputTag>("electronTag"))),
triggerType1_(iConfig.template getParameter<int>("triggerType1")),
triggerType2_(iConfig.template getParameter<int>("triggerType2")),
minDR_ (iConfig.template getParameter<double>("MinDR")),
maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
min_N_ (iConfig.template getParameter<int>("MinN")),
checkSC_ (iConfig.template getParameter<bool>("checkSC")),
same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
Expand All @@ -120,6 +144,7 @@ HLTDoubletDZ<T1,T2>::fillDescriptions(edm::ConfigurationDescriptions& descriptio
desc.add<int>("triggerType2",0);
desc.add<double>("MinDR",-1.0);
desc.add<double>("MaxDZ",0.2);
desc.add<int>("MinPixHitsForDZ",0);
desc.add<bool>("checkSC",false);
desc.add<int>("MinN",1);
descriptions.add(defaultModuleLabel<HLTDoubletDZ<T1,T2>>(), desc);
Expand All @@ -141,6 +166,7 @@ HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>::fillDescriptions
desc.add<int>("triggerType2",0);
desc.add<double>("MinDR",-1.0);
desc.add<double>("MaxDZ",0.2);
desc.add<int>("MinPixHitsForDZ",0);
desc.add<bool>("checkSC",false);
desc.add<int>("MinN",1);
descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>>(), desc);
Expand All @@ -162,6 +188,7 @@ HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>::fillDescripti
desc.add<int>("triggerType2",0);
desc.add<double>("MinDR",-1.0);
desc.add<double>("MaxDZ",0.2);
desc.add<int>("MinPixHitsForDZ",0);
desc.add<bool>("checkSC",false);
desc.add<int>("MinN",1);
descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>>(), desc);
Expand All @@ -183,11 +210,33 @@ HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>::fillDescripti
desc.add<int>("triggerType2",0);
desc.add<double>("MinDR",-1.0);
desc.add<double>("MaxDZ",0.2);
desc.add<int>("MinPixHitsForDZ",0);
desc.add<bool>("checkSC",false);
desc.add<int>("MinN",1);
descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>>(), desc);
}

template<>
void
HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
makeHLTFilterDescription(desc);
std::vector<edm::InputTag> originTag1(1,edm::InputTag("hltOriginal1"));
std::vector<edm::InputTag> originTag2(1,edm::InputTag("hltOriginal2"));
desc.add<std::vector<edm::InputTag> >("originTag1",originTag1);
desc.add<std::vector<edm::InputTag> >("originTag2",originTag2);
desc.add<edm::InputTag>("inputTag1",edm::InputTag("hltFiltered1"));
desc.add<edm::InputTag>("inputTag2",edm::InputTag("hltFiltered2"));
desc.add<int>("triggerType1",0);
desc.add<int>("triggerType2",0);
desc.add<double>("MinDR",-1.0);
desc.add<double>("MaxDZ",0.2);
desc.add<int>("MinPixHitsForDZ",0);
desc.add<bool>("checkSC",false);
desc.add<int>("MinN",1);
descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>>(), desc);
}

template<typename T1, typename T2>
bool
//HLTDoubletDZ<T1, T2>::getCollections(edm::Event& iEvent, std::vector<T1Ref>& coll1, std::vector<T2Ref>& coll2) const {
Expand Down Expand Up @@ -273,8 +322,13 @@ HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>::computeDZ(edm
e1 = *(eleIt);
}

const reco::Candidate& candidate2(*r2);
if ( std::abs(e1.vz()-candidate2.vz()) > maxDZ_ )
const reco::RecoChargedCandidate& candidate2(*r2);
bool skipDZ = false;
if ( minPixHitsForDZ_ > 0 &&
(e1.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
|| candidate2.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
skipDZ = true;
if ( !skipDZ && std::abs(e1.vz()-candidate2.vz()) > maxDZ_ )
return false;

return true;
Expand All @@ -297,8 +351,13 @@ HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>::computeDZ(edm
e2 = *(eleIt);
}

const reco::Candidate& candidate1(*r1);
if ( std::abs(e2.vz()-candidate1.vz()) > maxDZ_ )
const reco::RecoChargedCandidate& candidate1(*r1);
bool skipDZ = false;
if ( minPixHitsForDZ_ > 0 &&
(candidate1.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
|| e2.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
skipDZ = true;
if ( !skipDZ && std::abs(e2.vz()-candidate1.vz()) > maxDZ_ )
return false;

return true;
Expand All @@ -323,7 +382,32 @@ HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>::computeDZ(edm::E
e1 = *(eleIt);
}

if ( std::abs(e2.vz()-e1.vz()) > maxDZ_ )
bool skipDZ = false;
if ( minPixHitsForDZ_ > 0 &&
(e1.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
|| e2.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
skipDZ = true;
if ( !skipDZ && std::abs(e2.vz()-e1.vz()) > maxDZ_ )
return false;

return true;
}

template<>
bool
HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>::computeDZ(edm::Event& iEvent, T1Ref& r1, T2Ref& r2) const
{

const reco::RecoChargedCandidate& candidate1(*r1);
const reco::RecoChargedCandidate& candidate2(*r2);
if ( reco::deltaR(candidate1, candidate2) < minDR_ )
return false;
bool skipDZ = false;
if ( minPixHitsForDZ_ > 0 &&
(candidate1.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
|| candidate2.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
skipDZ = true;
if ( !skipDZ && std::abs(candidate1.vz()-candidate2.vz()) > maxDZ_ )
return false;

return true;
Expand Down