allow GenKt in exclusive_jets

JuliaHEP · Feb 3, 2025 · 378f520 · 378f520
1 parent a904f6c
commit 378f520
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/src/ClusterSequence.jl b/src/ClusterSequence.jl
@@ -321,8 +321,10 @@ function exclusive_jets(clusterseq::ClusterSequence{U}; dcut = nothing, njets =
     end
 
     # Check that an algorithm was used that makes sense for exclusive jets
-    if !(clusterseq.algorithm ∈
-         (JetAlgorithm.CA, JetAlgorithm.Kt, JetAlgorithm.EEKt, JetAlgorithm.Durham))
+    if clusterseq.algorithm == JetAlgorithm.GenKt && clusterseq.power < 0
+        throw(ArgumentError("Algorithm GenKt requires power >= 0 (power=$(clusterseq.power))"))
+    elseif clusterseq.algorithm ∉
+           (JetAlgorithm.CA, JetAlgorithm.Kt, JetAlgorithm.EEKt, JetAlgorithm.Durham)
         throw(ArgumentError("Algorithm used is not suitable for exclusive jets ($(clusterseq.algorithm))"))
     end
 
@@ -378,8 +380,10 @@ n_exclusive_jets(clusterseq, dcut = 20.0)
 """
 function n_exclusive_jets(clusterseq::ClusterSequence; dcut::AbstractFloat)
     # Check that an algorithm was used that makes sense for exclusive jets
-    if !(clusterseq.algorithm ∈
-         (JetAlgorithm.CA, JetAlgorithm.Kt, JetAlgorithm.EEKt, JetAlgorithm.Durham))
+    if clusterseq.algorithm == JetAlgorithm.GenKt && clusterseq.power < 0
+        throw(ArgumentError("Algorithm GenKt requires power >= 0 (power=$(clusterseq.power))"))
+    elseif clusterseq.algorithm ∉
+           (JetAlgorithm.CA, JetAlgorithm.Kt, JetAlgorithm.EEKt, JetAlgorithm.Durham)
         throw(ArgumentError("Algorithm used is not suitable for exclusive jets ($(clusterseq.algorithm))"))
     end