[PWGJE] Setting the phi selection configs (#14888)

Author: nzardosh

PWGJE/Core/JetFinder.h

@@ -48,15 +48,15 @@ class JetFinder
   /// \return ClusterSequenceArea object needed to access constituents
   // fastjet::ClusterSequenceArea findJets(std::vector<fastjet::PseudoJet> &inputParticles, std::vector<fastjet::PseudoJet> &jets);
 
-  float phiMin = -2. * M_PI;
+  float phiMin = -1. * M_PI;
   float phiMax = 2. * M_PI;
   float etaMin = -.9;
   float etaMax = .9;
 
   float jetR = .4;
   float jetPtMin = 0.;
   float jetPtMax = 1000.;
-  float jetPhiMin = -2. * M_PI;
+  float jetPhiMin = -1. * M_PI;
   float jetPhiMax = 2. * M_PI;
   float jetEtaMin = -99.;
   float jetEtaMax = 99.;

PWGJE/JetFinders/jetFinder.h

@@ -137,9 +137,15 @@ struct JetFinderTask {
     jetFinder.etaMax = trackEtaMax;
     jetFinder.phiMin = trackPhiMin;
     jetFinder.phiMax = trackPhiMax;
+    if (trackPhiMin < -98.0) {
+      jetFinder.phiMin = -1.0 * M_PI;
+      jetFinder.phiMax = 2.0 * M_PI;
+    }
+    jetFinder.jetPhiMin = jetPhiMin;
+    jetFinder.jetPhiMax = jetPhiMax;
     if (jetPhiMin < -98.0) {
-      jetFinder.jetPhiMin = -2. * M_PI;
-      jetFinder.jetPhiMax = 2. * M_PI;
+      jetFinder.jetPhiMin = -1.0 * M_PI;
+      jetFinder.jetPhiMax = 2.0 * M_PI;
     }
     jetFinder.jetEtaMin = jetEtaMin;
     jetFinder.jetEtaMax = jetEtaMax;

PWGJE/JetFinders/jetFinderHF.h

@@ -142,9 +142,17 @@ struct JetFinderHFTask {
     jetFinder.etaMax = trackEtaMax;
     jetFinder.jetPtMin = jetPtMin;
     jetFinder.jetPtMax = jetPtMax;
+    jetFinder.phiMin = trackPhiMin;
+    jetFinder.phiMax = trackPhiMax;
+    if (trackPhiMin < -98.0) {
+      jetFinder.phiMin = -1.0 * M_PI;
+      jetFinder.phiMax = 2.0 * M_PI;
+    }
+    jetFinder.jetPhiMin = jetPhiMin;
+    jetFinder.jetPhiMax = jetPhiMax;
     if (jetPhiMin < -98.0) {
-      jetFinder.jetPhiMin = -2. * M_PI;
-      jetFinder.jetPhiMax = 2. * M_PI;
+      jetFinder.jetPhiMin = -1.0 * M_PI;
+      jetFinder.jetPhiMax = 2.0 * M_PI;
     }
     jetFinder.jetEtaMin = jetEtaMin;
     jetFinder.jetEtaMax = jetEtaMax;

PWGJE/JetFinders/jetFinderHFHFBar.h

@@ -14,8 +14,8 @@
 /// \author Nima Zardoshti <nima.zardoshti@cern.ch>
 /// \author Jochen Klein <jochen.klein@cern.ch>
 
-#ifndef PWGJE_JETFINDERS_JETFINDERHF_H_
-#define PWGJE_JETFINDERS_JETFINDERHF_H_
+#ifndef PWGJE_JETFINDERS_JETFINDERHFHFBAR_H_
+#define PWGJE_JETFINDERS_JETFINDERHFHFBAR_H_
 
 #include "PWGJE/Core/JetDerivedDataUtilities.h"
 #include "PWGJE/Core/JetFinder.h"
@@ -139,9 +139,17 @@ struct JetFinderHFHFBarTask {
     jetFinder.etaMax = trackEtaMax;
     jetFinder.jetPtMin = jetPtMin;
     jetFinder.jetPtMax = jetPtMax;
+    jetFinder.phiMin = trackPhiMin;
+    jetFinder.phiMax = trackPhiMax;
+    if (trackPhiMin < -98.0) {
+      jetFinder.phiMin = -1.0 * M_PI;
+      jetFinder.phiMax = 2.0 * M_PI;
+    }
+    jetFinder.jetPhiMin = jetPhiMin;
+    jetFinder.jetPhiMax = jetPhiMax;
     if (jetPhiMin < -98.0) {
-      jetFinder.jetPhiMin = -2. * M_PI;
-      jetFinder.jetPhiMax = 2. * M_PI;
+      jetFinder.jetPhiMin = -1.0 * M_PI;
+      jetFinder.jetPhiMax = 2.0 * M_PI;
     }
     jetFinder.jetEtaMin = jetEtaMin;
     jetFinder.jetEtaMax = jetEtaMax;
@@ -323,4 +331,4 @@ struct JetFinderHFHFBarTask {
   PROCESS_SWITCH(JetFinderHFHFBarTask, processChargedJetsMCP, "hf jet finding on MC particle level", false);
 };
 
-#endif // PWGJE_JETFINDERS_JETFINDERHF_H_
+#endif // PWGJE_JETFINDERS_JETFINDERHFHFBAR_H_

PWGJE/JetFinders/jetFinderV0.h

@@ -123,9 +123,17 @@ struct JetFinderV0Task {
     jetFinder.etaMax = trackEtaMax;
     jetFinder.jetPtMin = jetPtMin;
     jetFinder.jetPtMax = jetPtMax;
+    jetFinder.phiMin = trackPhiMin;
+    jetFinder.phiMax = trackPhiMax;
+    if (trackPhiMin < -98.0) {
+      jetFinder.phiMin = -1.0 * M_PI;
+      jetFinder.phiMax = 2.0 * M_PI;
+    }
+    jetFinder.jetPhiMin = jetPhiMin;
+    jetFinder.jetPhiMax = jetPhiMax;
     if (jetPhiMin < -98.0) {
-      jetFinder.jetPhiMin = -2. * M_PI;
-      jetFinder.jetPhiMax = 2. * M_PI;
+      jetFinder.jetPhiMin = -1.0 * M_PI;
+      jetFinder.jetPhiMax = 2.0 * M_PI;
     }
     jetFinder.jetEtaMin = jetEtaMin;
     jetFinder.jetEtaMax = jetEtaMax;