Skip to content

Commit 97baa8e

Browse files
authored
PWGHF: Revert use of partitions in skimming (#63)
1 parent 77e19f5 commit 97baa8e

1 file changed

Lines changed: 29 additions & 12 deletions

File tree

PWGHF/TableProducer/HFTrackIndexSkimsCreator.cxx

Lines changed: 29 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -773,8 +773,9 @@ struct HfTrackIndexSkimsCreator {
773773
using SelectedCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::HFSelCollision>>;
774774
using SelectedTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksCov, aod::TracksExtra, aod::HFSelTrack>>;
775775

776-
Partition<SelectedTracks> tracksPos = aod::track::signed1Pt > 0.f;
777-
Partition<SelectedTracks> tracksNeg = aod::track::signed1Pt < 0.f;
776+
// FIXME
777+
//Partition<SelectedTracks> tracksPos = aod::track::signed1Pt > 0.f;
778+
//Partition<SelectedTracks> tracksNeg = aod::track::signed1Pt < 0.f;
778779

779780
// int nColls{0}; //can be added to run over limited collisions per file - for tesing purposes
780781

@@ -833,12 +834,16 @@ struct HfTrackIndexSkimsCreator {
833834
auto nCand3 = rowTrackIndexProng3.lastIndex();
834835

835836
// if there isn't at least a positive and a negative track, continue immediately
836-
if (tracksPos.size() < 1 || tracksNeg.size() < 1) {
837-
return;
838-
}
837+
//if (tracksPos.size() < 1 || tracksNeg.size() < 1) {
838+
// return;
839+
//}
839840

840841
// first loop over positive tracks
841-
for (auto trackPos1 = tracksPos.begin(); trackPos1 != tracksPos.end(); ++trackPos1) {
842+
//for (auto trackPos1 = tracksPos.begin(); trackPos1 != tracksPos.end(); ++trackPos1) {
843+
for (auto trackPos1 = tracks.begin(); trackPos1 != tracks.end(); ++trackPos1) {
844+
if (trackPos1.signed1Pt() < 0) {
845+
continue;
846+
}
842847
bool sel2ProngStatusPos = TESTBIT(trackPos1.isSelProng(), CandidateType::Cand2Prong);
843848
bool sel3ProngStatusPos1 = TESTBIT(trackPos1.isSelProng(), CandidateType::Cand3Prong);
844849
if (!sel2ProngStatusPos && !sel3ProngStatusPos1) {
@@ -848,7 +853,11 @@ struct HfTrackIndexSkimsCreator {
848853
auto trackParVarPos1 = getTrackParCov(trackPos1);
849854

850855
// first loop over negative tracks
851-
for (auto trackNeg1 = tracksNeg.begin(); trackNeg1 != tracksNeg.end(); ++trackNeg1) {
856+
//for (auto trackNeg1 = tracksNeg.begin(); trackNeg1 != tracksNeg.end(); ++trackNeg1) {
857+
for (auto trackNeg1 = tracks.begin(); trackNeg1 != tracks.end(); ++trackNeg1) {
858+
if (trackNeg1.signed1Pt() > 0) {
859+
continue;
860+
}
852861
bool sel2ProngStatusNeg = TESTBIT(trackNeg1.isSelProng(), CandidateType::Cand2Prong);
853862
bool sel3ProngStatusNeg1 = TESTBIT(trackNeg1.isSelProng(), CandidateType::Cand3Prong);
854863
if (!sel2ProngStatusNeg && !sel3ProngStatusNeg1) {
@@ -945,11 +954,15 @@ struct HfTrackIndexSkimsCreator {
945954
continue;
946955
}
947956

948-
// if (tracks.size() < 2) {
949-
// continue;
950-
// }
957+
if (tracks.size() < 2) {
958+
continue;
959+
}
951960
// second loop over positive tracks
952-
for (auto trackPos2 = trackPos1 + 1; trackPos2 != tracksPos.end(); ++trackPos2) {
961+
//for (auto trackPos2 = trackPos1 + 1; trackPos2 != tracksPos.end(); ++trackPos2) {
962+
for (auto trackPos2 = trackPos1 + 1; trackPos2 != tracks.end(); ++trackPos2) {
963+
if (trackPos2.signed1Pt() < 0) {
964+
continue;
965+
}
953966
if (!TESTBIT(trackPos2.isSelProng(), CandidateType::Cand3Prong)) {
954967
continue;
955968
}
@@ -1055,7 +1068,11 @@ struct HfTrackIndexSkimsCreator {
10551068
}
10561069

10571070
// second loop over negative tracks
1058-
for (auto trackNeg2 = trackNeg1 + 1; trackNeg2 != tracksNeg.end(); ++trackNeg2) {
1071+
//for (auto trackNeg2 = trackNeg1 + 1; trackNeg2 != tracksNeg.end(); ++trackNeg2) {
1072+
for (auto trackNeg2 = trackNeg1 + 1; trackNeg2 != tracks.end(); ++trackNeg2) {
1073+
if (trackNeg2.signed1Pt() > 0) {
1074+
continue;
1075+
}
10591076
if (!TESTBIT(trackNeg2.isSelProng(), CandidateType::Cand3Prong)) {
10601077
continue;
10611078
}

0 commit comments

Comments
 (0)