Skip to content

Commit 4db093d

Browse files
stefanopolitanovkucera
stefanopolitano
and
vkucera
authored
PWGHF: Add Ds task and selector (#941)
* PWGHF: updating Ds task * PWGHF: updating cuts, candidate selector and candidate creator for Ds analysis * Fixing spacing in HFCandidateCreator3Prong * Reverting unwanted change in HFCandidateSelectionTables * Reverting unwanted change in HFCandidateSelectionTables * PWGHF: Addressing comments to Ds task * PWGHF: Updating Ds task and selector * Adjusting Dkkpi label * Moving Ds task in D2H folder * Fix clang-format * Updating inv. mass label accordingly to other patticles * Adding space in CMakeLists.txt * Adding subdirectories to CMakeLists.txt * Removing subdirectory from CMakeList * Removing Ds task from old folder * Update PWGHF/D2H/Tasks/taskDs.cxx Co-authored-by: Vít Kučera <vit.kucera@cern.ch> * Update PWGHF/TableProducer/HFDsToKKPiCandidateSelector.cxx Co-authored-by: Vít Kučera <vit.kucera@cern.ch> * Update PWGHF/TableProducer/HFDsToKKPiCandidateSelector.cxx Co-authored-by: Vít Kučera <vit.kucera@cern.ch> * Update PWGHF/D2H/Tasks/taskDs.cxx Co-authored-by: Vít Kučera <vit.kucera@cern.ch> * Updating maximum decay tree level for D+, Ds and Xic+ * Adding minimal documentation for Ds selection table * Update PWGHF/D2H/Tasks/taskDs.cxx Co-authored-by: Vít Kučera <vit.kucera@cern.ch> * Updating y binning definition, moving include line and fixing typos in Ds task * Removing documentation in HFCandidateSelectionTables * Updating Ds candidates selector with multiple PID possibilities * Renaming histos * Adjusting flag in HFFilter Co-authored-by: Vít Kučera <vit.kucera@cern.ch>
1 parent c3586b4 commit 4db093d

13 files changed

Lines changed: 554 additions & 29 deletions

‎EventFiltering/PWGHF/HFFilter.cxx‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -913,7 +913,7 @@ struct HfFilter { // Main struct for HF triggers
913913

914914
std::array<int8_t, kNCharmParticles - 1> is3Prong = {
915915
TESTBIT(cand3Prong.hfflag(), o2::aod::hf_cand_prong3::DecayType::DPlusToPiKPi),
916-
TESTBIT(cand3Prong.hfflag(), o2::aod::hf_cand_prong3::DecayType::DsToPiKK),
916+
TESTBIT(cand3Prong.hfflag(), o2::aod::hf_cand_prong3::DecayType::DsToKKPi),
917917
TESTBIT(cand3Prong.hfflag(), o2::aod::hf_cand_prong3::DecayType::LcToPKPi),
918918
TESTBIT(cand3Prong.hfflag(), o2::aod::hf_cand_prong3::DecayType::XicToPKPi)};
919919
if (!std::accumulate(is3Prong.begin(), is3Prong.end(), 0)) { // check if it's a D+, Ds+, Lc+ or Xic+
@@ -1169,7 +1169,7 @@ struct HfFilter { // Main struct for HF triggers
11691169

11701170
float deltaMassKKFirst = -1.f;
11711171
float deltaMassKKSecond = -1.f;
1172-
if (TESTBIT(cand3Prong.hfflag(), o2::aod::hf_cand_prong3::DecayType::DsToPiKK)) {
1172+
if (TESTBIT(cand3Prong.hfflag(), o2::aod::hf_cand_prong3::DecayType::DsToKKPi)) {
11731173
deltaMassKKFirst = std::abs(RecoDecay::m(std::array{pVecFirst, pVecSecond}, std::array{massK, massK}) - massPhi);
11741174
deltaMassKKSecond = std::abs(RecoDecay::m(std::array{pVecThird, pVecSecond}, std::array{massK, massK}) - massPhi);
11751175
}

‎PWGHF/Core/HFSelectorCuts.h‎

Lines changed: 54 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -27,6 +27,7 @@ enum Code {
2727
kD0 = 421,
2828
kD0bar = -421,
2929
kDPlus = 411,
30+
kDs = 431,
3031
kLambdaCPlus = 4122,
3132
kXiCPlus = 4232,
3233
kXiCCPlusPlus = 4422,
@@ -358,6 +359,59 @@ static const std::vector<std::string> pTBinLabels = {
358359
static const std::vector<std::string> cutVarLabels = {"deltaM", "pT Pi", "pT K", "decay length", "normalized decay length XY", "cos pointing angle", "cos pointing angle XY", "max normalized deltaIP"};
359360
} // namespace hf_cuts_dplus_topikpi
360361

362+
namespace hf_cuts_ds_tokkpi
363+
{
364+
static const int npTBins = 12;
365+
static const int nCutVars = 8;
366+
// momentary cuts
367+
constexpr double pTBins[npTBins + 1] = {
368+
1.,
369+
2.,
370+
3.,
371+
4.,
372+
5.,
373+
6.,
374+
7.,
375+
8.,
376+
10.,
377+
12.,
378+
16.,
379+
24.,
380+
36.};
381+
auto pTBins_v = std::vector<double>{pTBins, pTBins + npTBins + 1};
382+
383+
constexpr double cuts[npTBins][nCutVars] = {{0.2, 0.3, 0.3, 0.07, 6., 0.96, 0.985, 2.5}, /* 1 < pT < 2 */
384+
{0.2, 0.3, 0.3, 0.07, 5., 0.96, 0.985, 2.5}, /* 2 < pT < 3 */
385+
{0.2, 0.3, 0.3, 0.10, 5., 0.96, 0.980, 2.5}, /* 3 < pT < 4 */
386+
{0.2, 0.3, 0.3, 0.10, 5., 0.96, 0.000, 2.5}, /* 4 < pT < 5 */
387+
{0.2, 0.3, 0.3, 0.10, 5., 0.96, 0.000, 2.5}, /* 5 < pT < 6 */
388+
{0.2, 0.3, 0.3, 0.10, 5., 0.96, 0.000, 2.5}, /* 6 < pT < 7 */
389+
{0.2, 0.3, 0.3, 0.10, 5., 0.96, 0.000, 2.5}, /* 7 < pT < 8 */
390+
{0.2, 0.3, 0.3, 0.12, 5., 0.96, 0.000, 2.5}, /* 8 < pT < 10 */
391+
{0.2, 0.3, 0.3, 0.12, 5., 0.96, 0.000, 2.5}, /* 10 < pT < 12 */
392+
{0.2, 0.3, 0.3, 0.12, 5., 0.96, 0.000, 2.5}, /* 12 < pT < 16 */
393+
{0.2, 0.3, 0.3, 0.12, 5., 0.96, 0.000, 2.5}, /* 16 < pT < 24 */
394+
{0.2, 0.3, 0.3, 0.20, 5., 0.94, 0.000, 2.5}}; /* 24 < pT < 36 */
395+
396+
// row labels
397+
static const std::vector<std::string> pTBinLabels = {
398+
"pT bin 0",
399+
"pT bin 1",
400+
"pT bin 2",
401+
"pT bin 3",
402+
"pT bin 4",
403+
"pT bin 5",
404+
"pT bin 6",
405+
"pT bin 7",
406+
"pT bin 8",
407+
"pT bin 9",
408+
"pT bin 10",
409+
"pT bin 11"};
410+
411+
// column labels
412+
static const std::vector<std::string> cutVarLabels = {"m", "pT Pi", "pT K", "decay length", "normalized decay length XY", "cos pointing angle", "cos pointing angle XY", "max normalized deltaIP"};
413+
} // namespace hf_cuts_ds_tokkpi
414+
361415
namespace hf_cuts_xic_topkpi
362416
{
363417
static const int npTBins = 10;

‎PWGHF/D2H/CMakeLists.txt‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -8,3 +8,5 @@
88
# In applying this license CERN does not waive the privileges and immunities
99
# granted to it by virtue of its status as an Intergovernmental Organization
1010
# or submit itself to any jurisdiction.
11+
12+
add_subdirectory(Tasks)

‎PWGHF/D2H/Tasks/CMakeLists.txt‎

Lines changed: 15 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,15 @@
1+
# Copyright 2019-2020 CERN and copyright holders of ALICE O2.
2+
# See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
3+
# All rights not expressly granted are reserved.
4+
#
5+
# This software is distributed under the terms of the GNU General Public
6+
# License v3 (GPL Version 3), copied verbatim in the file "COPYING".
7+
#
8+
# In applying this license CERN does not waive the privileges and immunities
9+
# granted to it by virtue of its status as an Intergovernmental Organization
10+
# or submit itself to any jurisdiction.
11+
12+
o2physics_add_dpl_workflow(task-ds
13+
SOURCES taskDs.cxx
14+
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2::DetectorsVertexing
15+
COMPONENT_NAME Analysis)

0 commit comments

Comments
 (0)