Skip to content

Commit aa958e0

Browse files
committed
fix for filling MC particle info in run 3 + toggle additional ev sels
1 parent 71b5d25 commit aa958e0

File tree

1 file changed

+6
-12
lines changed

1 file changed

+6
-12
lines changed

PWGLF/TableProducer/Nuspex/ebyeMaker.cxx

Lines changed: 6 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -230,7 +230,7 @@ struct EbyeMaker {
230230
Configurable<float> etaMaxV0dau{"etaMaxV0dau", 0.8f, "maximum eta V0 daughters"};
231231
Configurable<float> outerPIDMin{"outerPIDMin", -4.f, "minimum outer PID"};
232232

233-
Configurable<std::string> genName{"genname", "", "Genearator name: HIJING, PYTHIA8, ... Default: \"\""};
233+
Configurable<bool> useAllEvSel{"useAllEvSel", false, "use additional event selections fo run 3 analyses"};
234234
Configurable<uint8_t> triggerCut{"triggerCut", 0x0, "trigger cut to select"};
235235
Configurable<bool> kINT7Intervals{"kINT7Intervals", false, "toggle kINT7 trigger selection in the 10-30% and 50-90% centrality intervals (2018 Pb-Pb)"};
236236
Configurable<bool> kUsePileUpCut{"kUsePileUpCut", false, "toggle strong correlation cuts (Run 2)"};
@@ -816,7 +816,7 @@ struct EbyeMaker {
816816
auto mcTrack = mcLab.template mcParticle_as<aod::McParticles>();
817817
if (std::abs(mcTrack.pdgCode()) != kPartPdg[iP])
818818
continue;
819-
if (((mcTrack.flags() & 0x8) && (doprocessMcRun2 || doprocessMiniMcRun2)) || (mcTrack.flags() & 0x2) || ((mcTrack.flags() & 0x1) && !doprocessMiniMcRun2))
819+
if ((((mcTrack.flags() & 0x8) || (mcTrack.flags() & 0x2)) && (doprocessMcRun2 || doprocessMiniMcRun2)) || ((mcTrack.flags() & 0x1) && !doprocessMiniMcRun2))
820820
continue;
821821

822822
if (!mcTrack.isPhysicalPrimary() && !doprocessMiniMcRun2)
@@ -854,7 +854,7 @@ struct EbyeMaker {
854854
}
855855
if (!posMother.isPhysicalPrimary() && !posMother.has_mothers())
856856
continue;
857-
if (((posMother.flags() & 0x8) && (doprocessMcRun2 || doprocessMiniMcRun2)) || (posMother.flags() & 0x2) || (posMother.flags() & 0x1))
857+
if (((posMother.flags() & 0x8) || (posMother.flags() & 0x2) || (posMother.flags() & 0x1)) && (doprocessMcRun2 || doprocessMiniMcRun2))
858858
continue;
859859

860860
auto genPt = std::hypot(posMother.px(), posMother.py());
@@ -877,7 +877,7 @@ struct EbyeMaker {
877877
if (std::abs(genEta) > etaMax) {
878878
continue;
879879
}
880-
if (((mcPart.flags() & 0x8) && (doprocessMcRun2 || doprocessMiniMcRun2)) || (mcPart.flags() & 0x2) || ((mcPart.flags() & 0x1) && !doprocessMiniMcRun2))
880+
if ((((mcPart.flags() & 0x8) || (mcPart.flags() & 0x2)) && (doprocessMcRun2 || doprocessMiniMcRun2)) || ((mcPart.flags() & 0x1) && !doprocessMiniMcRun2))
881881
continue;
882882
auto pdgCode = mcPart.pdgCode();
883883
if (std::abs(pdgCode) == PDG_t::kLambda0) {
@@ -938,10 +938,7 @@ struct EbyeMaker {
938938
auto bc = collision.bc_as<aod::BCsWithTimestamps>();
939939
initCCDB(bc);
940940

941-
if (!collision.sel8() || std::abs(collision.posZ()) > zVtxMax)
942-
continue;
943-
944-
if (!collision.selection_bit(aod::evsel::kNoITSROFrameBorder) || !collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))
941+
if (std::abs(collision.posZ()) > zVtxMax || !collision.selection_bit(aod::evsel::kNoITSROFrameBorder) || !collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kIsTriggerTVX) || ((!collision.selection_bit(aod::evsel::kIsGoodITSLayersAll) || !collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV)) && useAllEvSel))
945942
continue;
946943

947944
histos.fill(HIST("QA/zVtx"), collision.posZ());
@@ -1119,10 +1116,7 @@ struct EbyeMaker {
11191116
auto bc = collision.bc_as<aod::BCsWithTimestamps>();
11201117
initCCDB(bc);
11211118

1122-
if (!collision.sel8() || std::abs(collision.posZ()) > zVtxMax)
1123-
continue;
1124-
1125-
if (!collision.selection_bit(aod::evsel::kNoITSROFrameBorder) || !collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))
1119+
if (std::abs(collision.posZ()) > zVtxMax || !collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kIsTriggerTVX) || ((!collision.selection_bit(aod::evsel::kIsGoodITSLayersAll) || !collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV)) && useAllEvSel))
11261120
continue;
11271121

11281122
auto centrality = collision.centFT0C();

0 commit comments

Comments
 (0)