Skip to content

Commit 003c9ac

Browse files
abylinkinalibuild
andauthored
[PWGUD] Fixed truegap function (#5947)
* Add files via upload * Add files via upload * Please consider the following formatting changes * Update SGTrackSelector.h * Add files via upload * Please consider the following formatting changes * Add files via upload * Please consider the following formatting changes * Update SGSelector.h * Please consider the following formatting changes --------- Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent b19aba5 commit 003c9ac

File tree

1 file changed

+15
-5
lines changed

1 file changed

+15
-5
lines changed

PWGUD/Core/SGSelector.h

Lines changed: 15 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -119,19 +119,29 @@ class SGSelector
119119
float fit_cut[3] = {fv0, ft0a, ft0c};
120120
int gap = collision.gapSide();
121121
int true_gap = gap;
122+
float FV0A, FT0A, FT0C, ZNA, ZNC;
123+
FV0A = collision.totalFV0AmplitudeA();
124+
FT0A = collision.totalFT0AmplitudeA();
125+
FT0C = collision.totalFT0AmplitudeC();
126+
ZNA = collision.energyCommonZNA();
127+
ZNC = collision.energyCommonZNC();
122128
if (gap == 0) {
123-
if (collision.totalFV0AmplitudeA() > fit_cut[0] || collision.totalFT0AmplitudeA() > fit_cut[1] || collision.energyCommonZNA() > zdc_cut)
129+
if (FV0A > fit_cut[0] || FT0A > fit_cut[1] || ZNA > zdc_cut)
124130
true_gap = -1;
125131
} else if (gap == 1) {
126-
if (collision.totalFT0AmplitudeC() > fit_cut[2] || collision.energyCommonZNC() > zdc_cut)
132+
if (FT0C > fit_cut[2] || ZNC > zdc_cut)
127133
true_gap = -1;
128134
} else if (gap == 2) {
129-
if ((collision.totalFV0AmplitudeA() > fit_cut[0] || collision.totalFT0AmplitudeA() > fit_cut[1] || collision.energyCommonZNA() > zdc_cut) && (collision.totalFT0AmplitudeC() > fit_cut[2] || collision.energyCommonZNC() > zdc_cut))
135+
if ((FV0A > fit_cut[0] || FT0A > fit_cut[1] || ZNA > zdc_cut) && (FT0C > fit_cut[2] || ZNC > zdc_cut))
130136
true_gap = -1;
131-
else if (collision.totalFV0AmplitudeA() > fit_cut[0] || collision.totalFT0AmplitudeA() > fit_cut[1] || collision.energyCommonZNA() > zdc_cut)
137+
else if ((FV0A > fit_cut[0] || FT0A > fit_cut[1] || ZNA > zdc_cut) && (FT0C <= fit_cut[2] && ZNC <= zdc_cut))
132138
true_gap = 1;
133-
else if (collision.totalFT0AmplitudeC() > fit_cut[2] || collision.energyCommonZNC() > zdc_cut)
139+
else if ((FV0A <= fit_cut[0] && FT0A <= fit_cut[1] && ZNA <= zdc_cut) && (FT0C > fit_cut[2] || ZNC > zdc_cut))
134140
true_gap = 0;
141+
else if (FV0A <= fit_cut[0] && FT0A <= fit_cut[1] && ZNA <= zdc_cut && FT0C <= fit_cut[2] && ZNC <= zdc_cut)
142+
true_gap = 2;
143+
else
144+
std::cout << "Something wrong with DG" << std::endl;
135145
}
136146
return true_gap;
137147
}

0 commit comments

Comments
 (0)