Skip to content

[ALICE3] Add DCAFitterN QA #15164

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
jesgum wants to merge 2 commits into
base: master
Choose a base branch
from
+73 −0
Open

[ALICE3] Add DCAFitterN QA #15164

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
9f9a243
Add DCAFitterN QA
jesgum Feb 25, 2026
2dcc53a
megalinter fix
jesgum Feb 25, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
48 changes: 48 additions & 0 deletions ALICE3/TableProducer/alice3MulticharmFinder.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -227,6 +227,9 @@ struct Alice3MulticharmFinder {
if (nCand == 0) {
return false;
}

const u_int8_t fitterStatusCode = fitter.getFitStatus();
histos.fill(HIST("hFitterStatusCode"), fitterStatusCode);
//}-{}-{}-{}-{}-{}-{}-{}-{}-{}

o2::track::TrackParCov t0new = fitter.getTrack(0);
Expand Down Expand Up @@ -298,6 +301,9 @@ struct Alice3MulticharmFinder {
if (nCand == 0) {
return false;
}

const u_int8_t fitter3StatusCode = fitter3.getFitStatus();
histos.fill(HIST("hFitter3StatusCode"), fitter3StatusCode);
//}-{}-{}-{}-{}-{}-{}-{}-{}-{}

t0 = fitter3.getTrack(0);
Expand Down Expand Up @@ -433,6 +439,48 @@ struct Alice3MulticharmFinder {
fitter3.setBz(magneticField);
fitter3.setMatCorrType(o2::base::Propagator::MatCorrType::USEMatCorrNONE);

auto hFitterStatusCode = histos.add<TH1>("hFitterStatusCode", "hFitterStatusCode", kTH1D, {{15, -0.5, 14.5}});
hFitterStatusCode->GetXaxis()->SetBinLabel(1, "None"); // no status set (should not be possible!)

/* Good Conditions */
hFitterStatusCode->GetXaxis()->SetBinLabel(2, "Converged"); // fit converged
hFitterStatusCode->GetXaxis()->SetBinLabel(3, "MaxIter"); // max iterations reached before fit convergence

/* Error Conditions */
hFitterStatusCode->GetXaxis()->SetBinLabel(4, "NoCrossing"); // no reasaonable crossing was found
hFitterStatusCode->GetXaxis()->SetBinLabel(5, "RejRadius"); // radius of crossing was not acceptable
hFitterStatusCode->GetXaxis()->SetBinLabel(6, "RejTrackX"); // one candidate track x was below the mimimum required radius
hFitterStatusCode->GetXaxis()->SetBinLabel(7, "RejTrackRoughZ"); // rejected by rough cut on tracks Z difference
hFitterStatusCode->GetXaxis()->SetBinLabel(8, "RejChi2Max"); // rejected by maximum chi2 cut
hFitterStatusCode->GetXaxis()->SetBinLabel(9, "FailProp"); // propagation of at least prong to PCA failed
hFitterStatusCode->GetXaxis()->SetBinLabel(10, "FailInvCov"); // inversion of cov.-matrix failed
hFitterStatusCode->GetXaxis()->SetBinLabel(11, "FailInvWeight"); // inversion of Ti weight matrix failed
hFitterStatusCode->GetXaxis()->SetBinLabel(12, "FailInv2ndDeriv"); // inversion of 2nd derivatives failed
hFitterStatusCode->GetXaxis()->SetBinLabel(13, "FailCorrTracks"); // correction of tracks to updated x failed
hFitterStatusCode->GetXaxis()->SetBinLabel(14, "FailCloserAlt"); // alternative PCA is closer
hFitterStatusCode->GetXaxis()->SetBinLabel(15, "NStatusesDefined");

auto hFitter3StatusCode = histos.add<TH1>("hFitter3StatusCode", "hFitter3StatusCode", kTH1D, {{15, -0.5, 14.5}});
hFitter3StatusCode->GetXaxis()->SetBinLabel(1, "None"); // no status set (should not be possible!)

/* Good Conditions */
hFitter3StatusCode->GetXaxis()->SetBinLabel(2, "Converged"); // fit converged
hFitter3StatusCode->GetXaxis()->SetBinLabel(3, "MaxIter"); // max iterations reached before fit convergence

/* Error Conditions */
hFitter3StatusCode->GetXaxis()->SetBinLabel(4, "NoCrossing"); // no reasaonable crossing was found
hFitter3StatusCode->GetXaxis()->SetBinLabel(5, "RejRadius"); // radius of crossing was not acceptable
hFitter3StatusCode->GetXaxis()->SetBinLabel(6, "RejTrackX"); // one candidate track x was below the mimimum required radius
hFitter3StatusCode->GetXaxis()->SetBinLabel(7, "RejTrackRoughZ"); // rejected by rough cut on tracks Z difference
hFitter3StatusCode->GetXaxis()->SetBinLabel(8, "RejChi2Max"); // rejected by maximum chi2 cut
hFitter3StatusCode->GetXaxis()->SetBinLabel(9, "FailProp"); // propagation of at least prong to PCA failed
hFitter3StatusCode->GetXaxis()->SetBinLabel(10, "FailInvCov"); // inversion of cov.-matrix failed
hFitter3StatusCode->GetXaxis()->SetBinLabel(11, "FailInvWeight"); // inversion of Ti weight matrix failed
hFitter3StatusCode->GetXaxis()->SetBinLabel(12, "FailInv2ndDeriv"); // inversion of 2nd derivatives failed
hFitter3StatusCode->GetXaxis()->SetBinLabel(13, "FailCorrTracks"); // correction of tracks to updated x failed
hFitter3StatusCode->GetXaxis()->SetBinLabel(14, "FailCloserAlt"); // alternative PCA is closer
hFitter3StatusCode->GetXaxis()->SetBinLabel(15, "NStatusesDefined");

INSERT_HIST(std::string("h2dGenXi"), "h2dGenXi", {kTH2D, {{axisPt, axisEta}}});
INSERT_HIST(std::string("h2dGenXiC"), "h2dGenXiC", {kTH2D, {{axisPt, axisEta}}});
INSERT_HIST(std::string("h2dGenXiCC"), "h2dGenXiCC", {kTH2D, {{axisPt, axisEta}}});
Expand Down
25 changes: 25 additions & 0 deletions ALICE3/TableProducer/alice3strangenessFinder.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,6 +44,7 @@
#include <TLorentzVector.h>

#include <cstdlib>
#include <vector>

using namespace o2;
// using namespace o2::analysis;
Expand Down Expand Up @@ -153,6 +154,27 @@ struct Alice3strangenessFinder {
fitter.setMatCorrType(o2::base::Propagator::MatCorrType::USEMatCorrNONE);

histos.add("hFitterQA", "", kTH1D, {{10, 0, 10}}); // For QA reasons, counting found candidates at different stages
auto hFitterStatusCode = histos.add<TH1>("hFitterStatusCode", "hFitterStatusCode", kTH1D, {{15, -0.5, 14.5}});
hFitterStatusCode->GetXaxis()->SetBinLabel(1, "None"); // no status set (should not be possible!)

/* Good Conditions */
hFitterStatusCode->GetXaxis()->SetBinLabel(2, "Converged"); // fit converged
hFitterStatusCode->GetXaxis()->SetBinLabel(3, "MaxIter"); // max iterations reached before fit convergence

/* Error Conditions */
hFitterStatusCode->GetXaxis()->SetBinLabel(4, "NoCrossing"); // no reasaonable crossing was found
hFitterStatusCode->GetXaxis()->SetBinLabel(5, "RejRadius"); // radius of crossing was not acceptable
hFitterStatusCode->GetXaxis()->SetBinLabel(6, "RejTrackX"); // one candidate track x was below the mimimum required radius
hFitterStatusCode->GetXaxis()->SetBinLabel(7, "RejTrackRoughZ"); // rejected by rough cut on tracks Z difference
hFitterStatusCode->GetXaxis()->SetBinLabel(8, "RejChi2Max"); // rejected by maximum chi2 cut
hFitterStatusCode->GetXaxis()->SetBinLabel(9, "FailProp"); // propagation of at least prong to PCA failed
hFitterStatusCode->GetXaxis()->SetBinLabel(10, "FailInvCov"); // inversion of cov.-matrix failed
hFitterStatusCode->GetXaxis()->SetBinLabel(11, "FailInvWeight"); // inversion of Ti weight matrix failed
hFitterStatusCode->GetXaxis()->SetBinLabel(12, "FailInv2ndDeriv"); // inversion of 2nd derivatives failed
hFitterStatusCode->GetXaxis()->SetBinLabel(13, "FailCorrTracks"); // correction of tracks to updated x failed
hFitterStatusCode->GetXaxis()->SetBinLabel(14, "FailCloserAlt"); // alternative PCA is closer
hFitterStatusCode->GetXaxis()->SetBinLabel(15, "NStatusesDefined");

histos.add("hPtPosDau", "", kTH1D, {axisPt});
histos.add("hPtNegDau", "", kTH1D, {axisPt});
histos.add("hPtPosDauAfterV0Finding", "", kTH2D, {axisPt, axisPt});
Expand Down Expand Up @@ -220,6 +242,9 @@ struct Alice3strangenessFinder {
} catch (...) {
return false;
}

const u_int8_t fitterStatusCode = fitter.getFitStatus();
histos.fill(HIST("hFitterStatusCode"), fitterStatusCode);
histos.fill(HIST("hFitterQA"), 1.5);
if (nCand == 0) {
LOG(info) << "0 candidates found by fitter";
Expand Down
Loading