Skip to content

Commit

Permalink
Merge pull request #43027 from llunerti/from-CMSSW_13_2_6
Browse files Browse the repository at this point in the history
fixed dimuon invariant mass plots range
  • Loading branch information
cmsbuild authored Feb 21, 2024
2 parents 522751f + 6f31d60 commit 9b282b2
Show file tree
Hide file tree
Showing 3 changed files with 46 additions and 46 deletions.
8 changes: 4 additions & 4 deletions DQMOffline/Muon/interface/DiMuonHistograms.h
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
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
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],
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

0 comments on commit 9b282b2

Please sign in to comment.