@@ -152,27 +152,6 @@ struct HfTaskCharmPolarisation {
152152 Configurable<int > minCent{" minCent" , 30 , " Minimum centrality (0-100) to be considered in the analysis" };
153153 Configurable<int > maxCent{" maxCent" , 50 , " Maximum centrality (0-100) to be considered in the analysis" };
154154
155- ConfigurableAxis configThnAxisInvMass{" configThnAxisInvMass" , {200 , 0 .139f , 0 .179f }, " #it{M} (GeV/#it{c}^{2})" };
156- ConfigurableAxis configThnAxisPt{" configThnAxisPt" , {100 , 0 .f , 100 .f }, " #it{p}_{T} (GeV/#it{c})" };
157- ConfigurableAxis configThnAxisY{" configThnAxisY" , {20 , -1 .f , 1 .f }, " #it{y}" };
158- ConfigurableAxis configThnAxisCosThetaStarHelicity{" configThnAxisCosThetaStarHelicity" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{helicity})" };
159- ConfigurableAxis configThnAxisCosThetaStarProduction{" configThnAxisCosThetaStarProduction" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{production})" };
160- ConfigurableAxis configThnAxisCosThetaStarRandom{" configThnAxisCosThetaStarRandom" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{random})" };
161- ConfigurableAxis configThnAxisCosThetaStarBeam{" configThnAxisCosThetaStarBeam" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{beam})" };
162- ConfigurableAxis configThnAxisCosThetaStarEP{" configThnAxisCosThetaStarEP" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{EP})" };
163- ConfigurableAxis configThnAxisMlBkg{" configThnAxisMlBkg" , {100 , 0 .f , 1 .f }, " ML bkg" };
164- ConfigurableAxis configThnAxisInvMassD0{" configThnAxisInvMassD0" , {250 , 1 .65f , 2 .15f }, " #it{M}(D^{0}) (GeV/#it{c}^{2})" }; // only for D*+
165- ConfigurableAxis configThnAxisInvMassKPiLc{" configThnAxisInvMassKPiLc" , {120 , 0 .65f , 1 .25f }, " #it{M}(K#pi) from #Lambda_{c}^{+} (GeV/#it{c}^{2})" }; // only for Lc+->pKpi
166- // ConfigurableAxis configThnAxisMlPrompt{"configThnAxisMlPrompt", {100, 0.f, 1.f}, "ML prompt"};
167- ConfigurableAxis configThnAxisMlNonPrompt{" configThnAxisMlNonPrompt" , {100 , 0 .f , 1 .f }, " ML non-prompt" };
168- ConfigurableAxis configThnAxisCent{" configThnAxisCent" , {102 , -1 .f , 101 .f }, " centrality (%)" };
169- ConfigurableAxis configThnAxisNumPvContributors{" configThnAxisNumPvContributors" , {300 , -0 .5f , 299 .5f }, " num PV contributors" };
170- ConfigurableAxis configThnAxisPtB{" configThnAxisPtB" , {3000 , 0 .f , 300 .f }, " #it{p}_{T}(B mother) (GeV/#it{c})" };
171- ConfigurableAxis configThnAxisAbsEtaTrackMin{" configThnAxisEtaTrackMin" , {3 , 0 .f , 0 .3f }, " min |#it{#eta_{track}}|" };
172- ConfigurableAxis configThnAxisNumItsClsMin{" configThnAxisNumItsClsMin" , {4 , 3 .5f , 7 .5f }, " min #it{N}_{cls ITS}" };
173- ConfigurableAxis configThnAxisNumTpcClsMin{" configThnAxisNumTpcClsMin" , {3 , 79 .5f , 140 .5f }, " min #it{N}_{cls TPC}" };
174- ConfigurableAxis configThnAxisCharge{" configThnAxisCharge" , {2 , -2 .f , 2 .f }, " electric charge" };
175-
176155 // / activate rotational background
177156 Configurable<int > nBkgRotations{" nBkgRotations" , 0 , " Number of rotated copies (background) per each original candidate" };
178157 Configurable<float > minRotAngleMultByPi{" minRotAngleMultByPi" , 5 . / 6 , " Minimum angle rotation for track rotation, to be multiplied by pi" };
@@ -213,7 +192,7 @@ struct HfTaskCharmPolarisation {
213192 struct : ConfigurableGroup {
214193 // / monitoring histograms (Dalitz plot)
215194 Configurable<bool > activateTHnLcChannelMonitor{" activateTHnLcChannelMonitor" , false , " Flag to switch on the monitoring THnSparse of M2(Kpi), M2(pK), M2(ppi), pt correlation for Lc -> pKpi" };
216- ConfigurableAxis configThnAxisInvMass2KPiLcMonitoring{" configThnAxisInvMassKPiLcMonitoring " , {200 , 0 .3f , 2 .3f }, " #it{M}^{2}(K#pi) from #Lambda_{c}^{+} (GeV/#it{c}^{2})" };
195+ ConfigurableAxis configThnAxisInvMass2KPiLcMonitoring{" configThnAxisInvMass2KPiLcMonitoring " , {200 , 0 .3f , 2 .3f }, " #it{M}^{2}(K#pi) from #Lambda_{c}^{+} (GeV/#it{c}^{2})" };
217196 ConfigurableAxis configThnAxisInvMass2PKLcMonitoring{" configThnAxisInvMass2PKLcMonitoring" , {320 , 2 .f , 5 .2f }, " #it{M}^{2}(pK) from #Lambda_{c}^{+} (GeV/#it{c}^{2})" };
218197 ConfigurableAxis configThnAxisInvMass2PPiLcMonitoring{" configThnAxisInvMass2PPiLcMonitoring" , {400 , 1 .f , 5 .f }, " #it{M}^{2}(p#pi) from #Lambda_{c}^{+} (GeV/#it{c}^{2})" };
219198
@@ -226,13 +205,13 @@ struct HfTaskCharmPolarisation {
226205
227206 // / Monitoring of phi Euler angle
228207 Configurable<bool > activateTHnEulerPhiMonitor{" activateTHnEulerPhiMonitor" , false , " Flag to switch on the monitoring THnSparse vs. Euler angle phi (Lc -> pKpi)" };
229- ConfigurableAxis configTHnAxisEulerPhi{" configTHnAxisEulerPhi" , {24 , -o2::constants::math::PI, o2::constants::math::PI}, " Euler polar angle #phi" };
230208
231209 // / Application of rapidity cut for reconstructed candidates
232210 Configurable<float > rapidityCut{" rapidityCut" , 999 .f , " Max. value of reconstructed candidate rapidity (abs. value)" };
233211
234- Filter filterSelectDstarCandidates = aod::hf_sel_candidate_dstar::isSelDstarToD0Pi == selectionFlagDstarToD0Pi;
235- Filter filterSelectLcToPKPiCandidates = (aod::hf_sel_candidate_lc::isSelLcToPKPi >= selectionFlagLcToPKPi) || (aod::hf_sel_candidate_lc::isSelLcToPiKP >= selectionFlagLcToPKPi);
212+ HfHelper hfHelper;
213+ SliceCache cache;
214+ EventPlaneHelper epHelper;
236215
237216 using CollisionsWithMcLabels = soa::SmallGroups<soa::Join<aod::Collisions, aod::McCollisionLabels>>;
238217 using CollsWithQvecs = soa::Join<aod::Collisions, aod::EvSels, aod::QvectorFT0Cs, aod::QvectorFT0Ms, aod::QvectorFV0As, aod::CentFT0Ms, aod::CentFT0Cs>;
@@ -257,7 +236,9 @@ struct HfTaskCharmPolarisation {
257236 using FilteredCandLcToPKPiWSelFlagAndMc = soa::Filtered<soa::Join<CandLcToPKPiWSelFlag, aod::HfCand3ProngMcRec>>;
258237 using FilteredCandLcToPKPiWSelFlagAndMcAndMl = soa::Filtered<soa::Join<CandLcToPKPiWSelFlag, aod::HfCand3ProngMcRec, aod::HfMlLcToPKPi>>;
259238
260- SliceCache cache;
239+ Filter filterSelectDstarCandidates = aod::hf_sel_candidate_dstar::isSelDstarToD0Pi == selectionFlagDstarToD0Pi;
240+ Filter filterSelectLcToPKPiCandidates = (aod::hf_sel_candidate_lc::isSelLcToPKPi >= selectionFlagLcToPKPi) || (aod::hf_sel_candidate_lc::isSelLcToPiKP >= selectionFlagLcToPKPi);
241+
261242 Preslice<FilteredCandDstarWSelFlag> dstarPerCollision = aod::hf_cand::collisionId;
262243 Preslice<FilteredCandDstarWSelFlagAndMl> dstarWithMlPerCollision = aod::hf_cand::collisionId;
263244 Preslice<FilteredCandDstarWSelFlagAndMc> dstarWithMcPerCollision = aod::hf_cand::collisionId;
@@ -268,8 +249,28 @@ struct HfTaskCharmPolarisation {
268249 Preslice<FilteredCandLcToPKPiWSelFlagAndMc> lcToPKPiWithMcPerCollision = aod::hf_cand::collisionId;
269250 Preslice<FilteredCandLcToPKPiWSelFlagAndMcAndMl> lcToPKPiWithMcAndMlPerCollision = aod::hf_cand::collisionId;
270251
271- HfHelper hfHelper;
272- EventPlaneHelper epHelper;
252+ ConfigurableAxis configTHnAxisEulerPhi{" configTHnAxisEulerPhi" , {24 , -o2::constants::math::PI, o2::constants::math::PI}, " Euler polar angle #phi" };
253+ ConfigurableAxis configThnAxisInvMass{" configThnAxisInvMass" , {200 , 0 .139f , 0 .179f }, " #it{M} (GeV/#it{c}^{2})" };
254+ ConfigurableAxis configThnAxisPt{" configThnAxisPt" , {100 , 0 .f , 100 .f }, " #it{p}_{T} (GeV/#it{c})" };
255+ ConfigurableAxis configThnAxisY{" configThnAxisY" , {20 , -1 .f , 1 .f }, " #it{y}" };
256+ ConfigurableAxis configThnAxisCosThetaStarHelicity{" configThnAxisCosThetaStarHelicity" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{helicity})" };
257+ ConfigurableAxis configThnAxisCosThetaStarProduction{" configThnAxisCosThetaStarProduction" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{production})" };
258+ ConfigurableAxis configThnAxisCosThetaStarRandom{" configThnAxisCosThetaStarRandom" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{random})" };
259+ ConfigurableAxis configThnAxisCosThetaStarBeam{" configThnAxisCosThetaStarBeam" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{beam})" };
260+ ConfigurableAxis configThnAxisCosThetaStarEP{" configThnAxisCosThetaStarEP" , {20 , -1 .f , 1 .f }, " cos(#vartheta_{EP})" };
261+ ConfigurableAxis configThnAxisMlBkg{" configThnAxisMlBkg" , {100 , 0 .f , 1 .f }, " ML bkg" };
262+ ConfigurableAxis configThnAxisInvMassD0{" configThnAxisInvMassD0" , {250 , 1 .65f , 2 .15f }, " #it{M}(D^{0}) (GeV/#it{c}^{2})" }; // only for D*+
263+ ConfigurableAxis configThnAxisInvMassKPiLc{" configThnAxisInvMassKPiLc" , {120 , 0 .65f , 1 .25f }, " #it{M}(K#pi) from #Lambda_{c}^{+} (GeV/#it{c}^{2})" }; // only for Lc+->pKpi
264+ // ConfigurableAxis configThnAxisMlPrompt{"configThnAxisMlPrompt", {100, 0.f, 1.f}, "ML prompt"};
265+ ConfigurableAxis configThnAxisMlNonPrompt{" configThnAxisMlNonPrompt" , {100 , 0 .f , 1 .f }, " ML non-prompt" };
266+ ConfigurableAxis configThnAxisCent{" configThnAxisCent" , {102 , -1 .f , 101 .f }, " centrality (%)" };
267+ ConfigurableAxis configThnAxisNumPvContributors{" configThnAxisNumPvContributors" , {300 , -0 .5f , 299 .5f }, " num PV contributors" };
268+ ConfigurableAxis configThnAxisPtB{" configThnAxisPtB" , {3000 , 0 .f , 300 .f }, " #it{p}_{T}(B mother) (GeV/#it{c})" };
269+ ConfigurableAxis configThnAxisAbsEtaTrackMin{" configThnAxisAbsEtaTrackMin" , {3 , 0 .f , 0 .3f }, " min |#it{#eta_{track}}|" };
270+ ConfigurableAxis configThnAxisNumItsClsMin{" configThnAxisNumItsClsMin" , {4 , 3 .5f , 7 .5f }, " min #it{N}_{cls ITS}" };
271+ ConfigurableAxis configThnAxisNumTpcClsMin{" configThnAxisNumTpcClsMin" , {3 , 79 .5f , 140 .5f }, " min #it{N}_{cls TPC}" };
272+ ConfigurableAxis configThnAxisCharge{" configThnAxisCharge" , {2 , -2 .f , 2 .f }, " electric charge" };
273+
273274 HistogramRegistry registry{" registry" , {}};
274275
275276 void init (InitContext&)
0 commit comments