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

80X backport: Added Z mass histograms in Online DQM to monitor L1/L2/L3 dimuon triggers for pPb 2016 #15913

Merged
merged 1 commit into from Oct 4, 2016
Merged

80X backport: Added Z mass histograms in Online DQM to monitor L1/L2/L3 dimuon triggers for pPb 2016 #15913

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
90 changes: 90 additions & 0 deletions DQM/HLTEvF/plugins/HLTObjectMonitor.cc
View file
Open in desktop
Expand Up @@ -162,6 +162,9 @@ class HLTObjectMonitor : public DQMEDAnalyzer {
edm::ParameterSet bJetCSVCalo_pset;
edm::ParameterSet bJetCSVPF_pset;
edm::ParameterSet diMuonMass_pset;
edm::ParameterSet pAL1DoubleMuZMass_pset;
edm::ParameterSet pAL2DoubleMuZMass_pset;
edm::ParameterSet pAL3DoubleMuZMass_pset;
edm::ParameterSet diElecMass_pset;
edm::ParameterSet muonDxy_pset;
edm::ParameterSet wallTime_pset;
Expand Down Expand Up @@ -202,6 +205,9 @@ class HLTObjectMonitor : public DQMEDAnalyzer {
hltPlot bJetCSVCalo_;
hltPlot bJetCSVPF_;
hltPlot diMuonMass_;
hltPlot pAL1DoubleMuZMass_;
hltPlot pAL2DoubleMuZMass_;
hltPlot pAL3DoubleMuZMass_;
hltPlot diElecMass_;
hltPlot muonDxy_;
hltPlot wallTime_;
Expand Down Expand Up @@ -296,6 +302,12 @@ HLTObjectMonitor::HLTObjectMonitor(const edm::ParameterSet& iConfig)
plotMap[&bJetCSVPF_] = &bJetCSVPF_pset;
diMuonMass_pset = iConfig.getParameter<edm::ParameterSet>("diMuonMass");
plotMap[&diMuonMass_] = &diMuonMass_pset;
pAL1DoubleMuZMass_pset = iConfig.getParameter<edm::ParameterSet>("pAL1DoubleMuZMass");
plotMap[&pAL1DoubleMuZMass_] = &pAL1DoubleMuZMass_pset;
pAL2DoubleMuZMass_pset = iConfig.getParameter<edm::ParameterSet>("pAL2DoubleMuZMass");
plotMap[&pAL2DoubleMuZMass_] = &pAL2DoubleMuZMass_pset;
pAL3DoubleMuZMass_pset = iConfig.getParameter<edm::ParameterSet>("pAL3DoubleMuZMass");
plotMap[&pAL3DoubleMuZMass_] = &pAL3DoubleMuZMass_pset;
diElecMass_pset = iConfig.getParameter<edm::ParameterSet>("diElecMass");
plotMap[&diElecMass_] = &diElecMass_pset;
muonDxy_pset = iConfig.getParameter<edm::ParameterSet>("muonDxy");
Expand Down Expand Up @@ -669,6 +681,84 @@ HLTObjectMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetu
}
}

else if (pathName == pAL1DoubleMuZMass_.pathName)
{
const double mu_mass(.105658);
unsigned int kCnt0 = 0;
for (const auto & key0: keys)
{
unsigned int kCnt1 = 0;
for (const auto & key1: keys)
{
if (key0 != key1 && kCnt1 > kCnt0) // avoid filling hists with same objs && avoid double counting separate objs
{
// if (abs(objects[key0].id()) == 13 && (objects[key0].id()+objects[key1].id()==0)) // id is not filled for l1 stage2 muons
// {
TLorentzVector mu1, mu2, dimu;
mu1.SetPtEtaPhiM(objects[key0].pt(), objects[key0].eta(), objects[key0].phi(), mu_mass);
mu2.SetPtEtaPhiM(objects[key1].pt(), objects[key1].eta(), objects[key1].phi(), mu_mass);
dimu = mu1+mu2;
if(dimu.M()>pAL1DoubleMuZMass_.xMin && dimu.M()<pAL1DoubleMuZMass_.xMax) pAL1DoubleMuZMass_.ME->Fill(dimu.M());
// }
}
kCnt1 +=1;
}
kCnt0 +=1;
}
}

else if (pathName == pAL2DoubleMuZMass_.pathName)
{
const double mu_mass(.105658);
unsigned int kCnt0 = 0;
for (const auto & key0: keys)
{
unsigned int kCnt1 = 0;
for (const auto & key1: keys)
{
if (key0 != key1 && kCnt1 > kCnt0) // avoid filling hists with same objs && avoid double counting separate objs
{
if (abs(objects[key0].id()) == 13 && (objects[key0].id()+objects[key1].id()==0)) // check muon id and dimuon charge
{
TLorentzVector mu1, mu2, dimu;
mu1.SetPtEtaPhiM(objects[key0].pt(), objects[key0].eta(), objects[key0].phi(), mu_mass);
mu2.SetPtEtaPhiM(objects[key1].pt(), objects[key1].eta(), objects[key1].phi(), mu_mass);
dimu = mu1+mu2;
if(dimu.M()>pAL2DoubleMuZMass_.xMin && dimu.M()<pAL2DoubleMuZMass_.xMax) pAL2DoubleMuZMass_.ME->Fill(dimu.M());
}
}
kCnt1 +=1;
}
kCnt0 +=1;
}
}

else if (pathName == pAL3DoubleMuZMass_.pathName)
{
const double mu_mass(.105658);
unsigned int kCnt0 = 0;
for (const auto & key0: keys)
{
unsigned int kCnt1 = 0;
for (const auto & key1: keys)
{
if (key0 != key1 && kCnt1 > kCnt0) // avoid filling hists with same objs && avoid double counting separate objs
{
if (abs(objects[key0].id()) == 13 && (objects[key0].id()+objects[key1].id()==0)) // check muon id and dimuon charge
{
TLorentzVector mu1, mu2, dimu;
mu1.SetPtEtaPhiM(objects[key0].pt(), objects[key0].eta(), objects[key0].phi(), mu_mass);
mu2.SetPtEtaPhiM(objects[key1].pt(), objects[key1].eta(), objects[key1].phi(), mu_mass);
dimu = mu1+mu2;
if(dimu.M()>pAL3DoubleMuZMass_.xMin && dimu.M()<pAL3DoubleMuZMass_.xMax) pAL3DoubleMuZMass_.ME->Fill(dimu.M());
}
}
kCnt1 +=1;
}
kCnt0 +=1;
}
}

else if (pathName == diElecMass_.pathName)
{
unsigned int kCnt0 = 0;
Expand Down
30 changes: 30 additions & 0 deletions DQM/HLTEvF/python/HLTObjectMonitor_cfi.py
View file
Open in desktop
Expand Up @@ -384,6 +384,36 @@
Xmin = cms.double(60),
Xmax = cms.double(160)
),
pAL1DoubleMuZMass = cms.PSet(
pathName = cms.string("HLT_PAL1DoubleMu10"),
moduleName = cms.string("hltL1fL1sDoubleMu10BptxANDL1Filtered0"),
plotLabel = cms.string("PAL1DoubleMu10_ZMass"),
axisLabel = cms.string("dimuon mass [GeV]"),
mainWorkspace = cms.bool(True),
NbinsX = cms.int32(50),
Xmin = cms.double(60.0),
Xmax = cms.double(160.0)
),
pAL2DoubleMuZMass = cms.PSet(
pathName = cms.string("HLT_PAL2DoubleMu10"),
moduleName = cms.string("hltL2fL1sDoubleMu10BptxANDL1f0L2Filtered10"),
plotLabel = cms.string("PAL2DoubleMu10_ZMass"),
axisLabel = cms.string("dimuon mass [GeV]"),
mainWorkspace = cms.bool(True),
NbinsX = cms.int32(50),
Xmin = cms.double(60.0),
Xmax = cms.double(160.0)
),
pAL3DoubleMuZMass = cms.PSet(
pathName = cms.string("HLT_PAL3DoubleMu10"),
moduleName = cms.string("hltL3fL1sDoubleMu10BptxANDL1f0L2f10L3Filtered10"),
plotLabel = cms.string("PAL3DoubleMu10_ZMass"),
axisLabel = cms.string("dimuon mass [GeV]"),
mainWorkspace = cms.bool(True),
NbinsX = cms.int32(50),
Xmin = cms.double(60.0),
Xmax = cms.double(160.0)
),
diElecMass = cms.PSet(
pathName = cms.string("HLT_Ele23_Ele12_CaloIdL_TrackIdL_IsoVL_DZ"),
moduleName = cms.string("hltEle23Ele12CaloIdLTrackIdLIsoVLDZFilter"),
Expand Down