Skip to content

Commit 8e91e96

Browse files
committed
Add flag to disable ITS based track selection cuts
1 parent 5e3822f commit 8e91e96

File tree

1 file changed

+16
-9
lines changed

1 file changed

+16
-9
lines changed

PWGCF/Flow/Tasks/flowPtEfficiency.cxx

Lines changed: 16 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -60,6 +60,7 @@ struct FlowPtEfficiency {
6060
O2_DEFINE_CONFIGURABLE(cfgCutDCAz, float, 2.0f, "DCAz cut for tracks")
6161
O2_DEFINE_CONFIGURABLE(cfgCutDCAxyppPass3Enabled, bool, false, "switch of ppPass3 DCAxy pt dependent cut")
6262
O2_DEFINE_CONFIGURABLE(cfgCutDCAzPtDepEnabled, bool, false, "switch of DCAz pt dependent cut")
63+
O2_DEFINE_CONFIGURABLE(cfgEnableITSCuts, bool, true, "switch of enabling ITS based track selection cuts")
6364
O2_DEFINE_CONFIGURABLE(cfgSelRunNumberEnabled, bool, false, "switch of run number selection")
6465
O2_DEFINE_CONFIGURABLE(cfgFlowEnabled, bool, false, "switch of calculating flow")
6566
O2_DEFINE_CONFIGURABLE(cfgFlowNbootstrap, int, 30, "Number of subsamples")
@@ -219,21 +220,25 @@ struct FlowPtEfficiency {
219220
fGFWReco->CreateRegions();
220221
}
221222

222-
if (cfgTrkSelRun3ITSMatch) {
223-
myTrackSel = getGlobalTrackSelectionRun3ITSMatch(TrackSelection::GlobalTrackRun3ITSMatching::Run3ITSall7Layers, TrackSelection::GlobalTrackRun3DCAxyCut::Default);
224-
} else {
225-
myTrackSel = getGlobalTrackSelectionRun3ITSMatch(TrackSelection::GlobalTrackRun3ITSMatching::Run3ITSibAny, TrackSelection::GlobalTrackRun3DCAxyCut::Default);
223+
if(cfgEnableITSCuts){
224+
if (cfgTrkSelRun3ITSMatch) {
225+
myTrackSel = getGlobalTrackSelectionRun3ITSMatch(TrackSelection::GlobalTrackRun3ITSMatching::Run3ITSall7Layers, TrackSelection::GlobalTrackRun3DCAxyCut::Default);
226+
} else {
227+
myTrackSel = getGlobalTrackSelectionRun3ITSMatch(TrackSelection::GlobalTrackRun3ITSMatching::Run3ITSibAny, TrackSelection::GlobalTrackRun3DCAxyCut::Default);
228+
}
226229
}
227230
if (cfgCutDCAxyppPass3Enabled) {
228231
myTrackSel.SetMaxDcaXYPtDep([](float pt) { return 0.004f + 0.013f / pt; });
229232
} else {
230233
myTrackSel.SetMaxDcaXY(cfgCutDCAxy);
231234
}
232235
myTrackSel.SetMinNClustersTPC(cfgCutTPCclu);
233-
myTrackSel.SetMinNClustersITS(cfgCutITSclu);
234236
myTrackSel.SetMinNCrossedRowsTPC(cfgCutTPCcrossedrows);
237+
if(cfgEnableITSCuts) myTrackSel.SetMinNClustersITS(cfgCutITSclu);
235238
if (!cfgCutDCAzPtDepEnabled)
236239
myTrackSel.SetMaxDcaZ(cfgCutDCAz);
240+
241+
237242
}
238243

239244
template <char... chars>
@@ -358,10 +363,14 @@ struct FlowPtEfficiency {
358363
template <typename TTrack>
359364
bool trackSelected(TTrack track)
360365
{
361-
if (cfgkIsTrackGlobal && !track.isGlobalTrack())
366+
if (cfgkIsTrackGlobal && !track.isGlobalTrack()){
367+
LOGF(info,"Track is not global");
362368
return false;
363-
if (cfgCutDCAzPtDepEnabled && (track.dcaZ() > (0.004f + 0.013f / track.pt())))
369+
}
370+
if (cfgCutDCAzPtDepEnabled && (track.dcaZ() > (0.004f + 0.013f / track.pt()))){
371+
LOGF(info,"pt-dep DCAz not valid");
364372
return false;
373+
}
365374
return myTrackSel.IsSelected(track);
366375
}
367376

@@ -378,7 +387,6 @@ struct FlowPtEfficiency {
378387
if (!std::count(cfgRunNumberList.value.begin(), cfgRunNumberList.value.end(), runNumber))
379388
return;
380389
}
381-
382390
float imp = 0;
383391
bool impFetched = false;
384392
float evPhi = 0;
@@ -393,7 +401,6 @@ struct FlowPtEfficiency {
393401

394402
fGFWReco->Clear();
395403
}
396-
397404
for (const auto& track : tracks) {
398405
if (!trackSelected(track))
399406
continue;

0 commit comments

Comments
 (0)