Skip to content

Commit 89b4d21

Browse files
SuJeong-Jialibuild
andauthored
[PWGLF] Added background histogram (#12401)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 522e496 commit 89b4d21

File tree

1 file changed

+12
-8
lines changed

1 file changed

+12
-8
lines changed

PWGLF/Tasks/Resonances/chk892pp.cxx

Lines changed: 12 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -24,7 +24,7 @@
2424
#include <TObjArray.h>
2525
#include <TFile.h>
2626
#include <TH2F.h>
27-
#include <TDatabasePDG.h> // FIXME
27+
// #include <TDatabasePDG.h> // FIXME
2828
#include <TPDGCode.h> // FIXME
2929

3030
#include <vector>
@@ -289,6 +289,10 @@ struct Chk892pp {
289289
histos.add("QA/before/VtxZ", "Centrality distribution", {HistType::kTH1D, {vtxzAxis}});
290290
histos.add("QA/before/hEvent", "Number of Events", HistType::kTH1F, {{1, 0.5, 1.5}});
291291

292+
if (BkgEstimationConfig.cfgFillRotBkg) {
293+
histos.add("QA/RotBkg/hRotBkg", "Rotated angle of rotated background", HistType::kTH1F, {{360, 0.0, o2::constants::math::TwoPI}});
294+
}
295+
292296
// Bachelor pion
293297
histos.add("QA/before/trkbpionDCAxy", "DCAxy distribution of bachelor pion candidates", HistType::kTH1D, {dcaxyAxis});
294298
histos.add("QA/before/trkbpionDCAz", "DCAz distribution of bachelor pion candidates", HistType::kTH1D, {dcazAxis});
@@ -666,16 +670,16 @@ struct Chk892pp {
666670
if (std::abs(motherbTrack.pdgCode()) != kKstarPlus) // Are you charged Kstar's daughter?
667671
return false; // Apply first since it's more restrictive
668672

669-
if (std::abs(motherkV0.pdgCode()) != 310) // Is it K0s?
673+
if (std::abs(motherkV0.pdgCode()) != kPDGK0s) // Is it K0s?
670674
return false;
671675
// Check if K0s's mother is K0 (311)
672676
auto motherK0 = motherkV0.template mothers_as<aod::McParticles>();
673-
if (std::abs(motherK0.pdgCode()) != 311)
677+
if (std::abs(motherK0.pdgCode()) != kPDGK0)
674678
return false;
675679

676680
// Check if K0's mother is Kstar (323)
677681
auto motherKstar = motherK0.template mothers_as<aod::McParticles>();
678-
if (std::abs(motherKstar.pdgCode()) != 323)
682+
if (std::abs(motherKstar.pdgCode()) != kKstarPlus)
679683
return false;
680684

681685
// Check if bTrack and K0 have the same mother (global index)
@@ -696,7 +700,7 @@ struct Chk892pp {
696700
std::vector<int> trackIndicies = {};
697701
std::vector<int> k0sIndicies = {};
698702

699-
for (auto& bTrack : dTracks1) {
703+
for (const auto& bTrack : dTracks1) {
700704
auto trkbpt = bTrack.pt();
701705
auto istrkbhasTOF = bTrack.hasTOF();
702706
auto trkbNSigmaPiTPC = bTrack.tpcNSigmaPi();
@@ -733,7 +737,7 @@ struct Chk892pp {
733737
trackIndicies.push_back(bTrack.index());
734738
}
735739

736-
for (auto& k0sCand : dTracks2) {
740+
for (const auto& k0sCand : dTracks2) {
737741
auto posDauTrack = k0sCand.template posTrack_as<TrackCandidates>();
738742
auto negDauTrack = k0sCand.template negTrack_as<TrackCandidates>();
739743

@@ -838,8 +842,8 @@ struct Chk892pp {
838842
k0sIndicies.push_back(k0sCand.index());
839843
}
840844

841-
for (auto& trackIndex : trackIndicies) {
842-
for (auto& k0sIndex : k0sIndicies) {
845+
for (const auto& trackIndex : trackIndicies) {
846+
for (const auto& k0sIndex : k0sIndicies) {
843847
auto bTrack = dTracks1.rawIteratorAt(trackIndex);
844848
auto k0sCand = dTracks2.rawIteratorAt(k0sIndex);
845849
auto trkkMass = k0sCand.mK0Short();

0 commit comments

Comments
 (0)