Skip to content

[PWGCF] add occupancy cut for MC data #18428

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
abilandz merged 1 commit into from
Oct 23, 2026
+20 −5
Merged

[PWGCF] add occupancy cut for MC data #18428

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
25 changes: 20 additions & 5 deletions PWGCF/MultiparticleCorrelations/Tasks/multiparticleCumulants.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -230,6 +230,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
Configurable<bool> cfCentCutSwitch{"cfCentCutSwitch", true, "centrality cut switch"};
Configurable<bool> cfNumContribCutSwitch{"cfNumContribCutSwitch", true, "NContribution cut switch"};
Configurable<bool> cfNumTrackCutSwitch{"cfNumTrackCutSwitch", true, "NTrack cut switch"};
Configurable<bool> cfOccupancyCutSwitch{"cfOccupancyCutSwitch", true, "occupancy cut switch"};
Configurable<bool> cfCentCorrCutSwitch{"cfCentCorrCutSwitch", true, "centrality correlation outlier cut switch"};
Configurable<bool> cfMultCorrCutSwitch{"cfMultCorrCutSwitch", true, "multiplicity correlation outlier cut switch"};
Configurable<bool> cfNumContribMultCorrCutSwitch{"cfNumContribMultCorrCutSwitch", true, "NContribution vs. multiplicity correlation outlier cut switch"};
Expand All @@ -248,6 +249,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
Configurable<std::vector<float>> cfCentCut{"cfCentCut", {10., 20.}, "centrality range: {min, max}[%]"};
Configurable<std::vector<float>> cfNumContribCut{"cfNumContribCut", {0, 3000.}, "NContribution range: {min, max}"};
Configurable<std::vector<float>> cfNumTrackCut{"cfNumTrackCut", {0, 3000.}, "NTrack range: {min, max}"};
Configurable<std::vector<float>> cfOccupancyCut{"cfOccupancyCut", {0, 6000.}, "Occupancy range: {min, max}"};
Configurable<std::vector<float>> cfCentABCorrCut{"cfCentABCorrCut", {1.4, 300., 1.4, 300.}, "parameters of limits in FT0C vs. FT0M centrality 2D histograms, cx+d < y < ax+b: {a,b,c,d}"};
Configurable<std::vector<float>> cfCentACCorrCut{"cfCentACCorrCut", {1.4, 300., 1.4, 300.}, "parameters of limits in FT0C vs. FT0M centrality 2D histograms, cx+d < y < ax+b: {a,b,c,d}"};
Configurable<std::vector<float>> cfCentBCCorrCut{"cfCentBCCorrCut", {1.4, 300., 1.4, 300.}, "parameters of limits in FT0C vs. FT0M centrality 2D histograms, cx+d < y < ax+b: {a,b,c,d}"};
Expand Down Expand Up @@ -302,6 +304,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
bool fCentCutSwitch = true;
bool fNumContribCutSwitch = true;
bool fNumTrackCutSwitch = true;
bool fOccupancyCutSwitch = true;
bool fCentCorrCutSwitch = true;
bool fMultCorrCutSwitch = true;
bool fNumContribMultCorrCutSwitch = true;
Expand All @@ -318,6 +321,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
std::vector<float> fCentCut = {10., 20.};
std::vector<float> fNumContribCut = {0, 3000.};
std::vector<float> fNumTrackCut = {0, 350.};
std::vector<float> fOccupancyCut = {0, 6000.};
std::vector<std::vector<std::vector<float>>> fCentAllCorrCut = {{{0., 0., 0., 0.}, {0., 0., 0., 0.}, {0., 0., 0., 0.}}, {{0., 0., 0., 0.}, {0., 0., 0., 0.}, {0., 0., 0., 0.}}, {{0., 0., 0., 0.}, {0., 0., 0., 0.}, {0., 0., 0., 0.}}};
std::vector<std::vector<std::vector<float>>> fMultAllCorrCut = {{{0., 0., 0., 0.}, {0., 0., 0., 0.}, {0., 0., 0., 0.}}, {{0., 0., 0., 0.}, {0., 0., 0., 0.}, {0., 0., 0., 0.}}, {{0., 0., 0., 0.}, {0., 0., 0., 0.}, {0., 0., 0., 0.}}};
std::vector<float> fNumContribMultCorrCut = {0.024, -600., 0.035, 100.};
Expand Down Expand Up @@ -455,6 +459,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
bool bSel8Cut = true;
bool bCentCut = true;
bool bNumContribCut = true;
bool bNumOccupancyCut = true;
bool bCentCorrCut = true;
bool bMultCorrCut = true;
bool bNumContribMultCorrCut = true;
Expand All @@ -481,6 +486,9 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
// *) NumContribCut
bNumContribCut = ebye.fNumContrib < tc.fNumContribCut[1] &&
ebye.fNumContrib > tc.fNumContribCut[0];
// *) NumOccupancyCut
bNumOccupancyCut = collision.trackOccupancyInTimeRange() < tc.fOccupancyCut[1] &&
collision.trackOccupancyInTimeRange() > tc.fOccupancyCut[0];
// *) CentCorrCut
float iCent = 0.;
float jCent = 0.;
Expand Down Expand Up @@ -538,6 +546,9 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
if (tc.fNumContribCutSwitch) {
pass &= bNumContribCut;
}
if (tc.fOccupancyCutSwitch) {
pass &= bNumOccupancyCut;
}
if (tc.fCentCorrCutSwitch) {
pass &= bCentCorrCut;
}
Expand Down Expand Up @@ -1405,11 +1416,13 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
double diffPtMax = (tc.fPtBins[2] - tc.fPtBins[1]) / tc.fPtBins[0];
if (diffPt < diffPtMax) {
wt.fPtMCByRunMap.at(ebye.fRunNumber)->Fill(mcparticle.pt());
} else {
if (tc.fPrintSwitch) {
LOGF(info, "|RecPt - SimPt| = %e > %e", diffPt, diffPtMax);
}
}
// else {
// if (tc.fPrintSwitch) {
// LOGF(info, "|RecPt - SimPt| = %e > %e", diffPt, diffPtMax);
// }
// }

// Fill MC particle histograms after cut:
pc.fParticleHistograms[ePt][eSim][eAfter]->Fill(mcparticle.pt());
pc.fParticleHistograms[ePhi][eSim][eAfter]->Fill(mcparticle.phi());
Expand All @@ -1423,6 +1436,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
ev.fEventHistograms[eNumTrack][eRec][eBefore]->Fill(nTracksBefore);
ev.fEventHistograms[eNumTrack][eRec][eAfter]->Fill(nTracksAfter);
bool bNumTrackCut = nTracksAfter > tc.fNumTrackCut[0] && nTracksAfter < tc.fNumTrackCut[1];
// LOGF(info, "nTracksAfter = %d", nTracksAfter);

// 2-p correlations:
for (int i = 0; i < NumTwoPCorrBins; i++) {
Expand Down Expand Up @@ -1560,7 +1574,6 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
double wFourRecursionAfter = mccRecursion(4, harmonicsFourDen, eAfter).Re();
ebye.fFourParticleCorrelationEbye[eAfter][i] = fourRecursionAfter.Re();

// LOGF(info, "nTracksAfter = %d", nTracksAfter);
if (nTracksAfter > k - 1 && wFourRecursionAfter > 0. && bNumTrackCut) {
// Fill vn^2vm^2 profile and histogram:
mc.fFourParticleCorrelationProfiles[eAfter]->Fill(mc.fFourParticleCorrelationProfiles[eAfter]->GetXaxis()->GetBinCenter(i + 1), ebye.fFourParticleCorrelationEbye[eAfter][i], wFourRecursionAfter);
Expand Down Expand Up @@ -1971,6 +1984,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
tc.fCentCutSwitch = cfCentCutSwitch;
tc.fNumContribCutSwitch = cfNumContribCutSwitch;
tc.fNumTrackCutSwitch = cfNumTrackCutSwitch;
tc.fOccupancyCutSwitch = cfOccupancyCutSwitch;
tc.fCentCorrCutSwitch = cfCentCorrCutSwitch;
tc.fMultCorrCutSwitch = cfMultCorrCutSwitch;
tc.fNumContribMultCorrCutSwitch = cfNumContribMultCorrCutSwitch;
Expand All @@ -1987,6 +2001,7 @@ struct MultiparticleCumulants { // this name is used in lower-case format to nam
tc.fCentCut = cfCentCut;
tc.fNumContribCut = cfNumContribCut;
tc.fNumTrackCut = cfNumTrackCut;
tc.fOccupancyCut = cfOccupancyCut;
tc.fCentAllCorrCut[0][1] = cfCentABCorrCut;
tc.fCentAllCorrCut[0][2] = cfCentACCorrCut;
tc.fCentAllCorrCut[1][2] = cfCentBCCorrCut;
Expand Down
Loading