From 97046e335f227eee45938a26afcc8095ed7777dc Mon Sep 17 00:00:00 2001 From: Graeme A Stewart Date: Fri, 6 Dec 2024 16:50:32 +0100 Subject: [PATCH] Fix for jetreco.jl jet types We explicitly use Float64 here, because this is deliberately the simple case. --- examples/jetreco.jl | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/examples/jetreco.jl b/examples/jetreco.jl index 82a6a11..ae64327 100644 --- a/examples/jetreco.jl +++ b/examples/jetreco.jl @@ -140,17 +140,18 @@ function main() logger = ConsoleLogger(stdout, Logging.Info) global_logger(logger) # Try to read events into the correct type! - if JetReconstruction.is_ee(args[:algorithm]) - jet_type = EEjet + # If we don't have an algorithm we default to PseudoJet + if !isnothing(args[:algorithm]) JetReconstruction.is_ee(args[:algorithm]) + jet_type = EEjet{Float64} else - jet_type = PseudoJet + jet_type = PseudoJet{Float64} end events::Vector{Vector{jet_type}} = read_final_state_particles(args[:file], maxevents = args[:maxevents], skipevents = args[:skip], T = jet_type) if isnothing(args[:algorithm]) && isnothing(args[:power]) - @warn "Neither algorithm nor power specified, defaulting to AntiKt" + @warn "Neither algorithm nor power specified, defaulting to pp event AntiKt" args[:algorithm] = JetAlgorithm.AntiKt end jet_process(events, distance = args[:distance], algorithm = args[:algorithm],