Skip to content

Commit 75e6ca4

Browse files
author
Christopher Klumm
committed
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 bf75c17 commit 75e6ca4

4 files changed

Lines changed: 20 additions & 436 deletions

File tree

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)