Skip to content

[PWGCF] Please consider the following formatting changes to #13050 #1

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

Closed
alibuild wants to merge 1 commit into from
Closed

[PWGCF] Please consider the following formatting changes to #13050 #1

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
84 changes: 40 additions & 44 deletions PWGCF/Femto/Core/kinkBuilder.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -25,9 +25,9 @@
#include "PWGCF/Femto/DataModel/FemtoTables.h"

#include "Common/Core/RecoDecay.h"
#include "CommonConstants/PhysicsConstants.h"

#include "CommonConstants/MathConstants.h"
#include "CommonConstants/PhysicsConstants.h"
#include "Framework/AnalysisHelpers.h"
#include "Framework/Configurable.h"

Expand Down Expand Up @@ -64,7 +64,7 @@
o2::framework::Configurable<std::vector<float>> kinkTopoDcaMax{"kinkTopoDcaMax", {2.0f}, "Maximum kink topological DCA"}; \
o2::framework::Configurable<std::vector<float>> transRadMin{"transRadMin", {0.2f}, "Minimum transverse radius (cm)"}; \
o2::framework::Configurable<std::vector<float>> transRadMax{"transRadMax", {100.f}, "Maximum transverse radius (cm)"}; \
o2::framework::Configurable<std::vector<float>> dauAbsEtaMax{"DauAbsEtaMax", {0.8f}, "Maximum absolute pseudorapidity for daughter track"}; \

Check failure on line 67 in PWGCF/Femto/Core/kinkBuilder.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/configurable] 

Use lowerCamelCase for names of configurables and use the same name for the struct member as for the JSON string. (Declare the type and names on the same line.)
o2::framework::Configurable<std::vector<float>> dauDcaPvMin{"dauDcaPvMin", {0.0f}, "Minimum DCA of daughter from primary vertex (cm)"}; \
o2::framework::Configurable<std::vector<float>> mothDcaPvMax{"mothDcaPvMax", {1.0f}, "Maximum DCA of mother from primary vertex (cm)"}; \
o2::framework::Configurable<std::vector<float>> alphaAPMax{"alphaAPMax", {0.0f}, "Maximum Alpha_AP for Sigma candidates"}; \
Expand All @@ -82,16 +82,16 @@
#undef KINK_DEFAULT_BITS

// base selection for analysis task for kinks
#define KINK_DEFAULT_SELECTIONS(defaultMassMin, defaultMassMax, defaultPdgCode) \
o2::framework::Configurable<int> pdgCode{"pdgCode", defaultPdgCode, "Kink PDG code"}; \
o2::framework::Configurable<float> ptMin{"ptMin", 0.f, "Minimum pT"}; \
o2::framework::Configurable<float> ptMax{"ptMax", 999.f, "Maximum pT"}; \
o2::framework::Configurable<float> etaMin{"etaMin", -10.f, "Minimum eta"}; \
o2::framework::Configurable<float> etaMax{"etaMax", 10.f, "Maximum eta"}; \
o2::framework::Configurable<float> phiMin{"phiMin", 0.f, "Minimum phi"}; \
o2::framework::Configurable<float> phiMax{"phiMax", 1.f * o2::constants::math::TwoPI, "Maximum phi"}; \
o2::framework::Configurable<float> massMin{"massMin", defaultMassMin, "Minimum invariant mass for Sigma"}; \
o2::framework::Configurable<float> massMax{"massMax", defaultMassMax, "Maximum invariant mass for Sigma"}; \
#define KINK_DEFAULT_SELECTIONS(defaultMassMin, defaultMassMax, defaultPdgCode) \
o2::framework::Configurable<int> pdgCode{"pdgCode", defaultPdgCode, "Kink PDG code"}; \
o2::framework::Configurable<float> ptMin{"ptMin", 0.f, "Minimum pT"}; \
o2::framework::Configurable<float> ptMax{"ptMax", 999.f, "Maximum pT"}; \
o2::framework::Configurable<float> etaMin{"etaMin", -10.f, "Minimum eta"}; \
o2::framework::Configurable<float> etaMax{"etaMax", 10.f, "Maximum eta"}; \
o2::framework::Configurable<float> phiMin{"phiMin", 0.f, "Minimum phi"}; \
o2::framework::Configurable<float> phiMax{"phiMax", 1.f * o2::constants::math::TwoPI, "Maximum phi"}; \
o2::framework::Configurable<float> massMin{"massMin", defaultMassMin, "Minimum invariant mass for Sigma"}; \
o2::framework::Configurable<float> massMax{"massMax", defaultMassMax, "Maximum invariant mass for Sigma"}; \
o2::framework::Configurable<o2::aod::femtodatatypes::KinkMaskType> mask{"mask", 0, "Bitmask for kink selection"};

// base selection for analysis task for sigmas
Expand Down Expand Up @@ -139,8 +139,7 @@
{kAlphaAPMax, "alphaAPMax"},
{kQtAPMin, "qtAPMin"},
{kQtAPMax, "qtAPMax"},
{kCosPointingAngleMin, "cosPointingAngleMin"}
};
{kCosPointingAngleMin, "cosPointingAngleMin"}};

/// \class KinkCuts
/// \brief Cut class to contain and execute all cuts applied to kinks
Expand All @@ -155,12 +154,12 @@
void configure(T1& config, T2& filter)
{
mPtMin = filter.ptMin.value;
mPtMax = filter.ptMax.value;
mPtMax = filter.ptMax.value;
mEtaMin = filter.etaMin.value;
mEtaMax = filter.etaMax.value;
mPhiMin = filter.phiMin.value;
mPhiMax = filter.phiMax.value;

if constexpr (modes::isEqual(kinkType, modes::Kink::kSigma)) {
mMassSigmaLowerLimit = filter.massMinSigma.value;
mMassSigmaUpperLimit = filter.massMaxSigma.value;
Expand Down Expand Up @@ -206,23 +205,22 @@
std::array<float, 3> vMother = {kinkCand.xDecVtx(), kinkCand.yDecVtx(), kinkCand.zDecVtx()};
float pMother = std::sqrt(std::inner_product(momMother.begin(), momMother.end(), momMother.begin(), 0.f));
float vMotherNorm = std::sqrt(std::inner_product(vMother.begin(), vMother.end(), vMother.begin(), 0.f));
float cosPointingAngle = (vMotherNorm > 0.f && pMother > 0.f) ?
(std::inner_product(momMother.begin(), momMother.end(), vMother.begin(), 0.f)) / (pMother * vMotherNorm) : 0.f;
float cosPointingAngle = (vMotherNorm > 0.f && pMother > 0.f) ? (std::inner_product(momMother.begin(), momMother.end(), vMother.begin(), 0.f)) / (pMother * vMotherNorm) : 0.f;
this->evaluateObservable(kCosPointingAngleMin, cosPointingAngle);

this->evaluateObservable(kKinkTopoDcaMax, kinkCand.dcaKinkTopo());

// Compute transRadius
float transRadius = std::hypot(kinkCand.xDecVtx(), kinkCand.yDecVtx());
this->evaluateObservable(kTransRadMin, transRadius);
this->evaluateObservable(kTransRadMax, transRadius);

// Compute daughter eta
float px_daug = kinkCand.pxDaug();

Check failure on line 219 in PWGCF/Femto/Core/kinkBuilder.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable] 

Use lowerCamelCase for names of functions and variables.
float py_daug = kinkCand.pyDaug();

Check failure on line 220 in PWGCF/Femto/Core/kinkBuilder.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable] 

Use lowerCamelCase for names of functions and variables.
float pz_daug = kinkCand.pzDaug();

Check failure on line 221 in PWGCF/Femto/Core/kinkBuilder.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable] 

Use lowerCamelCase for names of functions and variables.
float p_daug = std::sqrt(px_daug*px_daug + py_daug*py_daug + pz_daug*pz_daug);
float p_daug = std::sqrt(px_daug * px_daug + py_daug * py_daug + pz_daug * pz_daug);

Check failure on line 222 in PWGCF/Femto/Core/kinkBuilder.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable] 

Use lowerCamelCase for names of functions and variables.
float eta_daug = (p_daug > 0.f) ? 0.5f * std::log((p_daug + pz_daug) / (p_daug - pz_daug)) : 0.f;

Check failure on line 223 in PWGCF/Femto/Core/kinkBuilder.h

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable] 

Use lowerCamelCase for names of functions and variables.
this->evaluateObservable(kDauAbsEtaMax, std::fabs(eta_daug));

this->evaluateObservable(kDauDcaPvMin, std::abs(kinkCand.dcaDaugPv()));
Expand All @@ -242,10 +240,10 @@
float py = kink.pyMoth();
float pz = kink.pzMoth();
float pt = std::hypot(px, py);
float p = std::sqrt(px*px + py*py + pz*pz);
float p = std::sqrt(px * px + py * py + pz * pz);
float eta = (p > 0.f) ? 0.5f * std::log((p + pz) / (p - pz)) : 0.f;
float phi = std::atan2(py, px);

return ((pt > mPtMin && pt < mPtMax) &&
(eta > mEtaMin && eta < mEtaMax) &&
(phi > mPhiMin && phi < mPhiMax));
Expand All @@ -265,17 +263,16 @@
float pxneut = pxmoth - pxch;
float pyneut = pymoth - pych;
float pzneut = pzmoth - pzch;

float sigmaMass = RecoDecay::m(
std::array{std::array{pxch, pych, pzch}, std::array{pxneut, pyneut, pzneut}},
std::array{o2::constants::physics::MassPionCharged, o2::constants::physics::MassNeutron}
);

std::array{std::array{pxch, pych, pzch}, std::array{pxneut, pyneut, pzneut}},
std::array{o2::constants::physics::MassPionCharged, o2::constants::physics::MassNeutron});

return (sigmaMass > mMassSigmaLowerLimit && sigmaMass < mMassSigmaUpperLimit);
}
return false;
}

protected:
float mMassSigmaLowerLimit = 1.15f;
float mMassSigmaUpperLimit = 1.25f;
Expand Down Expand Up @@ -361,46 +358,45 @@
float pxneut = pxmoth - pxch;
float pyneut = pymoth - pych;
float pzneut = pzmoth - pzch;

float mass = RecoDecay::m(
std::array{std::array{pxch, pych, pzch}, std::array{pxneut, pyneut, pzneut}},
std::array{o2::constants::physics::MassPionCharged, o2::constants::physics::MassNeutron}
);

std::array{std::array{pxch, pych, pzch}, std::array{pxneut, pyneut, pzneut}},
std::array{o2::constants::physics::MassPionCharged, o2::constants::physics::MassNeutron});

if (mProduceSigmas) {
// Compute mother eta and phi
float px = kink.pxMoth();
float py = kink.pyMoth();
float pz = kink.pzMoth();
float pt = std::hypot(px, py);
float p = std::sqrt(px*px + py*py + pz*pz);
float p = std::sqrt(px * px + py * py + pz * pz);
float eta = (p > 0.f) ? 0.5f * std::log((p + pz) / (p - pz)) : 0.f;
float phi = std::atan2(py, px);

kinkProducts.producedSigmas(collisionProducts.producedCollision.lastIndex(),
kink.mothSign() * pt,
eta,
phi,
mass,
daughterIndex);
kink.mothSign() * pt,
eta,
phi,
mass,
daughterIndex);
}
if (mProduceSigmaMasks) {
kinkProducts.producedSigmaMasks(mKinkSelection.getBitmask());
}
if (mProduceSigmaExtras) {
// Compute kink angle and transRadius
float pMoth = std::sqrt(pxmoth*pxmoth + pymoth*pymoth + pzmoth*pzmoth);
float pDaug = std::sqrt(pxch*pxch + pych*pych + pzch*pzch);
float pMoth = std::sqrt(pxmoth * pxmoth + pymoth * pymoth + pzmoth * pzmoth);
float pDaug = std::sqrt(pxch * pxch + pych * pych + pzch * pzch);
float kinkAngle = 0.f;
if (pMoth > 0.f && pDaug > 0.f) {
float dotProduct = pxmoth*pxch + pymoth*pych + pzmoth*pzch;
float dotProduct = pxmoth * pxch + pymoth * pych + pzmoth * pzch;
float cosAngle = dotProduct / (pMoth * pDaug);
cosAngle = std::max(-1.0f, std::min(1.0f, cosAngle)); // Clamp
kinkAngle = std::acos(cosAngle);
}

float transRadius = std::hypot(kink.xDecVtx(), kink.yDecVtx());

kinkProducts.producedSigmaExtras(
kinkAngle,
kink.dcaDaugPv(),
Expand Down
50 changes: 25 additions & 25 deletions PWGCF/Femto/Core/kinkHistManager.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -62,7 +62,7 @@ enum KinkHist {
kKinkHistLast
};

#define KINK_DEFAULT_BINNING(defaultMassMin, defaultMassMax) \
#define KINK_DEFAULT_BINNING(defaultMassMin, defaultMassMax) \
o2::framework::ConfigurableAxis pt{"pt", {{600, 0, 6}}, "Pt"}; \
o2::framework::ConfigurableAxis eta{"eta", {{300, -1.5, 1.5}}, "Eta"}; \
o2::framework::ConfigurableAxis phi{"phi", {{720, 0, 1.f * o2::constants::math::TwoPI}}, "Phi"}; \
Expand Down Expand Up @@ -96,23 +96,23 @@ using ConfSigmaQaBinning1 = ConfKinkQaBinning<PrefixSigmaQaBinning1>;
// the enum gives the correct index in the array
constexpr std::array<histmanager::HistInfo<KinkHist>, kKinkHistLast> HistTable = {
{{kPt, o2::framework::kTH1F, "hPt", "Transverse Momentum; p_{T} (GeV/#it{c}); Entries"},
{kEta, o2::framework::kTH1F, "hEta", "Pseudorapidity; #eta; Entries"},
{kPhi, o2::framework::kTH1F, "hPhi", "Azimuthal angle; #varphi; Entries"},
{kMass, o2::framework::kTH1F, "hMass", "Invariant Mass; m_{Inv} (GeV/#it{c}^{2}); Entries"},
{kSign, o2::framework::kTH1F, "hSign", "Sign; sign; Entries"},
{kKinkAngle, o2::framework::kTH1F, "hKinkAngle", "Kink Angle; Angle (rad); Entries"},
{kDcaMothToPV, o2::framework::kTH1F, "hDcaMothToPV", "Mother DCA to PV; DCA (cm); Entries"},
{kDcaDaugToPV, o2::framework::kTH1F, "hDcaDaugToPV", "Daughter DCA to PV; DCA (cm); Entries"},
{kDecayVtxX, o2::framework::kTH1F, "hDecayVtxX", "Decay Vertex X; x (cm); Entries"},
{kDecayVtxY, o2::framework::kTH1F, "hDecayVtxY", "Decay Vertex Y; y (cm); Entries"},
{kDecayVtxZ, o2::framework::kTH1F, "hDecayVtxZ", "Decay Vertex Z; z (cm); Entries"},
{kDecayVtx, o2::framework::kTH1F, "hDecayVtx", "Decay Distance from PV; r (cm); Entries"},
{kTransRadius, o2::framework::kTH1F, "hTransRadius", "Transverse Decay Radius; r_{xy} (cm); Entries"},
{kPtVsEta, o2::framework::kTH2F, "hPtVsEta", "p_{T} vs #eta; p_{T} (GeV/#it{c}); #eta"},
{kPtVsPhi, o2::framework::kTH2F, "hPtVsPhi", "p_{T} vs #varphi; p_{T} (GeV/#it{c}); #varphi"},
{kPhiVsEta, o2::framework::kTH2F, "hPhiVsEta", "#varphi vs #eta; #varphi; #eta"},
{kPtVsKinkAngle, o2::framework::kTH2F, "hPtVsKinkAngle", "p_{T} vs kink angle; p_{T} (GeV/#it{c}); kink angle (rad)"},
{kPtVsDecayRadius, o2::framework::kTH2F, "hPtVsDecayRadius", "p_{T} vs transverse decay radius; p_{T} (GeV/#it{c}); r_{xy} (cm)"}}};
{kEta, o2::framework::kTH1F, "hEta", "Pseudorapidity; #eta; Entries"},
{kPhi, o2::framework::kTH1F, "hPhi", "Azimuthal angle; #varphi; Entries"},
{kMass, o2::framework::kTH1F, "hMass", "Invariant Mass; m_{Inv} (GeV/#it{c}^{2}); Entries"},
{kSign, o2::framework::kTH1F, "hSign", "Sign; sign; Entries"},
{kKinkAngle, o2::framework::kTH1F, "hKinkAngle", "Kink Angle; Angle (rad); Entries"},
{kDcaMothToPV, o2::framework::kTH1F, "hDcaMothToPV", "Mother DCA to PV; DCA (cm); Entries"},
{kDcaDaugToPV, o2::framework::kTH1F, "hDcaDaugToPV", "Daughter DCA to PV; DCA (cm); Entries"},
{kDecayVtxX, o2::framework::kTH1F, "hDecayVtxX", "Decay Vertex X; x (cm); Entries"},
{kDecayVtxY, o2::framework::kTH1F, "hDecayVtxY", "Decay Vertex Y; y (cm); Entries"},
{kDecayVtxZ, o2::framework::kTH1F, "hDecayVtxZ", "Decay Vertex Z; z (cm); Entries"},
{kDecayVtx, o2::framework::kTH1F, "hDecayVtx", "Decay Distance from PV; r (cm); Entries"},
{kTransRadius, o2::framework::kTH1F, "hTransRadius", "Transverse Decay Radius; r_{xy} (cm); Entries"},
{kPtVsEta, o2::framework::kTH2F, "hPtVsEta", "p_{T} vs #eta; p_{T} (GeV/#it{c}); #eta"},
{kPtVsPhi, o2::framework::kTH2F, "hPtVsPhi", "p_{T} vs #varphi; p_{T} (GeV/#it{c}); #varphi"},
{kPhiVsEta, o2::framework::kTH2F, "hPhiVsEta", "#varphi vs #eta; #varphi; #eta"},
{kPtVsKinkAngle, o2::framework::kTH2F, "hPtVsKinkAngle", "p_{T} vs kink angle; p_{T} (GeV/#it{c}); kink angle (rad)"},
{kPtVsDecayRadius, o2::framework::kTH2F, "hPtVsDecayRadius", "p_{T} vs transverse decay radius; p_{T} (GeV/#it{c}); r_{xy} (cm)"}}};

template <typename T>
auto makeKinkHistSpecMap(const T& confBinningAnalysis)
Expand Down Expand Up @@ -169,7 +169,7 @@ class KinkHistManager
public:
/// Destructor
virtual ~KinkHistManager() = default;

/// Initializes histograms for the task
/// \param registry Histogram registry to be passed
///
Expand Down Expand Up @@ -240,7 +240,7 @@ class KinkHistManager
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kEta, HistTable)), kinkcandidate.eta());
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kPhi, HistTable)), kinkcandidate.phi());
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kMass, HistTable)), kinkcandidate.mass());

if constexpr (isEqual(kink, modes::Kink::kSigma)) {
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kSign, HistTable)), kinkcandidate.sign());
}
Expand All @@ -252,11 +252,11 @@ class KinkHistManager
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kDecayVtxX, HistTable)), kinkcandidate.decayVtxX());
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kDecayVtxY, HistTable)), kinkcandidate.decayVtxY());
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kDecayVtxZ, HistTable)), kinkcandidate.decayVtxZ());

// Calculate decay distance from PV
float decayDistance = std::sqrt(kinkcandidate.decayVtxX() * kinkcandidate.decayVtxX() +
kinkcandidate.decayVtxY() * kinkcandidate.decayVtxY() +
kinkcandidate.decayVtxZ() * kinkcandidate.decayVtxZ());
float decayDistance = std::sqrt(kinkcandidate.decayVtxX() * kinkcandidate.decayVtxX() +
kinkcandidate.decayVtxY() * kinkcandidate.decayVtxY() +
kinkcandidate.decayVtxZ() * kinkcandidate.decayVtxZ());
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kDecayVtx, HistTable)), decayDistance);
mHistogramRegistry->fill(HIST(kinkPrefix) + HIST(QaDir) + HIST(GetHistName(kTransRadius, HistTable)), kinkcandidate.transRadius());

Expand All @@ -270,7 +270,7 @@ class KinkHistManager
}

/// Fill histograms for kink candidates - overload with track table argument
/// \param kinkcandidate Kink candidate to fill histograms for
/// \param kinkcandidate Kink candidate to fill histograms for
/// \param tracks Track table for daughter access
template <typename T1, typename T2>
void fill(T1 const& kinkcandidate, T2 const& /*tracks*/)
Expand Down
44 changes: 22 additions & 22 deletions PWGCF/Femto/DataModel/FemtoTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -508,34 +508,34 @@ DECLARE_SOA_INDEX_COLUMN_FULL(ChaDau, chaDau, int32_t, FTracks, "_ChaDau"); //!

// table for basic sigma minus information
DECLARE_SOA_TABLE_STAGED_VERSIONED(FSigmas_001, "FSIGMA", 1,
o2::soa::Index<>,
femtobase::stored::CollisionId, // use sign to differentiate between sigma minus (-1) and anti sigma minus (+1)
femtobase::stored::SignedPt,
femtobase::stored::Eta,
femtobase::stored::Phi,
femtobase::stored::Mass,
femtokinks::ChaDauId,
femtobase::dynamic::Sign<femtobase::stored::SignedPt>,
femtobase::dynamic::Pt<femtobase::stored::SignedPt>,
femtobase::dynamic::P<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Px<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Py<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Pz<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Theta<femtobase::stored::Eta>);
o2::soa::Index<>,
femtobase::stored::CollisionId, // use sign to differentiate between sigma minus (-1) and anti sigma minus (+1)
femtobase::stored::SignedPt,
femtobase::stored::Eta,
femtobase::stored::Phi,
femtobase::stored::Mass,
femtokinks::ChaDauId,
femtobase::dynamic::Sign<femtobase::stored::SignedPt>,
femtobase::dynamic::Pt<femtobase::stored::SignedPt>,
femtobase::dynamic::P<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Px<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Py<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Pz<femtobase::stored::SignedPt, femtobase::stored::Eta>,
femtobase::dynamic::Theta<femtobase::stored::Eta>);
using FSigmas = FSigmas_001;

DECLARE_SOA_TABLE_STAGED_VERSIONED(FSigmaMasks_001, "FSIGMAMASKS", 1,
femtokinks::Mask);
femtokinks::Mask);
using FSigmaMasks = FSigmaMasks_001;

DECLARE_SOA_TABLE_STAGED_VERSIONED(FSigmaExtras_001, "FSIGMAEXTRAS", 1,
femtokinks::KinkAngle,
femtokinks::DcaDaugToPV,
femtokinks::DcaMothToPV,
femtokinks::DecayVtxX,
femtokinks::DecayVtxY,
femtokinks::DecayVtxZ,
femtokinks::TransRadius);
femtokinks::KinkAngle,
femtokinks::DcaDaugToPV,
femtokinks::DcaMothToPV,
femtokinks::DecayVtxX,
femtokinks::DecayVtxY,
femtokinks::DecayVtxZ,
femtokinks::TransRadius);

using FSigmaExtras = FSigmaExtras_001;

Expand Down
Loading
Loading