diff --git a/PWGJE/Core/JetFinder.h b/PWGJE/Core/JetFinder.h index af78c69eea2..8f727bbdbbd 100644 --- a/PWGJE/Core/JetFinder.h +++ b/PWGJE/Core/JetFinder.h @@ -48,7 +48,7 @@ class JetFinder /// \return ClusterSequenceArea object needed to access constituents // fastjet::ClusterSequenceArea findJets(std::vector &inputParticles, std::vector &jets); - float phiMin = -2. * M_PI; + float phiMin = -1. * M_PI; float phiMax = 2. * M_PI; float etaMin = -.9; float etaMax = .9; @@ -56,7 +56,7 @@ class JetFinder 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.; diff --git a/PWGJE/JetFinders/jetFinder.h b/PWGJE/JetFinders/jetFinder.h index 123d0e1b4a4..72c370e44ac 100644 --- a/PWGJE/JetFinders/jetFinder.h +++ b/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; diff --git a/PWGJE/JetFinders/jetFinderHF.h b/PWGJE/JetFinders/jetFinderHF.h index b44b4d497f2..cb60dae3e09 100644 --- a/PWGJE/JetFinders/jetFinderHF.h +++ b/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; diff --git a/PWGJE/JetFinders/jetFinderHFHFBar.h b/PWGJE/JetFinders/jetFinderHFHFBar.h index ee29349a9cf..46200d87356 100644 --- a/PWGJE/JetFinders/jetFinderHFHFBar.h +++ b/PWGJE/JetFinders/jetFinderHFHFBar.h @@ -14,8 +14,8 @@ /// \author Nima Zardoshti /// \author Jochen Klein -#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_ diff --git a/PWGJE/JetFinders/jetFinderV0.h b/PWGJE/JetFinders/jetFinderV0.h index 89ea38d3d96..82272837509 100644 --- a/PWGJE/JetFinders/jetFinderV0.h +++ b/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;