Skip to content

Commit 6d5351d

Browse files
[PWGCF] added removed workflow in CMakeList; electron rejection cut added in PID (#9647)
1 parent fb4e6b5 commit 6d5351d

File tree

2 files changed

+13
-0
lines changed

2 files changed

+13
-0
lines changed

PWGCF/EbyEFluctuations/Tasks/CMakeLists.txt

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -24,6 +24,11 @@ o2physics_add_dpl_workflow(netproton-cumulants
2424
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::PWGCFCore
2525
COMPONENT_NAME Analysis)
2626

27+
o2physics_add_dpl_workflow(netproton-cumulants-mc
28+
SOURCES netprotonCumulantsMc.cxx
29+
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::PWGCFCore
30+
COMPONENT_NAME Analysis)
31+
2732
o2physics_add_dpl_workflow(identified-meanpt-fluctuations
2833
SOURCES IdentifiedMeanPtFluctuations.cxx
2934
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::PWGCFCore

PWGCF/EbyEFluctuations/Tasks/netprotonCumulantsMc.cxx

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -816,6 +816,11 @@ struct NetprotonCumulantsMc {
816816
template <typename T>
817817
bool selectionPIDnew(const T& candidate)
818818
{
819+
// electron rejection
820+
if (candidate.tpcNSigmaEl() > -3.0f && candidate.tpcNSigmaEl() < 5.0f && std::abs(candidate.tpcNSigmaPi()) > 3.0f && std::abs(candidate.tpcNSigmaKa()) > 3.0f && std::abs(candidate.tpcNSigmaPr()) > 3.0f) {
821+
return false;
822+
}
823+
819824
//! if pt < threshold
820825
if (candidate.pt() > 0.2f && candidate.pt() <= cfgCutPtUpperTPC) {
821826
if (!candidate.hasTOF() && std::abs(candidate.tpcNSigmaPr()) < cfgnSigmaCutTPC && std::abs(candidate.tpcNSigmaPi()) > cfgnSigmaCutTPC && std::abs(candidate.tpcNSigmaKa()) > cfgnSigmaCutTPC) {
@@ -1070,6 +1075,9 @@ struct NetprotonCumulantsMc {
10701075
{
10711076
continue;
10721077
}
1078+
if ((track.pt() < cfgCutPtLower) || (track.pt() > 5.0f) || (std::abs(track.eta()) > 0.8f)) {
1079+
continue;
1080+
}
10731081

10741082
histos.fill(HIST("hrecPtAll"), track.pt());
10751083
histos.fill(HIST("hrecEtaAll"), track.eta());

0 commit comments

Comments
 (0)