Skip to content

[PWGLF] fixing bb implementation and adding cfg for tof mass cut #12183

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
maciacco merged 11 commits into from
Jul 24, 2025
Merged

[PWGLF] fixing bb implementation and adding cfg for tof mass cut #12183

Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
6416658
Update trHeAnalysis.cxx
mherzer28 Jul 20, 2025
785163b
Merge branch 'AliceO2Group:master' into master
mherzer28 Jul 20, 2025
e1204d1
add cfg for enabling/disabling tof mass cut
mherzer28 Jul 22, 2025
23a3081
Merge branch 'AliceO2Group:master' into master
mherzer28 Jul 22, 2025
e4dfcff
Please consider the following formatting changes
alibuild Jul 22, 2025
843322c
Merge pull request #13 from alibuild/alibot-cleanup-12183
mherzer28 Jul 22, 2025
abcd489
Update trHeAnalysis.cxx
mherzer28 Jul 22, 2025
4fa92ca
fix1
mherzer28 Jul 23, 2025
019535f
fix2
mherzer28 Jul 23, 2025
f00d366
Update trHeAnalysis.cxx
mherzer28 Jul 23, 2025
4d71e2e
Update trHeAnalysis.cxx
mherzer28 Jul 23, 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
44 changes: 30 additions & 14 deletions PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -194,6 +194,8 @@
Configurable<bool> cfgTPCPidMethod{"cfgTPCPidMethod", false, "Using own or built in bethe parametrization"}; // false for built in
Configurable<int> cfgMassMethod{"cfgMassMethod", 0, "0: Using built in 1: mass calculated with beta 2: mass calculated with the event time"};
Configurable<bool> cfgEnableItsClusterSizeCut{"cfgEnableItsClusterSizeCut", false, "Enable ITS cluster size cut"};
Configurable<bool> cfgEnableTofMassCut{"cfgEnableTofMassCut", false, "Enable TOF mass cut"};
Configurable<float> cfgTofMassCutPt{"cfgTofMassCutPt", 1.6f, "Pt value for which the TOF-cut starts to be used"};
// Set the multiplity event limits
Configurable<float> cfgLowMultCut{"cfgLowMultCut", 0.0f, "Accepted multiplicity percentage lower limit"};
Configurable<float> cfgHighMultCut{"cfgHighMultCut", 100.0f, "Accepted multiplicity percentage higher limit"};
Expand All @@ -219,8 +221,11 @@
Configurable<float> cfgCutMinItsClusterSizeHe{"cfgCutMinItsClusterSizeHe", 1.f, "Minimum ITS Cluster Size for He"};
Configurable<float> cfgCutMaxItsClusterSizeH3{"cfgCutMaxItsClusterSizeH3", 4.f, "Maximum ITS Cluster Size for Tr"};
Configurable<float> cfgCutMinItsClusterSizeH3{"cfgCutMinItsClusterSizeH3", 1.f, "Minimum ITS Cluster Size for Tr"};
Configurable<float> cfgCutMinTofMassH3{"cfgCutMinTofMassH3", 2.24f, "Minimum Tof mass H3"};
Configurable<float> cfgCutMaxTofMassH3{"cfgCutMaxTofMassH3", 3.32f, "Maximum TOF mass H3"};
Configurable<float> cfgCutMinTofMassH3{"cfgCutMinTofMassH3", 5.f, "Minimum Tof mass H3"};
Configurable<float> cfgCutMaxTofMassH3{"cfgCutMaxTofMassH3", 11.f, "Maximum TOF mass H3"};
Configurable<float> cfgMaxRigidity{"cfgMaxRigidity", 10.f, "Maximum rigidity value"};
Configurable<float> cfgMaxPt{"cfgMaxPt", 10.f, "Maximum pT value"};

// Set the kinematic and PID cuts for tracks
struct : ConfigurableGroup {
Configurable<float> pCut{"pCut", 0.6f, "Value of the p selection for spectra (default 0.3)"};
Expand Down Expand Up @@ -345,6 +350,9 @@
histos.fill(HIST("histogram/cuts"), 2);
continue;
}
if (track.pt() > cfgMaxPt || getRigidity(track) > cfgMaxRigidity) {
continue;
}
if (track.tpcNClsFound() < cfgCutTpcClusters) {
histos.fill(HIST("histogram/cuts"), 3);
continue;
Expand Down Expand Up @@ -399,9 +407,11 @@
continue;
}
}
if (getMass(track) < cfgCutMinTofMassH3 || getMass(track) > cfgCutMaxTofMassH3) {
histos.fill(HIST("histogram/cuts"), 13);
continue;
if (cfgEnableTofMassCut && track.pt() > cfgTofMassCutPt) {
if (getMass(track) < cfgCutMinTofMassH3 || getMass(track) > cfgCutMaxTofMassH3) {
histos.fill(HIST("histogram/cuts"), 13);
continue;
}
}
histos.fill(HIST("histogram/H3/H3-TPCsignVsTPCmomentum"),
getRigidity(track) * track.sign(),
Expand Down Expand Up @@ -510,6 +520,9 @@
histos.fill(HIST("histogram/cuts"), 2);
continue;
}
if (track.pt() > cfgMaxPt || getRigidity(track) > cfgMaxRigidity) {
continue;
}
if (track.tpcNClsFound() < cfgCutTpcClusters) {
histos.fill(HIST("histogram/cuts"), 3);
continue;
Expand Down Expand Up @@ -562,9 +575,11 @@
continue;
}
}
if (getMass(track) < cfgCutMinTofMassH3 || getMass(track) > cfgCutMaxTofMassH3) {
histos.fill(HIST("histogram/cuts"), 13);
continue;
if (cfgEnableTofMassCut && track.pt() > cfgTofMassCutPt) {
if (getMass(track) < cfgCutMinTofMassH3 || getMass(track) > cfgCutMaxTofMassH3) {
histos.fill(HIST("histogram/cuts"), 13);
continue;
}
}
histos.fill(HIST("histogram/H3/H3-TPCsignVsTPCmomentum"),
getRigidity(track) * (1.f * track.sign()),
Expand Down Expand Up @@ -709,16 +724,17 @@
float getMass(const T& track)
{
if (cfgMassMethod == 0) {
return track.mass();
float m = track.mass();
return m * m;
}
if (cfgMassMethod == 1) {
const float beta = track.beta();
const float rigidity = getRigidity(track);
float gamma = 1 / std::sqrt(1 - beta * beta);
float mass = (rigidity / std::sqrt(gamma * gamma - 1.f));
return mass;
float gamma = 1.f / std::sqrt(1.f - beta * beta);
float mass = rigidity / std::sqrt(gamma * gamma - 1.f);
return mass * mass;
}
if (cfgMassMethod == 2) {

Check failure on line 737 in PWGLF/TableProducer/Nuspex/trHeAnalysis.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.
const float rigidity = getRigidity(track);
float tofStartTime = track.evTimeForTrack();
float tofTime = track.tofSignal();
Expand All @@ -727,9 +743,9 @@
float time = tofTime - tofStartTime;
if (time > 0.f && length > 0.f) {
float beta = length / (CInCmPs * time);
float gamma = 1 / std::sqrt(1 - beta * beta);
float gamma = 1.f / std::sqrt(1.f - beta * beta);
float mass = rigidity / std::sqrt(gamma * gamma - 1.f);
return mass;
return mass * mass;
}
return -1.f;
}
Expand Down
Loading