Skip to content

Commit d00fd34

Browse files
mtorresc15mtorresc15
committed
Fix Omegac QA selections and DCA rapidity
1 parent 988a58a commit d00fd34

6 files changed

Lines changed: 26 additions & 11 deletions

File tree

‎PWGHF/Core/HfMlResponseOmegacToOmegaPiQa.h‎

100755100644
Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -10,7 +10,7 @@
1010
// or submit itself to any jurisdiction.
1111

1212
/// \file HfMlResponseOmegacToOmegaPiQa.h
13-
/// \brief Class to compute the ML response for Ωc± → Ω∓ π± analysis selections
13+
/// \brief Class to compute the ML response for Ωc0 → Ω- π+ (and charge-conjugate) analysis selections
1414
/// \author Yunfan Liu <yunfan.liu@cern.ch>, China University of Geosciences
1515
/// \author Maria Fernanda Torres Cabrera <maria.fernanda.torres.cabrera@cern.ch>, University of Houston
1616

‎PWGHF/D2H/Tasks/CMakeLists.txt‎

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -139,6 +139,11 @@ o2physics_add_dpl_workflow(task-omegac0-to-omega-pi
139139
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
140140
COMPONENT_NAME Analysis)
141141

142+
o2physics_add_dpl_workflow(task-omegac0-to-omega-pi-qa
143+
SOURCES taskOmegac0ToOmegaPiQa.cxx
144+
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
145+
COMPONENT_NAME Analysis)
146+
142147
o2physics_add_dpl_workflow(task-pt-fluc-charm-hadrons
143148
SOURCES taskPtFlucCharmHadrons.cxx
144149
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::EventFilteringUtils

‎PWGHF/D2H/Tasks/taskOmegac0ToOmegaPiQa.cxx‎

100755100644
Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -175,7 +175,7 @@ struct HfTaskOmegac0ToOmegaPiQa {
175175
if constexpr (UseKfParticle) {
176176
yOmegac = candidate.kfRapOmegac();
177177
} else {
178-
yOmegac = candidate.y(o2::constants::physics::MassOmegaC0);
178+
yOmegac = RecoDecay::y(std::array{candidate.pxCharmBaryon(), candidate.pyCharmBaryon(), candidate.pzCharmBaryon()}, o2::constants::physics::MassOmegaC0);
179179
}
180180
if (yCandRecoMax >= 0. && std::abs(yOmegac) > yCandRecoMax) {
181181
continue;
@@ -200,7 +200,7 @@ struct HfTaskOmegac0ToOmegaPiQa {
200200
if constexpr (UseKfParticle) {
201201
yOmegac = candidate.kfRapOmegac();
202202
} else {
203-
yOmegac = candidate.y(o2::constants::physics::MassOmegaC0);
203+
yOmegac = RecoDecay::y(std::array{candidate.pxCharmBaryon(), candidate.pyCharmBaryon(), candidate.pzCharmBaryon()}, o2::constants::physics::MassOmegaC0);
204204
}
205205

206206
if (yCandRecoMax >= 0. && std::abs(yOmegac) > yCandRecoMax) {
@@ -239,7 +239,7 @@ struct HfTaskOmegac0ToOmegaPiQa {
239239
if constexpr (UseKfParticle) {
240240
yOmegac = candidate.kfRapOmegac();
241241
} else {
242-
yOmegac = candidate.y(o2::constants::physics::MassOmegaC0);
242+
yOmegac = RecoDecay::y(std::array{candidate.pxCharmBaryon(), candidate.pyCharmBaryon(), candidate.pzCharmBaryon()}, o2::constants::physics::MassOmegaC0);
243243
}
244244

245245
if (yCandRecoMax >= 0. && std::abs(yOmegac) > yCandRecoMax) {
@@ -285,7 +285,7 @@ struct HfTaskOmegac0ToOmegaPiQa {
285285
if constexpr (UseKfParticle) {
286286
yOmegac = candidate.kfRapOmegac();
287287
} else {
288-
yOmegac = candidate.y(o2::constants::physics::MassOmegaC0);
288+
yOmegac = RecoDecay::y(std::array{candidate.pxCharmBaryon(), candidate.pyCharmBaryon(), candidate.pzCharmBaryon()}, o2::constants::physics::MassOmegaC0);
289289
}
290290

291291
if (yCandRecoMax >= 0. && std::abs(yOmegac) > yCandRecoMax) {

‎PWGHF/TableProducer/CMakeLists.txt‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -172,6 +172,11 @@ o2physics_add_dpl_workflow(candidate-selector-omegac0-to-omega-pi
172172
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::MLCore
173173
COMPONENT_NAME Analysis)
174174

175+
o2physics_add_dpl_workflow(candidate-selector-omegac0-to-omega-pi-qa
176+
SOURCES candidateSelectorOmegac0ToOmegaPiQa.cxx
177+
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::MLCore
178+
COMPONENT_NAME Analysis)
179+
175180
o2physics_add_dpl_workflow(candidate-selector-omegac0-xic0-to-omega-ka
176181
SOURCES candidateSelectorOmegac0Xic0ToOmegaKa.cxx
177182
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::MLCore
@@ -279,6 +284,11 @@ o2physics_add_dpl_workflow(tree-creator-omegac0-to-omega-pi
279284
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
280285
COMPONENT_NAME Analysis)
281286

287+
o2physics_add_dpl_workflow(tree-creator-omegac0-to-omega-pi-qa
288+
SOURCES treeCreatorOmegac0ToOmegaPiQa.cxx
289+
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
290+
COMPONENT_NAME Analysis)
291+
282292
o2physics_add_dpl_workflow(tree-creator-to-xi-pi
283293
SOURCES treeCreatorToXiPi.cxx
284294
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore

‎PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPiQa.cxx‎

100755100644
Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -343,8 +343,8 @@ struct HfCandidateSelectorToOmegaPiQa {
343343
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(13, "dcaXYToPvCascDau");
344344
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(14, "ptKaFromCasc");
345345
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(15, "cosPaV0ToCasc");
346-
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(16, "chi2GeoV0");
347-
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(17, "chi2GeoCasc");
346+
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(16, "v0Chi2OverNdf");
347+
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(17, "cascChi2OverNdf");
348348
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(18, "chi2TopoCascToPv");
349349
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(19, "chi2TopoV0ToCasc");
350350
registry.get<TH1>(HIST("hSelStatusLf"))->GetXaxis()->SetBinLabel(20, "v0ldl");
@@ -359,7 +359,7 @@ struct HfCandidateSelectorToOmegaPiQa {
359359
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(4, "ptPiFromCharmBaryon");
360360
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(5, "kfptOmegac");
361361
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(6, "cosPaCascToOmegac");
362-
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(7, "chi2GeoOmegac");
362+
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(7, "omegacChi2OverNdf");
363363
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(8, "chi2TopoOmegacToPv");
364364
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(9, "chi2TopoCascToOmegac");
365365
registry.get<TH1>(HIST("hSelStatusHf"))->GetXaxis()->SetBinLabel(10, "decayLenXYOmegac");
@@ -513,11 +513,11 @@ struct HfCandidateSelectorToOmegaPiQa {
513513
registry.fill(HIST("hSelStatusLf"), 14.0);
514514

515515
// Chi2
516-
if (candidate.chi2GeoV0() < 0 || candidate.chi2GeoV0() > kfConfigurableGroup.v0Chi2OverNdfMax) {
516+
if (candidate.v0Chi2OverNdf() < 0 || candidate.v0Chi2OverNdf() > kfConfigurableGroup.v0Chi2OverNdfMax) {
517517
return false;
518518
}
519519
registry.fill(HIST("hSelStatusLf"), 15.0);
520-
if (candidate.chi2GeoCasc() < 0 || candidate.chi2GeoCasc() > kfConfigurableGroup.cascChi2OverNdfMax) {
520+
if (candidate.cascChi2OverNdf() < 0 || candidate.cascChi2OverNdf() > kfConfigurableGroup.cascChi2OverNdfMax) {
521521
return false;
522522
}
523523
registry.fill(HIST("hSelStatusLf"), 16.0);
@@ -621,7 +621,7 @@ struct HfCandidateSelectorToOmegaPiQa {
621621
registry.fill(HIST("hSelStatusHf"), 5.0);
622622

623623
// Chi2
624-
if (candidate.chi2GeoOmegac() < 0 || candidate.chi2GeoOmegac() > kfConfigurableGroup.omegacChi2OverNdfMax) {
624+
if (candidate.omegacChi2OverNdf() < 0 || candidate.omegacChi2OverNdf() > kfConfigurableGroup.omegacChi2OverNdfMax) {
625625
return false;
626626
}
627627
registry.fill(HIST("hSelStatusHf"), 6.0);

‎PWGHF/TableProducer/treeCreatorOmegac0ToOmegaPiQa.cxx‎

100755100644
File mode changed.

0 commit comments

Comments
 (0)