Skip to content

[PWGCF] Add MultTPC vs cent cut in flowGfwOmegaXi.cxx #16591

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
wenyaCern merged 1 commit into from
Jun 9, 2026
+14 −5
Merged

[PWGCF] Add MultTPC vs cent cut in flowGfwOmegaXi.cxx #16591

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
19 changes: 14 additions & 5 deletions PWGCF/Flow/Tasks/flowGfwOmegaXi.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -187,7 +187,8 @@ struct FlowGfwOmegaXi {
O2_DEFINE_CONFIGURABLE(cfgDoIsGoodITSLayersAll, bool, true, "check kIsGoodITSLayersAll")
O2_DEFINE_CONFIGURABLE(cfgCutOccupancyHigh, int, 500, "High cut on TPC occupancy")
O2_DEFINE_CONFIGURABLE(cfgDoZResAndNumContribCut, bool, true, "check kNoTimeFrameBorder")
O2_DEFINE_CONFIGURABLE(cfgDoMultPVCut, bool, true, "do multNTracksPV vs cent cut")
O2_DEFINE_CONFIGURABLE(cfgDoMultPVCut, bool, false, "do multNTracksPV vs cent cut")
O2_DEFINE_CONFIGURABLE(cfgDoMultTPCCut, bool, false, "do multNTracksTPC vs cent cut")
O2_DEFINE_CONFIGURABLE(cfgMultPVCut, std::vector<float>, (std::vector<float>{3074.43, -106.192, 1.46176, -0.00968364, 2.61923e-05, 182.128, -7.43492, 0.193901, -0.00256715, 1.22594e-05}), "Used MultPVCut function parameter")
O2_DEFINE_CONFIGURABLE(cfgDoV0AT0Acut, bool, true, "do V0A-T0A cut")
O2_DEFINE_CONFIGURABLE(cfgCutminIR, float, -1, "cut min IR")
Expand Down Expand Up @@ -396,7 +397,7 @@ struct FlowGfwOmegaXi {
registry.add("hEta", "", {HistType::kTH1D, {cfgaxisEta}});
registry.add("hVtxZ", "", {HistType::kTH1D, {cfgaxisVertex}});
registry.add("hMult", "", {HistType::kTH1D, {cfgaxisNch}});
registry.add("hMultTPC", "", {HistType::kTH1D, {cfgaxisNch}});
registry.add("hCentvsMultTPC", "", {HistType::kTH2D, {{100, 0, 100}, cfgaxisNch}});
registry.add("hCent", "", {HistType::kTH1D, {{90, 0, 90}}});
registry.add("hNTracksPVvsCentrality", "", {HistType::kTH2D, {{5000, 0, 5000}, {100, 0, 100}}});
registry.add("hmultFV0AvsmultFT0A", "", {HistType::kTH2D, {{5000, 0, 5000}, {5000, 0, 5000}}});
Expand Down Expand Up @@ -449,7 +450,7 @@ struct FlowGfwOmegaXi {
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(7, "after kIsGoodZvtxFT0vsPV");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(8, "after kNoCollInTimeRangeStandard");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(9, "after kIsGoodITSLayersAll");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(10, "after MultPVCut");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(10, "after MultvsCentCut");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(11, "after TPC occupancy cut");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(12, "after V0AT0Acut");
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(13, "after IRmincut");
Expand Down Expand Up @@ -1062,6 +1063,15 @@ struct FlowGfwOmegaXi {
return false;
}

int nMultTPC = collision.multTPC();
registry.fill(HIST("hCentvsMultTPC"), centrality, nMultTPC);
if (evtSeleOpts.cfgDoMultTPCCut.value) {
if (nMultTPC < fMultPVCutLow->Eval(centrality))
return false;
if (nMultTPC > fMultPVCutLow->Eval(centrality))
return false;
}

registry.fill(HIST("hEventCount"), 9.5);

if (occupancy > evtSeleOpts.cfgCutOccupancyHigh.value)
Expand Down Expand Up @@ -1092,7 +1102,7 @@ struct FlowGfwOmegaXi {
{
o2::aod::ITSResponse itsResponse;
int nTot = tracks.size();
float nMultTPC = collision.multTPC();
int nMultTPC = collision.multTPC();
auto bc = collision.bc_as<aod::BCsWithTimestamps>();
int runNumber = bc.runNumber();
double interactionRate = rateFetcher.fetch(ccdb.service, bc.timestamp(), runNumber, "ZNC hadronic") * 1.e-3;
Expand Down Expand Up @@ -1128,7 +1138,6 @@ struct FlowGfwOmegaXi {
float vtxz = collision.posZ();
registry.fill(HIST("hVtxZ"), vtxz);
registry.fill(HIST("hMult"), nTot);
registry.fill(HIST("hMultTPC"), nMultTPC);
registry.fill(HIST("hCent"), cent);

float weff = 1;
Expand Down
Loading