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

fixed dimuon invariant mass plots range #43027

Merged
merged 3 commits into from
Feb 21, 2024
Merged

fixed dimuon invariant mass plots range #43027

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
b374509
fixed dimuon invariant mass plots range
Oct 16, 2023
2fb45d2
variables name fix
Dec 11, 2023
6f31d60
missing fix on cfg file
Dec 11, 2023
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
8 changes: 4 additions & 4 deletions DQMOffline/Muon/interface/DiMuonHistograms.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -64,10 +64,10 @@ class DiMuonHistograms : public DQMEDAnalyzer {
double etaECMax;

//Defining the relevant invariant mass regions
double LowMassMin;
double LowMassMax;
double HighMassMin;
double HighMassMax;
double lowMassMin;
double lowMassMax;
double highMassMin;
double highMassMax;

std::vector<MonitorElement*> GlbGlbMuon_LM;
std::vector<MonitorElement*> GlbGlbMuon_HM;
Expand Down
24 changes: 12 additions & 12 deletions DQMOffline/Muon/python/diMuonHistograms_cfi.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -17,10 +17,10 @@
etaECMin = cms.double(0.9),
etaECMax = cms.double(2.4),

LowMassMin = cms.double(2.0),
LowMassMax = cms.double(12.0),
HighMassMin = cms.double(70.0),
HighMassMax = cms.double(110.0),
lowMassMin = cms.double(2.0),
lowMassMax = cms.double(12.0),
highMassMin = cms.double(70.0),
highMassMax = cms.double(110.0),
folder = cms.string("Muons/diMuonHistograms")
)
diMuonHistos_miniAOD = DQMEDAnalyzer('DiMuonHistograms',
Expand All @@ -37,10 +37,10 @@
etaECMin = cms.double(0.9),
etaECMax = cms.double(2.4),

LowMassMin = cms.double(2.0),
LowMassMax = cms.double(12.0),
HighMassMin = cms.double(70.0),
HighMassMax = cms.double(110.0),
lowMassMin = cms.double(2.0),
lowMassMax = cms.double(12.0),
highMassMin = cms.double(70.0),
highMassMax = cms.double(110.0),
folder = cms.string("Muons_miniAOD/diMuonHistograms")
)

Expand All @@ -50,8 +50,8 @@
etaBBin = 350,
etaEBin = 350,

LowMassMin = 2.0,
LowMassMax = 51.0,
HighMassMin = 55.0,
HighMassMax = 125.0
lowMassMin = 2.0,
lowMassMax = 51.0,
highMassMin = 55.0,
highMassMax = 125.0
)
60 changes: 30 additions & 30 deletions DQMOffline/Muon/src/DiMuonHistograms.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -53,10 +53,10 @@ DiMuonHistograms::DiMuonHistograms(const edm::ParameterSet& pSet) {
etaECMin = parameters.getParameter<double>("etaECMin");
etaECMax = parameters.getParameter<double>("etaECMax");

LowMassMin = parameters.getParameter<double>("LowMassMin");
LowMassMax = parameters.getParameter<double>("LowMassMax");
HighMassMin = parameters.getParameter<double>("HighMassMin");
HighMassMax = parameters.getParameter<double>("HighMassMax");
lowMassMin = parameters.getParameter<double>("lowMassMin");
lowMassMax = parameters.getParameter<double>("lowMassMax");
highMassMin = parameters.getParameter<double>("highMassMin");
highMassMax = parameters.getParameter<double>("highMassMax");

theFolder = parameters.getParameter<string>("folder");
}
Expand All @@ -78,51 +78,51 @@ void DiMuonHistograms::bookHistograms(DQMStore::IBooker& ibooker,
GlbGlbMuon_LM.push_back(ibooker.book1D("GlbGlbMuon_LM" + EtaName[iEtaRegion],
"InvMass_{GLB,GLB}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
LowMassMin,
LowMassMax));
lowMassMin,
lowMassMax));
TrkTrkMuon_LM.push_back(ibooker.book1D("TrkTrkMuon_LM" + EtaName[iEtaRegion],
"InvMass_{TRK,TRK}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
LowMassMin,
LowMassMax));
lowMassMin,
lowMassMax));
StaTrkMuon_LM.push_back(ibooker.book1D("StaTrkMuon_LM" + EtaName[iEtaRegion],
"InvMass_{STA,TRK}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
LowMassMin,
LowMassMax));
lowMassMin,
lowMassMax));

GlbGlbMuon_HM.push_back(ibooker.book1D("GlbGlbMuon_HM" + EtaName[iEtaRegion],
"InvMass_{GLB,GLB}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
HighMassMin,
HighMassMax));
highMassMin,
highMassMax));
TrkTrkMuon_HM.push_back(ibooker.book1D("TrkTrkMuon_HM" + EtaName[iEtaRegion],
"InvMass_{TRK,TRK}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
HighMassMin,
HighMassMax));
highMassMin,
highMassMax));
StaTrkMuon_HM.push_back(ibooker.book1D("StaTrkMuon_HM" + EtaName[iEtaRegion],
"InvMass_{STA,TRK}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
HighMassMin,
HighMassMax));
highMassMin,
highMassMax));

// arround the Z peak
TightTightMuon.push_back(ibooker.book1D("TightTightMuon" + EtaName[iEtaRegion],
"InvMass_{Tight,Tight}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
LowMassMin,
LowMassMax));
highMassMin,
highMassMax));
MediumMediumMuon.push_back(ibooker.book1D("MediumMediumMuon" + EtaName[iEtaRegion],
"InvMass_{Medium,Medium}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
rappoccio marked this conversation as resolved.
Show resolved Hide resolved
LowMassMin,
LowMassMax));
highMassMin,
highMassMax));
LooseLooseMuon.push_back(ibooker.book1D("LooseLooseMuon" + EtaName[iEtaRegion],
"InvMass_{Loose,Loose}" + EtaName[iEtaRegion],
nBin[iEtaRegion],
LowMassMin,
LowMassMax));
highMassMin,
highMassMax));
//Fraction of bad hits in the tracker track to the total
TightTightMuonBadFrac.push_back(ibooker.book1D(
"TightTightMuonBadFrac" + EtaName[iEtaRegion], "BadFrac_{Tight,Tight}" + EtaName[iEtaRegion], 10, 0, 0.4));
Expand Down Expand Up @@ -222,9 +222,9 @@ void DiMuonHistograms::analyze(const edm::Event& iEvent, const edm::EventSetup&
fabs(recoCombinedGlbTrack1->eta()) < EtaCutMax[iEtaRegion] &&
fabs(recoCombinedGlbTrack2->eta()) > EtaCutMin[iEtaRegion] &&
fabs(recoCombinedGlbTrack2->eta()) < EtaCutMax[iEtaRegion]) {
if (InvMass < LowMassMax)
if (InvMass < lowMassMax)
GlbGlbMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin)
if (InvMass > highMassMin)
GlbGlbMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
Expand Down Expand Up @@ -300,9 +300,9 @@ void DiMuonHistograms::analyze(const edm::Event& iEvent, const edm::EventSetup&
if (fabs(recoStaTrack->eta()) > EtaCutMin[iEtaRegion] &&
fabs(recoStaTrack->eta()) < EtaCutMax[iEtaRegion] && fabs(recoTrack->eta()) > EtaCutMin[iEtaRegion] &&
fabs(recoTrack->eta()) < EtaCutMax[iEtaRegion]) {
if (InvMass < LowMassMax)
if (InvMass < lowMassMax)
StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin)
if (InvMass > highMassMin)
StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
Expand All @@ -322,9 +322,9 @@ void DiMuonHistograms::analyze(const edm::Event& iEvent, const edm::EventSetup&
if (fabs(recoStaTrack->eta()) > EtaCutMin[iEtaRegion] &&
fabs(recoStaTrack->eta()) < EtaCutMax[iEtaRegion] && fabs(recoTrack->eta()) > EtaCutMin[iEtaRegion] &&
fabs(recoTrack->eta()) < EtaCutMax[iEtaRegion]) {
if (InvMass < LowMassMax)
if (InvMass < lowMassMax)
StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin)
if (InvMass > highMassMin)
StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
Expand All @@ -345,9 +345,9 @@ void DiMuonHistograms::analyze(const edm::Event& iEvent, const edm::EventSetup&
for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
if (fabs(recoTrack1->eta()) > EtaCutMin[iEtaRegion] && fabs(recoTrack1->eta()) < EtaCutMax[iEtaRegion] &&
fabs(recoTrack2->eta()) > EtaCutMin[iEtaRegion] && fabs(recoTrack2->eta()) < EtaCutMax[iEtaRegion]) {
if (InvMass < LowMassMax)
if (InvMass < lowMassMax)
TrkTrkMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin)
if (InvMass > highMassMin)
TrkTrkMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
Expand Down