Skip to content

Commit d1101e2

Browse files
authored
[PWGLF] [PWGMM] Modification to the event selection criteria (#8672)
1 parent d3dd8a4 commit d1101e2

File tree

1 file changed

+9
-17
lines changed

1 file changed

+9
-17
lines changed

PWGMM/UE/Tasks/uecharged.cxx

Lines changed: 9 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -53,7 +53,7 @@ struct ueCharged {
5353
selectedTracks.SetRequireTPCRefit(true);
5454
// selectedTracks.SetRequireGoldenChi2(true);
5555
selectedTracks.SetMinNCrossedRowsTPC(70);
56-
selectedTracks.SetMinNCrossedRowsOverFindableClustersTPC(0.4f);
56+
selectedTracks.SetMinNCrossedRowsOverFindableClustersTPC(0.8f);
5757
selectedTracks.SetMaxChi2PerClusterTPC(4.f);
5858
selectedTracks.SetRequireHitsInITSLayers(1, {0, 1}); // one hit in any SPD layer
5959
selectedTracks.SetMaxChi2PerClusterITS(36.f);
@@ -70,7 +70,7 @@ struct ueCharged {
7070
selectedTracks.SetRequireTPCRefit(true);
7171
// selectedTracks.SetRequireGoldenChi2(true);
7272
selectedTracks.SetMinNCrossedRowsTPC(70);
73-
selectedTracks.SetMinNCrossedRowsOverFindableClustersTPC(0.4f);
73+
selectedTracks.SetMinNCrossedRowsOverFindableClustersTPC(0.8f);
7474
selectedTracks.SetMaxChi2PerClusterTPC(4.f);
7575
selectedTracks.SetRequireHitsInITSLayers(1, {0, 1}); // one hit in any SPD layer
7676
selectedTracks.SetMaxChi2PerClusterITS(36.f);
@@ -86,7 +86,6 @@ struct ueCharged {
8686
float DeltaPhi(float phia, float phib, float rangeMin, float rangeMax);
8787
// Configurable for event selection
8888
Configurable<bool> isRun3{"isRun3", true, "is Run3 dataset"};
89-
Configurable<bool> timeEvsel{"timeEvsel", true, "TPC Time frame boundary cut"};
9089
Configurable<bool> piluprejection{"piluprejection", true, "Pileup rejection"};
9190
Configurable<bool> goodzvertex{"goodzvertex", true, "removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference."};
9291

@@ -431,21 +430,16 @@ void ueCharged::processMeas(const C& collision, const T& tracks)
431430

432431
ue.fill(HIST("hCounter"), 1);
433432

434-
if (timeEvsel && (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoITSROFrameBorder))) {
435-
return;
436-
}
437-
438-
ue.fill(HIST("hCounter"), 2);
439433
if (piluprejection && !collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
440434
return;
441435
}
442436

443-
ue.fill(HIST("hCounter"), 3);
437+
ue.fill(HIST("hCounter"), 2);
444438
if (goodzvertex && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
445439
return;
446440
}
447441

448-
ue.fill(HIST("hCounter"), 4);
442+
ue.fill(HIST("hCounter"), 3);
449443

450444
ue.fill(HIST("hStat"), collision.size());
451445
auto vtxZ = collision.posZ();
@@ -454,7 +448,7 @@ void ueCharged::processMeas(const C& collision, const T& tracks)
454448
return;
455449
}
456450

457-
ue.fill(HIST("hCounter"), 5);
451+
ue.fill(HIST("hCounter"), 4);
458452

459453
ue.fill(HIST("hvtxZ"), vtxZ);
460454

@@ -726,28 +720,26 @@ void ueCharged::processMeasMC(const C& collision, const T& tracks, const P& part
726720
phiArrayTrue.clear();
727721
indexArrayTrue.clear();
728722

729-
if (!collision.sel8()) {
723+
if (!collision.selection_bit(o2::aod::evsel::kIsTriggerTVX) || !collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
730724
return;
731725
}
732726
ue.fill(HIST("hCounter"), 1);
733-
// TODO:Implement time frame selection (only if MC includes this effect)
734-
ue.fill(HIST("hCounter"), 2);
735727

736728
if (piluprejection && !collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
737729
return;
738730
}
731+
ue.fill(HIST("hCounter"), 2);
739732

740-
ue.fill(HIST("hCounter"), 3);
741733
if (goodzvertex && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
742734
return;
743735
}
744-
ue.fill(HIST("hCounter"), 4);
736+
ue.fill(HIST("hCounter"), 3);
745737

746738
// only PS
747739
if ((std::abs(collision.posZ()) >= 10.f)) {
748740
return;
749741
}
750-
ue.fill(HIST("hCounter"), 5);
742+
ue.fill(HIST("hCounter"), 4);
751743

752744
ue.fill(HIST(pNumDenTruePS[0]), flPtTrue, ue_true[0]);
753745
ue.fill(HIST(pSumPtTruePS[0]), flPtTrue, ue_true[3]);

0 commit comments

Comments
 (0)