Skip to content

[PWGMM] : Lumi. New definition of leading BC #14743

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
jgcn merged 23 commits into from
Feb 2, 2026
+30 −68
Merged

[PWGMM] : Lumi. New definition of leading BC #14743

Changes from all commits
Commits
Show all changes
23 commits
Select commit Hold shift + click to select a range
9e17f61
new changes
JosueMtzGar May 20, 2024
3670d73
new changes
JosueMtzGar Aug 15, 2024
8255956
Ordering the kind of trigger to be used in the histogram entries
JosueMtzGar Oct 2, 2024
ad8e520
correcting the histo of time for fdd
JosueMtzGar Oct 7, 2024
92f239d
adding the flags to 2022 data
JosueMtzGar Oct 8, 2024
6d014e7
Issue with access histos
JosueMtzGar Oct 31, 2024
c0e65a3
resolving issues
JosueMtzGar Oct 31, 2024
513fe23
resolving issues
JosueMtzGar Oct 31, 2024
8c59dcc
a solution for histo acces
JosueMtzGar Oct 31, 2024
42fd6f6
new changes,
JosueMtzGar Dec 10, 2024
b1df0b7
new changes, histo for rate calculus added
JosueMtzGar Dec 17, 2024
10d8803
new changes
JosueMtzGar Dec 20, 2024
ea1cbb5
new changes
JosueMtzGar Dec 24, 2024
67d02a5
new changes
JosueMtzGar Dec 24, 2024
3871ebd
new changes
JosueMtzGar Dec 24, 2024
8bc6ee2
adding leading BC pattern and corrections for rate
JosueMtzGar Apr 1, 2025
60345a9
new changes
JosueMtzGar Apr 11, 2025
0712249
doing a rebase from master branch
JosueMtzGar Nov 13, 2025
cbb7559
adding histo for rate correction by filled jobs in alihyperloop
JosueMtzGar Nov 18, 2025
4d3f704
fixing Megalinter issues
JosueMtzGar Nov 18, 2025
36d532b
fixing typo and logical filling of histos
JosueMtzGar Nov 21, 2025
9606ecf
new def od leading BC and addition of configurable for bc offset
JosueMtzGar Feb 1, 2026
173d8cc
remove of cout line
JosueMtzGar Feb 1, 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
98 changes: 30 additions & 68 deletions PWGMM/Lumi/Tasks/lumiStability.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -55,7 +55,8 @@ struct LumiStabilityTask {
Configurable<int> nOrbitsPerTF{"nOrbitsPerTF", 128, "number of orbits per time frame"};
Configurable<double> minOrbitConf{"minOrbitConf", 0, "minimum orbit"};
Configurable<bool> is2022Data{"is2022Data", true, "To 2022 data"};
Configurable<int> minEmpty{"minEmpty", 5, "number of BCs empty for leading BC"};
Configurable<int> minEmpty{"minEmpty", 13, "number of BCs empty for leading BC"};
Configurable<int> nBCsOffSet{"nBCsOffSet", 7, "number of BCs offset for FDD"};

Service<o2::ccdb::BasicCCDBManager> ccdb;
parameters::GRPLHCIFData* grplhcif = nullptr;
Expand All @@ -74,6 +75,7 @@ struct LumiStabilityTask {
std::bitset<o2::constants::lhc::LHCMaxBunches> bcPatternA;
std::bitset<o2::constants::lhc::LHCMaxBunches> bcPatternC;
std::bitset<o2::constants::lhc::LHCMaxBunches> bcPatternB;
std::bitset<o2::constants::lhc::LHCMaxBunches> bcPatternL;
std::bitset<o2::constants::lhc::LHCMaxBunches> bcPatternE;

void init(InitContext const&)
Expand Down Expand Up @@ -287,6 +289,7 @@ struct LumiStabilityTask {
bcPatternC = ~beamPatternA & beamPatternC;
bcPatternB = beamPatternA & beamPatternC;
bcPatternE = ~beamPatternA & ~beamPatternC;
bcPatternL = beamPatternA & beamPatternC;

for (int i = 0; i < nBCsPerOrbit; i++) {
if (bcPatternA[i]) {
Expand All @@ -297,18 +300,23 @@ struct LumiStabilityTask {
}
if (bcPatternB[i]) {
histos.fill(HIST("hBcB"), i);
}
if (bcPatternL[i]) {
bool isLeadBC = true;
for (int jbit = i - minEmpty; jbit < i; jbit++) {
int kbit = jbit;
if (kbit < 0)
kbit += nbin;
if (bcPatternB[kbit]) {
if (!bcPatternE[kbit]) {
bcPatternL[i] = false;
isLeadBC = false;
break;
}
}
if (isLeadBC)
if (isLeadBC) {
// bcPatternL[i] = true;
histos.fill(HIST("hBcBL"), i);
}
}
if (bcPatternE[i]) {
histos.fill(HIST("hBcE"), i);
Expand Down Expand Up @@ -379,8 +387,8 @@ struct LumiStabilityTask {
}

if (trgFDD) {
histos.fill(HIST("FDD/bcVertexTriggerCTP"), localBC + 7);
if (bcPatternB[localBC]) {
histos.fill(HIST("FDD/bcVertexTriggerCTP"), localBC + nBCsOffSet);
if (bcPatternB[localBC + nBCsOffSet]) {
histos.fill(HIST("FDD/nBCsVsTime"), timeSinceSOF);
histos.fill(HIST("FDD/hTimeForRateCTP"), (bc.timestamp() - tsSOR) * 1.e-3); // Converting ms into seconds
}
Expand All @@ -398,30 +406,21 @@ struct LumiStabilityTask {
histos.fill(HIST("FV0/hTimeForRateCTP"), (bc.timestamp() - tsSOR) * 1.e-3); // Converting ms into seconds
}
}
bool isLeadBC = true;
for (int jbit = localBC - minEmpty; jbit < localBC; jbit++) {
int kbit = jbit;
if (kbit < 0)
kbit += nbin;
if (bcPatternB[kbit]) {
isLeadBC = false;
break;
}
}
if (isLeadBC) {
if (bcPatternL[localBC + nBCsOffSet]) {
if (trgFDD) {
histos.fill(HIST("FDD/nBCsVsTimeLeadingBC"), timeSinceSOF);
histos.fill(HIST("FDD/hTimeForRateLeadingBCCTP"), (bc.timestamp() - tsSOR) * 1.e-3);
}
}
if (bcPatternL[localBC]) {
if (trgFT0) {
histos.fill(HIST("FT0/nBCsVsTimeLeadingBC"), timeSinceSOF);
histos.fill(HIST("FT0/hTimeForRateLeadingBCCTP"), (bc.timestamp() - tsSOR) * 1.e-3);
}
if (trgFV0) {
histos.fill(HIST("FV0/hTimeForRateLeadingBCCTP"), (bc.timestamp() - tsSOR) * 1.e-3);
}
}
// }
} // if over leading pattern
} // loop over bcs

for (auto const& fdd : fdds) {
Expand Down Expand Up @@ -465,18 +464,9 @@ struct LumiStabilityTask {

if (bcPatternB[localBC]) {
histos.fill(HIST("FDD/hTimeForRate"), (bc.timestamp() - tsSOR) * 1.e-3); // Converting ms into seconds
bool isLeadBC = true;
for (int jbit = localBC - minEmpty; jbit < localBC; jbit++) {
int kbit = jbit;
if (kbit < 0)
kbit += nbin;
if (bcPatternB[kbit]) {
isLeadBC = false;
break;
}
}
if (isLeadBC)
histos.fill(HIST("FDD/hTimeForRateLeadingBC"), (bc.timestamp() - tsSOR) * 1.e-3);
}
if (bcPatternL[localBC]) {
histos.fill(HIST("FDD/hTimeForRateLeadingBC"), (bc.timestamp() - tsSOR) * 1.e-3);
}

int deltaIndex = 0; // backward move counts
Expand Down Expand Up @@ -514,18 +504,6 @@ struct LumiStabilityTask {
if (bcPatternB[localBC]) {
histos.fill(HIST("FDD/timeACbcBVertex"), fdd.timeA(), fdd.timeC());
histos.fill(HIST("FDD/hBcBVertex"), localBC);
bool isLeadBC = true;
for (int jbit = localBC - minEmpty; jbit < localBC; jbit++) {
int kbit = jbit;
if (kbit < 0)
kbit += nbin;
if (bcPatternB[kbit]) {
isLeadBC = false;
break;
}
}
if (isLeadBC)
histos.fill(HIST("FDD/hBcBVertexL"), localBC);
histos.fill(HIST("FDD/hTimeAVertex"), fdd.timeA());
histos.fill(HIST("FDD/hTimeCVertex"), fdd.timeC());
if (is2022Data) {
Expand Down Expand Up @@ -555,6 +533,9 @@ struct LumiStabilityTask {
}
}
}
if (bcPatternL[localBC]) {
histos.fill(HIST("FDD/hBcBVertexL"), localBC);
}
if (bcPatternE[localBC]) {
histos.fill(HIST("FDD/timeACbcEVertex"), fdd.timeA(), fdd.timeC());
histos.fill(HIST("FDD/hBcEVertex"), localBC);
Expand All @@ -577,18 +558,6 @@ struct LumiStabilityTask {
if (bcPatternB[localBC]) {
histos.fill(HIST("FDD/timeACbcB"), fdd.timeA(), fdd.timeC());
histos.fill(HIST("FDD/hBcB"), localBC);
bool isLeadBC = true;
for (int jbit = localBC - minEmpty; jbit < localBC; jbit++) {
int kbit = jbit;
if (kbit < 0)
kbit += nbin;
if (bcPatternB[kbit]) {
isLeadBC = false;
break;
}
}
if (isLeadBC)
histos.fill(HIST("FDD/hBcBL"), localBC);
histos.fill(HIST("FDD/hTimeACoinc"), fdd.timeA());
histos.fill(HIST("FDD/hTimeCCoinc"), fdd.timeC());
if (!is2022Data) {
Expand Down Expand Up @@ -618,6 +587,9 @@ struct LumiStabilityTask {
}
}
}
if (bcPatternL[localBC]) {
histos.fill(HIST("FDD/hBcBL"), localBC);
}
if (bcPatternE[localBC]) {
histos.fill(HIST("FDD/timeACbcE"), fdd.timeA(), fdd.timeC());
histos.fill(HIST("FDD/hBcE"), localBC);
Expand Down Expand Up @@ -738,20 +710,6 @@ struct LumiStabilityTask {
histos.fill(HIST("FT0/timeACbcB"), ft0.timeA(), ft0.timeC());
histos.fill(HIST("FT0/hBcB"), localBC);
histos.fill(HIST("FT0/hTimeForRate"), (bc.timestamp() - tsSOR) * 1.e-3); // Converting ms into seconds
bool isLeadBC = true;
for (int jbit = localBC - minEmpty; jbit < localBC; jbit++) {
int kbit = jbit;
if (kbit < 0)
kbit += nbin;
if (bcPatternB[kbit]) {
isLeadBC = false;
break;
}
}
if (isLeadBC) {
histos.fill(HIST("FT0/hTimeForRateLeadingBC"), (bc.timestamp() - tsSOR) * 1.e-3); // Converting ms into seconds
histos.fill(HIST("FT0/hBcBL"), localBC);
}
histos.fill(HIST("FT0/hTimeA"), ft0.timeA());
histos.fill(HIST("FT0/hTimeC"), ft0.timeC());

Expand Down Expand Up @@ -780,6 +738,10 @@ struct LumiStabilityTask {
histos.fill(HIST("FT0/hValidTimevsBC"), localBC);
}
}
if (bcPatternL[localBC]) {
histos.fill(HIST("FT0/hTimeForRateLeadingBC"), (bc.timestamp() - tsSOR) * 1.e-3); // Converting ms into seconds
histos.fill(HIST("FT0/hBcBL"), localBC);
}
if (bcPatternE[localBC]) {
histos.fill(HIST("FT0/timeACbcE"), ft0.timeA(), ft0.timeC());
histos.fill(HIST("FT0/hBcE"), localBC);
Expand Down
Loading