Skip to content

Commit 1be224f

Browse files
authored
[PWGCF] fix dpt selection bug (#13756)
1 parent 286e1c9 commit 1be224f

File tree

1 file changed

+5
-4
lines changed

1 file changed

+5
-4
lines changed

PWGCF/Flow/Tasks/flowTask.cxx

Lines changed: 5 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -156,7 +156,7 @@ struct FlowTask {
156156
TF1* fPhiCutLow = nullptr;
157157
TF1* fPhiCutHigh = nullptr;
158158
// for deltaPt/<pT> vs centrality
159-
O2_DEFINE_CONFIGURABLE(cfgDptDisEnable, bool, false, "Produce deltaPt/meanPt vs c22 vs centrality")
159+
O2_DEFINE_CONFIGURABLE(cfgDptDisEnable, bool, false, "Produce deltaPt/meanPt vs centrality")
160160
O2_DEFINE_CONFIGURABLE(cfgDptDisSelectionSwitch, int, 0, "0: disable, 1: use low cut, 2:use high cut")
161161
TH1D* hEvAvgMeanPt = nullptr;
162162
TH1D* fDptDisCutLow = nullptr;
@@ -1108,11 +1108,12 @@ struct FlowTask {
11081108
if (cfgFuncParas.cfgDptDisEnable) {
11091109
double meanPt = ptSum / weffEvent;
11101110
double deltaPt = meanPt - cfgFuncParas.hEvAvgMeanPt->GetBinContent(cfgFuncParas.hEvAvgMeanPt->FindBin(independent));
1111-
registry.fill(HIST("hNormDeltaPt_X"), meanPt, independent, weffEvent);
1112-
if (cfgFuncParas.cfgDptDisSelectionSwitch == kLowDptCut && deltaPt > cfgFuncParas.fDptDisCutLow->GetBinContent(cfgFuncParas.fDptDisCutLow->FindBin(independent))) {
1111+
double normDeltaPt = deltaPt / meanPt;
1112+
registry.fill(HIST("hNormDeltaPt_X"), normDeltaPt, independent, weffEvent);
1113+
if (cfgFuncParas.cfgDptDisSelectionSwitch == kLowDptCut && normDeltaPt > cfgFuncParas.fDptDisCutLow->GetBinContent(cfgFuncParas.fDptDisCutLow->FindBin(independent))) {
11131114
// only keep low 10% dpt event
11141115
return;
1115-
} else if (cfgFuncParas.cfgDptDisSelectionSwitch == kHighDptCut && deltaPt < cfgFuncParas.fDptDisCutHigh->GetBinContent(cfgFuncParas.fDptDisCutHigh->FindBin(independent))) {
1116+
} else if (cfgFuncParas.cfgDptDisSelectionSwitch == kHighDptCut && normDeltaPt < cfgFuncParas.fDptDisCutHigh->GetBinContent(cfgFuncParas.fDptDisCutHigh->FindBin(independent))) {
11161117
// only keep high 10% dpt event
11171118
return;
11181119
}

0 commit comments

Comments
 (0)