diff --git a/PWGLF/Tasks/Resonances/phianalysisrun3.cxx b/PWGLF/Tasks/Resonances/phianalysisrun3.cxx
index a865482393a..7cd608e955c 100644
--- a/PWGLF/Tasks/Resonances/phianalysisrun3.cxx
+++ b/PWGLF/Tasks/Resonances/phianalysisrun3.cxx
@@ -56,6 +56,7 @@ struct phianalysisrun3 {
   SliceCache cache;
   HistogramRegistry histos{"histos", {}, OutputObjHandlingPolicy::AnalysisObject};
   // events
+  Configurable<bool> applyEvsel{"applyEvsel", false, "applyEvsel"};
   Configurable<float> cfgCutVertex{"cfgCutVertex", 10.0f, "Accepted z-vertex range"};
   // track
   Configurable<float> cfgCutPT{"cfgCutPT", 0.2, "PT cut on daughter track"};
@@ -252,7 +253,7 @@ struct phianalysisrun3 {
 
   void processSameEvent(EventCandidates::iterator const& collision, TrackCandidates const& tracks, aod::BCs const&)
   {
-    if (!collision.sel8()) {
+    if (applyEvsel && !collision.sel8()) {
       return;
     }
     if (timFrameEvsel && (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoITSROFrameBorder))) {
@@ -309,10 +310,10 @@ struct phianalysisrun3 {
     BinningTypeVertexContributor binningOnPositions{{axisVertex, axisMultiplicity}, true};
     SameKindPair<EventCandidates, TrackCandidates, BinningTypeVertexContributor> pair{binningOnPositions, cfgNoMixedEvents, -1, collisions, tracksTuple, &cache};
     for (auto& [c1, tracks1, c2, tracks2] : pair) {
-      if (!c1.sel8()) {
+      if (applyEvsel && !c1.sel8()) {
         continue;
       }
-      if (!c2.sel8()) {
+      if (applyEvsel && !c2.sel8()) {
         continue;
       }
       if (timFrameEvsel && (!c1.selection_bit(aod::evsel::kNoTimeFrameBorder) || !c2.selection_bit(aod::evsel::kNoTimeFrameBorder) || !c1.selection_bit(aod::evsel::kNoITSROFrameBorder) || !c2.selection_bit(aod::evsel::kNoITSROFrameBorder))) {