Skip to content

Commit 2f61fbf

Browse files
5ub-Man5ub-Man
authored
[PWGLF] Minor changes (#12760)
1 parent 552ad72 commit 2f61fbf

File tree

1 file changed

+34
-14
lines changed

1 file changed

+34
-14
lines changed

‎PWGLF/Tasks/Resonances/kstarqa.cxx‎

Lines changed: 34 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -251,8 +251,8 @@ struct Kstarqa {
251251
hOthers.add("hCRFC_after", "CRFC after distribution", kTH1F, {{100, 0.0f, 10.0f}});
252252
hOthers.add("hCRFC_before", "CRFC before distribution", kTH1F, {{100, 0.0f, 10.0f}});
253253

254-
hOthers.add("hKstar_Rap", "Pair rapidity distribution; y; Counts", kTH1F, {{1000, -5.0f, 5.0f}});
255-
hOthers.add("hKstar_Eta", "Pair eta distribution; #eta; Counts", kTH1F, {{1000, -5.0f, 5.0f}});
254+
hOthers.add("hKstar_rap_pt", "Pair rapidity distribution; y; p_{T}; Counts", kTH2F, {{400, -2.0f, 2.0f}, ptAxis});
255+
hOthers.add("hKstar_eta_pt", "Pair eta distribution; #eta; p_{T}; Counts", kTH2F, {{400, -2.0f, 2.0f}, ptAxis});
256256

257257
hPID.add("Before/hNsigmaTPC_Ka_before", "N #sigma Kaon TPC before", kTH2F, {{50, 0.0f, 10.0f}, {100, -10.0f, 10.0f}});
258258
hPID.add("Before/hNsigmaTOF_Ka_before", "N #sigma Kaon TOF before", kTH2F, {{50, 0.0f, 10.0f}, {100, -10.0f, 10.0f}});
@@ -387,7 +387,7 @@ struct Kstarqa {
387387
double massKa = o2::constants::physics::MassKPlus;
388388

389389
template <typename Coll>
390-
bool selectionEvent(const Coll& collision, bool fillHist = true)
390+
bool selectionEvent(const Coll& collision, bool fillHist = false) // default to false
391391
{
392392
if (fillHist)
393393
rEventSelection.fill(HIST("hEventCut"), 0);
@@ -937,7 +937,7 @@ struct Kstarqa {
937937
int occupancy = collision.trackOccupancyInTimeRange();
938938
rEventSelection.fill(HIST("hOccupancy"), occupancy);
939939

940-
if (!selectionEvent(collision, true)) {
940+
if (!selectionEvent(collision, true)) { // fill event cut histogram
941941
return;
942942
}
943943

@@ -1095,8 +1095,8 @@ struct Kstarqa {
10951095
continue;
10961096
}
10971097

1098-
hOthers.fill(HIST("hKstar_Rap"), mother.Rapidity());
1099-
hOthers.fill(HIST("hKstar_Eta"), mother.Eta());
1098+
hOthers.fill(HIST("hKstar_rap_pt"), mother.Rapidity(), mother.Pt());
1099+
hOthers.fill(HIST("hKstar_eta_pt"), mother.Eta(), mother.Pt());
11001100

11011101
isMix = false;
11021102
fillInvMass(daughter1, daughter2, mother, multiplicity, isMix, track1, track2);
@@ -1148,7 +1148,7 @@ struct Kstarqa {
11481148
// if (!c1.sel8() || !c2.sel8())
11491149
// continue;
11501150

1151-
if (!selectionEvent(c1, false) || !selectionEvent(c2, false)) {
1151+
if (!selectionEvent(c1, false) || !selectionEvent(c2, false)) { // don't fill event cut histogram
11521152
continue;
11531153
}
11541154

@@ -1214,7 +1214,7 @@ struct Kstarqa {
12141214
auto runMixing = [&](auto& pair, auto multiplicityGetter) {
12151215
for (const auto& [c1, tracks1, c2, tracks2] : pair) {
12161216

1217-
if (!selectionEvent(c1, false) || !selectionEvent(c2, false)) {
1217+
if (!selectionEvent(c1, false) || !selectionEvent(c2, false)) { // don't fill event cut histogram
12181218
continue;
12191219
}
12201220

@@ -1297,7 +1297,7 @@ struct Kstarqa {
12971297
int occupancy = collision.trackOccupancyInTimeRange();
12981298
rEventSelection.fill(HIST("hOccupancy"), occupancy);
12991299

1300-
if (!selectionEvent(collision, false)) {
1300+
if (!selectionEvent(collision, false)) { // don't fill event cut histogram
13011301
return;
13021302
}
13031303

@@ -1533,7 +1533,7 @@ struct Kstarqa {
15331533
rEventSelection.fill(HIST("eventsCheckGen"), 2.5);
15341534

15351535
for (const auto& collision : collisions) {
1536-
if (!selectionEvent(collision, true)) {
1536+
if (!selectionEvent(collision, false)) { // don't fill event cut histogram
15371537
continue;
15381538
}
15391539
multiplicity = collision.centFT0M();
@@ -1657,7 +1657,7 @@ struct Kstarqa {
16571657
for (const auto& RecCollision : recCollisions) {
16581658
if (!RecCollision.has_mcCollision())
16591659
continue;
1660-
if (!selectionEvent(RecCollision, false))
1660+
if (!selectionEvent(RecCollision, false)) // don't fill event cut histogram
16611661
continue;
16621662
// multiplicity1 = RecCollision.centFT0M();
16631663
const auto& mcCollisionRec = RecCollision.mcCollision_as<EventMCGenerated>();
@@ -1732,13 +1732,17 @@ struct Kstarqa {
17321732
hInvMass.fill(HIST("hAllRecCollisions"), multiplicity);
17331733
hInvMass.fill(HIST("hAllRecCollisionsCalib"), multiplicityRec);
17341734

1735-
if (!selectionEvent(collision, false)) {
1735+
if (!selectionEvent(collision, false)) { // don't fill event cut histogram
17361736
return;
17371737
}
17381738

17391739
hInvMass.fill(HIST("h1RecMult"), multiplicity);
17401740
hInvMass.fill(HIST("h1RecMult2"), multiplicityRec);
17411741

1742+
if (cQAevents) {
1743+
rEventSelection.fill(HIST("hVertexZRec"), collision.posZ());
1744+
}
1745+
17421746
auto oldindex = -999;
17431747
for (const auto& track1 : tracks) {
17441748
if (!selectionTrack(track1)) {
@@ -1749,6 +1753,11 @@ struct Kstarqa {
17491753
continue;
17501754
}
17511755

1756+
if (cQAevents) {
1757+
rEventSelection.fill(HIST("hDcaxy"), track1.dcaXY());
1758+
rEventSelection.fill(HIST("hDcaz"), track1.dcaZ());
1759+
}
1760+
17521761
auto track1ID = track1.index();
17531762
for (const auto& track2 : tracks) {
17541763
rEventSelection.fill(HIST("recMCparticles"), 0.5);
@@ -1777,7 +1786,12 @@ struct Kstarqa {
17771786
const auto mctrack2 = track2.mcParticle();
17781787
int track1PDG = std::abs(mctrack1.pdgCode());
17791788
int track2PDG = std::abs(mctrack2.pdgCode());
1780-
1789+
if (cQAplots) {
1790+
hPID.fill(HIST("Before/hTPCnsigKa_mult_pt"), track1.tpcNSigmaKa(), multiplicity, track1.pt());
1791+
hPID.fill(HIST("Before/hTPCnsigPi_mult_pt"), track2.tpcNSigmaPi(), multiplicity, track2.pt());
1792+
hPID.fill(HIST("Before/hTOFnsigKa_mult_pt"), track1.tofNSigmaKa(), multiplicity, track1.pt());
1793+
hPID.fill(HIST("Before/hTOFnsigPi_mult_pt"), track2.tofNSigmaKa(), multiplicity, track2.pt());
1794+
}
17811795
if (cQAplots && (mctrack2.pdgCode() == PDG_t::kPiPlus)) { // pion
17821796
hPID.fill(HIST("Before/h1PID_TPC_pos_pion"), track2.tpcNSigmaPi());
17831797
hPID.fill(HIST("Before/h1PID_TOF_pos_pion"), track2.tofNSigmaPi());
@@ -1906,6 +1920,12 @@ struct Kstarqa {
19061920
continue;
19071921
rEventSelection.fill(HIST("recMCparticles"), 13.5);
19081922
}
1923+
if (cQAplots) {
1924+
hPID.fill(HIST("After/hTPCnsigKa_mult_pt"), track1.tpcNSigmaKa(), multiplicity, track1.pt());
1925+
hPID.fill(HIST("After/hTPCnsigPi_mult_pt"), track2.tpcNSigmaPi(), multiplicity, track2.pt());
1926+
hPID.fill(HIST("After/hTOFnsigKa_mult_pt"), track1.tofNSigmaKa(), multiplicity, track1.pt());
1927+
hPID.fill(HIST("After/hTOFnsigPi_mult_pt"), track2.tofNSigmaKa(), multiplicity, track2.pt());
1928+
}
19091929

19101930
if (selectionConfig.isApplyCutsOnMother) {
19111931
if (mothertrack1.pt() >= selectionConfig.cMaxPtMotherCut) // excluding candidates in overflow
@@ -1975,7 +1995,7 @@ struct Kstarqa {
19751995
hInvMass.fill(HIST("hAllRecCollisions"), multiplicity);
19761996
hInvMass.fill(HIST("hAllRecCollisionsCalib"), multiplicityRec);
19771997

1978-
if (!selectionEvent(collision, false)) {
1998+
if (!selectionEvent(collision, false)) { // don't fill event cut histogram
19791999
return;
19802000
}
19812001

0 commit comments

Comments
 (0)