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

BTV DQM Updates #37832

Merged
merged 2 commits into from
May 21, 2022
Merged

BTV DQM Updates #37832

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
d4aef3d
fixed rebase
terrill37 May 15, 2022
31cea07
code-checks and -format
terrill37 May 15, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
242 changes: 240 additions & 2 deletions DQMOffline/Trigger/plugins/BTVHLTOfflineSource.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -172,6 +172,37 @@ class BTVHLTOfflineSource : public DQMEDAnalyzer {
MonitorElement* h_3d_ip_distance = nullptr;
MonitorElement* h_3d_ip_error = nullptr;
MonitorElement* h_3d_ip_sig = nullptr;

//NEW
MonitorElement* h_jetNSecondaryVertices = nullptr;
MonitorElement* h_jet_pt = nullptr;
MonitorElement* h_jet_eta = nullptr;
MonitorElement* h_trackSumJetEtRatio = nullptr;
MonitorElement* h_trackSip2dValAboveCharm = nullptr;
MonitorElement* h_trackSip2dSigAboveCharm = nullptr;
MonitorElement* h_trackSip3dValAboveCharm = nullptr;
MonitorElement* h_trackSip3dSigAboveCharm = nullptr;
MonitorElement* h_jetNSelectedTracks = nullptr;
MonitorElement* h_jetNTracksEtaRel = nullptr;
MonitorElement* h_vertexCategory = nullptr;
MonitorElement* h_trackSumJetDeltaR = nullptr;

MonitorElement* h_trackJetDistVal = nullptr;
MonitorElement* h_trackPtRel = nullptr;
MonitorElement* h_trackDeltaR = nullptr;
MonitorElement* h_trackPtRatio = nullptr;
MonitorElement* h_trackSip3dSig = nullptr;
MonitorElement* h_trackSip2dSig = nullptr;
MonitorElement* h_trackDecayLenVal = nullptr;
MonitorElement* h_trackEtaRel = nullptr;

MonitorElement* h_vertexEnergyRatio = nullptr;
MonitorElement* h_vertexJetDeltaR = nullptr;
MonitorElement* h_flightDistance2dVal = nullptr;
MonitorElement* h_flightDistance2dSig = nullptr;
MonitorElement* h_flightDistance3dVal = nullptr;
MonitorElement* h_flightDistance3dSig = nullptr;

ObjME OnlineTrkEff_Pt;
ObjME OnlineTrkEff_Eta;
ObjME OnlineTrkEff_3d_ip_distance;
Expand Down Expand Up @@ -484,7 +515,7 @@ void BTVHLTOfflineSource::analyze(const edm::Event& iEvent, const edm::EventSetu
v.PVz_HLTMinusRECO->Fill(VertexHandler->begin()->z() - offlineVertexHandler->begin()->z());
}
}
} // caloTagsValid
} // caloTagsValid or PFTagsValid

// specific to Calo b-tagging
if (caloTags.isValid() && v.getTriggerType() == "Calo" && !caloTags->empty()) {
Expand Down Expand Up @@ -521,21 +552,99 @@ void BTVHLTOfflineSource::analyze(const edm::Event& iEvent, const edm::EventSetu

// n secondary vertices and n selected tracks
for (const auto& tagVar : tagVars.getList(reco::btau::jetNSecondaryVertices, false)) {
v.h_jetNSecondaryVertices->Fill(tagVar);
v.n_vtx->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::jetNSelectedTracks, false)) {
v.n_sel_tracks->Fill(tagVar);
v.h_jetNSelectedTracks->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::jetPt, false)) {
v.h_jet_pt->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::jetEta, false)) {
v.h_jet_eta->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::trackSumJetEtRatio, false)) {
v.h_trackSumJetEtRatio->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::trackSumJetDeltaR, false)) {
v.h_trackSumJetDeltaR->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::vertexCategory, false)) {
v.h_vertexCategory->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::trackSip2dValAboveCharm, false)) {
v.h_trackSip2dValAboveCharm->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::trackSip2dSigAboveCharm, false)) {
v.h_trackSip2dSigAboveCharm->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::trackSip3dValAboveCharm, false)) {
v.h_trackSip3dValAboveCharm->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::trackSip3dSigAboveCharm, false)) {
v.h_trackSip3dSigAboveCharm->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::jetNTracksEtaRel, false)) {
v.h_jetNTracksEtaRel->Fill(tagVar);
}

// impact parameter
// and new info
const auto& trackSip3dVal = tagVars.getList(reco::btau::trackSip3dVal, false);
const auto& trackSip3dSig = tagVars.getList(reco::btau::trackSip3dSig, false);
const auto& trackJetDistVal = tagVars.getList(reco::btau::trackJetDistVal, false);
const auto& trackPtRel = tagVars.getList(reco::btau::trackPtRel, false);
const auto& trackSip2dSig = tagVars.getList(reco::btau::trackSip2dSig, false);
const auto& trackDeltaR = tagVars.getList(reco::btau::trackDeltaR, false);
const auto& trackPtRatio = tagVars.getList(reco::btau::trackPtRatio, false);
const auto& trackDecayLenVal = tagVars.getList(reco::btau::trackDecayLenVal, false);
const auto& trackEtaRel = tagVars.getList(reco::btau::trackEtaRel, false);

for (unsigned i_trk = 0; i_trk < trackEtaRel.size(); i_trk++) {
v.h_trackEtaRel->Fill(trackEtaRel[i_trk]);
}

for (unsigned i_trk = 0; i_trk < trackJetDistVal.size(); i_trk++) {
v.h_trackJetDistVal->Fill(trackJetDistVal[i_trk]);
}

for (unsigned i_trk = 0; i_trk < trackPtRel.size(); i_trk++) {
v.h_trackPtRel->Fill(trackPtRel[i_trk]);
}

for (unsigned i_trk = 0; i_trk < trackDeltaR.size(); i_trk++) {
v.h_trackDeltaR->Fill(trackDeltaR[i_trk]);
}

for (unsigned i_trk = 0; i_trk < trackPtRatio.size(); i_trk++) {
v.h_trackPtRatio->Fill(trackPtRatio[i_trk]);
}

for (unsigned i_trk = 0; i_trk < trackDecayLenVal.size(); i_trk++) {
v.h_trackDecayLenVal->Fill(trackDecayLenVal[i_trk]);
}

for (unsigned i_trk = 0; i_trk < trackSip3dVal.size(); i_trk++) {
float val = trackSip3dVal[i_trk];
float sig = trackSip3dSig[i_trk];
v.h_3d_ip_distance->Fill(val);
v.h_3d_ip_error->Fill(val / sig);
v.h_3d_ip_sig->Fill(sig);

v.h_trackSip2dSig->Fill(trackSip2dSig[i_trk]);
}

// vertex mass and tracks per vertex
Expand All @@ -546,6 +655,30 @@ void BTVHLTOfflineSource::analyze(const edm::Event& iEvent, const edm::EventSetu
v.n_vtx_trks->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::vertexEnergyRatio, false)) {
v.h_vertexEnergyRatio->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::vertexJetDeltaR, false)) {
v.h_vertexJetDeltaR->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::flightDistance2dVal, false)) {
v.h_flightDistance2dVal->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::flightDistance2dSig, false)) {
v.h_flightDistance2dSig->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::flightDistance3dVal, false)) {
v.h_flightDistance3dVal->Fill(tagVar);
}

for (const auto& tagVar : tagVars.getList(reco::btau::flightDistance3dSig, false)) {
v.h_flightDistance3dSig->Fill(tagVar);
}

// // track N total/pixel hits
// for (const auto & tagVar : tagVars.getList(reco::btau::trackNPixelHits, false)) {
// v.n_pixel_hits->Fill(tagVar);}
Expand Down Expand Up @@ -589,7 +722,7 @@ void BTVHLTOfflineSource::analyze(const edm::Event& iEvent, const edm::EventSetu
// }
// }
// }
}
} //end paths loop
}

std::vector<const reco::Track*> BTVHLTOfflineSource::getOfflineBTagTracks(float hltJetEta,
Expand Down Expand Up @@ -759,6 +892,111 @@ void BTVHLTOfflineSource::bookHistograms(DQMStore::IBooker& iBooker, edm::Run co
title = "3D IP significance of tracks (cm)" + trigPath;
v.h_3d_ip_sig = iBooker.book1D(histoname.c_str(), title.c_str(), 40, -40, 40);

//new
histoname = "jetNSecondaryVertices";
title = "jet N Secondary Vertices" + trigPath;
v.h_jetNSecondaryVertices = iBooker.book1D(histoname.c_str(), title.c_str(), 10, -0.5, 9.5);

histoname = "jet_pt";
title = "jet pt" + trigPath;
v.h_jet_pt = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 100);

histoname = "jet_eta";
title = "jet eta" + trigPath;
v.h_jet_eta = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -2.5, 2.5);

histoname = "trackSumJetEtRatio";
title = "trackSumJetEtRatio" + trigPath;
v.h_trackSumJetEtRatio = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 1.5);

histoname = "trackSip2dValAboveCharm";
title = "trackSip2dSigAboveCharm" + trigPath;
v.h_trackSip2dSigAboveCharm = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.2, 0.2);

histoname = "trackSip2dSigAboveCharm";
title = "trackSip2dSigAboveCharm" + trigPath;
v.h_trackSip2dValAboveCharm = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -50, 50);

histoname = "trackSip3dValAboveCharm";
title = "trackSip3dValAboveCharm" + trigPath;
v.h_trackSip3dValAboveCharm = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.2, 0.2);

histoname = "trackSip3dSigAboveCharm";
title = "trackSip3dSigAboveCharm" + trigPath;
v.h_trackSip3dSigAboveCharm = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -50, 50);

histoname = "jetNSelectedTracks";
title = "jet N Selected Tracks" + trigPath;
v.h_jetNSelectedTracks = iBooker.book1D(histoname.c_str(), title.c_str(), 42, -1.5, 40.5);

histoname = "jetNTracksEtaRel";
title = "jetNTracksEtaRel" + trigPath;
v.h_jetNTracksEtaRel = iBooker.book1D(histoname.c_str(), title.c_str(), 42, -1.5, 40.5);

histoname = "vertexCategory";
title = "vertex category" + trigPath;
v.h_vertexCategory = iBooker.book1D(histoname.c_str(), title.c_str(), 4, -1.5, 2.5);

histoname = "trackSumJetDeltaR";
title = "trackSumJetDeltaR" + trigPath;
v.h_trackSumJetDeltaR = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 0.35);

//new 2 below
histoname = "trackJetDistVal";
title = "trackJetDistVal" + trigPath;
v.h_trackJetDistVal = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -1, 0.01);

histoname = "trackPtRel";
title = "track pt rel" + trigPath;
v.h_trackPtRel = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 7);

histoname = "trackDeltaR";
title = "trackDeltaR" + trigPath;
v.h_trackDeltaR = iBooker.book1D(histoname.c_str(), title.c_str(), 160, -0.05, .47);

histoname = "trackPtRatio";
title = "trackPtRatio" + trigPath;
v.h_trackPtRatio = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.01, 0.3);

histoname = "trackSip2dSig";
title = "trackSip2dSig" + trigPath;
v.h_trackSip2dSig = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -55, 55);

histoname = "trackDecayLenVal";
title = "trackDecayLenVal" + trigPath;
v.h_trackDecayLenVal = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 22);

histoname = "trackEtaRel";
title = "trackEtaRel" + trigPath;
v.h_trackEtaRel = iBooker.book1D(histoname.c_str(), title.c_str(), 31, 0, 30);

//new 3 below
histoname = "vertexEnergyRatio";
title = "vertexEnergyRatio" + trigPath;
v.h_vertexEnergyRatio = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 3);

histoname = "vertexJetDeltaR";
title = "vertexJetDeltaR" + trigPath;
v.h_vertexJetDeltaR = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.01, 0.4);

histoname = "flightDistance2dVal";
title = "flightDistance2dVal" + trigPath;
v.h_flightDistance2dVal = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 5);

histoname = "flightDistance2dSig";
title = "flightDistance2dSig" + trigPath;
v.h_flightDistance2dSig = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -10, 150);

histoname = "flightDistance3dVal";
title = "flightDistance3dVal" + trigPath;
v.h_flightDistance3dVal = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -0.1, 5);

histoname = "flightDistance3dSig";
title = "flightDistance3dSig" + trigPath;
v.h_flightDistance3dSig = iBooker.book1D(histoname.c_str(), title.c_str(), 100, -10, 150);

//end new

histoname = "OnlineTrkEff_Pt";
title = "Relative Online Trk Efficiency vs Pt " + trigPath;
v.bookME(iBooker, v.OnlineTrkEff_Pt, histoname, title, 50, -0.5, 20.);
Expand Down
Loading