Skip to content

Commit

Permalink
[PWGLF] Change V0 Cos PA to regular column for filterability (AliceO2…
Browse files Browse the repository at this point in the history 
…Group#3834)

* Change V0 Cos PA to regular for filter-ability

* Please consider the following formatting changes (#183)

---------

Co-authored-by: ALICE Builder <[email protected]>
  • Loading branch information
@ddobrigk @alibuild
ddobrigk and alibuild authored Nov 20, 2023
1 parent 9e00cf5 commit 0f0165c
Show file tree
Hide file tree
Showing 47 changed files with 148 additions and 127 deletions.
2 changes: 1 addition & 1 deletion DPG/Tasks/AOTTrack/PID/qaPIDWithV0s.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -124,7 +124,7 @@ struct pidQaWithV0s {
template <typename T1, typename C>
bool acceptV0s(const T1& v0, const C& collision) // Apply general selections on V0
{
float cospa = v0.v0cosPA(collision.posX(), collision.posY(), collision.posZ());
float cospa = v0.v0cosPA();
fillHistogram("beforeselections/cospa", cospa);
fillHistogram("beforeselections/decayradius", v0.v0radius());
if (cospa < v0cospa) {
Expand Down
2 changes: 1 addition & 1 deletion DPG/Tasks/AOTTrack/perfK0sResolution.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -70,7 +70,7 @@ struct perfK0sResolution {
// Apply selections on V0
if (TMath::Abs(v0.yK0Short()) > rapidity)
return kFALSE;
if (v0.v0cosPA(collision.posX(), collision.posY(), collision.posZ()) < v0setting_cospa)
if (v0.v0cosPA() < v0setting_cospa)
return kFALSE;
if (v0.v0radius() < v0setting_radius)
return kFALSE;
Expand Down
2 changes: 1 addition & 1 deletion DPG/Tasks/AOTTrack/qaK0sTrackingEfficiency.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -89,7 +89,7 @@ struct qaK0sTrackingEfficiency {
bool acceptV0(const T1& v0, const T2& ntrack, const T2& ptrack, const C& collision)
{
// Apply selections on V0
if (v0.v0cosPA(collision.posX(), collision.posY(), collision.posZ()) < v0cospa)
if (v0.v0cosPA() < v0cospa)
return kFALSE;
if (TMath::Abs(v0.yK0Short()) > rapidity)
return kFALSE;
Expand Down
2 changes: 1 addition & 1 deletion DPG/Tasks/TPC/tpcSkimsTableCreator.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -82,7 +82,7 @@ struct TreeWriterTpcV0 {

const float alpha = v0.alpha();
const float qt = v0.qtarm();
const float cosPA = v0.v0cosPA(collision.posX(), collision.posY(), collision.posZ());
const float cosPA = v0.v0cosPA();
const float pT = v0.pt();
const float v0radius = v0.v0radius();
const float gammapsipair = v0.psipair();
Expand Down
8 changes: 4 additions & 4 deletions EventFiltering/PWGCF/CFFilterAll.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1080,7 +1080,7 @@ struct CFFilter {
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();
const float invMassAntiLambda = v0.mAntiLambda();
Expand Down Expand Up @@ -1506,7 +1506,7 @@ struct CFFilter {
registry.fill(HIST("TrackCuts/V0Before/fInvMassLambdavsAntiLambda"), v0.mLambda(), v0.mAntiLambda());
registry.fill(HIST("TrackCuts/V0Before/fInvMassV0BeforeKaonvsV0Before"), v0.mK0Short(), v0.mLambda());
registry.fill(HIST("TrackCuts/V0Before/fV0DCADaugh"), v0.dcaV0daughters());
registry.fill(HIST("TrackCuts/V0Before/fV0CPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
registry.fill(HIST("TrackCuts/V0Before/fV0CPA"), v0.v0cosPA());
registry.fill(HIST("TrackCuts/V0Before/fV0TranRad"), v0.v0radius());
registry.fill(HIST("TrackCuts/V0Before/f0DecVtxX"), v0.x());
registry.fill(HIST("TrackCuts/V0Before/f0DecVtxY"), v0.y());
Expand Down Expand Up @@ -1540,7 +1540,7 @@ struct CFFilter {
registry.fill(HIST("TrackCuts/Lambda/fInvMassLambda"), v0.mLambda());
registry.fill(HIST("TrackCuts/Lambda/fInvMassLambdaKaonvsLambda"), v0.mK0Short(), v0.mLambda());
registry.fill(HIST("TrackCuts/Lambda/fV0DCADaugh"), v0.dcaV0daughters());
registry.fill(HIST("TrackCuts/Lambda/fV0CPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
registry.fill(HIST("TrackCuts/Lambda/fV0CPA"), v0.v0cosPA());
registry.fill(HIST("TrackCuts/Lambda/fV0TranRad"), v0.v0radius());
registry.fill(HIST("TrackCuts/Lambda/f0DecVtxX"), v0.x());
registry.fill(HIST("TrackCuts/Lambda/f0DecVtxY"), v0.y());
Expand All @@ -1566,7 +1566,7 @@ struct CFFilter {
registry.fill(HIST("TrackCuts/AntiLambda/fInvMassAntiLambda"), v0.mAntiLambda());
registry.fill(HIST("TrackCuts/AntiLambda/fInvMassAntiLambdaKaonvsAntiLambda"), v0.mK0Short(), v0.mAntiLambda());
registry.fill(HIST("TrackCuts/AntiLambda/fV0DCADaugh"), v0.dcaV0daughters());
registry.fill(HIST("TrackCuts/AntiLambda/fV0CPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
registry.fill(HIST("TrackCuts/AntiLambda/fV0CPA"), v0.v0cosPA());
registry.fill(HIST("TrackCuts/AntiLambda/fV0TranRad"), v0.v0radius());
registry.fill(HIST("TrackCuts/AntiLambda/f0DecVtxX"), v0.x());
registry.fill(HIST("TrackCuts/AntiLambda/f0DecVtxY"), v0.y());
Expand Down
8 changes: 4 additions & 4 deletions EventFiltering/PWGCF/CFFilterQA.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -934,7 +934,7 @@ struct CFFilterQA {
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();
const float invMassAntiLambda = v0.mAntiLambda();
Expand Down Expand Up @@ -1293,7 +1293,7 @@ struct CFFilterQA {
registry.fill(HIST("TrackCuts/V0Before/fInvMassLambdavsAntiLambda"), v0.mLambda(), v0.mAntiLambda());
registry.fill(HIST("TrackCuts/V0Before/fInvMassV0BeforeKaonvsV0Before"), v0.mK0Short(), v0.mLambda());
registry.fill(HIST("TrackCuts/V0Before/fV0DCADaugh"), v0.dcaV0daughters());
registry.fill(HIST("TrackCuts/V0Before/fV0CPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
registry.fill(HIST("TrackCuts/V0Before/fV0CPA"), v0.v0cosPA());
registry.fill(HIST("TrackCuts/V0Before/fV0TranRad"), v0.v0radius());
registry.fill(HIST("TrackCuts/V0Before/f0DecVtxX"), v0.x());
registry.fill(HIST("TrackCuts/V0Before/f0DecVtxY"), v0.y());
Expand Down Expand Up @@ -1327,7 +1327,7 @@ struct CFFilterQA {
registry.fill(HIST("TrackCuts/Lambda/fInvMassLambda"), v0.mLambda());
registry.fill(HIST("TrackCuts/Lambda/fInvMassLambdaKaonvsLambda"), v0.mK0Short(), v0.mLambda());
registry.fill(HIST("TrackCuts/Lambda/fV0DCADaugh"), v0.dcaV0daughters());
registry.fill(HIST("TrackCuts/Lambda/fV0CPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
registry.fill(HIST("TrackCuts/Lambda/fV0CPA"), v0.v0cosPA());
registry.fill(HIST("TrackCuts/Lambda/fV0TranRad"), v0.v0radius());
registry.fill(HIST("TrackCuts/Lambda/f0DecVtxX"), v0.x());
registry.fill(HIST("TrackCuts/Lambda/f0DecVtxY"), v0.y());
Expand All @@ -1354,7 +1354,7 @@ struct CFFilterQA {
registry.fill(HIST("TrackCuts/AntiLambda/fInvMassAntiLambda"), v0.mAntiLambda());
registry.fill(HIST("TrackCuts/AntiLambda/fInvMassAntiLambdaKaonvsAntiLambda"), v0.mK0Short(), v0.mAntiLambda());
registry.fill(HIST("TrackCuts/AntiLambda/fV0DCADaugh"), v0.dcaV0daughters());
registry.fill(HIST("TrackCuts/AntiLambda/fV0CPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
registry.fill(HIST("TrackCuts/AntiLambda/fV0CPA"), v0.v0cosPA());
registry.fill(HIST("TrackCuts/AntiLambda/fV0TranRad"), v0.v0radius());
registry.fill(HIST("TrackCuts/AntiLambda/f0DecVtxX"), v0.x());
registry.fill(HIST("TrackCuts/AntiLambda/f0DecVtxY"), v0.y());
Expand Down
4 changes: 2 additions & 2 deletions EventFiltering/PWGHF/HFFilterHelpers.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -845,7 +845,7 @@ inline int8_t HfFilterHelper::isSelectedV0(const V0& v0, const std::array<T, 2>&
}
}

auto v0CosinePa = v0.v0cosPA(collision.posX(), collision.posY(), collision.posZ());
auto v0CosinePa = v0.v0cosPA();
// cosine of pointing angle
if (TESTBIT(isSelected, kPhoton) && v0CosinePa < mMinGammaCosinePa) {
CLRBIT(isSelected, kPhoton);
Expand Down Expand Up @@ -890,7 +890,7 @@ inline int8_t HfFilterHelper::isSelectedV0(const V0& v0, const std::array<T, 2>&

// DCA V0 and V0 daughters
for (int iV0{kK0S}; iV0 < kNV0; ++iV0) {
if (TESTBIT(isSelected, iV0) && v0.dcav0topv(collision.posX(), collision.posY(), collision.posZ()) > 0.1f) { // we want only primary V0s
if (TESTBIT(isSelected, iV0) && v0.dcav0topv() > 0.1f) { // we want only primary V0s
CLRBIT(isSelected, iV0);
if (activateQA > 1) {
hV0Selected->Fill(5., iV0);
Expand Down
4 changes: 2 additions & 2 deletions EventFiltering/PWGLF/filterf1proton.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -306,15 +306,15 @@ struct filterf1proton {
template <typename Collision, typename V0>
bool SelectionV0(Collision const& collision, V0 const& candidate)
{
if (fabs(candidate.dcav0topv(collision.posX(), collision.posY(), collision.posZ())) > cMaxV0DCA) {
if (fabs(candidate.dcav0topv()) > cMaxV0DCA) {
return false;
}

const float pT = candidate.pt();
const std::vector<float> decVtx = {candidate.x(), candidate.y(), candidate.z()};
const float tranRad = candidate.v0radius();
const float dcaDaughv0 = candidate.dcaV0daughters();
const float cpav0 = candidate.v0cosPA(collision.posX(), collision.posY(), collision.posZ());
const float cpav0 = candidate.v0cosPA();

float CtauK0s = candidate.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * o2::analysis::pdg::MassK0Short;
float lowmasscutks0 = 0.497 - 2.0 * cSigmaMassKs0;
Expand Down
10 changes: 5 additions & 5 deletions PWGCF/FemtoDream/FemtoDreamV0Selection.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -415,7 +415,7 @@ bool FemtoDreamV0Selection::isSelectedMinimal(C const& col, V const& v0,
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();
const float invMassAntiLambda = v0.mAntiLambda();
Expand Down Expand Up @@ -498,7 +498,7 @@ void FemtoDreamV0Selection::fillLambdaQA(C const& col, V const& v0,
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();

Expand Down Expand Up @@ -587,7 +587,7 @@ std::array<cutContainerType, 5>
const auto eta = v0.eta();
const auto tranRad = v0.v0radius();
const auto dcaDaughv0 = v0.dcaV0daughters();
const auto cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const auto cpav0 = v0.v0cosPA();
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};

float observable = 0.;
Expand Down Expand Up @@ -680,11 +680,11 @@ void FemtoDreamV0Selection::fillQA(C const& col, V const& v0, T const& posTrack,
mQAHistogramRegistry->fill(
HIST(o2::aod::femtodreamparticle::ParticleTypeName[part]) +
HIST("/hCPA"),
v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
v0.v0cosPA());
mQAHistogramRegistry->fill(
HIST(o2::aod::femtodreamparticle::ParticleTypeName[part]) +
HIST("/hCPAvsPt"),
v0.pt(), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
v0.pt(), v0.v0cosPA());
mQAHistogramRegistry->fill(
HIST(o2::aod::femtodreamparticle::ParticleTypeName[part]) +
HIST("/hInvMassLambda"),
Expand Down
2 changes: 1 addition & 1 deletion PWGCF/FemtoDream/femtoDreamProducerTask.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -483,7 +483,7 @@ struct femtoDreamProducerTask {
aod::femtodreamparticle::ParticleType::kV0,
cutContainerV0.at(femtoDreamV0Selection::V0ContainerPosition::kV0),
0,
v0.v0cosPA(col.posX(), col.posY(), col.posZ()),
v0.v0cosPA(),
indexChildID,
v0.mLambda(),
v0.mAntiLambda());
Expand Down
2 changes: 1 addition & 1 deletion PWGCF/FemtoDream/femtoDreamProducerTaskV0Only.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -442,7 +442,7 @@ struct femtoDreamProducerTaskV0Only {
aod::femtodreamparticle::ParticleType::kV0,
cutContainerV0.at(
femtoDreamV0Selection::V0ContainerPosition::kV0),
0, v0.v0cosPA(col.posX(), col.posY(), col.posZ()),
0, v0.v0cosPA(),
indexChildID, v0.mLambda(), v0.mAntiLambda());
if (ConfDebugOutput) {
outputDebugParts(
Expand Down
10 changes: 5 additions & 5 deletions PWGCF/FemtoUniverse/Core/FemtoUniverseV0Selection.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -415,7 +415,7 @@ bool FemtoUniverseV0Selection::isSelectedMinimal(C const& col, V const& v0,
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();
const float invMassAntiLambda = v0.mAntiLambda();
Expand Down Expand Up @@ -498,7 +498,7 @@ void FemtoUniverseV0Selection::fillLambdaQA(C const& col, V const& v0,
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();

Expand Down Expand Up @@ -587,7 +587,7 @@ std::array<cutContainerType, 5>
const auto eta = v0.eta();
const auto tranRad = v0.v0radius();
const auto dcaDaughv0 = v0.dcaV0daughters();
const auto cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const auto cpav0 = v0.v0cosPA();
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};

float observable = 0.;
Expand Down Expand Up @@ -680,11 +680,11 @@ void FemtoUniverseV0Selection::fillQA(C const& col, V const& v0, T const& posTra
mHistogramRegistry->fill(
HIST(o2::aod::femtouniverseparticle::ParticleTypeName[part]) +
HIST("/hCPA"),
v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
v0.v0cosPA());
mHistogramRegistry->fill(
HIST(o2::aod::femtouniverseparticle::ParticleTypeName[part]) +
HIST("/hCPAvsPt"),
v0.pt(), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
v0.pt(), v0.v0cosPA());
mHistogramRegistry->fill(
HIST(o2::aod::femtouniverseparticle::ParticleTypeName[part]) +
HIST("/hInvMassLambda"),
Expand Down
4 changes: 2 additions & 2 deletions PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTask.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -664,7 +664,7 @@ struct femtoUniverseProducerTask {
aod::femtouniverseparticle::ParticleType::kV0,
cutContainerV0.at(femtoUniverseV0Selection::V0ContainerPosition::kV0),
0,
v0.v0cosPA(col.posX(), col.posY(), col.posZ()),
v0.v0cosPA(),
indexChildID,
v0.mLambda(),
v0.mAntiLambda());
Expand Down Expand Up @@ -780,7 +780,7 @@ struct femtoUniverseProducerTask {
aod::femtouniverseparticle::ParticleType::kPhi,
-999, // cutContainerV0.at(femtoUniverseV0Selection::V0ContainerPosition::kV0),
0,
phiM, // v0.v0cosPA(col.posX(), col.posY(), col.posZ()),
phiM, // v0.v0cosPA(),
indexChildID,
phiM, // phi.mLambda(), //for now it will have a mLambda getter, maybe we will change it in the future so it's more logical
-999); // v0.mAntiLambda()
Expand Down
2 changes: 1 addition & 1 deletion PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTaskV0Only.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -459,7 +459,7 @@ struct femtoUniverseProducerTaskV0Only {
aod::femtouniverseparticle::ParticleType::kV0,
cutContainerV0.at(
femtoUniverseV0Selection::V0ContainerPosition::kV0),
0, v0.v0cosPA(col.posX(), col.posY(), col.posZ()),
0, v0.v0cosPA(),
indexChildID, v0.mLambda(), v0.mAntiLambda());
if (ConfDebugOutput) {
outputDebugParts(
Expand Down
10 changes: 5 additions & 5 deletions PWGCF/FemtoWorld/Core/FemtoWorldPhiSelection.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -400,7 +400,7 @@ bool FemtoWorldPhiSelection::isSelectedMinimal(C const& col, V const& v0, T cons
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaPhidaughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassPhi = v0.mPhi();
const float invMassAntiPhi = v0.mAntiPhi();
Expand Down Expand Up @@ -472,7 +472,7 @@ void FemtoWorldPhiSelection::fillPhiQA(C const& col, V const& v0, T const& posTr
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaPhidaughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassPhi = v0.mPhi();

Expand Down Expand Up @@ -570,7 +570,7 @@ std::array<cutContainerType, 5> FemtoWorldPhiSelection::getCutContainer(C const&
const auto pT = sumVec.Pt();
// const auto tranRad = v0.v0radius();
// const auto dcaDaughv0 = v0.dcaPhidaughters();
// const auto cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
// const auto cpav0 = v0.v0cosPA();
const std::vector<float> decVtx = {posTrack.x(), posTrack.y(), posTrack.z()};

float observable = 0.;
Expand Down Expand Up @@ -663,8 +663,8 @@ void FemtoWorldPhiSelection::fillQAPhi(C const& col, V const& v0, T const& posTr
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hDecayVtxX"), v0.x());
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hDecayVtxY"), v0.y());
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hDecayVtxZ"), v0.z());
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPAvsPt"), v0.pt(), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPA"), v0.v0cosPA());
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPAvsPt"), v0.pt(), v0.v0cosPA());
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hInvMassPhi"), v0.mPhi());
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hInvMassAntiPhi"), v0.mAntiPhi());
// mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hInvMassPhiAntiPhi"), v0.mPhi(), v0.mAntiPhi());
Expand Down
10 changes: 5 additions & 5 deletions PWGCF/FemtoWorld/Core/FemtoWorldV0Selection.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -328,7 +328,7 @@ bool FemtoWorldV0Selection::isSelectedMinimal(C const& col, V const& v0, T const
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();
const float invMassAntiLambda = v0.mAntiLambda();
Expand Down Expand Up @@ -400,7 +400,7 @@ void FemtoWorldV0Selection::fillLambdaQA(C const& col, V const& v0, T const& pos
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};
const float tranRad = v0.v0radius();
const float dcaDaughv0 = v0.dcaV0daughters();
const float cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const float cpav0 = v0.v0cosPA();

const float invMassLambda = v0.mLambda();

Expand Down Expand Up @@ -476,7 +476,7 @@ std::array<cutContainerType, 5> FemtoWorldV0Selection::getCutContainer(C const&
const auto pT = v0.pt();
const auto tranRad = v0.v0radius();
const auto dcaDaughv0 = v0.dcaV0daughters();
const auto cpav0 = v0.v0cosPA(col.posX(), col.posY(), col.posZ());
const auto cpav0 = v0.v0cosPA();
const std::vector<float> decVtx = {v0.x(), v0.y(), v0.z()};

float observable = 0.;
Expand Down Expand Up @@ -531,8 +531,8 @@ void FemtoWorldV0Selection::fillQA(C const& col, V const& v0, T const& posTrack,
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hDecayVtxX"), v0.x());
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hDecayVtxY"), v0.y());
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hDecayVtxZ"), v0.z());
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPA"), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPAvsPt"), v0.pt(), v0.v0cosPA(col.posX(), col.posY(), col.posZ()));
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPA"), v0.v0cosPA());
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hCPAvsPt"), v0.pt(), v0.v0cosPA());
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hInvMassLambda"), v0.mLambda());
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hInvMassAntiLambda"), v0.mAntiLambda());
mHistogramRegistry->fill(HIST(o2::aod::femtoworldparticle::ParticleTypeName[part]) + HIST("/hInvMassLambdaAntiLambda"), v0.mLambda(), v0.mAntiLambda());
Expand Down
Loading

0 comments on commit 0f0165c

Please sign in to comment.