Skip to content

Commit d056c27

Browse files
sdudi123sandeep dudi
sdudi123
and
sandeep dudi
authored
[PWGUD] RCT flag added (#12194)
Co-authored-by: sandeep dudi <sandeep.dudi@cern.ch>
1 parent 981fbd5 commit d056c27

File tree

1 file changed

+138
-24
lines changed

1 file changed

+138
-24
lines changed

‎PWGUD/Tasks/sginclusivePhiKstarSD.cxx‎

Lines changed: 138 additions & 24 deletions
Original file line numberDiff line numberDiff line change
@@ -51,6 +51,7 @@ struct SginclusivePhiKstarSD {
5151

5252
HistogramRegistry registry{"registry", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
5353

54+
Configurable<int> cutRCTflag{"cutRCTflag", 0, {"0 = off, 1 = CBT, 2 = CBT+ZDC, 3 = CBThadron, 4 = CBThadron+ZDC"}};
5455
Configurable<float> fv0Cut{"fv0Cut", 50., "FV0A threshold"};
5556
Configurable<float> ft0aCut{"ft0aCut", 100., "FT0A threshold"};
5657
Configurable<float> ft0cCut{"ft0cCut", 50., "FT0C threshold"};
@@ -68,6 +69,7 @@ struct SginclusivePhiKstarSD {
6869
Configurable<int> useSbp{"useSbp", -1, "kNoSameBunchPileup cut"};
6970
Configurable<int> useZvtxftovpv{"useZvtxftovpv", -1, "kIsGoodZvtxFT0vsPV cut"};
7071
Configurable<int> useVtxItsTpc{"useVtxItsTpc", -1, "kIsVertexITSTPC cut"};
72+
Configurable<int> upcflag{"upcflag", -1, "upc run selection, 0 = std, 1= upc"};
7173

7274
// Track Selections
7375
Configurable<float> pvCut{"pvCut", 1.0, "Use Only PV tracks"};
@@ -89,7 +91,6 @@ struct SginclusivePhiKstarSD {
8991
Configurable<float> nsigmaTpcCut1{"nsigmaTpcCut1", 3.0, "nsigma tpc cut1"};
9092
Configurable<float> nsigmaTpcCut2{"nsigmaTpcCut2", 3.0, "nsigma tpc cut2"};
9193
Configurable<float> nsigmaTpcCut3{"nsigmaTpcCut3", 3.0, "nsigma tpc cut3"};
92-
Configurable<float> nsigmaTpcCut4{"nsigmaTpcCut4", 3.0, "nsigma tpc cut4"};
9394
Configurable<float> nsigmaTpcCut{"nsigmaTpcCut", 3.0, "nsigma tpc cut"};
9495
Configurable<float> nsigmaTofCut{"nsigmaTofCut", 9.0, "nsigma tpc+tof cut"};
9596
Configurable<float> nsigmaTofCut1{"nsigmaTofCut1", 3.0, "nsigma tof cut"};
@@ -120,8 +121,24 @@ struct SginclusivePhiKstarSD {
120121
Configurable<bool> reconstruction{"reconstruction", true, ""};
121122
Configurable<int> generatedId{"generatedId", 31, ""};
122123

124+
// Configurable axes for histogram
125+
ConfigurableAxis dcaAxisConfig{"dcaAxisConfig", {600, -0.3f, 0.3f}, "DCAxy & DCAz axis"};
126+
ConfigurableAxis etaAxisConfig{"etaAxisConfig", {400, -1.0f, 1.0f}, "Pseudorapidity & Rapidity axis"};
127+
ConfigurableAxis VrtxXAxisConfig{"VrtxXAxisConfig", {400, -0.1f, 0.1f}, "Vertex X axis"};
128+
ConfigurableAxis VrtxYAxisConfig{"VrtxYAxisConfig", {200, -0.05f, 0.05f}, "Vertex Y axis"};
129+
ConfigurableAxis VrtxZAxisConfig{"VrtxZAxisConfig", {600, -15.0f, 15.0f}, "Vertex Z axis"};
130+
123131
void init(InitContext const& context)
124132
{
133+
// Axes
134+
AxisSpec dcaxyAxis = {dcaAxisConfig, "DCAxy (cm)"};
135+
AxisSpec dcazAxis = {dcaAxisConfig, "DCAz (cm)"};
136+
AxisSpec etaAxis = {etaAxisConfig, "#eta"};
137+
AxisSpec rapAxis = {etaAxisConfig, "y"};
138+
AxisSpec VrtxXAxis = {VrtxXAxisConfig, "Vertex X (cm)"};
139+
AxisSpec VrtxYAxis = {VrtxYAxisConfig, "Vertex Y (cm)"};
140+
AxisSpec VrtxZAxis = {VrtxZAxisConfig, "Vertex Z (cm)"};
141+
125142
registry.add("GapSide", "Gap Side; Entries", kTH1F, {{4, -1.5, 2.5}});
126143
registry.add("TrueGapSide", "Gap Side; Entries", kTH1F, {{4, -1.5, 2.5}});
127144
registry.add("nPVContributors_data", "Multiplicity_dist_before track cut gap A", kTH1F, {{110, 0, 110}});
@@ -130,56 +147,56 @@ struct SginclusivePhiKstarSD {
130147
if (phi) {
131148
registry.add("os_KK_pT_0", "pt kaon pair", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
132149
registry.add("os_KK_pT_1", "pt kaon pair", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
133-
registry.add("os_KK_pT_2", "pt kaon pair", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
150+
registry.add("os_KK_pT_2", "pt kaon pair", kTH3F, {{305, 0.98, 2.2}, {80, -2.0, 2.0}, {100, 0, 10}});
134151
registry.add("os_KK_ls_pT_0", "kaon pair like sign", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
135152
registry.add("os_KK_ls_pT_1", "kaon pair like sign", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
136-
registry.add("os_KK_ls_pT_2", "kaon pair like sign", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
153+
registry.add("os_KK_ls_pT_2", "kaon pair like sign", kTH3F, {{305, 0.98, 2.2}, {80, -2.0, 2.0}, {100, 0, 10}});
137154

138155
registry.add("os_KK_mix_pT_0", "kaon pair mix event", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
139156
registry.add("os_KK_mix_pT_1", "kaon pair mix event", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
140-
registry.add("os_KK_mix_pT_2", "kaon pair mix event", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
157+
registry.add("os_KK_mix_pT_2", "kaon pair mix event", kTH3F, {{305, 0.98, 2.2}, {80, -2.0, 2.0}, {100, 0, 10}});
141158

142159
registry.add("os_KK_rot_pT_0", "kaon pair mix event", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
143160
registry.add("os_KK_rot_pT_1", "kaon pair mix event", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
144-
registry.add("os_KK_rot_pT_2", "kaon pair mix event", kTH3F, {{220, 0.98, 1.2}, {80, -2.0, 2.0}, {100, 0, 10}});
161+
registry.add("os_KK_rot_pT_2", "kaon pair mix event", kTH3F, {{305, 0.98, 2.2}, {80, -2.0, 2.0}, {100, 0, 10}});
145162
}
146163
if (rho) {
147-
registry.add("os_pp_pT_0", "pt pion pair", kTH3F, {{120, 1.44, 2.04}, {80, -2.0, 2.0}, {100, 0, 10}});
148-
registry.add("os_pp_pT_1", "pt pion pair", kTH3F, {{120, 1.44, 2.04}, {80, -2.0, 2.0}, {100, 0, 10}});
149-
registry.add("os_pp_pT_2", "pt pion pair", kTH3F, {{120, 1.44, 2.04}, {80, -2.0, 2.0}, {100, 0, 10}});
150-
registry.add("os_pp_ls_pT_0", "pion pair like sign", kTH3F, {{120, 1.44, 2.04}, {80, -2.0, 2.0}, {100, 0, 10}});
151-
registry.add("os_pp_ls_pT_1", "pion pair like sign", kTH3F, {{120, 1.44, 2.04}, {80, -2.0, 2.0}, {100, 0, 10}});
152-
registry.add("os_pp_ls_pT_2", "pion pair like sign", kTH3F, {{120, 1.44, 2.04}, {80, -2.0, 2.0}, {100, 0, 10}});
164+
registry.add("os_pp_pT_0", "pt pion pair", kTH3F, {{200, 1.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
165+
registry.add("os_pp_pT_1", "pt pion pair", kTH3F, {{200, 1.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
166+
registry.add("os_pp_pT_2", "pt pion pair", kTH3F, {{200, 1.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
167+
registry.add("os_pp_ls_pT_0", "pion pair like sign", kTH3F, {{200, 1.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
168+
registry.add("os_pp_ls_pT_1", "pion pair like sign", kTH3F, {{200, 1.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
169+
registry.add("os_pp_ls_pT_2", "pion pair like sign", kTH3F, {{200, 1.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
153170
}
154171
if (kstar) {
155172
registry.add("os_pk_pT_0", "pion-kaon pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
156173
registry.add("os_pk_pT_1", "pion-kaon pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
157-
registry.add("os_pk_pT_2", "pion-kaon pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
174+
registry.add("os_pk_pT_2", "pion-kaon pair", kTH3F, {{600, 0.0, 3.0}, {80, -2.0, 2.0}, {1000, 0, 10}});
158175

159176
registry.add("os_pk_mix_pT_0", "pion-kaon mix pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
160177
registry.add("os_pk_mix_pT_1", "pion-kaon mix pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
161-
registry.add("os_pk_mix_pT_2", "pion-kaon mix pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
178+
registry.add("os_pk_mix_pT_2", "pion-kaon mix pair", kTH3F, {{600, 0.0, 3.0}, {80, -2.0, 2.0}, {1000, 0, 10}});
162179

163180
registry.add("os_pk_rot_pT_0", "pion-kaon rotional pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
164181
registry.add("os_pk_rot_pT_1", "pion-kaon rotional pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
165-
registry.add("os_pk_rot_pT_2", "pion-kaon rotional pair", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
182+
registry.add("os_pk_rot_pT_2", "pion-kaon rotional pair", kTH3F, {{600, 0.0, 3.0}, {80, -2.0, 2.0}, {1000, 0, 10}});
166183

167184
registry.add("os_pk_ls_pT_0", "pion-kaon pair like sign", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
168185
registry.add("os_pk_ls_pT_1", "pion-kaon like sign", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
169-
registry.add("os_pk_ls_pT_2", "pion-kaon like sign", kTH3F, {{400, 0.0, 2.0}, {80, -2.0, 2.0}, {100, 0, 10}});
186+
registry.add("os_pk_ls_pT_2", "pion-kaon like sign", kTH3F, {{600, 0.0, 3.0}, {80, -2.0, 2.0}, {1000, 0, 10}});
170187

171188
registry.add("hRotation", "hRotation", kTH1F, {{360, 0.0, o2::constants::math::TwoPI}});
172189
}
173190
// qa plots
174191
if (qa) {
175-
registry.add("tpc_dedx", "p vs dE/dx", kTH2F, {{100, 0.0, 10.0}, {10000, 0.0, 2500.0}});
176-
registry.add("tof_beta", "p vs beta", kTH2F, {{100, 0.0, 10.0}, {100, 0.0, 5.0}});
177-
178-
registry.add("tpc_dedx_kaon", "p#k dE/dx", kTH2F, {{100, 0.0, 10.0}, {10000, 0.0, 2500.0}});
179-
registry.add("tpc_dedx_pion", "p#pi dE/dx", kTH2F, {{100, 0.0, 10.0}, {10000, 0.0, 2500.0}});
180-
registry.add("tpc_dedx_kaon_1", "tpc+tof pid cut p#k dE/dx", kTH2F, {{100, 0.0, 10.0}, {10000, 0.0, 2500.0}});
181-
registry.add("tpc_dedx_kaon_2", "tpc+tof pid cut1 p#k dE/dx", kTH2F, {{100, 0.0, 10.0}, {10000, 0.0, 2500.0}});
182-
registry.add("tpc_dedx_pion_1", "tpc+tof pid cut p#pi dE/dx", kTH2F, {{100, 0.0, 10.0}, {10000, 0.0, 25000.0}});
192+
registry.add("tpc_dedx", "p vs dE/dx", kTH2F, {{500, 0.0, 10.0}, {5000, 0.0, 5000.0}});
193+
registry.add("tof_beta", "p vs beta", kTH2F, {{500, 0.0, 10.0}, {500, 0.0, 1.0}});
194+
195+
registry.add("tpc_dedx_kaon", "p#k dE/dx", kTH2F, {{500, 0.0, 10.0}, {5000, 0.0, 5000.0}});
196+
registry.add("tpc_dedx_pion", "p#pi dE/dx", kTH2F, {{500, 0.0, 10.0}, {5000, 0.0, 5000.0}});
197+
registry.add("tpc_dedx_kaon_1", "tpc+tof pid cut p#k dE/dx", kTH2F, {{500, 0.0, 10.0}, {5000, 0.0, 5000.0}});
198+
registry.add("tpc_dedx_kaon_2", "tpc+tof pid cut1 p#k dE/dx", kTH2F, {{500, 0.0, 10.0}, {5000, 0.0, 5000.0}});
199+
registry.add("tpc_dedx_pion_1", "tpc+tof pid cut p#pi dE/dx", kTH2F, {{500, 0.0, 10.0}, {5000, 0.0, 5000.0}});
183200
registry.add("tpc_nsigma_kaon", "p#k n#sigma", kTH2F, {{100, 0.0, 10.0}, {100, -10.0, 10.0}});
184201
registry.add("tpc_nsigma_pion", "p#pi n#sigma", kTH2F, {{100, 0.0, 10.0}, {100, -10.0, 10.0}});
185202
registry.add("tpc_tof_nsigma_kaon", "p#k n#sigma TPC vs TOF", kTH2F, {{100, -10.0, 10.0}, {100, -10.0, 10.0}});
@@ -208,6 +225,31 @@ struct SginclusivePhiKstarSD {
208225
registry.add("V0A_0", "V0A amplitude", kTH1F, {{1000, 0.0, 1000.0}});
209226
registry.add("V0A_1", "V0A amplitude", kTH1F, {{1000, 0.0, 1000.0}});
210227

228+
registry.add("hDcaxy_all_before", "DCAxy Distribution of all tracks before track selection; DCAxy (cm); Counts", kTH1F, {dcaxyAxis});
229+
registry.add("hDcaz_all_before", "DCAz Distribution of all tracks before track selection; DCAz (cm); Counts", kTH1F, {dcazAxis});
230+
231+
registry.add("hDcaxy_all_after", "DCAxy Distribution of all tracks after track selection; DCAxy (cm); Counts", kTH1F, {dcaxyAxis});
232+
registry.add("hDcaz_all_after", "DCAz Distribution of all tracks after track selection; DCAz (cm); Counts", kTH1F, {dcazAxis});
233+
234+
registry.add("hDcaxy_pi", "DCAxy Distribution of selected pions; DCAxy (cm); Counts", kTH1F, {dcaxyAxis});
235+
registry.add("hDcaz_pi", "DCAz Distribution of selected pions; DCAz (cm); Counts", kTH1F, {dcazAxis});
236+
237+
registry.add("hDcaxy_ka", "DCAxy Distribution of selected kaons; DCAxy (cm); Counts", kTH1F, {dcaxyAxis});
238+
registry.add("hDcaz_ka", "DCAz Distribution of selected kaons; DCAz (cm); Counts", kTH1F, {dcazAxis});
239+
240+
registry.add("hVertexX", "Vertex X distribution; Vertex X (cm); Counts", kTH1F, {VrtxXAxis});
241+
registry.add("hVertexY", "Vertex Y distribution; Vertex Y (cm); Counts", kTH1F, {VrtxYAxis});
242+
registry.add("hVertexZ", "VertexZ distribution; Vertex Z (cm); Counts", kTH1F, {VrtxZAxis});
243+
244+
registry.add("hEta_all_after", "Pseudorapidity of all tracks after track selection; #eta; Counts", kTH1F, {etaAxis});
245+
registry.add("hRap_all_after", "Rapidity of all tracks after track selection; y; Counts", kTH1F, {rapAxis});
246+
247+
registry.add("hEta_pi", "Pseudorapidity of selected Pions; #eta; Counts", kTH1F, {etaAxis});
248+
registry.add("hRap_pi", "Rapidity of selected Pions; y; Counts", kTH1F, {rapAxis});
249+
250+
registry.add("hEta_ka", "Pseudorapidity of selected Kaons; #eta; Counts", kTH1F, {etaAxis});
251+
registry.add("hRap_ka", "Rapidity of selected Kaons; y; Counts", kTH1F, {rapAxis});
252+
211253
if (rapidityGap) {
212254
registry.add("mult_0", "mult0", kTH1F, {{150, 0, 150}});
213255
registry.add("mult_1", "mult1", kTH1F, {{150, 0, 150}});
@@ -432,6 +474,23 @@ struct SginclusivePhiKstarSD {
432474
return cosThetaCs;
433475
}
434476

477+
template <typename C>
478+
bool isGoodRCTflag(C const& coll)
479+
{
480+
switch (cutRCTflag) {
481+
case 1:
482+
return sgSelector.isCBTOk(coll);
483+
case 2:
484+
return sgSelector.isCBTZdcOk(coll);
485+
case 3:
486+
return sgSelector.isCBTHadronOk(coll);
487+
case 4:
488+
return sgSelector.isCBTHadronZdcOk(coll);
489+
default:
490+
return true;
491+
}
492+
}
493+
435494
template <typename T>
436495
bool selectionPIDKaon1(const T& candidate)
437496
{
@@ -575,7 +634,10 @@ struct SginclusivePhiKstarSD {
575634
return;
576635
if (useVtxItsTpc != -1 && collision.vtxITSTPC() != useVtxItsTpc)
577636
return;
578-
637+
if (!isGoodRCTflag(collision))
638+
return;
639+
if (upcflag != -1 && collision.flags() != upcflag)
640+
return;
579641
int mult = collision.numContrib();
580642
if (gapSide == 0) {
581643
registry.fill(HIST("gap_mult0"), mult);
@@ -597,6 +659,12 @@ struct SginclusivePhiKstarSD {
597659
int trackextra = 0;
598660
int trackextraDG = 0;
599661

662+
if (qa) {
663+
registry.fill(HIST("hVertexX"), collision.posX());
664+
registry.fill(HIST("hVertexY"), collision.posY());
665+
registry.fill(HIST("hVertexZ"), collision.posZ());
666+
}
667+
600668
/* Partition<UDtracksfull> pvContributors1 = aod::udtrack::isPVContributor == true;
601669
pvContributors1.bindTable(tracks);
602670
if (gapSide == 0) {
@@ -607,9 +675,24 @@ struct SginclusivePhiKstarSD {
607675
}
608676
*/
609677
for (const auto& track1 : tracks) {
678+
679+
if (qa) {
680+
registry.fill(HIST("hDcaxy_all_before"), track1.dcaXY());
681+
registry.fill(HIST("hDcaz_all_before"), track1.dcaZ());
682+
}
683+
610684
if (!trackselector(track1, parameters))
611685
continue;
686+
612687
v0.SetCoordinates(track1.px(), track1.py(), track1.pz(), o2::constants::physics::MassPionCharged);
688+
689+
if (qa) {
690+
registry.fill(HIST("hDcaxy_all_after"), track1.dcaXY());
691+
registry.fill(HIST("hDcaz_all_after"), track1.dcaZ());
692+
registry.fill(HIST("hEta_all_after"), v0.Eta());
693+
registry.fill(HIST("hRap_all_after"), v0.Rapidity());
694+
}
695+
613696
if (selectionPIDPion1(track1)) {
614697
onlyPionTrackspm.push_back(v0);
615698
rawPionTrackspm.push_back(track1);
@@ -660,13 +743,21 @@ struct SginclusivePhiKstarSD {
660743
registry.fill(HIST("tpc_nsigma_kaon"), v0.Pt(), track1.tpcNSigmaKa());
661744
registry.fill(HIST("tof_nsigma_kaon"), v0.Pt(), track1.tofNSigmaKa());
662745
registry.fill(HIST("tpc_tof_nsigma_kaon"), track1.tpcNSigmaKa(), track1.tofNSigmaKa());
746+
registry.fill(HIST("hEta_ka"), v0.Eta());
747+
registry.fill(HIST("hRap_ka"), v0.Rapidity());
748+
registry.fill(HIST("hDcaxy_ka"), track1.dcaXY());
749+
registry.fill(HIST("hDcaz_ka"), track1.dcaZ());
663750
}
664751

665752
if (selectionPIDPion1(track1)) {
666753
registry.fill(HIST("tpc_dedx_pion_1"), v0.P(), track1.tpcSignal());
667754
registry.fill(HIST("tpc_nsigma_pion"), v0.Pt(), track1.tpcNSigmaPi());
668755
registry.fill(HIST("tof_nsigma_pion"), v0.Pt(), track1.tofNSigmaPi());
669756
registry.fill(HIST("tpc_tof_nsigma_pion"), track1.tpcNSigmaPi(), track1.tofNSigmaPi());
757+
registry.fill(HIST("hEta_pi"), v0.Eta());
758+
registry.fill(HIST("hRap_pi"), v0.Rapidity());
759+
registry.fill(HIST("hDcaxy_pi"), track1.dcaXY());
760+
registry.fill(HIST("hDcaz_pi"), track1.dcaZ());
670761
}
671762
}
672763
}
@@ -1384,6 +1475,29 @@ struct SginclusivePhiKstarSD {
13841475
return;
13851476
if (std::abs(collision.occupancyInTime()) > occCut)
13861477
return;
1478+
if (std::abs(collision.hadronicRate()) > hadronicRate)
1479+
return;
1480+
if (useTrs != -1 && collision.trs() != useTrs)
1481+
return;
1482+
if (useTrofs != -1 && collision.trofs() != useTrofs)
1483+
return;
1484+
if (useHmpr != -1 && collision.hmpr() != useHmpr)
1485+
return;
1486+
if (useTfb != -1 && collision.tfb() != useTfb)
1487+
return;
1488+
if (useItsrofb != -1 && collision.itsROFb() != useItsrofb)
1489+
return;
1490+
if (useSbp != -1 && collision.sbp() != useSbp)
1491+
return;
1492+
if (useZvtxftovpv != -1 && collision.zVtxFT0vPV() != useZvtxftovpv)
1493+
return;
1494+
if (useVtxItsTpc != -1 && collision.vtxITSTPC() != useVtxItsTpc)
1495+
return;
1496+
if (!isGoodRCTflag(collision))
1497+
return;
1498+
if (upcflag != -1 && collision.flags() != upcflag)
1499+
return;
1500+
13871501
registry.get<TH1>(HIST("Reco/Stat"))->Fill(truegapSide, 1.);
13881502
if (truegapSide != gapsideMC)
13891503
return;

0 commit comments

Comments
 (0)