Skip to content

Commit 39b53ee

Browse files
rashiguptrashi234
andauthored
[PWGHF] Add histograms in Non-HFE analysis (#15491)
Co-authored-by: rashi.gupta@cern.ch <rashi.gupta@cern.ch>
1 parent 344bfe5 commit 39b53ee

File tree

1 file changed

+26
-11
lines changed

1 file changed

+26
-11
lines changed

PWGHF/HFL/TableProducer/electronSelectionWithTpcEmcal.cxx

Lines changed: 26 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -189,14 +189,20 @@ struct HfElectronSelectionWithTpcEmcal {
189189
registry.add("hZvertex", "z vertex", {HistType::kTH1D, {axisPosZ}});
190190
registry.add("hNeventsAfterPassEmcal", "No of events pass the Emcal", {HistType::kTH1D, {{3, 1, 4}}});
191191
registry.add("hNevents", "No of events", {HistType::kTH1D, {{3, 1, 4}}});
192-
registry.add("hLikeMass", "Like mass", {HistType::kTH1D, {{axisMass}}});
193-
registry.add("hUnLikeMass", "unLike mass", {HistType::kTH1D, {{axisMass}}});
194-
registry.add("hLikeSignPt", "Like sign Momentum ", {HistType::kTH1D, {{axisPt}}});
195-
registry.add("hUnLikeSignPt", "UnLike sign Momentum", {HistType::kTH1D, {{axisPt}}});
192+
registry.add("hLikeMass_EMCAL", "Like mass Emcal", {HistType::kTH1D, {{axisMass}}});
193+
registry.add("hUnLikeMass_EMCAL", "unLike mass Emcal", {HistType::kTH1D, {{axisMass}}});
194+
registry.add("hLikeSignPt_EMCAL", "Like sign Momentum Emcal ", {HistType::kTH1D, {{axisPt}}});
195+
registry.add("hUnLikeSignPt_EMCAL", "UnLike sign Momentum Emcal", {HistType::kTH1D, {{axisPt}}});
196+
registry.add("hLikeMass_NoEMCAL", "Like mass NoEMCAL", {HistType::kTH1D, {{axisMass}}});
197+
registry.add("hUnLikeMass_NoEMCAL", "unLike mass NoEMCAL", {HistType::kTH1D, {{axisMass}}});
198+
registry.add("hLikeSignPt_NoEMCAL", "Like sign Momentum NoEMCAL ", {HistType::kTH1D, {{axisPt}}});
199+
registry.add("hUnLikeSignPt_NoEMCAL", "UnLike sign Momentum NoEMCAL", {HistType::kTH1D, {{axisPt}}});
200+
196201
registry.add("hMcgenInElectron", "Mc Gen Inclusive Electron", {HistType::kTH1D, {{axisPt}}});
197202
registry.add("hMcRecInElectron", "Mc Rec Inclusive Electron", {HistType::kTH1D, {{axisPt}}});
198203
registry.add("hMcRecwithoutEMCalInElectron", "Mc Rec Inclusive Electron without Emcal", {HistType::kTH1D, {{axisPt}}});
199-
204+
registry.add("hphiElectron", "hphiElectron", {HistType::kTH1D, {axisPhi}});
205+
registry.add("hphiElectronPassEmcal", "hphiElectron pass Emcal", {HistType::kTH1D, {axisPhi}});
200206
registry.add("hMcgenAllNonHfeElectron", "Mc Gen All NonHf Electron", {HistType::kTH1D, {{axisPt}}});
201207
registry.add("hMcgenNonHfeElectron", "Mc Gen NonHf Electron with mother", {HistType::kTH1D, {{axisPt}}});
202208
registry.add("hPi0eEmbTrkPt", "Mc Gen Pi0 mother NonHf Electron", {HistType::kTH1D, {{axisPt}}});
@@ -353,7 +359,9 @@ struct HfElectronSelectionWithTpcEmcal {
353359
vecLSMass.push_back(massLike);
354360
isLSElectron = true;
355361
if (isEMcal) {
356-
registry.fill(HIST("hLikeMass"), massLike);
362+
registry.fill(HIST("hLikeMass_EMCAL"), massLike);
363+
} else {
364+
registry.fill(HIST("hLikeMass_NoEMCAL"), massLike);
357365
}
358366
}
359367
// for unlike charge
@@ -362,7 +370,9 @@ struct HfElectronSelectionWithTpcEmcal {
362370
vecULSMass.push_back(massUnLike);
363371
isULSElectron = true;
364372
if (isEMcal) {
365-
registry.fill(HIST("hUnLikeMass"), massUnLike);
373+
registry.fill(HIST("hUnLikeMass_EMCAL"), massUnLike);
374+
} else {
375+
registry.fill(HIST("hUnLikeMass_NoEMCAL"), massUnLike);
366376
}
367377
}
368378

@@ -371,15 +381,19 @@ struct HfElectronSelectionWithTpcEmcal {
371381
massLike = invMassElectron;
372382
++nElPairsLS;
373383
if (isEMcal) {
374-
registry.fill(HIST("hLikeSignPt"), electron.pt());
384+
registry.fill(HIST("hLikeSignPt_EMCAL"), electron.pt());
385+
} else {
386+
registry.fill(HIST("hLikeSignPt_NoEMCAL"), electron.pt());
375387
}
376388
}
377389
// for unlike charge
378390
if (isULSElectron && (invMassElectron <= invariantMass)) {
379391
massUnLike = invMassElectron;
380392
++nElPairsUS;
381393
if (isEMcal) {
382-
registry.fill(HIST("hUnLikeSignPt"), electron.pt());
394+
registry.fill(HIST("hUnLikeSignPt_EMCAL"), electron.pt());
395+
} else {
396+
registry.fill(HIST("hUnLikeSignPt_NoEMCAL"), electron.pt());
383397
}
384398
}
385399
}
@@ -445,7 +459,7 @@ struct HfElectronSelectionWithTpcEmcal {
445459
}
446460

447461
if (fillTrackInfo) {
448-
registry.fill(HIST("hTrackEtaPhi"), etaTrack, phiTrack, passEMCal); // track etaphi infor after filter bit
462+
// track etaphi infor after filter bit
449463
registry.fill(HIST("hTrackEnergyLossVsP"), track.tpcSignal(), pTrack, passEMCal); // track etaphi infor after filter bit
450464
registry.fill(HIST("hTrackEnergyLossVsPt"), track.tpcSignal(), ptTrack, passEMCal); // track etaphi infor after filter bit
451465
registry.fill(HIST("hTracknSigmaVsP"), tpcNsigmaTrack, pTrack, passEMCal); // track etaphi infor after filter bit
@@ -530,7 +544,7 @@ struct HfElectronSelectionWithTpcEmcal {
530544
if (eop < eopElectronMin || eop > eopElectronMax) {
531545
continue;
532546
}
533-
547+
registry.fill(HIST("hphiElectronPassEmcal"), track.phi());
534548
///////////////// NonHf electron Selection with Emcal ////////////////////////
535549
if constexpr (IsMc) {
536550
if (matchTrack.has_mcParticle()) {
@@ -648,6 +662,7 @@ struct HfElectronSelectionWithTpcEmcal {
648662
if ((track.tpcNSigmaEl() < tpcNsigmaElectronMin || track.tpcNSigmaEl() > tpcNsigmaElectronMax)) {
649663
continue;
650664
}
665+
registry.fill(HIST("hphiElectron"), track.phi());
651666
if constexpr (IsMc) {
652667
if (track.has_mcParticle()) {
653668
auto mcParticle = track.template mcParticle_as<aod::McParticles>();

0 commit comments

Comments
 (0)