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

Razor run2 hlt 7e33 #7349

Merged
merged 3 commits into from Jan 26, 2015
Merged

Razor run2 hlt 7e33 #7349

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
8629055
DQM update for 74X for 7e33 triggers
Jan 23, 2015
cf8fcb4
DQM update for 74X for 7e33 triggers
Jan 23, 2015
8c0df56
DQM update for 74X for 7e33 triggers
Jan 23, 2015
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
3 changes: 3 additions & 0 deletions HLTriggerOffline/SUSYBSM/interface/SUSY_HLT_Razor.h
View file
Open in desktop
Expand Up @@ -36,6 +36,7 @@ class SUSY_HLT_Razor: public DQMEDAnalyzer{
static double CalcMR(TLorentzVector ja,TLorentzVector jb);
static double CalcR(double MR, TLorentzVector ja,TLorentzVector jb, edm::Handle<edm::View<reco::MET> > met, const std::vector<math::XYZTLorentzVector>& muons);
virtual ~SUSY_HLT_Razor();
static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

protected:
void dqmBeginRun(edm::Run const &, edm::EventSetup const &) override;
Expand All @@ -57,6 +58,8 @@ class SUSY_HLT_Razor: public DQMEDAnalyzer{
edm::EDGetTokenT<std::vector<math::XYZTLorentzVector> > theHemispheres_;

std::string triggerPath_;
std::string denomPath_; // denominator path for 1.4e34 lumi
std::string denomPathLoose_; // denominator path for 7e33 lumi
edm::InputTag triggerFilter_;
edm::InputTag caloFilter_;

Expand Down
16 changes: 8 additions & 8 deletions HLTriggerOffline/SUSYBSM/python/SUSYBSM_Razor_cff.py
View file
Open in desktop
Expand Up @@ -8,7 +8,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR300_Rsq0p09_MR200_v'),
TriggerFilter = cms.InputTag('hltRsqMR300Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR200Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand All @@ -19,7 +19,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR300_Rsq0p09_MR200_4jet_v'),
TriggerFilter = cms.InputTag('hltRsqMR300Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR200Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand All @@ -41,7 +41,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR260_Rsq0p09_MR200_v'),
TriggerFilter = cms.InputTag('hltRsqMR260Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR160Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand All @@ -52,7 +52,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR260_Rsq0p09_MR200_4jet_v'),
TriggerFilter = cms.InputTag('hltRsqMR260Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR160Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand All @@ -65,7 +65,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR300_Rsq0p09_MR200_v'),
TriggerFilter = cms.InputTag('hltRsqMR300Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR200Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand All @@ -76,7 +76,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR300_Rsq0p09_MR200_4jet_v'),
TriggerFilter = cms.InputTag('hltRsqMR300Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR200Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand All @@ -98,7 +98,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR260_Rsq0p09_MR200_v'),
TriggerFilter = cms.InputTag('hltRsqMR260Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR160Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand All @@ -109,7 +109,7 @@
TriggerResults = cms.InputTag('TriggerResults','','HLT'),
TriggerPath = cms.string('HLT_RsqMR260_Rsq0p09_MR200_4jet_v'),
TriggerFilter = cms.InputTag('hltRsqMR260Rsq0p09MR200', '', 'HLT'), #the last filter in the path
CaloFilter = cms.InputTag('hltRsqMR160Rsq0p01MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
CaloFilter = cms.InputTag('hltRsqMR240Rsq0p0196MR100Calo', '', 'HLT'), #filter implementing cuts on calo MR and Rsq
hemispheres = cms.InputTag('hemispheres')
)

Expand Down
4 changes: 4 additions & 0 deletions HLTriggerOffline/SUSYBSM/python/SusyExoValidation_cff.py
View file
Open in desktop
Expand Up @@ -50,6 +50,8 @@
SUSY_HLT_Razor_Main +
SUSY_HLT_Razor_QuadJet +
SUSY_HLT_Razor_DM +
SUSY_HLT_Razor_Main_7e33 +
SUSY_HLT_Razor_QuadJet_7e33 +
SUSY_HLT_CaloHT200 +
SUSY_HLT_CaloHT250 +
SUSY_HLT_CaloHT300 +
Expand Down Expand Up @@ -91,6 +93,8 @@
SUSY_HLT_Razor_Main_FASTSIM +
SUSY_HLT_Razor_QuadJet_FASTSIM +
SUSY_HLT_Razor_DM_FASTSIM +
SUSY_HLT_Razor_Main_7e33_FASTSIM +
SUSY_HLT_Razor_QuadJet_7e33_FASTSIM +
SUSY_HLT_CaloHT200_FASTSIM +
SUSY_HLT_CaloHT250_FASTSIM +
SUSY_HLT_CaloHT300_FASTSIM +
Expand Down
5 changes: 1 addition & 4 deletions HLTriggerOffline/SUSYBSM/python/razorHemispheres_cff.py
View file
Open in desktop
@@ -1,8 +1,5 @@
import FWCore.ParameterSet.Config as cms

from RecoJets.JetProducers.ak4PFJets_cfi import ak4PFJets
from RecoJets.Configuration.RecoPFJets_cff import ak4PFJetsCHS

hemispheres = cms.EDFilter(
"HLTRHemisphere",
inputTag = cms.InputTag("ak4PFJetsCHS"),
Expand All @@ -11,4 +8,4 @@
maxNJ = cms.int32(9)
)

hemisphereSequence = cms.Sequence(ak4PFJets*ak4PFJetsCHS*hemispheres)
hemisphereSequence = cms.Sequence(hemispheres)
31 changes: 28 additions & 3 deletions HLTriggerOffline/SUSYBSM/src/SUSY_HLT_Razor.cc
View file
Open in desktop
Expand Up @@ -15,6 +15,8 @@ SUSY_HLT_Razor::SUSY_HLT_Razor(const edm::ParameterSet& ps)
theHemispheres_ = consumes<std::vector<math::XYZTLorentzVector> >(ps.getParameter<edm::InputTag>("hemispheres"));
triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
triggerPath_ = ps.getParameter<std::string>("TriggerPath");
denomPath_ = ps.getParameter<std::string>("DenomPath");
denomPathLoose_ = ps.getParameter<std::string>("DenomPathLoose");
triggerFilter_ = ps.getParameter<edm::InputTag>("TriggerFilter");
caloFilter_ = ps.getParameter<edm::InputTag>("CaloFilter");
thePfJetCollection_ = consumes<reco::PFJetCollection>(ps.getParameter<edm::InputTag>("pfJetCollection"));
Expand Down Expand Up @@ -43,6 +45,21 @@ void SUSY_HLT_Razor::beginLuminosityBlock(edm::LuminosityBlock const& lumiSeg,
edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::beginLuminosityBlock" << std::endl;
}

void SUSY_HLT_Razor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
desc.add<edm::InputTag>("pfJetCollection",edm::InputTag("ak4PFJetsCHS"));
desc.add<edm::InputTag>("CaloFilter",edm::InputTag("hltRsqMR200Rsq0p01MR100Calo","","HLT"))->setComment("Calo HLT filter module used to save razor variable objects");
desc.add<edm::InputTag>("pfMETCollection",edm::InputTag("pfMet"));
desc.add<edm::InputTag>("hemispheres",edm::InputTag("hemispheres"))->setComment("hemisphere jets used to compute razor variables");
desc.add<std::string>("TriggerPath","HLT_RsqMR300_Rsq0p09_MR200_v")->setComment("trigger path name");
desc.add<std::string>("DenomPath","HLT_Ele40_eta2p1_WP85_Gsf_v")->setComment("tight denominator trigger path name (for 1.4e34 lumi)");
desc.add<std::string>("DenomPathLoose","HLT_Ele35_eta2p1_WP85_Gsf_v")->setComment("loose denominator trigger path name (for 7e33 lumi)");
desc.add<edm::InputTag>("TriggerFilter",edm::InputTag("hltRsqMR300Rsq0p09MR200","","HLT"))->setComment("PF HLT filter module used to save razor variable objects");
desc.add<edm::InputTag>("TriggerResults",edm::InputTag("TriggerResults","","HLT"));
desc.add<edm::InputTag>("trigSummary",edm::InputTag("hltTriggerSummaryAOD"));
descriptions.add("SUSY_HLT_Razor_Main",desc);
}

void SUSY_HLT_Razor::analyze(edm::Event const& e, edm::EventSetup const& eSetup){

edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::analyze" << std::endl;
Expand Down Expand Up @@ -117,13 +134,21 @@ void SUSY_HLT_Razor::analyze(edm::Event const& e, edm::EventSetup const& eSetup)
}

bool hasFired = false;
std::string denomPath = "HLT_Ele27_eta2p1_WP85_Gsf_v1"; //trigger path used as a reference for computing turn-ons and efficiencies
bool denomFired = false;
const edm::TriggerNames& trigNames = e.triggerNames(*hltresults);
unsigned int numTriggers = trigNames.size();
for( unsigned int hltIndex=0; hltIndex<numTriggers; ++hltIndex ){
if (trigNames.triggerName(hltIndex).find(triggerPath_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) hasFired = true;
if (trigNames.triggerName(hltIndex).find(denomPath) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) denomFired = true;
if (trigNames.triggerName(hltIndex).find(triggerPath_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex))
{
hasFired = true;
}
if (trigNames.triggerName(hltIndex).find(denomPath_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) {
denomFired = true;
}
if (trigNames.triggerName(hltIndex).find(denomPathLoose_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex))
{
denomFired = true;
}
}

float pfHT = 0.0;
Expand Down