Skip to content

[DPG] Update lightIonsEvSelQa.cxx #13447

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
altsybee merged 2 commits into from
Oct 19, 2025
Merged

[DPG] Update lightIonsEvSelQa.cxx #13447

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
3648f73
Update lightIonsEvSelQa.cxx
altsybee Oct 18, 2025
9146328
Update lightIonsEvSelQa.cxx
altsybee Oct 18, 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
146 changes: 142 additions & 4 deletions DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -56,7 +56,7 @@
Configurable<int> nBinsTracks{"nBinsTracks", 450, "N bins in n tracks histo"}; // o2-linter: disable=name/configurable (temporary fix)
Configurable<int> nMaxTracks{"nMaxTracks", 450, "N max in n tracks histo"}; // o2-linter: disable=name/configurable (temporary fix)
Configurable<int> nMaxGlobalTracks{"nMaxGlobalTracks", 450, "N max in n tracks histo"}; // o2-linter: disable=name/configurable (temporary fix)
Configurable<int> nBinsMultFwd{"nBinsMultFwd", 1000, "N bins in mult fwd histo"}; // o2-linter: disable=name/configurable (temporary fix)
Configurable<int> nBinsMultFwd{"nBinsMultFwd", 800, "N bins in mult fwd histo"}; // o2-linter: disable=name/configurable (temporary fix)
Configurable<float> nMaxMultFwd{"nMaxMultFwd", 100000, "N max in mult fwd histo"}; // o2-linter: disable=name/configurable (temporary fix)
Configurable<float> timeBinWidthInSec{"TimeBinWidthInSec", 10, "Width of time bins in seconds"}; // o2-linter: disable=name/configurable (temporary fix)

Expand All @@ -65,6 +65,8 @@

Configurable<int> confUseDiffVzCutFromEvSel{"UseDiffVzCutFromEvSel", 0, "0 - custom diffVz cut from this task, 1 - cut from event selection"}; // o2-linter: disable=name/configurable (temporary fix)

Configurable<bool> isMC{"isMC", false, "Run MC"};

uint64_t minGlobalBC = 0;
Service<o2::ccdb::BasicCCDBManager> ccdb;
HistogramRegistry histos{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject};
Expand Down Expand Up @@ -137,12 +139,16 @@
const AxisSpec axisVtxZ{800, -20., 20., ""};
const AxisSpec axisBcDiff{600, -300., 300., "bc difference"};

const AxisSpec axisNcontrib{801, -0.5, 800.5, "n contributors"};
const AxisSpec axisNcontrib{601, -0.5, 600.5, "n contributors"};
const AxisSpec axisColTimeRes{1500, 0., 1500., "collision time resolution (ns)"};

AxisSpec axisVertexChi2{100, 0, 500, "Chi2 of vertex fit"};
AxisSpec axisVertexChi2perContrib{100, 0, 10, "Chi2 of vertex fit"};

const AxisSpec axisTimeDiff{100, -10., 10., ""};
const AxisSpec axisTimeSum{100, -10., 10., ""};
const AxisSpec axisZNampl{200, 0, 3000, ""};

histos.add("noSpecSelections/hBcColNoSel8", "", kTH1F, {axisBCs});
histos.add("noSpecSelections/hBcOrigNoSel8", "", kTH1F, {axisBCs});
// histos.add("noSpecSelections/hBcColNoSel8TOF", "", kTH1F, {axisBCs});
Expand All @@ -165,6 +171,7 @@
histos.add("noSpecSelections/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("noSpecSelections/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
histos.add("noSpecSelections/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
histos.add("noSpecSelections/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});

histos.add("noPU/hBcColNoSel8", "", kTH1F, {axisBCs});
histos.add("noPU/hBcOrigNoSel8", "", kTH1F, {axisBCs});
Expand All @@ -187,6 +194,8 @@
histos.add("noPU/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("noPU/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
histos.add("noPU/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
histos.add("noPU/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
histos.add("noPU/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});

histos.add("noPU_pvTOFmatched/hBcColNoSel8", "", kTH1F, {axisBCs});
histos.add("noPU_pvTOFmatched/hBcTVX", "", kTH1F, {axisBCs});
Expand Down Expand Up @@ -394,6 +403,7 @@
histos.add("noPU_lowMultCut/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("noPU_lowMultCut/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
histos.add("noPU_lowMultCut/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
histos.add("noPU_lowMultCut/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});

histos.add("highMultCloudCut/hBcColNoSel8", "", kTH1F, {axisBCs});
histos.add("highMultCloudCut/hBcOrigNoSel8", "", kTH1F, {axisBCs});
Expand All @@ -420,6 +430,7 @@
histos.add("noPU_highMultCloudCut/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("noPU_highMultCloudCut/hVertexChi2VsNcontrib", "", kTH2F, {axisNcontrib, axisVertexChi2perContrib});
histos.add("noPU_highMultCloudCut/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
histos.add("noPU_highMultCloudCut/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});

histos.add("noPU_badVzDiff/hBcColNoSel8", "", kTH1F, {axisBCs});
histos.add("noPU_badVzDiff/hBcTVX", "", kTH1F, {axisBCs});
Expand Down Expand Up @@ -559,6 +570,7 @@
histos.add("noPU_cutByVzDiff_pvTOF/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
histos.add("noPU_cutByVzDiff_pvTOF/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("noPU_cutByVzDiff_pvTOF/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
histos.add("noPU_cutByVzDiff_pvTOF/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});

histos.add("noPU_cutByVzDiff_noFT0activityNearby/hBcColNoSel8", "", kTH1F, {axisBCs});
histos.add("noPU_cutByVzDiff_noFT0activityNearby/hBcTVX", "", kTH1F, {axisBCs});
Expand All @@ -579,8 +591,46 @@
histos.add("noPU_cutByVzDiff_noFT0activityNearby/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("noPU_cutByVzDiff_noFT0activityNearby/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});

histos.add("noPU_CutOnZNACtime/hBcColNoSel8", "", kTH1F, {axisBCs});
histos.add("noPU_CutOnZNACtime/hBcTVX", "", kTH1F, {axisBCs});
histos.add("noPU_CutOnZNACtime/hBcOrig", "", kTH1F, {axisBCs});
histos.add("noPU_CutOnZNACtime/hBcFT0", "", kTH1F, {axisBCs});
histos.add("noPU_CutOnZNACtime/hBcFV0", "", kTH1F, {axisBCs});
histos.add("noPU_CutOnZNACtime/hBcZDC", "", kTH1F, {axisBCs});
histos.add("noPU_CutOnZNACtime/hVtxFT0VsVtxCol", "", kTH2F, {axisVtxZ, axisVtxZ});
histos.add("noPU_CutOnZNACtime/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
histos.add("noPU_CutOnZNACtime/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
histos.add("noPU_CutOnZNACtime/nTracksPV_vs_T0A", "", kTH2F, {axisMultFT0A, axisNtracks});
histos.add("noPU_CutOnZNACtime/nTracksPV_vs_T0C", "", kTH2F, {axisMultFT0C, axisNtracks});
histos.add("noPU_CutOnZNACtime/hTVXvsBcDiffwrtOrigBc", "", kTH1F, {axisBcDiff});
histos.add("noPU_CutOnZNACtime/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
histos.add("noPU_CutOnZNACtime/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("noPU_CutOnZNACtime/hNPVvsNch", "", kTH2F, {axisNcontrib, axisNcontrib});
histos.add("noPU_CutOnZNACtime/hTimeZAC", "ZNC-ZNA time (ns); ZNC+ZNA time (ns)", kTH2F, {axisTimeDiff, axisTimeSum});
histos.add("noPU_CutOnZNACtime/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});

histos.add("noPU_AntiCutOnZNACtime/hBcFV0", "", kTH1F, {axisBCs});
histos.add("noPU_AntiCutOnZNACtime/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
histos.add("noPU_AntiCutOnZNACtime/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
histos.add("noPU_AntiCutOnZNACtime/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});

histos.add("noPU_AntiCutOnZNAampl/hBcFV0", "", kTH1F, {axisBCs});
histos.add("noPU_AntiCutOnZNAampl/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
histos.add("noPU_AntiCutOnZNAampl/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
histos.add("noPU_AntiCutOnZNAampl/hAmplZAC", "ZNC vs ZNA", kTH2F, {axisZNampl, axisZNampl});

//
histos.add("hNcontribColFromData", "", kTH1F, {axisNcontrib});
histos.add("hNcontribAccFromData", "", kTH1F, {axisNcontrib});

// MC QA
const AxisSpec axisVtxZdiff{400, -10., 10., ""};
histos.add("MC/hMCdataVzDiff", "", kTH2F, {axisNcontrib, axisVtxZdiff});
histos.add("MC/hMCdataBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("MC/hMCdataFoundBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("MCsel8/hMCdataVzDiff", "", kTH2F, {axisNcontrib, axisVtxZdiff});
histos.add("MCsel8/hMCdataBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
histos.add("MCsel8/hMCdataFoundBcDiffVsMult", "", kTH2F, {axisNcontrib, axisBcDiff});
}

Preslice<FullTracksIU> perCollision = aod::track::collisionId;
Expand All @@ -604,6 +654,7 @@
ColEvSels const& cols,
FullTracksIU const& tracks,
BCsRun3 const& bcs,
aod::Zdcs const&,
aod::FT0s const&)
{
int runNumber = bcs.iteratorAt(0).runNumber();
Expand All @@ -613,7 +664,7 @@
int64_t tsSOR = 0; // dummy start-of-run timestamp
int64_t tsEOR = 1; // dummy end-of-run timestamp

if (runNumber >= 500000) {

Check failure on line 667 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
auto runInfo = o2::parameters::AggregatedRunInfo::buildAggregatedRunInfo(o2::ccdb::BasicCCDBManager::instance(), runNumber);
// first bc of the first orbit
bcSOR = runInfo.orbitSOR * o2::constants::lhc::LHCMaxBunches;
Expand Down Expand Up @@ -710,7 +761,7 @@
pastActivityFV0 |= bcPast.has_fv0a();
// pastActivityFDD |= bcPast.has_fdd();
}
if (deltaBC < 2) {

Check failure on line 764 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
if (bcPast.has_ft0()) {
std::bitset<8> triggers = bcPast.ft0().triggerMask();
nearbyFT0activity |= (triggers[o2::ft0::RecPoints::ETriggerBits::kIsActiveSideA] || triggers[o2::ft0::RecPoints::ETriggerBits::kIsActiveSideC]);
Expand Down Expand Up @@ -739,7 +790,7 @@
futureActivityFV0 |= bcFuture.has_fv0a();
futureActivityFDD |= bcFuture.has_fdd();
}
if (deltaBC < 2) {

Check failure on line 793 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
if (bcFuture.has_ft0()) {
std::bitset<8> triggers = bcFuture.ft0().triggerMask();
nearbyFT0activity |= (triggers[o2::ft0::RecPoints::ETriggerBits::kIsActiveSideA] || triggers[o2::ft0::RecPoints::ETriggerBits::kIsActiveSideC]);
Expand Down Expand Up @@ -817,7 +868,7 @@

// ### collision loop
for (const auto& col : cols) {
if (std::abs(col.posZ()) > 10)

Check failure on line 871 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
continue;

const auto& foundBC = col.foundBC_as<BCsRun3>();
Expand Down Expand Up @@ -853,21 +904,21 @@
// int nTOFtracks = 0;
auto tracksGrouped = tracks.sliceBy(perCollision, col.globalIndex());
for (const auto& track : tracksGrouped) {
if (track.itsNCls() < 5)

Check failure on line 907 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
continue;
if (track.pt() < 0.2 || track.pt() > 10)

Check failure on line 909 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
continue;
if (std::abs(track.eta()) > 0.8)

Check failure on line 911 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
continue;

if (track.hasITS() && track.hasTPC() && track.tpcNClsFound() > 50 && track.tpcNClsCrossedRows() > 50 && track.tpcChi2NCl() < 4)

Check failure on line 914 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
nGlobalTracksAll++;

if (!track.isPVContributor()) {
continue;
}
nPVtracks++;
if (track.hasITS() && track.hasTPC() && track.tpcNClsFound() > 50 && track.tpcNClsCrossedRows() > 50 && track.tpcChi2NCl() < 4)

Check failure on line 921 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.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.
nGlobalTracksPV++;
} // end of track loop

Expand Down Expand Up @@ -904,6 +955,17 @@
histos.fill(HIST("multT0M_vs_multSumFT0"), multT0M, multT0M_uncorr);
}

float timeZNA = foundBC.has_zdc() ? foundBC.zdc().timeZNA() : -999.f;
float timeZNC = foundBC.has_zdc() ? foundBC.zdc().timeZNC() : -999.f;
float znSum = timeZNA + timeZNC;
float znDiff = timeZNA - timeZNC;
// bool goodZNACtime = fabs(znSum) < 2 && fabs(znDiff) < 2;
bool goodZNACtime = (timeZNA > -6 && timeZNA < 4) && (timeZNC > -6 && timeZNC < 4);

float multZNA = foundBC.has_zdc() ? foundBC.zdc().energyCommonZNA() : -999;
float multZNC = foundBC.has_zdc() ? foundBC.zdc().energyCommonZNC() : -999;
bool cutZNACampl = multZNA < 400 && multZNC < 400;

// vZ diff (FT0 vs by tracks)
bool badVzDiff = 0;
if (confUseDiffVzCutFromEvSel)
Expand All @@ -925,7 +987,7 @@
if (lastRunNumber == 564468)
meanDiff = -0.60;

float stdDev = (multT0M > 10) ? 0.144723 + 13.5345 / sqrt(multT0M) : 1.5; // cm

Check failure on line 990 in DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix] 

Use std:: prefix for names from the std namespace.
if (multT0M > 5000)
stdDev = stdDev > 0.2 ? stdDev : 0.2; // 0.35; // cm
badVzDiff = diffVz < (meanDiff - stdDev * nSigmaForVzDiff - safetyDiffVzMargin) || diffVz > (meanDiff + stdDev * nSigmaForVzDiff + safetyDiffVzMargin);
Expand Down Expand Up @@ -968,11 +1030,14 @@
histos.fill(HIST("noSpecSelections/hBcColNoSel8"), localBC);
histos.fill(HIST("noSpecSelections/hBcOrigNoSel8"), bcOriginal);
histos.fill(HIST("noSpecSelections/Vz"), vZ);
histos.fill(HIST("noSpecSelections/hTimeZAC"), znDiff, znSum);

if (noPU) {
histos.fill(HIST("noPU/hBcColNoSel8"), localBC);
histos.fill(HIST("noPU/hBcOrigNoSel8"), bcOriginal);
histos.fill(HIST("noPU/Vz"), vZ);
histos.fill(HIST("noPU/hTimeZAC"), znDiff, znSum);
histos.fill(HIST("noPU/hAmplZAC"), multZNA, multZNC);
}
if (noPU && pvTOFmatched) {
histos.fill(HIST("noPU_pvTOFmatched/hBcColNoSel8"), localBC);
Expand Down Expand Up @@ -1028,6 +1093,7 @@
if (noPU && underLine) {
histos.fill(HIST("noPU_lowMultCut/hBcColNoSel8"), localBC);
histos.fill(HIST("noPU_lowMultCut/hBcOrigNoSel8"), bcOriginal);
histos.fill(HIST("noPU_lowMultCut/hTimeZAC"), znDiff, znSum);
}
if (grassOnTheRight) {
histos.fill(HIST("highMultCloudCut/hBcColNoSel8"), localBC);
Expand All @@ -1036,9 +1102,11 @@
if (noPU && grassOnTheRight) {
histos.fill(HIST("noPU_highMultCloudCut/hBcColNoSel8"), localBC);
histos.fill(HIST("noPU_highMultCloudCut/hBcOrigNoSel8"), bcOriginal);
histos.fill(HIST("noPU_highMultCloudCut/hTimeZAC"), znDiff, znSum);
}
if (noPU && !badVzDiff && pvTOFmatched) { // noPileup_cutByVzDiff_pvTOF_noFT0act
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hBcColNoSel8"), localBC);
histos.fill(HIST("noPU_cutByVzDiff_pvTOF/hTimeZAC"), znDiff, znSum);
}

// only here cut on sel8:
Expand Down Expand Up @@ -1224,6 +1292,14 @@
histos.fill(HIST("noPU_cutByVzDiff_noFT0activityNearby/hColTimeResVsNcontrib"), nContributors, timeRes);
histos.fill(HIST("noPU_cutByVzDiff_noFT0activityNearby/hNPVvsNch"), nPVtracks, nGlobalTracksAll);
}
if (noPU && goodZNACtime) {
histos.fill(HIST("noPU_CutOnZNACtime/hBcTVX"), localBC);
histos.fill(HIST("noPU_CutOnZNACtime/hBcOrig"), bcOriginal);
histos.fill(HIST("noPU_CutOnZNACtime/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
histos.fill(HIST("noPU_CutOnZNACtime/hColTimeResVsNcontrib"), nContributors, timeRes);
histos.fill(HIST("noPU_CutOnZNACtime/hNPVvsNch"), nPVtracks, nGlobalTracksAll);
histos.fill(HIST("noPU_CutOnZNACtime/hTimeZAC"), znDiff, znSum);
}

if (foundBC.has_ft0()) {
// float multT0A = foundBC.has_ft0() ? foundBC.ft0().sumAmpA() : -999.f;
Expand Down Expand Up @@ -1417,6 +1493,14 @@
histos.fill(HIST("noPU_cutByVzDiff_noFT0activityNearby/nTracksPV_vs_T0C"), multT0C, nPVtracks);
// histos.fill(HIST("noPU_cutByVzDiff_noFT0activityNearby/nTracksGlobal_vs_T0C"), multT0C, nGlobalTracksPV);
}
if (noPU && goodZNACtime) {
histos.fill(HIST("noPU_CutOnZNACtime/hBcFT0"), localBC);
histos.fill(HIST("noPU_CutOnZNACtime/hVtxFT0VsVtxCol"), vZft0, vZ);
histos.fill(HIST("noPU_CutOnZNACtime/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
histos.fill(HIST("noPU_CutOnZNACtime/nTracksPV_vs_T0A"), multT0A, nPVtracks);
histos.fill(HIST("noPU_CutOnZNACtime/nTracksPV_vs_T0C"), multT0C, nPVtracks);
histos.fill(HIST("noPU_CutOnZNACtime/hAmplZAC"), multZNA, multZNC);
}
}

if (foundBC.has_fv0a()) {
Expand Down Expand Up @@ -1534,6 +1618,23 @@
histos.fill(HIST("noPU_cutByVzDiff_noFT0activityNearby/nTracksPV_vs_V0A"), multV0A, nPVtracks);
// histos.fill(HIST("noPU_cutByVzDiff_noFT0activityNearby/nTracksGlobal_vs_V0A"), multV0A, nGlobalTracksPV);
}

if (noPU && goodZNACtime) {
histos.fill(HIST("noPU_CutOnZNACtime/hBcFV0"), localBC);
histos.fill(HIST("noPU_CutOnZNACtime/nTracksPV_vs_V0A"), multV0A, nPVtracks);
}
if (noPU && !goodZNACtime) {
histos.fill(HIST("noPU_AntiCutOnZNACtime/hBcFV0"), localBC);
histos.fill(HIST("noPU_AntiCutOnZNACtime/nTracksPV_vs_V0A"), multV0A, nPVtracks);
histos.fill(HIST("noPU_AntiCutOnZNACtime/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
histos.fill(HIST("noPU_AntiCutOnZNACtime/hAmplZAC"), multZNA, multZNC);
}
if (noPU && !cutZNACampl) {
histos.fill(HIST("noPU_AntiCutOnZNAampl/hBcFV0"), localBC);
histos.fill(HIST("noPU_AntiCutOnZNAampl/nTracksPV_vs_V0A"), multV0A, nPVtracks);
histos.fill(HIST("noPU_AntiCutOnZNAampl/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
histos.fill(HIST("noPU_AntiCutOnZNAampl/hAmplZAC"), multZNA, multZNC);
}
}
if (foundBC.has_zdc()) {
histos.fill(HIST("noSpecSelections/hBcZDC"), localBC);
Expand Down Expand Up @@ -1661,13 +1762,50 @@
if (noPU && !badVzDiff && noFT0activityNearby) {
histos.fill(HIST("noPU_cutByVzDiff_noFT0activityNearby/hTVXvsBcDiffwrtOrigBc"), bcDiffWrtOriginal);
}

if (noPU && goodZNACtime) {
histos.fill(HIST("noPU_CutOnZNACtime/hTVXvsBcDiffwrtOrigBc"), bcDiffWrtOriginal);
}
} // end of collisions loop
}
PROCESS_SWITCH(LightIonsEvSelQa, processRun3, "Process Run3 tracking vs detector occupancy QA", true);

// ### MC QA
using ColEvSelsWithMCLabels = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels>; //, aod::CentFT0Cs>;
using BCsInfo = soa::Join<aod::BCs, aod::Timestamps, aod::BcSels>;

void processMC(ColEvSelsWithMCLabels const& collisions,
BCsInfo const&,
aod::McCollisions const&)
{
for (const auto& col : collisions) {
bool isSel8 = col.sel8();
if (col.has_mcCollision()) {
const auto mcCollision = col.mcCollision();
LOGP(debug, "col.posZ() = {}, mcCollision.posZ() = {}", col.posZ(), mcCollision.posZ());
histos.fill(HIST("MC/hMCdataVzDiff"), col.numContrib(), col.posZ() - mcCollision.posZ());

auto bc = col.bc_as<BCsInfo>();
auto mcBc = mcCollision.bc_as<BCsInfo>();
auto foundBC = col.foundBC_as<BCsInfo>();

LOGP(debug, "col.bc().globalBC() = {}, mcCollision.bc().globalBC() = {}", col.bc().globalBC(), mcCollision.bc().globalBC());
histos.fill(HIST("MC/hMCdataBcDiffVsMult"), col.numContrib(), bc.globalBC() - mcBc.globalBC());
histos.fill(HIST("MC/hMCdataFoundBcDiffVsMult"), col.numContrib(), foundBC.globalBC() - mcBc.globalBC());

if (isSel8) {
histos.fill(HIST("MCsel8/hMCdataVzDiff"), col.numContrib(), col.posZ() - mcCollision.posZ());
histos.fill(HIST("MCsel8/hMCdataBcDiffVsMult"), col.numContrib(), bc.globalBC() - mcBc.globalBC());
histos.fill(HIST("MCsel8/hMCdataFoundBcDiffVsMult"), col.numContrib(), foundBC.globalBC() - mcBc.globalBC());
}
}
}
}

PROCESS_SWITCH(LightIonsEvSelQa, processMC, "Process MC", false);
};

WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
WorkflowSpec
defineDataProcessing(ConfigContext const& cfgc)
{
return WorkflowSpec{
adaptAnalysisTask<LightIonsEvSelQa>(cfgc)};
Expand Down
Loading