Skip to content

[DPG] Update aQCMFTTracks.cxx #12056

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

Merged
njacazio merged 10 commits into from
Jul 15, 2025
Merged

[DPG] Update aQCMFTTracks.cxx #12056

Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
f02fe25
Update aQCMFTTracks.cxx
njacazio Jul 15, 2025
f87d5ff
Update aQCMFTTracks.cxx
njacazio Jul 15, 2025
e0134ec
Update aQCMFTTracks.cxx
njacazio Jul 15, 2025
5d7f1da
Update aQCMFTTracks.cxx
njacazio Jul 15, 2025
bdf9dc3
Please consider the following formatting changes
alibuild Jul 15, 2025
0c2e3cd
Update aQCMFTTracks.cxx
njacazio Jul 15, 2025
741a348
Merge pull request #47 from alibuild/alibot-cleanup-12056
njacazio Jul 15, 2025
4e7309a
Please consider the following formatting changes
alibuild Jul 15, 2025
5669b09
Merge pull request #48 from alibuild/alibot-cleanup-12056
njacazio Jul 15, 2025
d7f821f
Update aQCMFTTracks.cxx
njacazio Jul 15, 2025
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
164 changes: 134 additions & 30 deletions DPG/Tasks/MFT/aQCMFTTracks.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.

Check failure on line 1 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/workflow-file] 

Name of a workflow file must match the name of the main struct in it (without the PWG prefix). (Class implementation files should be in "Core" directories.)
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
Expand All @@ -16,63 +16,124 @@
/// \author David Grund
/// \since

#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/Multiplicity.h"

#include "CCDB/BasicCCDBManager.h"
#include "Framework/runDataProcessing.h"
#include "Framework/AnalysisTask.h"
#include "Framework/AnalysisDataModel.h"
#include "CommonConstants/LHCConstants.h"
#include "DataFormatsITSMFT/ROFRecord.h"
#include "Framework/ASoAHelpers.h"

#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/DataTypes.h"
#include "Common/DataModel/Multiplicity.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/Centrality.h"
#include "CommonConstants/LHCConstants.h"
#include "Framework/TimingInfo.h"
#include "DataFormatsITSMFT/ROFRecord.h"
#include "Framework/runDataProcessing.h"

#include <TH1F.h>
#include <TH2F.h>

#include <algorithm>

using namespace o2;
using namespace o2::framework;
using namespace o2::aod;

struct CheckMFT {
HistogramRegistry registry{"registry",
{// 2d histograms
{"mMFTTrackEtaPhi_5_MinClusters", "Track #eta , #phi (NCls >= 5); #eta; #phi", {HistType::kTH2F, {{50, -4, -2}, {100, -3.2, 3.2}}}},
{"mMFTTrackXY_5_MinClusters", "Track Position (NCls >= 5); x; y", {HistType::kTH2F, {{320, -16, 16}, {320, -16, 16}}}},
{"mMFTTrackEtaPhi_7_MinClusters", "Track #eta , #phi (NCls >= 7); #eta; #phi", {HistType::kTH2F, {{50, -4, -2}, {100, -3.2, 3.2}}}},
{"mMFTTrackXY_7_MinClusters", "Track Position (NCls >= 7); x; y", {HistType::kTH2F, {{320, -16, 16}, {320, -16, 16}}}},
{"mMFTTrackEtaPhi_8_MinClusters", "Track #eta , #phi (NCls >= 8); #eta; #phi", {HistType::kTH2F, {{50, -4, -2}, {100, -3.2, 3.2}}}},
{"mMFTTrackXY_8_MinClusters", "Track Position (NCls >= 8); x; y", {HistType::kTH2F, {{320, -16, 16}, {320, -16, 16}}}},
// 1d histograms
{"mMFTTrackEta", "Track #eta; #eta; # entries", {HistType::kTH1F, {{50, -4, -2}}}},
{"mMFTTrackNumberOfClusters", "Number Of Clusters Per Track; # clusters; # entries", {HistType::kTH1F, {{10, 0.5, 10.5}}}},
{"mMFTTrackPhi", "Track #phi; #phi; # entries", {HistType::kTH1F, {{100, -3.2, 3.2}}}},
{"mMFTTrackTanl", "Track tan #lambda; tan #lambda; # entries", {HistType::kTH1F, {{100, -25, 0}}}},
{"mMFTTrackInvQPt", "Track q/p_{T}; q/p_{T} [1/GeV]; # entries", {HistType::kTH1F, {{250, -10, 10}}}}}};
HistogramRegistry registry{"registry"};
Configurable<bool> avClsPlots{"avClsPlots", false, "Enable average cluster plots"};

void init(o2::framework::InitContext&)
{

const AxisSpec etaAxis{50, -4, -2, "#eta"};
const AxisSpec phiAxis{100, -3.2, 3.2, "#phi"};
const AxisSpec xAxis{320, -16, 16, "x"};
const AxisSpec clsAxis{10, 0.5, 10.5, "# clusters"};
const AxisSpec yAxis{320, -16, 16, "y"};
const AxisSpec tanLamAxis{100, -25, 0, "tan #lambda"};
const AxisSpec invQPtAxis{250, -10, 10, "q/p_{T} [1/GeV]"};

registry.add("mMFTTrackPhi", "Track #phi", {HistType::kTH1F, {phiAxis}});
registry.add("mMFTTrackTanl", "Track tan #lambda", {HistType::kTH1F, {tanLamAxis}});
registry.add("mMFTTrackInvQPt", "Track q/p_{T}", {HistType::kTH1F, {invQPtAxis}});
registry.add("mMFTTrackEta", "Track #eta", {HistType::kTH1F, {etaAxis}});

registry.add("mMFTTrackEtaPhi_5_MinClusters", "Track Position (NCls >= 5)", {HistType::kTH2F, {etaAxis, phiAxis}});
registry.add("mMFTTrackEtaPhi_6_MinClusters", "Track Position (NCls >= 6)", {HistType::kTH2F, {etaAxis, phiAxis}});
registry.add("mMFTTrackEtaPhi_7_MinClusters", "Track Position (NCls >= 7)", {HistType::kTH2F, {etaAxis, phiAxis}});
registry.add("mMFTTrackEtaPhi_8_MinClusters", "Track Position (NCls >= 8)", {HistType::kTH2F, {etaAxis, phiAxis}});

registry.add("mMFTTrackXY_5_MinClusters", "Track Position (NCls >= 5)", {HistType::kTH2F, {xAxis, yAxis}});
registry.add("mMFTTrackXY_6_MinClusters", "Track Position (NCls >= 6)", {HistType::kTH2F, {xAxis, yAxis}});
registry.add("mMFTTrackXY_7_MinClusters", "Track Position (NCls >= 7)", {HistType::kTH2F, {xAxis, yAxis}});
registry.add("mMFTTrackXY_8_MinClusters", "Track Position (NCls >= 8)", {HistType::kTH2F, {xAxis, yAxis}});

registry.add("mMFTTrackNumberOfClusters", "Number Of Clusters Per Track", {HistType::kTH1F, {clsAxis}});

if (avClsPlots) {
registry.add("mMFTTrackAvgClusters", "Average number of clusters per track; p;# clusters; # entries", {HistType::kTH2F, {{100, 0, 100}, {100, 0, 100}}});
registry.add("mMFTTrackAvgClustersTru", "Average number of clusters per track; p;# clusters; # entries", {HistType::kTH2F, {{100, 0, 100}, {100, 0, 100}}});
if (doprocessMC) {
registry.add("mMFTTrackAvgClustersHe", "Average number of clusters per track; p;# clusters; # entries", {HistType::kTH2F, {{100, 0, 100}, {100, 0, 100}}});
registry.add("mMFTTrackAvgClustersTruHe", "Average number of clusters per track; p;# clusters; # entries", {HistType::kTH2F, {{100, 0, 100}, {100, 0, 100}}});
}
}
}
void process(aod::MFTTracks const& mfttracks)
{
for (auto& track : mfttracks) {
for (const auto& track : mfttracks) {
// 2d histograms
float x = track.x();
float y = track.y();
float eta = track.eta();
float phi = track.phi();
float nCls = track.nClusters();
if (nCls >= 5) {

Check failure on line 92 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
registry.fill(HIST("mMFTTrackXY_5_MinClusters"), x, y);
registry.fill(HIST("mMFTTrackEtaPhi_5_MinClusters"), eta, phi);
if (nCls >= 7) {
registry.fill(HIST("mMFTTrackXY_7_MinClusters"), x, y);
registry.fill(HIST("mMFTTrackEtaPhi_7_MinClusters"), eta, phi);
if (nCls >= 8) {
registry.fill(HIST("mMFTTrackXY_8_MinClusters"), x, y);
registry.fill(HIST("mMFTTrackEtaPhi_8_MinClusters"), eta, phi);
if (nCls >= 6) {

Check failure on line 95 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
registry.fill(HIST("mMFTTrackXY_6_MinClusters"), x, y);
registry.fill(HIST("mMFTTrackEtaPhi_6_MinClusters"), eta, phi);
if (nCls >= 7) {

Check failure on line 98 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
registry.fill(HIST("mMFTTrackXY_7_MinClusters"), x, y);
registry.fill(HIST("mMFTTrackEtaPhi_7_MinClusters"), eta, phi);
if (nCls >= 8) {

Check failure on line 101 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
registry.fill(HIST("mMFTTrackXY_8_MinClusters"), x, y);
registry.fill(HIST("mMFTTrackEtaPhi_8_MinClusters"), eta, phi);
}
}
}
}
if (avClsPlots) {
static constexpr int kNcls = 10;
std::array<float, kNcls> clsSize;
for (unsigned int layer = 0; layer < kNcls; layer++) {
clsSize[layer] = (track.mftClusterSizesAndTrackFlags() >> (layer * 6)) & 0x3f;
// LOG(info) << "Layer " << layer << ": " << clsSize[layer];
}
float avgCls = 0;
for (unsigned int layer = 0; layer < kNcls; layer++) {
avgCls += clsSize[layer];
}
avgCls /= track.nClusters();

std::sort(clsSize.begin(), clsSize.end());
float truncatedAvgCls = 0;
int ncls = 0;
for (unsigned int layer = 0; layer < kNcls; layer++) {
if (clsSize[layer] > 0) {
truncatedAvgCls += clsSize[layer];
ncls++;
if (ncls >= 3) {

Check failure on line 128 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
break; // we take the average of the first 5 non-zero clusters
}
}
}
truncatedAvgCls /= ncls;

registry.fill(HIST("mMFTTrackAvgClusters"), track.p(), avgCls);
registry.fill(HIST("mMFTTrackAvgClustersTru"), track.p(), truncatedAvgCls);
}
// 1d histograms
registry.fill(HIST("mMFTTrackEta"), eta);
Expand All @@ -82,6 +143,49 @@
registry.fill(HIST("mMFTTrackInvQPt"), track.signed1Pt());
}
}

void processMC(soa::Join<aod::MFTTracks, aod::McMFTTrackLabels> const& mfttracks,
aod::McParticles const&)
{
static constexpr int kNcls = 10;
for (const auto& track : mfttracks) {
if (avClsPlots) {
std::array<float, kNcls> clsSize;
for (unsigned int layer = 0; layer < kNcls; layer++) {
clsSize[layer] = (track.mftClusterSizesAndTrackFlags() >> (layer * 6)) & 0x3f;
// LOG(info) << "Layer " << layer << ": " << clsSize[layer];
}
float avgCls = 0;
for (unsigned int layer = 0; layer < kNcls; layer++) {
avgCls += clsSize[layer];
}
avgCls /= track.nClusters();

std::sort(clsSize.begin(), clsSize.end());
float truncatedAvgCls = 0;
int ncls = 0;
for (unsigned int layer = 0; layer < kNcls; layer++) {
if (clsSize[layer] > 0) {
truncatedAvgCls += clsSize[layer];
ncls++;
if (ncls >= 3) {

Check failure on line 171 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
break; // we take the average of the first 5 non-zero clusters
}
}
}
truncatedAvgCls /= ncls;

if (track.has_mcParticle()) {
const auto& mcParticle = track.mcParticle();
if (std::abs(mcParticle.pdgCode()) == 1000020040) { // He4

Check failure on line 180 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.

Check failure on line 180 in DPG/Tasks/MFT/aQCMFTTracks.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/explicit-code] 

Avoid hard-coded PDG codes. Use named values from PDG_t or o2::constants::physics::Pdg instead.
registry.fill(HIST("mMFTTrackAvgClustersHe"), track.p(), avgCls);
registry.fill(HIST("mMFTTrackAvgClustersTruHe"), track.p(), truncatedAvgCls);
}
}
}
}
}
PROCESS_SWITCH(CheckMFT, processMC, "Process MC", false);
};

WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
Expand Down
Loading