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

[12.4.X] Add Tracker Alignment di-muon monitoring to other resonances #39332

Merged
merged 4 commits into from Sep 8, 2022
Merged

[12.4.X] Add Tracker Alignment di-muon monitoring to other resonances #39332

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
be78f7a
add DiMuon monitioring to other TkAl X->mu mu producers
mmusich Aug 30, 2022
37d0d31
add di-muon mass bias harvesting sequences for charmonium and bottomo…
mmusich Aug 30, 2022
274affb
improve DiMuonVertexMonitor histogram titles: make them configurable …
mmusich Sep 2, 2022
4fc6e85
update unit tests for different resonances
mmusich Sep 2, 2022
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
5 changes: 4 additions & 1 deletion DQMOffline/Alignment/interface/DiMuonVertexMonitor.h
View file
Open in desktop
Expand Up @@ -50,8 +50,11 @@ class DiMuonVertexMonitor : public DQMEDAnalyzer {
//used to select what vertices to read from configuration file
const edm::EDGetTokenT<reco::VertexCollection> vertexToken_;

const std::string motherName_;
const std::string MEFolderName_; // Top-level folder name
const bool useClosestVertex_;

std::pair<float, float> massLimits_; /* for the mass plot x-range */
const float maxSVdist_;

// vertex quantities
Expand All @@ -73,7 +76,7 @@ class DiMuonVertexMonitor : public DQMEDAnalyzer {
MonitorElement *hCosPhi3D_;
MonitorElement *hCosPhiInv_;
MonitorElement *hCosPhiInv3D_;
MonitorElement *hTrackInvMass_;
MonitorElement *hInvMass_;
MonitorElement *hCutFlow_;

// impact parameters information
Expand Down
35 changes: 33 additions & 2 deletions DQMOffline/Alignment/python/ALCARECOTkAlDQM_cff.py
View file
Open in desktop
Expand Up @@ -186,6 +186,7 @@
TrackPtMax = 50
)


from Alignment.CommonAlignmentProducer.ALCARECOTkAlJpsiMuMu_cff import ALCARECOTkAlJpsiMuMuHLT
#ALCARECOTkAlJpsiMuMuHLTDQM = hltMonBitSummary.clone(
# directory = "AlCaReco/"+__selectionName+"/HLTSummary",
Expand All @@ -195,7 +196,22 @@
#)

#ALCARECOTkAlJpsiMuMuDQM = cms.Sequence( ALCARECOTkAlJpsiMuMuTrackingDQM + ALCARECOTkAlJpsiMuMuTkAlDQM + ALCARECOTkAlJpsiMuMuHLTDQM)
ALCARECOTkAlJpsiMuMuDQM = cms.Sequence( ALCARECOTkAlJpsiMuMuTrackingDQM + ALCARECOTkAlJpsiMuMuTkAlDQM )

ALCARECOTkAlJpsiMuMuVtxDQM = DQMOffline.Alignment.DiMuonVertexMonitor_cfi.DiMuonVertexMonitor.clone(
muonTracks = 'ALCARECO'+__selectionName,
decayMotherName = "J/#psi",
vertices = 'offlinePrimaryVertices',
FolderName = "AlCaReco/"+__selectionName,
maxSVdist = 50
)

ALCARECOTkAlJpsiMassBiasDQM = DQMOffline.Alignment.DiMuonMassBiasMonitor_cfi.DiMuonMassBiasMonitor.clone(
muonTracks = 'ALCARECO'+__selectionName,
FolderName = "AlCaReco/"+__selectionName,
decayMotherName = 'J/#psi',
DiMuMassConfig = dict(ymin = 2.7 ,ymax = 3.4))

ALCARECOTkAlJpsiMuMuDQM = cms.Sequence( ALCARECOTkAlJpsiMuMuTrackingDQM + ALCARECOTkAlJpsiMuMuTkAlDQM + ALCARECOTkAlJpsiMuMuVtxDQM + ALCARECOTkAlJpsiMassBiasDQM)

#########################################################
#############--- TkAlJpsiMuMuHI ---#####################
Expand Down Expand Up @@ -269,7 +285,22 @@
#)

#ALCARECOTkAlUpsilonMuMuDQM = cms.Sequence( ALCARECOTkAlUpsilonMuMuTrackingDQM + ALCARECOTkAlUpsilonMuMuTkAlDQM + ALCARECOTkAlUpsilonMuMuHLTDQM)
ALCARECOTkAlUpsilonMuMuDQM = cms.Sequence( ALCARECOTkAlUpsilonMuMuTrackingDQM + ALCARECOTkAlUpsilonMuMuTkAlDQM )

ALCARECOTkAlUpsilonMuMuVtxDQM = DQMOffline.Alignment.DiMuonVertexMonitor_cfi.DiMuonVertexMonitor.clone(
muonTracks = 'ALCARECO'+__selectionName,
decayMotherName = "#Upsilon",
vertices = 'offlinePrimaryVertices',
FolderName = "AlCaReco/"+__selectionName,
maxSVdist = 50
)

ALCARECOTkAlUpsilonMassBiasDQM = DQMOffline.Alignment.DiMuonMassBiasMonitor_cfi.DiMuonMassBiasMonitor.clone(
muonTracks = 'ALCARECO'+__selectionName,
FolderName = "AlCaReco/"+__selectionName,
decayMotherName = '#Upsilon',
DiMuMassConfig = dict(ymin = 8.9 ,ymax = 9.9))

ALCARECOTkAlUpsilonMuMuDQM = cms.Sequence( ALCARECOTkAlUpsilonMuMuTrackingDQM + ALCARECOTkAlUpsilonMuMuTkAlDQM + ALCARECOTkAlUpsilonMuMuVtxDQM + ALCARECOTkAlUpsilonMassBiasDQM)

############################################################
#############--- TkAlUpsilonMuMuHI ---#####################
Expand Down
25 changes: 25 additions & 0 deletions DQMOffline/Alignment/python/DiMuonMassBiasHarvesting_cff.py
View file
Open in desktop
Expand Up @@ -2,9 +2,34 @@

from DQMOffline.Alignment.DiMuonMassBiasClient_cfi import DiMuonMassBiasClient as diMuonMassBiasClient

# Z-> mm
__selectionName = 'TkAlDiMuonAndVertex'
ALCARECOTkAlZMuMuMassBiasClient = diMuonMassBiasClient.clone(
FolderName = "AlCaReco/"+__selectionName
)

alcaTkAlZMuMuBiasClients = cms.Sequence(ALCARECOTkAlZMuMuMassBiasClient)

# J/psi -> mm
__selectionName = 'TkAlJpsiMuMu'
ALCARECOTkAlJpsiMuMuMassBiasClient = diMuonMassBiasClient.clone(
FolderName = "AlCaReco/"+__selectionName,
fitBackground = True,
fit_par = dict(mean_par = [3.09, 2.7, 3.4],
width_par = [1.0, 0.0, 5.0],
sigma_par = [1.0, 0.0, 5.0])
)

alcaTkAlJpsiMuMuBiasClients = cms.Sequence(ALCARECOTkAlJpsiMuMuMassBiasClient)

# Upsilon -> mm
__selectionName = 'TkAlUpsilonMuMu'
ALCARECOTkAlUpsilonMuMuMassBiasClient = diMuonMassBiasClient.clone(
FolderName = "AlCaReco/"+__selectionName,
fitBackground = True,
fit_par = dict(mean_par = [9.46, 8.9, 9.9],
width_par = [1.0, 0.0, 5.0],
sigma_par = [1.0, 0.0, 5.0])
)

alcaTkAlUpsilonMuMuBiasClients = cms.Sequence(ALCARECOTkAlUpsilonMuMuMassBiasClient)
1 change: 1 addition & 0 deletions DQMOffline/Alignment/python/DiMuonVertexMonitor_cfi.py
View file
Open in desktop
Expand Up @@ -3,6 +3,7 @@
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
DiMuonVertexMonitor = DQMEDAnalyzer('DiMuonVertexMonitor',
muonTracks = cms.InputTag('ALCARECOTkAlDiMuon'),
decayMotherName = cms.string('Z'),
vertices = cms.InputTag('offlinePrimaryVertices'),
FolderName = cms.string('DiMuonVertexMonitor'),
maxSVdist = cms.double(50))
103 changes: 70 additions & 33 deletions DQMOffline/Alignment/src/DiMuonVertexMonitor.cc
View file
Open in desktop
Expand Up @@ -3,6 +3,8 @@
*
*/

#include <fmt/printf.h>

#include "DQMOffline/Alignment/interface/DiMuonVertexMonitor.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/Math/interface/deltaR.h"
Expand All @@ -28,43 +30,77 @@ DiMuonVertexMonitor::DiMuonVertexMonitor(const edm::ParameterSet& iConfig)
: ttbESToken_(esConsumes(edm::ESInputTag("", "TransientTrackBuilder"))),
tracksToken_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("muonTracks"))),
vertexToken_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertices"))),
motherName_(iConfig.getParameter<std::string>("decayMotherName")),
MEFolderName_(iConfig.getParameter<std::string>("FolderName")),
useClosestVertex_(iConfig.getParameter<bool>("useClosestVertex")),
maxSVdist_(iConfig.getParameter<double>("maxSVdist")) {}
maxSVdist_(iConfig.getParameter<double>("maxSVdist")) {
if (motherName_.find('Z') != std::string::npos) {
massLimits_ = std::make_pair(50., 120);
} else if (motherName_.find("J/#psi") != std::string::npos) {
massLimits_ = std::make_pair(2.7, 3.4);
} else if (motherName_.find("#Upsilon") != std::string::npos) {
massLimits_ = std::make_pair(8.9, 9.9);
} else {
edm::LogError("DiMuonVertexMonitor") << " unrecognized decay mother particle: " << motherName_
<< " setting the default for the Z->mm (50.,120.)" << std::endl;
massLimits_ = std::make_pair(50., 120);
}
}

void DiMuonVertexMonitor::bookHistograms(DQMStore::IBooker& iBooker, edm::Run const&, edm::EventSetup const&) {
iBooker.setCurrentFolder(MEFolderName_ + "/DiMuonVertexMonitor");
hSVProb_ = iBooker.book1D("VtxProb", ";ZV vertex probability;N(#mu#mu pairs)", 100, 0., 1.);
hSVChi2_ =
iBooker.book1D("VtxChi2", "#chi^{2} of the Z vertex; #chi^{2} of the Z vertex; N(#mu#mu pairs)", 200, 0., 200.);
hSVNormChi2_ = iBooker.book1D(
"VtxNormChi2", "#chi^{2}/ndf of the Z vertex; #chi^{2}/ndf of Z vertex; N(#mu#mu pairs)", 100, 0., 20.);
hSVDist_ = iBooker.book1D("VtxDist", ";PV-ZV xy distance [#mum];N(#mu#mu pairs)", 100, 0., 300.);
hSVDistErr_ = iBooker.book1D("VtxDistErr", ";PV-ZV xy distance error [#mum];N(#mu#mu pairs)", 100, 0., 1000.);
hSVDistSig_ = iBooker.book1D("VtxDistSig", ";PV-ZV xy distance signficance;N(#mu#mu pairs)", 100, 0., 5.);
hSVCompatibility_ = iBooker.book1D(
"VtxCompatibility", "compatibility of Z vertex; compatibility of Z vertex; N(#mu#mu pairs)", 100, 0., 100.);

hSVDist3D_ = iBooker.book1D("VtxDist3D", ";PV-ZV 3D distance [#mum];N(#mu#mu pairs)", 100, 0., 300.);
hSVDist3DErr_ = iBooker.book1D("VtxDist3DErr", ";PV-ZV 3D distance error [#mum];N(#mu#mu pairs)", 100, 0., 1000.);
hSVDist3DSig_ = iBooker.book1D("VtxDist3DSig", ";PV-ZV 3D distance signficance;N(#mu#mu pairs)", 100, 0., 5.);
hSVCompatibility3D_ = iBooker.book1D(
"VtxCompatibility3D", "3D compatibility of Z vertex;3D compatibility of Z vertex; N(#mu#mu pairs)", 100, 0., 100.);

hTrackInvMass_ = iBooker.book1D("TkTkInvMass", ";M(tk,tk) [GeV];N(tk tk pairs)", 70., 50., 120.);
hCosPhi_ = iBooker.book1D("CosPhi", ";cos(#phi_{xy});N(#mu#mu pairs)", 50, -1., 1.);
hCosPhi3D_ = iBooker.book1D("CosPhi3D", ";cos(#phi_{3D});N(#mu#mu pairs)", 50, -1., 1.);
hCosPhiInv_ = iBooker.book1D("CosPhiInv", ";inverted cos(#phi_{xy});N(#mu#mu pairs)", 50, -1., 1.);
hCosPhiInv3D_ = iBooker.book1D("CosPhiInv3D", ";inverted cos(#phi_{3D});N(#mu#mu pairs)", 50, -1., 1.);

hdxy_ = iBooker.book1D("dxy", ";muon track d_{xy}(PV) [#mum];muon tracks", 150, -300, 300);
hdz_ = iBooker.book1D("dz", ";muon track d_{z}(PV) [#mum];muon tracks", 150, -300, 300);
hdxyErr_ = iBooker.book1D("dxyErr", ";muon track err_{dxy} [#mum];muon tracks", 250, 0., 500.);
hdzErr_ = iBooker.book1D("dzErr", ";muon track err_{dz} [#mum];muon tracks", 250, 0., 500.);
hIP2d_ = iBooker.book1D("IP2d", ";muon track IP_{2D} [#mum];muon tracks", 150, -300, 300);
hIP3d_ = iBooker.book1D("IP3d", ";muon track IP_{3D} [#mum];muon tracks", 150, -300, 300);
hIP2dsig_ = iBooker.book1D("IP2Dsig", ";muon track IP_{2D} significance;muon tracks", 100, 0., 5.);
hIP3dsig_ = iBooker.book1D("IP3Dsig", ";muon track IP_{3D} significance;muon tracks", 100, 0., 5.);

// clang-format off
std::string ts = fmt::sprintf(";%s vertex probability;N(#mu#mu pairs)", motherName_);
std::string ps = "N(#mu#mu pairs)";
hSVProb_ = iBooker.book1D("VtxProb", ts, 100, 0., 1.);

ts = fmt::sprintf("#chi^{2} of the %s vertex; #chi^{2} of the %s vertex; %s", motherName_, motherName_, ps);
hSVChi2_ = iBooker.book1D("VtxChi2", ts, 200, 0., 200.);

ts = fmt::sprintf("#chi^{2}/ndf of the %s vertex; #chi^{2}/ndf of %s vertex; %s", motherName_, motherName_, ps);
hSVNormChi2_ = iBooker.book1D("VtxNormChi2", ts, 100, 0., 20.);

std::string histTit = motherName_ + " #rightarrow #mu^{+}#mu^{-}";
ts = fmt::sprintf("%s;PV- %sV xy distance [#mum];%s", histTit, motherName_, ps);
hSVDist_ = iBooker.book1D("VtxDist", ts, 100, 0., 300.);

ts = fmt::sprintf("%s;PV-%sV xy distance error [#mum];%s", histTit, motherName_, ps);
hSVDistErr_ = iBooker.book1D("VtxDistErr", ts, 100, 0., 1000.);

ts = fmt::sprintf("%s;PV-%sV xy distance signficance;%s", histTit, motherName_, ps);
hSVDistSig_ = iBooker.book1D("VtxDistSig", ts, 100, 0., 5.);

ts = fmt::sprintf("compatibility of %s vertex; compatibility of %s vertex; %s", motherName_, motherName_, ps);
hSVCompatibility_ = iBooker.book1D("VtxCompatibility", ts, 100, 0., 100.);

ts = fmt::sprintf("%s;PV-%sV 3D distance [#mum];%s", histTit, motherName_, ps);
hSVDist3D_ = iBooker.book1D("VtxDist3D", ts, 100, 0., 300.);

ts = fmt::sprintf("%s;PV-%sV 3D distance error [#mum];%s", histTit, motherName_, ps);
hSVDist3DErr_ = iBooker.book1D("VtxDist3DErr", ts, 100, 0., 1000.);

ts = fmt::sprintf("%s;PV-%sV 3D distance signficance;%s", histTit, motherName_, ps);
hSVDist3DSig_ = iBooker.book1D("VtxDist3DSig", ts, 100, 0., 5.);

ts = fmt::sprintf("3D compatibility of %s vertex;3D compatibility of %s vertex; %s", motherName_, motherName_, ps);
hSVCompatibility3D_ = iBooker.book1D("VtxCompatibility3D", ts, 100, 0., 100.);

hInvMass_ = iBooker.book1D("InvMass", fmt::sprintf("%s;M(#mu,#mu) [GeV];%s", histTit, ps), 70., massLimits_.first, massLimits_.second);
hCosPhi_ = iBooker.book1D("CosPhi", fmt::sprintf("%s;cos(#phi_{xy});%s", histTit, ps), 50, -1., 1.);
hCosPhi3D_ = iBooker.book1D("CosPhi3D", fmt::sprintf("%s;cos(#phi_{3D});%s", histTit, ps), 50, -1., 1.);
hCosPhiInv_ = iBooker.book1D("CosPhiInv", fmt::sprintf("%s;inverted cos(#phi_{xy});%s", histTit, ps), 50, -1., 1.);
hCosPhiInv3D_ = iBooker.book1D("CosPhiInv3D", fmt::sprintf("%s;inverted cos(#phi_{3D});%s", histTit, ps), 50, -1., 1.);

hdxy_ = iBooker.book1D("dxy", fmt::sprintf("%s;muon track d_{xy}(PV) [#mum];muon tracks", histTit), 150, -300, 300);
hdz_ = iBooker.book1D("dz", fmt::sprintf("%s;muon track d_{z}(PV) [#mum];muon tracks", histTit), 150, -300, 300);
hdxyErr_ = iBooker.book1D("dxyErr", fmt::sprintf("%s;muon track err_{dxy} [#mum];muon tracks", histTit), 250, 0., 500.);
hdzErr_ = iBooker.book1D("dzErr", fmt::sprintf("%s;muon track err_{dz} [#mum];muon tracks", histTit), 250, 0., 500.);
hIP2d_ = iBooker.book1D("IP2d", fmt::sprintf("%s;muon track IP_{2D} [#mum];muon tracks", histTit), 150, -300, 300);
hIP3d_ = iBooker.book1D("IP3d", fmt::sprintf("%s;muon track IP_{3D} [#mum];muon tracks", histTit), 150, -300, 300);
hIP2dsig_ = iBooker.book1D("IP2Dsig", fmt::sprintf("%s;muon track IP_{2D} significance;muon tracks", histTit), 100, 0., 5.);
hIP3dsig_ = iBooker.book1D("IP3Dsig", fmt::sprintf("%s;muon track IP_{3D} significance;muon tracks", histTit), 100, 0., 5.);
// clang-format on
}

void DiMuonVertexMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
Expand Down Expand Up @@ -100,7 +136,7 @@ void DiMuonVertexMonitor::analyze(const edm::Event& iEvent, const edm::EventSetu

const auto& Zp4 = p4_tplus + p4_tminus;
float track_invMass = Zp4.M();
hTrackInvMass_->Fill(track_invMass);
hInvMass_->Fill(track_invMass);

// creat the pair of TLorentVectors used to make the plos
std::pair<TLorentzVector, TLorentzVector> tktk_p4 = std::make_pair(p4_tplus, p4_tminus);
Expand Down Expand Up @@ -256,6 +292,7 @@ void DiMuonVertexMonitor::fillDescriptions(edm::ConfigurationDescriptions& descr
desc.add<edm::InputTag>("muonTracks", edm::InputTag("ALCARECOTkAlDiMuon"));
desc.add<edm::InputTag>("vertices", edm::InputTag("offlinePrimaryVertices"));
desc.add<std::string>("FolderName", "DiMuonVertexMonitor");
desc.add<std::string>("decayMotherName", "Z");
desc.add<bool>("useClosestVertex", true);
desc.add<double>("maxSVdist", 50.);
descriptions.addWithDefaultLabel(desc);
Expand Down
39 changes: 31 additions & 8 deletions ...Alignment/test/DiMuonVertex_HARVESTING.py → ...gnment/test/DiMuonTkAlDQMHarvester_cfg.py
View file
Open in desktop
Expand Up @@ -22,6 +22,11 @@
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"conditions")
options.register('resonance',
'Z',
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.string,
"resonance type")
options.parseArguments()

# import of standard configurations
Expand All @@ -43,7 +48,7 @@

# Input source
process.source = cms.Source("DQMRootSource",
fileNames = cms.untracked.vstring('file:step3_inDQM.root')
fileNames = cms.untracked.vstring('file:step3_inDQM_'+options.resonance+'.root')
)

process.options = cms.untracked.PSet(
Expand Down Expand Up @@ -93,16 +98,34 @@

process.dqmsave_step = cms.Path(process.DQMSaver)

if (options.resonance == 'Z'):
print('',30*"#",'\n # will harvest Z file \n',30*"#")
_folderName = cms.string('AlCaReco/TkAlDiMuonAndVertex')
_doBkgFit = cms.bool(False)
_fitPar = cms.PSet(mean_par = cms.vdouble(90.,60.,120.),
width_par = cms.vdouble(5.0,0.0,120.0),
sigma_par = cms.vdouble(5.0,0.0,120.0))
elif (options.resonance == 'Jpsi'):
print('',30*"#",'\n # will harvest J/psi file \n',30*"#")
_folderName = cms.string('AlCaReco/TkAlJpsiMuMu')
_doBkgFit = cms.bool(True)
_fitPar = cms.PSet(mean_par = cms.vdouble(3.09, 2.7, 3.4),
width_par = cms.vdouble(1.0, 0.0, 5.0),
sigma_par = cms.vdouble(1.0, 0.0, 5.0))
elif (options.resonance == 'Upsilon'):
print('',30*"#",'\n # will harvest Upsilon file \n',30*"#")
_folderName = cms.string('AlCaReco/TkAlUpsilonMuMu')
_doBkgFit = cms.bool(True)
_fitPar = cms.PSet(mean_par = cms.vdouble(9.46, 8.9, 9.9),
width_par = cms.vdouble(1.0, 0.0, 5.0),
sigma_par = cms.vdouble(1.0, 0.0, 5.0))

# the module to harvest
process.DiMuonMassBiasClient = cms.EDProducer("DiMuonMassBiasClient",
FolderName = cms.string('AlCaReco/TkAlDiMuonAndVertex'),
fitBackground = cms.bool(False),
FolderName = _folderName,
fitBackground = _doBkgFit,
debugMode = cms.bool(True),
fit_par = cms.PSet(
mean_par = cms.vdouble(90.,60.,120.),
width_par = cms.vdouble(5.0,0.0,120.0),
sigma_par = cms.vdouble(5.0,0.0,120.0)
),
fit_par = _fitPar,
MEtoHarvest = cms.vstring(
'DiMuMassVsMuMuPhi',
'DiMuMassVsMuMuEta',
Expand Down