Skip to content

Commit 21264c5

Browse files
Christopher Klummgmantzar
authored andcommitted
updated files for resonance selection
The partitions for particle selection in pair-task and debug-task were changed due to change in the way of how FDParticles is filled inside the producer task. The task was renamed and no longer needed files are removed.
1 parent 740c965 commit 21264c5

10 files changed

Lines changed: 60 additions & 2392 deletions

PWGCF/DataModel/FemtoDerived.h

Lines changed: 4 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -8,7 +8,6 @@
88
// In applying this license CERN does not waive the privileges and immunities
99
// granted to it by virtue of its status as an Intergovernmental Organization
1010
// or submit itself to any jurisdiction.
11-
1211
#ifndef PWGCF_DATAMODEL_FEMTODERIVED_H_
1312
#define PWGCF_DATAMODEL_FEMTODERIVED_H_
1413

@@ -98,15 +97,15 @@ enum ParticleType {
9897
kReso, //! Resonances (phi)
9998
kResoChild, // Child track of a Resonance
10099
kV0Child, //! Child track of a V0
101-
kPhiPosdaughTPC_NegdaughTPC, /// cases for Phi-daughters for TPC or TOF combinations
102-
kPhiPosdaughTPC_NegdaughTOF,
103-
kPhiPosdaughTOF_NegdaughTPC,
104-
kPhiPosdaughTOF_NegdaughTOF,
105100
kCascade, //! Cascade
106101
kCascadeV0,
107102
kCascadeV0Child,
108103
kCascadeBachelor, //! Bachelor track of a cascade
109104
kCharmHadron, //! Bachelor track of a cascade
105+
kPhiPosdaughTPC_NegdaughTPC, /// cases for Phi-daughters for TPC or TOF combinations
106+
kPhiPosdaughTPC_NegdaughTOF,
107+
kPhiPosdaughTOF_NegdaughTPC,
108+
kPhiPosdaughTOF_NegdaughTOF,
110109
kNParticleTypes //! Number of particle types
111110
};
112111

PWGCF/FemtoDream/Core/femtoDreamContainer.h

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -235,7 +235,7 @@ class FemtoDreamContainer
235235
if (extendedplots) {
236236
mHistogramRegistry->fill(HIST(mFolderSuffix[mEventType]) + HIST(o2::aod::femtodreamMCparticle::MCTypeName[mc]) + HIST("/relPairkstarmTPtPart1PtPart2MultPercentile"), femtoObs, mT, part1.pt(), part2.pt(), multPercentile);
237237

238-
if constexpr (std::is_same_v<T1, FDParticle> && std::is_same_v<T2, FDParticle>) {
238+
if constexpr (requires { part1.mLambda(); part2.mLambda(); }) {
239239
mHistogramRegistry->fill(HIST(mFolderSuffix[mEventType]) + HIST(o2::aod::femtodreamMCparticle::MCTypeName[mc]) + HIST("/pT1pT2kstarinvMassPart1invMassPart2"), part1.pt(), part2.pt(), femtoObs, part1.mLambda(), part2.mLambda());
240240
}
241241
}

PWGCF/FemtoDream/Core/femtoDreamResoSelection.h

Lines changed: 4 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -272,7 +272,7 @@ void FemtoDreamResoSelection::setDaughterCuts(femtoDreamResoSelection::Daughtert
272272
{
273273
if (daugh == femtoDreamResoSelection::kPosdaugh) {
274274
PosDaughTrack.setSelection(selVal, selVar, selType);
275-
};
275+
}
276276
if (daugh == femtoDreamResoSelection::kNegdaugh) {
277277
NegDaughTrack.setSelection(selVal, selVar, selType);
278278
}
@@ -346,25 +346,23 @@ template <typename cutContainerType, typename V>
346346
std::array<cutContainerType, 5> FemtoDreamResoSelection::getCutContainer(V const& track1, V const& track2, float sign)
347347
{
348348
cutContainerType outputSign = 0;
349+
cutContainerType outputPID = 0;
349350
size_t counter = 0;
350351
for (auto& sel : mSelections) { /// it should just be a 1D vector with sign
351352
const auto selVariable = sel.getSelectionVariable();
352353
if (selVariable == femtoDreamResoSelection::kResoSign) {
353354
sel.checkSelectionSetBit(sign, outputSign, counter, nullptr);
354355
}
355356
}
356-
cutContainerType outputPID = 0;
357357

358358
const auto DCA1 = std::sqrt(track1.dcaXY() * track1.dcaXY() + track1.dcaZ() * track1.dcaZ());
359359
const auto DCA2 = std::sqrt(track2.dcaXY() * track2.dcaXY() + track2.dcaZ() * track2.dcaZ());
360360

361361
auto outputPosTrack = PosDaughTrack.getCutContainer<false, cutContainerType>(track1, track1.pt(), track1.eta(), DCA1); // false for useItsPid
362362
auto outputNegTrack = NegDaughTrack.getCutContainer<false, cutContainerType>(track2, track2.pt(), track2.eta(), DCA2);
363363

364-
const auto shiftvalue = numBitsUsed(outputNegTrack.at(femtoDreamTrackSelection::TrackContainerPosition::kPID));
365-
outputPID = (outputPosTrack.at(femtoDreamTrackSelection::TrackContainerPosition::kPID) << shiftvalue) | outputNegTrack.at(femtoDreamTrackSelection::TrackContainerPosition::kPID);
366-
/// combine both outputs
367-
outputPID = (outputPID << counter) | outputSign;
364+
const auto shiftvalue = numBitsUsed(outputSign);
365+
outputPID = (outputNegTrack.at(femtoDreamTrackSelection::TrackContainerPosition::kPID) << shiftvalue) | outputSign;
368366

369367
std::array<cutContainerType, 5> bitmask = {outputPID,
370368
outputPosTrack.at(femtoDreamTrackSelection::TrackContainerPosition::kCuts),

0 commit comments

Comments
 (0)