Skip to content

[PWGLF,PWGMM] Adding bc idx and ZN amplitudes to vdm derived data #14283

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
mpuccio merged 3 commits into from
Dec 19, 2025
Merged

[PWGLF,PWGMM] Adding bc idx and ZN amplitudes to vdm derived data #14283

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
86ed94c
[PWGMM] Adding bc idx and ZN amplitudes to vdm derived data
Dec 18, 2025
795b788
Please consider the following formatting changes
alibuild Dec 18, 2025
6db3b0c
Merge pull request #4 from alibuild/alibot-cleanup-14283
lmattei01 Dec 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
20 changes: 14 additions & 6 deletions PWGMM/Lumi/Tasks/lumiStabilityLightIons.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,10 +43,13 @@
namespace myBc_aod
{
DECLARE_SOA_COLUMN(Timestamp, timestamp, uint64_t);
DECLARE_SOA_COLUMN(BCid, bcId, int);

Check failure on line 46 in PWGMM/Lumi/Tasks/lumiStabilityLightIons.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-column] 

Use UpperCamelCase for names of O2 columns and matching lowerCamelCase names for their getters.
DECLARE_SOA_COLUMN(TimeZNA, timeZNA, float);
DECLARE_SOA_COLUMN(TimeZNC, timeZNC, float);
DECLARE_SOA_COLUMN(AmplitudeZNA, amplitudeZNA, float);
DECLARE_SOA_COLUMN(AmplitudeZNC, amplitudeZNC, float);
} // namespace myBc_aod
DECLARE_SOA_TABLE(MyBCaod, "AOD", "MYBCAOD", myBc_aod::Timestamp, myBc_aod::TimeZNA, myBc_aod::TimeZNC);
DECLARE_SOA_TABLE(MyBCaod, "AOD", "MYBCAOD", myBc_aod::Timestamp, myBc_aod::BCid, myBc_aod::TimeZNA, myBc_aod::TimeZNC, myBc_aod::AmplitudeZNA, myBc_aod::AmplitudeZNC);
} // namespace o2::aod

using MyBCs = soa::Join<aod::BCs, aod::BcSels, aod::Timestamps, aod::Run3MatchedToBCSparse>;
Expand All @@ -66,8 +69,8 @@
Configurable<bool> cfgDoBCL{"cfgDoBCL", false, "Create and fill histograms for leading BCs of type B"};
Configurable<bool> cfgDoBCSL{"cfgDoBCSL", false, "Create and fill histograms for super-leading BCs (no preceding FT0/FDD activity) of type B"};

Configurable<bool> cfgRequireZDCTriggerForZDCQA{"cfgRequireZDCTriggerForZDCQA", false, "Require ZDC trigger (1ZNC or 1ZNA) for filling QA histograms"};
Configurable<bool> cfgRequireTVXTriggerForZDCQA{"cfgRequireTVXTriggerForZDCQA", false, "Require FT0 vertex trigger for filling ZDC QA histograms"};
Configurable<bool> cfgRequireZDCTriggerForZDCQA{"cfgRequireZDCTriggerForZDCQA", false, "Require ZDC trigger (1ZNC) for filling QA histograms"};
Configurable<bool> cfgRequireTVXTriggerForZDCQA{"cfgRequireTVXTriggerForZDCQA", false, "Require FT0 vertex trigger (MTVX) for filling ZDC QA histograms"};

Configurable<bool> cfgRequireNoT0ForSLBC{"cfgRequireNoT0ForSLBC", false, "Require no T0 signal for definition of super leading BC (otherwise only no FDD)"};

Expand Down Expand Up @@ -253,9 +256,9 @@
for (const auto& bc : bcs) {

std::bitset<64> ctpInputMask(bc.inputMask());
if (cfgRequireTVXTriggerForZDCQA && !(ctpInputMask.test(2)))
if (cfgRequireTVXTriggerForZDCQA && !(ctpInputMask.test(2))) // 2 = 3 - 1 -> MTVX
continue;
if (cfgRequireZDCTriggerForZDCQA && !(ctpInputMask.test(25)))
if (cfgRequireZDCTriggerForZDCQA && !(ctpInputMask.test(25))) // 25 = 26 - 1 -> 1ZNC
continue;

bool zdcHit = !bc.has_zdc() ? 0 : ((bc.zdc().energyCommonZNC() > -1 && std::abs(bc.zdc().timeZNC()) < 1E5) ? 1 : 0);
Expand All @@ -281,9 +284,14 @@
mHistManager.fill(HIST("ZDCQA/ZDCTimes"), timeZNA, timeZNC);

// For VdM analysis: fill timestamps and ZDC times in output tree, if enabled
// Fill BC idx and ZNA and ZNC amplitudes as well
uint64_t timestamp = bc.timestamp();
int64_t globalBC = bc.globalBC();
int localBC = globalBC % nBCsPerOrbit;
float amplitudeZNA = bc.zdc().amplitudeZNA();
float amplitudeZNC = bc.zdc().amplitudeZNC();
if (cfgFillBCao2d && timestamp >= cfgTstampStartFillingBCao2d && timestamp <= cfgTstampEndFillingBCao2d) {
BCaod(timestamp, timeZNA, timeZNC);
BCaod(timestamp, localBC, timeZNA, timeZNC, amplitudeZNA, amplitudeZNC);
}
}
}
Expand Down
Loading