Skip to content

[TOF] Add multiplicity in orbit and readout windows #2646

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

[TOF] Add multiplicity in orbit and readout windows #2646

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
b6d4d51
Add multiplicity in orbit
ercolessi Feb 27, 2026
43007a3
Increase bins for multiplicity orbit from 2000 to 6000
ercolessi Feb 27, 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
7 changes: 6 additions & 1 deletion Modules/TOF/include/TOF/TaskDigits.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -76,8 +76,11 @@ class TaskDigits final : public TaskInterface
static constexpr int mBinsOrbitPerTimeFrame = mRangeMaxOrbitPerTimeFrame * 3; /// Max range in the OrbitPerTimeFrame axis. 3 orbits are recorded per time frame
// Multiplicity
int mBinsMultiplicity = 2000; /// Number of bins in multiplicity plot
int mBinsMultiplicity2D = 2000; /// Number of bins in multiplicity plot
int mBinsMultiplicityOrbit = 6000; /// Number of bins in multiplicity plot
static constexpr int mRangeMinMultiplicity = 0; /// Min range in multiplicity plot
int mRangeMaxMultiplicity = mBinsMultiplicity; /// Max range in multiplicity plot
int mRangeMaxMultiplicityOrbit = mBinsMultiplicityOrbit; /// Max range in multiplicity plot
static constexpr int mBinsBCForMultiplicity = mRangeMaxBC; /// Number of bins for the BC axis in the multiplicity vs BC plot
// Time
int mBinsTime = 300; /// Number of bins in time plot
Expand Down Expand Up @@ -125,11 +128,13 @@ class TaskDigits final : public TaskInterface
std::shared_ptr<TH2S> mHistoNoisyChannels = nullptr; /// Channel flagged as noise (divided per flagged rate class)

// Multiplicity
std::shared_ptr<TH1I> mHistoMultiplicity = nullptr; /// TOF raw hit multiplicity per event
std::shared_ptr<TH1I> mHistoMultiplicity = nullptr; /// TOF raw hit multiplicity per event - all rw
std::shared_ptr<TH1I> mHistoMultiplicityIA = nullptr; /// TOF raw hit multiplicity per event - I/A side
std::shared_ptr<TH1I> mHistoMultiplicityOA = nullptr; /// TOF raw hit multiplicity per event - O/A side
std::shared_ptr<TH1I> mHistoMultiplicityIC = nullptr; /// TOF raw hit multiplicity per event - I/C side
std::shared_ptr<TH1I> mHistoMultiplicityOC = nullptr; /// TOF raw hit multiplicity per event - O/C side
std::shared_ptr<TH1I> mHistoMultiplicityOrbit = nullptr; /// TOF raw hit multiplicity per event - orbit
std::shared_ptr<TH1I> mHistoMultiplicityRW[3]; /// TOF raw hit multiplicity per event - RW1
std::shared_ptr<TH2F> mHitMultiplicityVsCrate = nullptr; /// TOF raw hit multiplicity per event vs Crate
std::shared_ptr<TProfile> mHitMultiplicityVsCratepro = nullptr; /// TOF raw hit multiplicity per event vs Crate (TProfile)
std::shared_ptr<TH2F> mHitMultiplicityVsBC = nullptr; /// TOF raw hit multiplicity per event vs BC
Expand Down
41 changes: 31 additions & 10 deletions Modules/TOF/src/TaskDigits.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -55,7 +55,10 @@ void TaskDigits::initialize(o2::framework::InitContext& /*ctx*/)
{
// Define parameters
utils::parseIntParameter(mCustomParameters, "NbinsMultiplicity", mBinsMultiplicity);
utils::parseIntParameter(mCustomParameters, "NbinsMultiplicity2D", mBinsMultiplicity2D);
utils::parseIntParameter(mCustomParameters, "RangeMaxMultiplicity", mRangeMaxMultiplicity);
utils::parseIntParameter(mCustomParameters, "NbinsMultiplicityOrbit", mBinsMultiplicityOrbit);
utils::parseIntParameter(mCustomParameters, "RangeMaxMultiplicityOrbit", mRangeMaxMultiplicityOrbit);

utils::parseIntParameter(mCustomParameters, "NbinsTime", mBinsTime);
utils::parseFloatParameter(mCustomParameters, "kNbinsWidthTime", fgkNbinsWidthTime);
Expand Down Expand Up @@ -144,17 +147,18 @@ void TaskDigits::initialize(o2::framework::InitContext& /*ctx*/)
}
getObjectsManager()->startPublishing(mHistoNoisyChannels.get());

// if mBinsMultiplicity > 1000 limit it in TH2F
int nBinsMultForTH2 = mBinsMultiplicity;
if (nBinsMultForTH2 > 1000) {
nBinsMultForTH2 = 1000;
ILOG(Info, Support) << "Requested Nbins in multiplicity is " << mBinsMultiplicity << " but limited to 1000 ONLY for TH2 " << ENDM;
}

// Multiplicity
mHistoMultiplicity = std::make_shared<TH1I>("Multiplicity/Integrated", "TOF hit multiplicity;TOF hits;Events ", mBinsMultiplicity, mRangeMinMultiplicity, mRangeMaxMultiplicity);
getObjectsManager()->startPublishing(mHistoMultiplicity.get());

for (int i = 0; i < 3; i++) {
mHistoMultiplicityRW[i] = std::make_shared<TH1I>(Form("Multiplicity/ReadoutWindow%i", i + 1), Form("TOF hit multiplicity in rw%i;TOF hits;Events ", i + 1), mBinsMultiplicity, mRangeMinMultiplicity, mRangeMaxMultiplicity);
getObjectsManager()->startPublishing(mHistoMultiplicityRW[i].get());
}

mHistoMultiplicityOrbit = std::make_shared<TH1I>("Multiplicity/IntegratedOrbit", "TOF hit multiplicity in orbit;TOF hits;Events ", mBinsMultiplicity, mRangeMinMultiplicity, mRangeMaxMultiplicityOrbit);
getObjectsManager()->startPublishing(mHistoMultiplicityOrbit.get());

mHistoMultiplicityIA = std::make_shared<TH1I>("Multiplicity/SectorIA", "TOF hit multiplicity - I/A side;TOF hits;Events ", mBinsMultiplicity, mRangeMinMultiplicity, mRangeMaxMultiplicity);
getObjectsManager()->startPublishing(mHistoMultiplicityIA.get());

Expand All @@ -167,13 +171,13 @@ void TaskDigits::initialize(o2::framework::InitContext& /*ctx*/)
mHistoMultiplicityOC = std::make_shared<TH1I>("Multiplicity/SectorOC", "TOF hit multiplicity - O/C side;TOF hits;Events ", mBinsMultiplicity, mRangeMinMultiplicity, mRangeMaxMultiplicity);
getObjectsManager()->startPublishing(mHistoMultiplicityOC.get());

mHitMultiplicityVsCrate = std::make_shared<TH2F>("Multiplicity/VsCrate", "TOF hit multiplicity vs Crate;Crate;TOF hits", RawDataDecoder::ncrates, 0, RawDataDecoder::ncrates, nBinsMultForTH2, mRangeMinMultiplicity, mRangeMaxMultiplicity);
mHitMultiplicityVsCrate = std::make_shared<TH2F>("Multiplicity/VsCrate", "TOF hit multiplicity vs Crate;Crate;TOF hits", RawDataDecoder::ncrates, 0, RawDataDecoder::ncrates, mBinsMultiplicity2D, mRangeMinMultiplicity, mRangeMaxMultiplicity);
getObjectsManager()->startPublishing(mHitMultiplicityVsCrate.get());

mHitMultiplicityVsCratepro = std::make_shared<TProfile>("Multiplicity/VsCratepro", "TOF hit multiplicity vs Crate;Crate;#LT TOF hits #GT", RawDataDecoder::ncrates, 0, RawDataDecoder::ncrates);
getObjectsManager()->startPublishing(mHitMultiplicityVsCratepro.get());

mHitMultiplicityVsBC = std::make_shared<TH2F>("Multiplicity/VsBC", "TOF hit multiplicity vs BC;BC;#TOF hits;Events", mBinsBCForMultiplicity, 0, mRangeMaxBC, nBinsMultForTH2, mRangeMinMultiplicity, mRangeMaxMultiplicity);
mHitMultiplicityVsBC = std::make_shared<TH2F>("Multiplicity/VsBC", "TOF hit multiplicity vs BC;BC;#TOF hits;Events", mBinsBCForMultiplicity, 0, mRangeMaxBC, mBinsMultiplicity2D, mRangeMinMultiplicity, mRangeMaxMultiplicity);
getObjectsManager()->startPublishing(mHitMultiplicityVsBC.get());

mHitMultiplicityVsBCpro = std::make_shared<TProfile>("Multiplicity/VsBCpro", "TOF hit multiplicity vs BC;BC;#TOF hits;Events", mBinsBCForMultiplicity, 0, mRangeMaxBC);
Expand Down Expand Up @@ -279,6 +283,7 @@ void TaskDigits::monitorData(o2::framework::ProcessingContext& ctx)

int currentReadoutWindow = 0;
int currentDiagnostics = 0;
int orbitMultiplicity = 0;
// Loop on readout windows
for (const auto& row : rows) {
const auto& digits_in_row = row.getBunchChannelData(digits); // Digits inside a readout window
Expand Down Expand Up @@ -345,7 +350,18 @@ void TaskDigits::monitorData(o2::framework::ProcessingContext& ctx)
currentDiagnostics += nDia;
}
}
currentReadoutWindow++;

int windowIndex = currentReadoutWindow % 3;
int ndigits_in_row = digits_in_row.size();

mHistoMultiplicityRW[windowIndex]->Fill(ndigits_in_row);
orbitMultiplicity += ndigits_in_row;

if (windowIndex == 2) { // End of readout windows, fill orbit multiplicity

mHistoMultiplicityOrbit->Fill(orbitMultiplicity);
orbitMultiplicity = 0; // Reset
}

// Loop on digits
for (auto const& digit : digits_in_row) {
Expand Down Expand Up @@ -467,6 +483,7 @@ void TaskDigits::monitorData(o2::framework::ProcessingContext& ctx)
ndigitsPerQuater[1] = 0;
ndigitsPerQuater[2] = 0;
ndigitsPerQuater[3] = 0;
currentReadoutWindow++;
}

for (int iorb = 0; iorb < nOrbits; iorb++) {
Expand Down Expand Up @@ -542,6 +559,10 @@ void TaskDigits::reset()
mHistoMultiplicityOA->Reset();
mHistoMultiplicityIC->Reset();
mHistoMultiplicityOC->Reset();
mHistoMultiplicityOrbit->Reset();
for (int i = 0; i < 3; i++) {
mHistoMultiplicityRW[i]->Reset();
}
mHitMultiplicityVsCrate->Reset();
mHitMultiplicityVsCratepro->Reset();
mHitMultiplicityVsBC->Reset();
Expand Down