Skip to content

[PWGLF] adding secondary vertex for the heptaquark from daughter lambda #14229

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 17, 2025
Merged

[PWGLF] adding secondary vertex for the heptaquark from daughter lambda #14229

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
dcdf3d1
adding vertex
jikim1290 Dec 16, 2025
3dc62cd
Please consider the following formatting changes
alibuild Dec 16, 2025
67a19d3
Merge pull request #172 from alibuild/alibot-cleanup-14229
jikim1290 Dec 16, 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
6 changes: 6 additions & 0 deletions PWGLF/DataModel/ReducedHeptaQuarkTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -48,16 +48,19 @@
namespace hqtrack
{
DECLARE_SOA_INDEX_COLUMN(RedHQEvent, redHQEvent);
DECLARE_SOA_COLUMN(HQId, hqId, int); //! HQ ID

Check failure on line 51 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQPx, hqPx, float); //! HQ Px

Check failure on line 52 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQPy, hqPy, float); //! HQ Py

Check failure on line 53 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQPz, hqPz, float); //! HQ Pz

Check failure on line 54 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQd1Px, hqd1Px, float); //! HQ d1 Px

Check failure on line 55 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQd1Py, hqd1Py, float); //! HQ d1 Py

Check failure on line 56 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQd1Pz, hqd1Pz, float); //! HQ d1 Pz

Check failure on line 57 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQd2Px, hqd2Px, float); //! HQ d2 Px

Check failure on line 58 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQd2Py, hqd2Py, float); //! HQ d2 Py

Check failure on line 59 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQd2Pz, hqd2Pz, float); //! HQ d2 Pz

Check failure on line 60 in PWGLF/DataModel/ReducedHeptaQuarkTables.h

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(HQx, hqx, float); //! HQ x
DECLARE_SOA_COLUMN(HQy, hqy, float); //! HQ y
DECLARE_SOA_COLUMN(HQz, hqz, float); //! HQ z
DECLARE_SOA_COLUMN(HQMass, hqMass, float); //! HQ Mass
DECLARE_SOA_COLUMN(HQd1Index, hqd1Index, int64_t); //! HQ d1 index
DECLARE_SOA_COLUMN(HQd2Index, hqd2Index, int64_t); //! HQ d2 index
Expand All @@ -84,6 +87,9 @@
hqtrack::HQd2Px,
hqtrack::HQd2Py,
hqtrack::HQd2Pz,
hqtrack::HQx,
hqtrack::HQy,
hqtrack::HQz,
hqtrack::HQMass,
hqtrack::HQd1Index,
hqtrack::HQd2Index,
Expand Down
12 changes: 11 additions & 1 deletion PWGLF/TableProducer/Resonances/HeptaQuarktable.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -38,6 +38,7 @@
#include <Framework/Configurable.h>

#include <Math/GenVector/Boost.h>
#include <Math/Vector3D.h>
#include <Math/Vector4D.h>
#include <TMath.h>

Expand Down Expand Up @@ -223,6 +224,7 @@ struct heptaquarktable {
}

ROOT::Math::PxPyPzMVector DauVec1, DauVec2, HQMesonMother, HQVectorDummy, HQd1dummy, HQd2dummy;
ROOT::Math::XYZVector HQPosVectorDummy;

void processHQReducedTable(EventCandidates::iterator const& collision, TrackCandidates const& /*tracks*/, aod::V0Datas const& V0s, aod::BCsWithTimestamps const&)
{
Expand Down Expand Up @@ -252,6 +254,7 @@ struct heptaquarktable {
std::vector<float> HQd2TOF = {};

std::vector<ROOT::Math::PtEtaPhiMVector> hqresonance, hqresonanced1, hqresonanced2;
std::vector<ROOT::Math::XYZVector> hqresonancePosition;

histos.fill(HIST("hEventstat"), 0.5);
if (!(collision.sel8() && collision.selection_bit(aod::evsel::kNoTimeFrameBorder) && collision.selection_bit(aod::evsel::kNoITSROFrameBorder) && collision.selection_bit(aod::evsel::kNoSameBunchPileup) && collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV)))
Expand Down Expand Up @@ -309,6 +312,9 @@ struct heptaquarktable {
hqresonanced2.push_back(temp2);
hqresonance.push_back(temp3);

ROOT::Math::XYZVector temppos(0, 0, 0);
hqresonancePosition.push_back(temppos);

HQId.push_back(333);

HQd1Index.push_back(track1.globalIndex());
Expand Down Expand Up @@ -388,6 +394,9 @@ struct heptaquarktable {
hqresonanced2.push_back(temp2);
hqresonance.push_back(temp3);

ROOT::Math::XYZVector temppos(v0.x(), v0.y(), v0.z());
hqresonancePosition.push_back(temppos);

HQd1Index.push_back(postrack_v0.globalIndex());
HQd2Index.push_back(negtrack_v0.globalIndex());

Expand Down Expand Up @@ -434,12 +443,13 @@ struct heptaquarktable {
//// Fill track table for HQ//////////////////
for (auto if1 = hqresonance.begin(); if1 != hqresonance.end(); ++if1) {
auto i5 = std::distance(hqresonance.begin(), if1);
HQPosVectorDummy = hqresonancePosition.at(i5);
HQVectorDummy = hqresonance.at(i5);
HQd1dummy = hqresonanced1.at(i5);
HQd2dummy = hqresonanced2.at(i5);
hqTrack(indexEvent, HQId.at(i5), HQVectorDummy.Px(), HQVectorDummy.Py(), HQVectorDummy.Pz(),
HQd1dummy.Px(), HQd1dummy.Py(), HQd1dummy.Pz(), HQd2dummy.Px(), HQd2dummy.Py(), HQd2dummy.Pz(),
HQVectorDummy.M(),
HQPosVectorDummy.X(), HQPosVectorDummy.Y(), HQPosVectorDummy.Z(), HQVectorDummy.M(),
HQd1Index.at(i5), HQd2Index.at(i5),
HQd1Charge.at(i5), HQd2Charge.at(i5), HQd1TPC.at(i5), HQd2TPC.at(i5),
HQd1TOFHit.at(i5), HQd2TOFHit.at(i5), HQd1TOF.at(i5), HQd2TOF.at(i5));
Expand Down
Loading