Merge pull request #5942 from prbbing/MinDxyCut

Backport from CMSSW_7_3_X, adding dxy cut in the HLTMuonL3PreFilter
cms-sw · Oct 24, 2014 · c14a1e1 · c14a1e1
2 parents 637aca9 + 9a77e27
commit c14a1e1
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/HLTrigger/Muon/interface/HLTMuonL3PreFilter.h b/HLTrigger/Muon/interface/HLTMuonL3PreFilter.h
@@ -50,6 +50,7 @@ class HLTMuonL3PreFilter : public HLTFilter {
       double nsigma_Pt_;        // pt uncertainty margin (in number of sigmas)
       double max_NormalizedChi2_; // cutoff in normalized chi2
       double max_DXYBeamSpot_; // cutoff in dxy from the beamspot
+      double min_DXYBeamSpot_; // minimum cut on dxy from the beamspot
   int min_NmuonHits_; // cutoff in minumum number of chi2 hits
   double max_PtDifference_; // cutoff in maximum different between global track and tracker track
   double min_TrackPt_; //cutoff in tracker track pt

diff --git a/HLTrigger/Muon/src/HLTMuonL3PreFilter.cc b/HLTrigger/Muon/src/HLTMuonL3PreFilter.cc
@@ -53,6 +53,7 @@ HLTMuonL3PreFilter::HLTMuonL3PreFilter(const ParameterSet& iConfig) : HLTFilter(
    nsigma_Pt_  (iConfig.getParameter<double> ("NSigmaPt")),
    max_NormalizedChi2_ (iConfig.getParameter<double> ("MaxNormalizedChi2")),
    max_DXYBeamSpot_ (iConfig.getParameter<double> ("MaxDXYBeamSpot")),
+   min_DXYBeamSpot_ (iConfig.getParameter<double> ("MinDXYBeamSpot")),
    min_NmuonHits_ (iConfig.getParameter<int> ("MinNmuonHits")),
    max_PtDifference_ (iConfig.getParameter<double> ("MaxPtDifference")),
    min_TrackPt_ (iConfig.getParameter<double> ("MinTrackPt")),
@@ -95,6 +96,7 @@ HLTMuonL3PreFilter::fillDescriptions(edm::ConfigurationDescriptions& description
   desc.add<double>("NSigmaPt",0.0);
   desc.add<double>("MaxNormalizedChi2",9999.0);
   desc.add<double>("MaxDXYBeamSpot",9999.0);
+  desc.add<double>("MinDXYBeamSpot",-1.0);
   desc.add<int>("MinNmuonHits",0);
   desc.add<double>("MaxPtDifference",9999.0);
   desc.add<double>("MinTrackPt",0.0);
@@ -183,7 +185,7 @@ HLTMuonL3PreFilter::hltFilter(Event& iEvent, const EventSetup& iSetup, trigger::
        if (tk->normalizedChi2() > max_NormalizedChi2_ ) continue;
 
        //dxy beamspot cut
-       if (fabs(tk->dxy(beamSpot.position())) > max_DXYBeamSpot_ ) continue;
+       if (std::abs(tk->dxy(beamSpot.position())) > max_DXYBeamSpot_ || std::abs(tk->dxy(beamSpot.position())) < min_DXYBeamSpot_ ) continue;
 
        //min muon hits cut
        reco::HitPattern trackHits = tk->hitPattern();