Skip to content

[PWGCF] Add MFT-MFT correlation task #18485

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
mintu1999 wants to merge 9 commits into
base: master
Choose a base branch
from
Open

[PWGCF] Add MFT-MFT correlation task #18485

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
686b5c5
[PWGCF] Add MFT-MFT correlation task
mintu1999 Oct 27, 2026
bd80453
Fix Whitespace
mintu1999 Oct 28, 2026
2dd63c4
Please consider the following formatting changes
alibuild Oct 28, 2026
ec7fb65
Merge pull request #3 from alibuild/alibot-cleanup-18485
mintu1999 Oct 28, 2026
5b603ea
Merge branch 'master' into mftcorrelation
mintu1999 Oct 29, 2026
c996701
Please consider the following formatting changes
alibuild Oct 29, 2026
f76c37d
Merge pull request #4 from alibuild/alibot-cleanup-18485
mintu1999 Oct 29, 2026
bd8dd6b
Fix Megalinter cpplint C-style casts in CorrelationContainer
mintu1999 Oct 29, 2026
aeaf83e
Fix Megalint cppcheck cstyleCast and uninitialized member warnings
mintu1999 Oct 30, 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
32 changes: 18 additions & 14 deletions PWGCF/Core/CorrelationContainer.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,11 @@
#include <TCanvas.h>
#include <TCollection.h>
#include <TF1.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TH3D.h>
#include <THn.h>
#include <THnSparse.h>
#include <TIterator.h>
#include <TList.h>
#include <TMath.h>
Expand Down Expand Up @@ -125,7 +129,7 @@ CorrelationContainer::CorrelationContainer(const char* name, const char* objTitl
triggerAxis.insert(triggerAxis.end(), userAxis.begin(), userAxis.end());
mTriggerHist = HistFactory::createHist<StepTHnF>({"mTriggerHist", "d^{2}N_{ch}/d#varphid#eta", {HistType::kStepTHnF, triggerAxis, fgkCFSteps}}).release();

mTrackHistEfficiency = HistFactory::createHist<StepTHnF>({"mTrackHistEfficiency", "Tracking efficiency", {HistType::kStepTHnF, {efficiencyAxis[0], efficiencyAxis[1], {5, -0.5, 4.5, "species"}, correlationAxis[3], efficiencyAxis[2]}, fgkCFSteps}}).release();
mTrackHistEfficiency = HistFactory::createHist<StepTHnF>({"mTrackHistEfficiency", "Tracking efficiency", {HistType::kStepTHnF, {efficiencyAxis[0], efficiencyAxis[1], {5, -0.5, 4.5, "species"}, correlationAxis[3], efficiencyAxis[2], efficiencyAxis[3]}, fgkCFSteps}}).release();

mEventCount = HistFactory::createHist<TH2F>({"mEventCount", ";step;centrality;count", {HistType::kTH2F, {{fgkCFSteps + 2, -2.5, -0.5 + fgkCFSteps, "step"}, correlationAxis[3]}}}).release();
}
Expand Down Expand Up @@ -155,7 +159,7 @@ CorrelationContainer::CorrelationContainer(const CorrelationContainer& c) : TNam
// CorrelationContainer copy constructor
//

((CorrelationContainer&)c).Copy(*this);
c.Copy(*this);
}

//____________________________________________________________________
Expand Down Expand Up @@ -195,7 +199,7 @@ CorrelationContainer& CorrelationContainer::operator=(const CorrelationContainer
// assigment operator

if (this != &c) {
((CorrelationContainer&)c).Copy(*this);
c.Copy(*this);
}

return *this;
Expand All @@ -206,7 +210,7 @@ void CorrelationContainer::Copy(TObject& c) const
{
// copy function

CorrelationContainer& target = (CorrelationContainer&)c;
CorrelationContainer& target = static_cast<CorrelationContainer&>(c);

if (mPairHist) {
target.mPairHist = dynamic_cast<StepTHn*>(mPairHist->Clone());
Expand Down Expand Up @@ -423,7 +427,7 @@ void CorrelationContainer::getHistsZVtxMult(CorrelationContainer::CFStep step, F

Int_t dimensions[] = {4, 0, 5, 3};
THnBase* tmpTrackHist = sparse->ProjectionND(4, dimensions, "E");
*eventHist = (TH2*)mTriggerHist->getTHn(step)->Projection(1, 2); // NOTE the syntax is Projection(Y, X) --> produces a histogram where x = axis 2 (vertex) and y = axis 1 (multiplicity)
*eventHist = static_cast<TH2*>(mTriggerHist->getTHn(step)->Projection(1, 2)); // NOTE the syntax is Projection(Y, X) --> produces a histogram where x = axis 2 (vertex) and y = axis 1 (multiplicity)
// convert to THn
*trackHist = changeToThn(tmpTrackHist);
delete tmpTrackHist;
Expand Down Expand Up @@ -728,7 +732,7 @@ TH2* CorrelationContainer::getSumOfRatios(CorrelationContainer* mixed, Correlati
*/

if (!totalTracks) {
totalTracks = (TH2*)tracksSame->Clone("totalTracks");
totalTracks = static_cast<TH2*>(tracksSame->Clone("totalTracks"));
} else {
totalTracks->Add(tracksSame);
}
Expand Down Expand Up @@ -836,7 +840,7 @@ THnBase* CorrelationContainer::getTrackEfficiencyND(EfficiencyStep step1, Effici
resetBinLimits(sourceContainer->getTHn(step1));
resetBinLimits(sourceContainer->getTHn(step2));

THnBase* clone = (THnBase*)measured->Clone();
THnBase* clone = static_cast<THnBase*>(measured->Clone());

clone->Divide(measured, generated, 1, 1, "B");

Expand Down Expand Up @@ -1010,7 +1014,7 @@ TH1* CorrelationContainer::getTrackEfficiency(EfficiencyStep step1, EfficiencySt
measured = new TH2D(Form("%s_rebinned", tmp->GetName()), tmp->GetTitle(), axis->GetNbins(), axis->GetXbins()->GetArray(), tmp->GetNbinsY(), tmp->GetYaxis()->GetXbins()->GetArray());
for (Int_t x = 1; x <= tmp->GetNbinsX(); x++) {
for (Int_t y = 1; y <= tmp->GetNbinsY(); y++) {
((TH2*)measured)->Fill(tmp->GetXaxis()->GetBinCenter(x), tmp->GetYaxis()->GetBinCenter(y), tmp->GetBinContent(x, y));
static_cast<TH2*>(measured)->Fill(tmp->GetXaxis()->GetBinCenter(x), tmp->GetYaxis()->GetBinCenter(y), tmp->GetBinContent(x, y));
measured->SetBinError(x, y, 0); // cannot trust bin error, set to 0
}
}
Expand All @@ -1020,7 +1024,7 @@ TH1* CorrelationContainer::getTrackEfficiency(EfficiencyStep step1, EfficiencySt
generated = new TH2D(Form("%s_rebinned", tmp->GetName()), tmp->GetTitle(), axis->GetNbins(), axis->GetXbins()->GetArray(), tmp->GetNbinsY(), tmp->GetYaxis()->GetXbins()->GetArray());
for (Int_t x = 1; x <= tmp->GetNbinsX(); x++) {
for (Int_t y = 1; y <= tmp->GetNbinsY(); y++) {
((TH2*)generated)->Fill(tmp->GetXaxis()->GetBinCenter(x), tmp->GetYaxis()->GetBinCenter(y), tmp->GetBinContent(x, y));
static_cast<TH2*>(generated)->Fill(tmp->GetXaxis()->GetBinCenter(x), tmp->GetYaxis()->GetBinCenter(y), tmp->GetBinContent(x, y));
generated->SetBinError(x, y, 0); // cannot trust bin error, set to 0
}
}
Expand Down Expand Up @@ -1170,14 +1174,14 @@ TH1* CorrelationContainer::getBias(EfficiencyStep step1, EfficiencyStep step2, c
resetBinLimits(mTriggerHist->getTHn(step2), 3);
setBinLimits(tmp->getTHn(step2));

TH1D* events1 = (TH1D*)mTriggerHist->getTHn(step1)->Projection(0);
TH3D* hist1 = (TH3D*)tmp->getTHn(step1)->Projection(0, tmp->getNVar() - 1, 2);
TH1D* events1 = static_cast<TH1D*>(mTriggerHist->getTHn(step1)->Projection(0));
TH3D* hist1 = static_cast<TH3D*>(tmp->getTHn(step1)->Projection(0, tmp->getNVar() - 1, 2));
if (weighting == 0) {
weightHistogram(hist1, events1);
}

TH1D* events2 = (TH1D*)mTriggerHist->getTHn(step2)->Projection(0);
TH3D* hist2 = (TH3D*)tmp->getTHn(step2)->Projection(0, tmp->getNVar() - 1, 2);
TH1D* events2 = static_cast<TH1D*>(mTriggerHist->getTHn(step2)->Projection(0));
TH3D* hist2 = static_cast<TH3D*>(tmp->getTHn(step2)->Projection(0, tmp->getNVar() - 1, 2));
if (weighting == 0) {
weightHistogram(hist2, events2);
}
Expand Down Expand Up @@ -1430,7 +1434,7 @@ void CorrelationContainer::symmetrizepTBins()
}

// for symmetric bins
THnBase* source = (THnBase*)target->Clone();
THnBase* source = static_cast<THnBase*>(target->Clone());

// axes: 0 delta eta; 1 pT,a; 2 pT,t; 3 centrality; 4 delta phi; 5 vtx-z
for (Int_t i3 = 1; i3 <= target->GetAxis(3)->GetNbins(); i3++) {
Expand Down
5 changes: 5 additions & 0 deletions PWGCF/TwoParticleCorrelations/Tasks/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -78,6 +78,11 @@ o2physics_add_dpl_workflow(long-range-dihadron-cor
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::PWGCFCore O2Physics::AnalysisCCDB O2Physics::GFWCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(mft-correlations
SOURCES mftcorrelations.cxx
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::PWGCFCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(two-particle-correlation-pp
SOURCES twoParticleCorrelationPp.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::PWGCFCore
Expand Down
Loading
Loading