Skip to content

[PWGCF] Add pT-dependent DCAz cut #16195

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 4 commits into from
May 11, 2026
+10 −3
Merged

[PWGCF] Add pT-dependent DCAz cut #16195

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
715624a
Commit after linter issue solution
May 11, 2026
93cd6e9
Fix linter issues
May 11, 2026
12ede94
Please consider the following formatting changes
alibuild May 11, 2026
5d6734f
Merge pull request #3 from alibuild/alibot-cleanup-16195
vojmach May 11, 2026
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
13 changes: 10 additions & 3 deletions PWGCF/Flow/Tasks/flowMc.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -89,6 +89,8 @@ struct FlowMc {
O2_DEFINE_CONFIGURABLE(cfgDCAxyFunction, std::string, "(0.0015+0.005/(x^1.1))", "Functional form of pt-dependent DCAxy cut");
O2_DEFINE_CONFIGURABLE(cfgCutDCAz, float, 2.0f, "DCAz cut for tracks")
O2_DEFINE_CONFIGURABLE(cfgCutDCAzPtDepEnabled, bool, false, "switch of DCAz pt dependent cut")
O2_DEFINE_CONFIGURABLE(cfgDCAzFunction, std::string, "(0.0015+0.005/(x^1.1))", "Functional form of pt-dependent DCAz cut");
O2_DEFINE_CONFIGURABLE(cfgDCAzNSigma, float, 7, "Cut on number of sigma deviations from expected DCA in the z direction");
O2_DEFINE_CONFIGURABLE(cfgEnableITSCuts, bool, true, "switch of enabling ITS based track selection cuts")
O2_DEFINE_CONFIGURABLE(cfgTrkSelRun3ITSMatch, bool, false, "GlobalTrackRun3ITSMatching::Run3ITSall7Layers selection")
O2_DEFINE_CONFIGURABLE(cfgFlowAcceptance, std::string, "", "CCDB path to acceptance object")
Expand Down Expand Up @@ -137,6 +139,7 @@ struct FlowMc {
// Additional filters for tracks
TrackSelection myTrackSel;
TF1* fPtDepDCAxy = nullptr;
TF1* fPtDepDCAz = nullptr;

// Cent vs IP
TH1D* mCentVsIPTruth = nullptr;
Expand Down Expand Up @@ -310,9 +313,13 @@ struct FlowMc {
myTrackSel.SetMinNCrossedRowsTPC(cfgCutTPCcrossedrows);
if (cfgEnableITSCuts)
myTrackSel.SetMinNClustersITS(cfgCutITSclu);
if (!cfgCutDCAzPtDepEnabled)
if (!cfgCutDCAzPtDepEnabled) {
myTrackSel.SetMaxDcaZ(cfgCutDCAz);

} else {
fPtDepDCAz = new TF1("ptDepDCAxy", Form("[0]*%s", cfgDCAzFunction->c_str()), 0.001, 100);
fPtDepDCAz->SetParameter(0, cfgDCAzNSigma);
LOGF(info, "DCAz pt-dependence function: %s", Form("[0]*%s", cfgDCAzFunction->c_str()));
}
if (cfgTrackDensityCorrUse) {
std::vector<double> pTEffBins = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.4, 1.8, 2.2, 2.6, 3.0};
hFindPtBin = new TH1D("hFindPtBin", "hFindPtBin", pTEffBins.size() - 1, &pTEffBins[0]);
Expand Down Expand Up @@ -472,7 +479,7 @@ struct FlowMc {
template <typename TTrack>
bool trackSelected(TTrack track)
{
if (cfgCutDCAzPtDepEnabled && (track.dcaZ() > (0.004f + 0.013f / track.pt()))) {
if (cfgCutDCAzPtDepEnabled && (std::fabs(track.dcaZ()) > fPtDepDCAz->Eval(track.pt()))) {
return false;
}
return myTrackSel.IsSelected(track);
Expand Down
Loading