You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
/// \brief task that creates the histograms necessary for computation of efficiency and purity functions in offline postprocess macros; also can make mcparticle and track QC histograms
15
15
16
16
#include<cmath>
17
17
#include<string>
@@ -85,6 +85,9 @@ struct TrackEfficiency {
85
85
86
86
std::vector<int> eventSelectionBits;
87
87
int trackSelection = -1;
88
+
int acceptSplitCollisionsCase0 = 0; // acceptSplitCollisions case: only look at mcCollisions that are not split
89
+
int acceptSplitCollisionsCase2 = 2; //acceptSplitCollisions case: accept split mcCollisions but only look at the first reco collision associated with it
90
+
// int acceptSplitCollisionsCase1 = 1; is in the else{} case
88
91
89
92
boolisChargedParticle(int code)
90
93
{
@@ -268,14 +271,14 @@ struct TrackEfficiency {
268
271
}
269
272
registry.fill(HIST("hMcCollCutsCounts"), 2.5); // mcCollisions with at least one reconstructed collision
270
273
271
-
if (acceptSplitCollisions == 0 && collisions.size() > 1) {
274
+
if (acceptSplitCollisions == acceptSplitCollisionsCase0 && collisions.size() > 1) {
if (acceptSplitCollisions == 2) { // check only that the first reconstructed collision passes the check
281
+
if (acceptSplitCollisions == acceptSplitCollisionsCase2) { // check only that the first reconstructed collision passes the check
279
282
if (jetderiveddatautilities::selectCollision(collisions.begin(), eventSelectionBits, skipMBGapEvents)) { // Skipping MC events that have not a single selected reconstructed collision ; effect unclear if mcColl is split
280
283
hasSel8Coll = true;
281
284
}
@@ -337,7 +340,7 @@ struct TrackEfficiency {
337
340
int splitCollCounter = 0;
338
341
for (autoconst& collision : collisions) {
339
342
splitCollCounter++;
340
-
if (acceptSplitCollisions == 2 && splitCollCounter > 1) {
343
+
if (acceptSplitCollisions == acceptSplitCollisionsCase2 && splitCollCounter > 1) {
341
344
return;
342
345
}
343
346
@@ -434,14 +437,14 @@ struct TrackEfficiency {
434
437
}
435
438
registry.fill(HIST("hMcCollCutsCounts"), 2.5); // mcCollisions with at least one reconstructed collision
436
439
437
-
if (acceptSplitCollisions == 0 && collisions.size() > 1) {
440
+
if (acceptSplitCollisions == acceptSplitCollisionsCase0 && collisions.size() > 1) {
if (acceptSplitCollisions == 2) { // check only that the first reconstructed collision passes the check
447
+
if (acceptSplitCollisions == acceptSplitCollisionsCase2) { // check only that the first reconstructed collision passes the check
445
448
if (jetderiveddatautilities::selectCollision(collisions.begin(), eventSelectionBits, skipMBGapEvents)) { // Skipping MC events that have not a single selected reconstructed collision ; effect unclear if mcColl is split
446
449
hasSel8Coll = true;
447
450
}
@@ -504,7 +507,7 @@ struct TrackEfficiency {
504
507
int splitCollCounter = 0;
505
508
for (autoconst& collision : collisions) {
506
509
splitCollCounter++;
507
-
if (acceptSplitCollisions == 2 && splitCollCounter > 1) {
510
+
if (acceptSplitCollisions == acceptSplitCollisionsCase2 && splitCollCounter > 1) {
508
511
return;
509
512
}
510
513
@@ -676,7 +679,7 @@ struct TrackEfficiency {
676
679
if (collisions.size() < 1) {
677
680
return;
678
681
}
679
-
if (acceptSplitCollisions == 0 && collisions.size() > 1) {
682
+
if (acceptSplitCollisions == acceptSplitCollisionsCase0 && collisions.size() > 1) {
680
683
return;
681
684
}
682
685
@@ -689,7 +692,7 @@ struct TrackEfficiency {
689
692
690
693
bool hasSel8Coll = false;
691
694
bool centralityCheck = false;
692
-
if (acceptSplitCollisions == 2) { // check only that the first reconstructed collision passes the check
695
+
if (acceptSplitCollisions == acceptSplitCollisionsCase2) { // check only that the first reconstructed collision passes the check
693
696
if (jetderiveddatautilities::selectCollision(collisions.begin(), eventSelectionBits, skipMBGapEvents)) { // Skipping MC events that have not a single selected reconstructed collision ; effect unclear if mcColl is split
694
697
hasSel8Coll = true;
695
698
}
@@ -737,7 +740,7 @@ struct TrackEfficiency {
737
740
if (collisions.size() < 1) {
738
741
return;
739
742
}
740
-
if (acceptSplitCollisions == 0 && collisions.size() > 1) {
743
+
if (acceptSplitCollisions == acceptSplitCollisionsCase0 && collisions.size() > 1) {
741
744
return;
742
745
}
743
746
@@ -750,7 +753,7 @@ struct TrackEfficiency {
750
753
751
754
bool hasSel8Coll = false;
752
755
bool centralityCheck = false;
753
-
if (acceptSplitCollisions == 2) { // check only that the first reconstructed collision passes the check
756
+
if (acceptSplitCollisions == acceptSplitCollisionsCase2) { // check only that the first reconstructed collision passes the check
754
757
if (jetderiveddatautilities::selectCollision(collisions.begin(), eventSelectionBits, skipMBGapEvents)) { // Skipping MC events that have not a single selected reconstructed collision ; effect unclear if mcColl is split
0 commit comments