Skip to content

Commit bbefec4

Browse files
authored
enum EventSelection, MuonSelection updated
1 parent f55d76d commit bbefec4

File tree

1 file changed

+15
-13
lines changed

1 file changed

+15
-13
lines changed

PWGHF/HFL/Tasks/taskSingleMuonMult.cxx

Lines changed: 15 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -53,19 +53,21 @@ using namespace o2::aod::fwdtrack;
5353
struct HfTaskSingleMuonMult {
5454

5555
// enum for event selection bins
56-
enum EventSelectionBinsEnum {
56+
enum EventSelection {
5757
AllEvents = 1,
5858
Sel8,
59-
VtxZAfterSel
59+
VtxZAfterSel,
60+
NEventSelection
6061
};
6162

6263
// enum for muon track selection bins
63-
enum MuonTrackSelectionBinsEnum {
64+
enum MuonSelection {
6465
NoCut = 1,
6566
EtaCut,
6667
RAbsorbCut,
6768
PDcaCut,
68-
Chi2Cut
69+
Chi2Cut,
70+
NMuonSelection
6971
};
7072

7173
Configurable<float> zVtxMax{"zVtxMax", 10., "maxium z of primary vertex [cm]"};
@@ -96,9 +98,9 @@ struct HfTaskSingleMuonMult {
9698
void init(InitContext&)
9799
{
98100
AxisSpec axisCent = {101, -0.5, 100.5, "centrality"};
99-
AxisSpec axisEvent{3, 0.5, 3.5, "Event Selection"};
101+
AxisSpec axisEvent{NEventSelection, 0.5, NEventSelection + 0.5, "Muon Selection"};
100102
AxisSpec axisVtxZ{80, -20., 20., "#it{z}_{vtx} (cm)"};
101-
AxisSpec axisMuTrk{5, 0.5, 5.5, "Muon Selection"};
103+
AxisSpec axisMuon{NMuonSelection, 0.5, NMuonSelection + 0.5, "Muon Selection"};
102104
AxisSpec axisNCh{500, 0.5, 500.5, "#it{N}_{ch}"};
103105
AxisSpec axisNMu{20, -0.5, 19.5, "#it{N}_{#mu}"};
104106
AxisSpec axisPt{1000, 0., 500., "#it{p}_{T} (GeV/#it{c})"};
@@ -121,7 +123,7 @@ struct HfTaskSingleMuonMult {
121123
registry.add("hVtxZBeforeSel", "Z-vertex distribution before zVtx Cut", {HistType::kTH1F, {axisVtxZ}});
122124
registry.add("hVtxZAfterSel", "Z-vertex distribution after zVtx Cut", {HistType::kTH1F, {axisVtxZ}});
123125

124-
registry.add("hMuTrkSel", "Selection of muon tracks at various kinematic cuts", {HistType::kTH1F, {axisMuTrk}});
126+
registry.add("hMuonSel", "Selection of muon tracks at various kinematic cuts", {HistType::kTH1F, {axisMuon}});
125127
registry.add("hMuBeforeMatchMFT", "Muon information before any Kinemeatic cuts applied", {HistType::kTHnSparseF, {axisCent, axisNCh, axisPt, axisEta, axisTheta, axisRAbsorb, axisDCA, axisPDca, axisChi2MatchMCHMFT, axisTrackType}, 10});
126128
registry.add("hMuBeforeAccCuts", "Muon information before applying Acceptance cuts", {HistType::kTHnSparseF, {axisCent, axisNCh, axisPt, axisEta, axisTheta, axisRAbsorb, axisDCA, axisPDca, axisChi2MatchMCHMFT, axisTrackType}, 10});
127129
registry.add("h3DCABeforeAccCuts", "DCAx,DCAy,DCAz information before Acceptance cuts", {HistType::kTH3F, {axisDCAx, axisDCAx, axisTrackType}});
@@ -143,7 +145,7 @@ struct HfTaskSingleMuonMult {
143145
xEv->SetBinLabel(Sel8, "sel8");
144146
xEv->SetBinLabel(VtxZAfterSel, "VtxZAfterSel");
145147

146-
auto hMustat = registry.get<TH1>(HIST("hMuTrkSel"));
148+
auto hMustat = registry.get<TH1>(HIST("hMuonSel"));
147149
auto* xMu = hMustat->GetXaxis();
148150
xMu->SetBinLabel(NoCut, "noCut");
149151
xMu->SetBinLabel(EtaCut, "etaCut");
@@ -203,7 +205,7 @@ struct HfTaskSingleMuonMult {
203205
registry.fill(HIST("hMuBeforeMatchMFT"), cent, nCh, pt, eta, theta, rAbsorb, dcaXY, pDca, chi2, muTrackType);
204206

205207
// histograms before the acceptance cuts
206-
registry.fill(HIST("hMuTrkSel"), 1.);
208+
registry.fill(HIST("hMuonSel"), 1.);
207209
registry.fill(HIST("hMuBeforeAccCuts"), cent, nCh, pt, eta, theta, rAbsorb, dcaXY, pDca, chi2, muTrackType);
208210
registry.fill(HIST("h3DCABeforeAccCuts"), muon.fwdDcaX(), muon.fwdDcaY(), muTrackType);
209211

@@ -220,14 +222,14 @@ struct HfTaskSingleMuonMult {
220222
if ((eta >= etaMax) || (eta < etaMin)) {
221223
continue;
222224
}
223-
registry.fill(HIST("hMuTrkSel"), 2.);
225+
registry.fill(HIST("hMuonSel"), 2.);
224226
registry.fill(HIST("hMuAfterEtaCuts"), cent, nCh, pt, eta, theta, rAbsorb, dcaXY, pDca, chi2, muTrackType);
225227

226228
// Rabsorb cuts
227229
if ((rAbsorb < rAbsorbMin) || (rAbsorb >= rAbsorbMax)) {
228230
continue;
229231
}
230-
registry.fill(HIST("hMuTrkSel"), 3.);
232+
registry.fill(HIST("hMuonSel"), 3.);
231233
registry.fill(HIST("hMuAfterRAbsorbCuts"), cent, nCh, pt, eta, theta, rAbsorb, dcaXY, pDca, chi2, muTrackType);
232234

233235
if ((rAbsorb < rAbsorbMid) && (pDca >= pDcaMin)) {
@@ -236,14 +238,14 @@ struct HfTaskSingleMuonMult {
236238
if ((rAbsorb >= rAbsorbMid) && (pDca >= pDcaMax)) {
237239
continue;
238240
}
239-
registry.fill(HIST("hMuTrkSel"), 4.);
241+
registry.fill(HIST("hMuonSel"), 4.);
240242
registry.fill(HIST("hMuAfterPdcaCuts"), cent, nCh, pt, eta, theta, rAbsorb, dcaXY, pDca, chi2, muTrackType);
241243

242244
// MCH-MFT matching chi2
243245
if (muon.chi2() >= 1e6) {
244246
continue;
245247
}
246-
registry.fill(HIST("hMuTrkSel"), 5.);
248+
registry.fill(HIST("hMuonSel"), 5.);
247249

248250
// histograms after acceptance cuts
249251
registry.fill(HIST("hMuAfterAccCuts"), cent, nCh, pt, eta, theta, rAbsorb, dcaXY, pDca, chi2, muTrackType);

0 commit comments

Comments
 (0)